Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,370,044 (T)9→8 intergenic (+41/‑172) pspE → / → ycjM thiosulfate:cyanide sulfurtransferase (rhodanese)/alpha amylase catalytic domain family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,370,0360T.100.0% 167.3 / NA 39intergenic (+33/‑180)pspE/ycjMthiosulfate:cyanide sulfurtransferase (rhodanese)/alpha amylase catalytic domain family protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (20/19);  total (20/19)

AAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTG  >  NC_000913/1369912‑1370169
                                                                                                                            |                                                                                                                                     
aagcaaAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTcac                                                                                                                        >  2:811389/1‑139 (MQ=255)
aagcaaAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTcac                                                                                                                        <  1:811389/139‑1 (MQ=255)
aagcaaAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATACACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTcac                                                                                                                        <  1:566636/139‑1 (MQ=255)
            tCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTATTTCTGTCACACCTTATTTACa                                                                                                            >  2:692826/1‑139 (MQ=255)
               ttAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATcc                                                                                                         <  1:51864/139‑1 (MQ=255)
                tAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTAc                                                                                                             >  1:159978/1‑134 (MQ=255)
                tAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTAc                                                                                                             <  2:159978/134‑1 (MQ=255)
                   cGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCAt                                                                                                     <  1:29849/139‑1 (MQ=255)
                                aCCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGt                                                                                        >  1:175090/1‑139 (MQ=255)
                                aCCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGt                                                                                        <  1:354193/139‑1 (MQ=255)
                                   cACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAg                                                                                     <  1:34286/139‑1 (MQ=255)
                                     cGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTca                                                                                                                         >  2:211866/1‑101 (MQ=255)
                                     cGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTca                                                                                                                         <  1:211866/101‑1 (MQ=255)
                                     cGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACAcc                                                                                                                     <  1:379958/105‑1 (MQ=255)
                                     cGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACAcc                                                                                                                     >  2:379958/1‑105 (MQ=255)
                                     cGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCt                                                                                   >  2:61811/1‑139 (MQ=255)
                                      gTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACa                                                                                                            >  2:709525/1‑113 (MQ=255)
                                      gTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACa                                                                                                            <  1:709525/113‑1 (MQ=255)
                                         gagaACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTcaca                                                                               >  2:768352/1‑139 (MQ=255)
                                          agaACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACAcc                                                                                                                     >  2:379534/1‑100 (MQ=255)
                                          agaACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACAcc                                                                                                                     <  1:379534/100‑1 (MQ=255)
                                          agaACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACaa                                                                              >  2:452753/1‑139 (MQ=255)
                                                    tGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTcaca                                                                                                                       <  1:611960/88‑1 (MQ=255)
                                                    tGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTcaca                                                                                                                       >  2:611960/1‑88 (MQ=255)
                                                    tGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAg                                                                                     >  1:612277/1‑122 (MQ=255)
                                                    tGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAg                                                                                     <  2:612277/122‑1 (MQ=255)
                                                      gCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTc                                                                  >  1:216748/1‑139 (MQ=255)
                                                                             aaGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGc                                           >  1:613278/1‑139 (MQ=255)
                                                                             aaGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGc                                           >  2:7016/1‑139 (MQ=255)
                                                                              aGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCg                                          <  2:407266/139‑1 (MQ=255)
                                                                                gTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTcac                                                                                                                        <  2:461517/59‑1 (MQ=255)
                                                                                gTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTcac                                                                                                                        >  1:461517/1‑59 (MQ=255)
                                                                                          aaaaGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAAccc                              <  2:234502/139‑1 (MQ=255)
                                                                                                 ttCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCt                       >  2:812005/1‑139 (MQ=255)
                                                                                                              gcgcTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGt                                                                                        >  2:303837/1‑61 (MQ=255)
                                                                                                              gcgcTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGt                                                                                        <  1:303837/61‑1 (MQ=255)
                                                                                                                gcTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACa                                                                                                            <  2:399587/39‑1 (MQ=255)
                                                                                                                gcTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACa                                                                                                            >  1:399587/1‑39 (MQ=255)
                                                                                                                      gcgcCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTg  <  2:216748/139‑1 (MQ=255)
                                                                                                                            |                                                                                                                                     
AAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTG  >  NC_000913/1369912‑1370169

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑