Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,193,141 G→A G211G (GGC→GGT)  mnmA ← tRNA‑specific 2‑thiouridylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,193,1410GA96.9% 87.5 / NA 32G211G (GGC→GGT) mnmAtRNA‑specific 2‑thiouridylase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (16/15);  minor base C (1/0);  total (17/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACTGCTGGGCAATCAACCCGACAGACATCAGCCGCGGGTGTTCATGGCCCTGAGCGACAACCAGAATGTTGTTT‑TCGACGTCTTTGTCCACCACATACCACGGTTCTTCGGTACCTTCTTTGGTGCCACCGATACCCAGACCTTTACGCTGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAA‑TTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGGCCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTGCCATCGACATCGGCCCGACGTACG  >  NC_000913/1192899‑1193387
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                      
aaCTGCTGGGCAATCAACCCGACAGTCATCAGCCGCGGGTGTTCAGGGCCCTGAGCGACAATCAGAATGT‑GTTTGTCGACGTCTTTGTCCACCACATACCACGGTTCTTCGGTACCGTCTTGGGTGCCACCGATACCCAGACCTTTACGCTGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAA‑TTTCATCGCCATCGACGGGAATGATTTGGCCCGGTTGCGCCGGGAGATAACggtgaagtgt                                                                                                                                                                                                                                                <  2:212146/251‑9 (MQ=255)
aaCTGCTGGGCAATCAACCCGACAGACATCAGCCGCGGGTGTTCATGGCCCTGAGCGACAACCAGAATGTTGTTT‑TCGACGTCTTTGTCCACCACATACCACGGTTCTTCGGTACCTTCTTAGGTGCCACCGATACCCAGACCTTTACGCTGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAA‑TTTCCTCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGAAAACGCCCCAGGa                                                                                                                                                                                                                                                 >  2:300618/1‑250 (MQ=255)
     gtgGGCAATCAACCCGACAGACATCAGCCGCGGGTGTTCATGGCCCTGAGCGACAACCAGAATGTTGTTT‑TCGACGTCTTTGTCCACCACATACCACGGTTCTTCGGTACCTTCTTTGGTGCCACCGATACCCAGACCTTTACGCTGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAA‑TTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTcgc                                                                                                                                                                                                                                           <  1:463741/250‑1 (MQ=255)
       gggCAATCAACCCGACAGACATCAGCCGCGGGTGTTCATGGCCCTGAGCGACAACCAGAATGTTGTTT‑TCGACGTCTTTGTCCACCACATACCACGGTTCTTCGGTACCTTCTTTGGTGCCACCGATACCCAGACCTTTACGCTGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAA‑TTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCgg                                                                                                                                                                                                                                         <  2:325349/251‑1 (MQ=255)
                       agacaTCAGCCGCGGGTGTTCATGGCCCTGAGCGACAACCAGAATGTTGTTT‑TCGACGTCTTTGTCCACCACATACCACGGTTCTTCGGTACCTTCTTTGGTGCCACCGATACCCAGACCTTTACGCTGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAA‑TTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCaa                                                                                                                                                                                                                         <  1:475660/251‑1 (MQ=255)
                             cAGCCGCGGGTGTTCATGGCCCTGAGCGACAACCAGAATGTTGTTT‑TCGACGTCTTTGTCCACCACATACCACGGTTCTTCGGTACCTTCTTTGGTGCCACCGATACCCAGACCTTTACGCTGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAA‑TTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAg                                                                                                                                                                                                                    <  1:65802/250‑1 (MQ=255)
                                ccGCGGGTGTTCATGGCCCTGAGCGACAACCAGAATGTTGTTT‑TCGACGTCTTTGTCCACCACATACCACGGTTCTTCGGTACCTTCTTTGGTGCCACCGATACCCAGACCTTTACGCTGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAA‑TTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCaaa                                                                                                                                                                                                                <  2:118469/251‑1 (MQ=255)
                                            agggccctgctggACAACCAGAAGGTTTTTT‑TCGCCGTCTTTTCCCACCACATACCAAGGTGCTTCGGTACATTCGTTGGGGCCACCGATACCCAGAACGTTAGGCTGACCGAGAGTGTGAGACATCAGCCCCTGGTGCT‑ACCAAATGTCATCGCCATCGACGGTAATGATTTTGCCCGGTAGCGCCGGGAGATAACGACACAGGAAATCGCTGAATTTACTTTCGCCAATGAAGCAAATGCCGGTagcgc                                                                                                                                                                                                    <  2:123999/239‑4 (MQ=255)
                                                 ccTGAGCGACAACCAGAATGTTGTTT‑TCGACGTCTTTGTACACCACATACCACGGTTCTTCGGTACTTTCTTTGGTGCCACCGATACCCAGATCTTTACGTTGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAA‑TTTCATTGCCATCGACGGTAATGATTTTGACCGGTTGCGCCGGGAGATCACGACCCAGGAACTCGCGGTATTTACGTTCGCCATTGAAGAAAATACCGGTAGAGTAtttt                                                                                                                                                                                               >  2:444025/1‑251 (MQ=255)
                                                        gACAACCAGAATGTTGTTT‑TCGACGTCTTTGTCCACCACATACCACGGTTCTTCGATACCTTCTTTGGTGCCACCGATACCCAGACCTTTACGCTGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAA‑TTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGACTTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTcg                                                                                                                                                                                        >  2:504861/1‑251 (MQ=255)
                                                                gAATGTTGTTT‑TCGACGTCTTTGTCCACCACATACCACGGTTCTTCGGTACCTTCTTTGGTGCCACCGATACCCAGACCTTTACGCTGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAA‑TTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTTGAGTCTTTTTTCTTCGCGGTgacc                                                                                                                                                                                >  2:166330/1‑251 (MQ=255)
                                                                gAATGTTGTTT‑TCGACGTCTTTGTCCACCACATACCACGGTTCTTCGGTACCTTCTTTGGTGCCACCGATACCCAGACCTTTACGCTGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAA‑TTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTgacc                                                                                                                                                                                >  1:230706/1‑251 (MQ=255)
                                                                            tCGACGTCTTTGGCCACCACATACCACGGTTCTTCGGTACCTTCTTTGGTGCCACCGATACCCAGACCTTTACGATGACCGAGAGTGTGATACATCAGCCACTGGTGCTCACCAC‑TTTCATCGCCAACGACGGTAATGATTTTGCCCGGGTGCGCCGGGAGCTAACGACCCAGGACATAGCGGATTTTACGTTCGCCAATGAAGCAAAAGCAGGAAGCGTCTTTTTTCTTCGCGTTAGCCAGACCAAGacc                                                                                                                                                                     >  2:303507/1‑249 (MQ=255)
                                                                                              acaTACCACGGTTCTTCGGTACCTTCTTTGGTGCCACCGATACCCAGACCTTTACGCTGACCGAGAGTGTGATACATCAGCCCCTGGTGATCACCAA‑TTTCATCGCCATCGACGGTACTGCTTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGCACTCGCGGAATTTACGTTCGCCAATGAAGCAAATACCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCCATCTTACGCAc                                                                                                                                                   >  1:254448/1‑251 (MQ=255)
                                                                                                       ggTTCTTCGGGACCTTCTTTGGTGCCACCGATACCCAGACCTTTACGCTGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAA‑TTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGCCCAGAACAAGGTCTTCAGCAATCTTACGCCCATGCgggtt                                                                                                                                          >  2:316326/1‑248 (MQ=255)
                                                                                                             tcGGTACCTTCTGTGGTGCCACCGATACCCAGACCTTTACGCTGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAA‑TTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCa                                                                                                                                     <  1:241243/250‑1 (MQ=255)
                                                                                                                      tcttTGGTGCCACCGATACCCAGACCTTTACGCTGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAA‑TTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGAc                                                                                                                           >  1:505240/1‑251 (MQ=255)
                                                                                                                                       aCCCAGACCTGTACGCTGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAA‑TTTCATCGCCATCGCCGGTAATGATTTTGCCCGGTTGCGCCGGGTGATAACGACCCAGGAACTCGCGGAATTTCGGTTGGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGTCGAACAGACCAAGCTCTTCAGCACTATTCAGCACCTGCGGGTTTTCCAGTTCGCCGACCGGGAACAGTCTTTg                                                                                                            >  2:76002/1‑249 (MQ=255)
                                                                                                                                                 ttACGCTGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAA‑TTTCATCGCCATCGACGGGAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCGTACGCACCGGCGGTTGTTCCAGTGCGCCGACCGGGAACAGGCTTTGCGCAATCTGc                                                                                                  <  2:435566/249‑1 (MQ=255)
                                                                                                                                                                            aGCCCCTGGTGCTCACCAA‑TTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGct                                                                     <  1:490211/251‑1 (MQ=255)
                                                                                                                                                                               cccTGGTGCTCACCAA‑TTTCATCGCCATCGCCGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGACGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTg                                                                  >  1:317482/1‑251 (MQ=255)
                                                                                                                                                                                   ggTGCTCACCAA‑TTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGt                                                                <  1:504861/249‑1 (MQ=255)
                                                                                                                                                                                           ccAA‑TTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCt                                                      >  1:386896/1‑251 (MQ=255)
                                                                                                                                                                                                           tCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGTATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCCGTTCGCCAACCGGGAACAGGCTTTGCGCAATCTGATCATGGCTGAGCGTATACAGGAAGTAGCTCTGGTCTTTTTTGCTGTCCAGACCACgc                                         >  2:442325/1‑249 (MQ=255)
                                                                                                                                                                                                                ggTAATGATTGTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAggc                                  <  1:481204/251‑1 (MQ=255)
                                                                                                                                                                                                                 gTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGCAGTAGCTCTGGTCTTTATTGCTGTCCAGCCCAc                                           >  2:220232/1‑241 (MQ=255)
                                                                                                                                                                                                                 gTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCAc                                           <  1:220232/241‑1 (MQ=255)
                                                                                                                                                                                                                  tAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAggcgg                                >  1:37174/1‑251 (MQ=255)
                                                                                                                                                                                                                                 gTTGCGCCGGGATATAACGACCCAGGAACACGCGGAACTTATGTTCGACAATGAAGCAAATGCCGGTAGATTCTTTTTTCTTCGCGTTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCATGTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCACATGGCTGAGCTTATAAAGGAAGTAGCTCTGGTCTTTATGGCTGTCCAGCCCACGCAGCAGGCGGCTCTTGCCATCGACa                 <  2:172076/251‑1 (MQ=255)
                                                                                                                                                                                                                                      gcCGGGAGATAACGACCCAGGAACTCGCTGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGACATAGCTCTGGTCTTTATTGCTGTCCAGACCCCGCAGCAGGCGGCTCTTGCCATCGCCATCGGc            >  1:241348/1‑251 (MQ=255)
                                                                                                                                                                                                                                           gagaTAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTGCCATCGACATCGGCCCGac       <  2:386896/251‑1 (MQ=255)
                                                                                                                                                                                                                                            agaTAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTGCCATCGACATCGGCCCGacg      <  1:383744/251‑1 (MQ=255)
                                                                                                                                                                                                                                                aaCGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGGCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTGCCATCGACATCGGCCCGacgtacg  >  2:385501/1‑251 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                      
AACTGCTGGGCAATCAACCCGACAGACATCAGCCGCGGGTGTTCATGGCCCTGAGCGACAACCAGAATGTTGTTT‑TCGACGTCTTTGTCCACCACATACCACGGTTCTTCGGTACCTTCTTTGGTGCCACCGATACCCAGACCTTTACGCTGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAA‑TTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGGCCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTGCCATCGACATCGGCCCGACGTACG  >  NC_000913/1192899‑1193387

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