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.4% 63.2 / NA 28G211G (GGC→GGT) mnmAtRNA‑specific 2‑thiouridylase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (10/17);  total (11/17)
Fisher's exact test for biased strand distribution p-value = 3.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTGCAACTGCTGGGCAATCAACCCGACAGACATCAGCCGCGGGTGTTCATGGCCCTGAGCGACAACCAGAATGTTGTTTTCGACGTCTTTGTCCACCACATACCACGGTTCTTCGGTACCTTCTTTGGTGCCAC‑CGATACCCAGACCTTTACGCTGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGGCCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAG‑ACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTGCCATCGACATCGGCCCGACGTACG  >  NC_000913/1192894‑1193387
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       
aGTGCAACTGCTGGGCGATCAACCCGACAGACATCAGCCGCGGGTGTTCATGGCCCTGAGCGACAACCAGAATGTTGTTTTAGACGTATTTGTCCACCACATACCACGGTTCTTCGGTACCTTCTTTGGTGCCAC‑CGATACCCAGACATTTAAGCTGACCGCGAGTGTGATACATCAGCCCCTGGGGCTCACCAATTTCATAGCCATCGACGGAAATGATTTTGCCCGGTTGCGCCGGGAGATAACGAccc                                                                                                                                                                                                                                                      >  1:312694/1‑251 (MQ=255)
     aaCTGCTGGGCAATCAACCCGACAGACATCAGCCGCGGGTGTTCATGGCCCTGAGCGACAACCAGAATGTTGTTTTCGACGTCTTTGTCCACCACATACCACGGTTCTTCGGTACCTTCTTTGGTGCCAC‑CGATACCCAGACCTTTACGCTGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGaa                                                                                                                                                                                                                                                 <  1:124501/251‑1 (MQ=255)
     aaCTGCTGGGCAATCAACCCGACAGACATCAGCCGCGGGTGTTCATGGCCCTGAGCGACAACCAGAATGTTGTTTTCGACGTCTTTGTCCACCACATACCACGGTTCTTCGGTACCTTCTTTGGTGCCAC‑CGATACCCAGACCTTTACGCTGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGaa                                                                                                                                                                                                                                                 <  1:50732/251‑1 (MQ=255)
         gctgGGCAATCAACCCGACAGACATCAGCCGCGGGTGTTCATGGCCCTGAGCGACAACCAGAATGTTGTTTTCGACGTCTTTGTCCACCACATACCACGGTTCTTCGGTACCTTCTTTGGTGCCAC‑CGATACCCAGACCTTTACGCTGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTc                                                                                                                                                                                                                                              <  1:146866/250‑1 (MQ=255)
                       ccGACAGTCATCAGCCGCGGGTGTTCATGGCCCTGAGCGACAACCAGAATGTTGTTTTCGACGTCATTTTCCACCACATACCACGGTTCTTCGGTACCTTCTGTGGTGCCAC‑CGATACCCAGACCCTTACGCGGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTGGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTc                                                                                                                                                                                                                               <  1:379375/251‑1 (MQ=255)
                          acagacagCTGCCGCGGGTGTTCATTGCCCTGAGCGACAGCCAGAATGTTGTTTTCGACGTCTTTGTCCACCACATACCACGTTTCGTCGGTTCCTTCTTTGGTGCCAC‑CGATACCCAGACCTTTACGCTGACCGAGAGTTTGATACATCAGCCCCTGGTGCTCACCAATTTCATCGCCATCGTCGGTAATGATTTTGCTCGGTTGCGTCGGGAGATAACGACCCAGTAACTCGCGGAATTTACGTTCGcc                                                                                                                                                                                                                            <  2:97606/251‑1 (MQ=255)
                           cagacaTCAGCCGCGGGTGTTCATGGCCCTGAGCGACAACCAGAATGTTGTTTTCGACGTCTTTGTCCACCACATACCACGGTTCTTCGGTACCTTCTTTGGTGCCAC‑CGATACCCAGACCTTTACGCTGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCa                                                                                                                                                                                                                           <  1:56263/251‑1 (MQ=255)
                                   aGCCGCGGGTGTTCAGGGCCCTGAGCGACAACCAGAATGTTGTTTTCGACGTCTTTGTCCACCACATACCACGGTGCTTCGGTACCTTCTGTGGTGCCAC‑CGATACCCAGACCTTTACGCTGACCGAGAGTTTGATACATCAGCCCCTGGTGCTCACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCa                                                                                                                                                                                                                   <  1:196116/251‑1 (MQ=255)
                                   aGCCGCGGGTGCTAATGGCCCTGAGCGACAACCAGAATGTTGTTTTCGACGTCTTTGTCCACCACATACCACGGTTCTTCGGTACCTTCTTTGGTGCCAC‑CGATACCCAGACCTTTACGCTGACCGAGAGTGTGCTACATCAGCCCCTGGTGCTCACCATTTTCATAGTCATCGACGGTTATGATTTTGCCCGCTTGCCCCGGGAGATAACGGCCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGc                                                                                                                                                                                                                    >  2:160865/1‑250 (MQ=255)
                                                    gCCCTGAGCGACAACCAGAATGTTGTTTTCGACGTCTTTGTCCACCACATACCACGGTTCTTCGGTACCTTCTTTGGTGCCAC‑CGATACCCAGACCTTTACGCTGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCtt                                                                                                                                                                                                  >  1:95709/1‑251 (MQ=255)
                                                             gACAACCAGAATGTTGTTTTCGACGTCTTTGTCCACCACATACCACGGTTCTTCGGTACCTTCTTTGGTGCCAC‑CGATACCCAGACCTTTACGCTGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTcg                                                                                                                                                                                         >  1:34733/1‑251 (MQ=255)
                                                                    aGAATGTTGTTTTCGACGTCTTTGTCCACCACATACCACGGTTCTTCGGTACCTTCTTTGGTGCCAC‑CGATACCCAGACCTTTACGCTGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTgac                                                                                                                                                                                  >  1:145178/1‑251 (MQ=255)
                                                                                    cGTCTTTGTCCACCACCTACCACGGTTCTTCGGTCCCTTCTTTGGTGACAC‑CGATACCCCGCCCTTAACGCTGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCACTTTCATCGCCAACGACGGTAATGATTTTGCCGGGTTGCGCCGGGAGATACCGACCCAGGAACCAGCGGAATTTCAGTTCGCCAATAAAGCAAATGCCGGAAGAGTCTTTTTTTTTCGCGGTAACCAG‑ACAAAGATATTcc                                                                                                                                                                 >  2:363258/1‑250 (MQ=255)
                                                                                     gtcgttGTCCACCACATACCACGATTCTTAGGTACCTTCTTTGGATACAC‑CGATACCCAGACCTTTACGATGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCCCCACTTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGAGCCGGGCGATACCGACCCAGGAGCTAGTGAAATTTACGTTCACCACTGGAGCACCTCCAGGCAGAGTTTTTGTTCTACGCGGTGACCAGGACCaattattagg                                                                                                                                                                 >  2:5115/5‑243 (MQ=255)
                                                                                          tGTCCACCTCATACCACGGTTCTTCGGTACCTCCTTTGGTGCCAC‑CGATACCCAGCCCTTTACGCTTACCGAGAGTGTGATACATCATCCCCGGGTGCTCACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGTTAGAGTCTTTTTTCTTCGCGGTGACCAG‑ACCAAGATCTTCAGCAAt                                                                                                                                                            <  1:63074/250‑1 (MQ=255)
                                                                                            tcacccataTACCACGGGTCTTCGGTACCTTCTTTAGTGCCAC‑CGATACCCAGAACTTTACGCTGACCGAGAGTGTTATACATCAGCACCTGATGCTCACAAATTTCATCGCCATCGACGGTAATGATTTTAACCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGATCACCCA‑CCCAAGATCTTCAGCAATCtt                                                                                                                                                         >  2:188175/9‑251 (MQ=255)
                                                                                              caccacATACCTCTGTTCTTCGGTACCTTCTTTGGTGC‑ACTCTCTACCCAGACCTTGCCTCAGACCGTGAGTTCGATACATCAGCCCCTGGTGCTCACCAATTTCATCGCCATCGACGGGAATGATTTTGCCCGGTTGCGCCGGTAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAATTGCCGGTAGAGTTTTTTTTCTTCGCGGTGACCAT‑TCCAAGATCTTTATCAATCTTAc                                                                                                                                                       <  2:260279/251‑1 (MQ=255)
                                                                                                        cACGGTTCTTCGGTACATTCTTTGGTGCCGC‑CGATACCCAGACCTTTACGCTGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAATTTCATCGCCATCGACGGTACTGATTTTGCCCGGTTGCGCCGGGCGATAACGACCCAGGACCTCGCGGATTTTACGTTCGCCAGTGAAGCAAATTCCTGTAGAGTCTTTTTTCTTCGCGGTGACCCG‑GCCAAGATTGTCAGCAATCTTTCGCACCTTCgg                                                                                                                                             >  2:145162/1‑251 (MQ=255)
                                                                                                                                          aTACCCAGACCTTTAAGCTGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAG‑ACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTg                                                                                                            >  1:9726/1‑251 (MQ=255)
                                                                                                                                                aGACCTTTACGCTTACCGTGATTGTGATTCATCAGCCCCTGGTGCTCACTAATTTCATCGCCCTCGACGGTAAGGATTTTGCCCGGTGGCGTCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCGTCTCGGTGACCAG‑ACCAAGATCTTCAGCAATCTTTCGCTCCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTTCGCAAt                                                                                                      <  2:249668/251‑1 (MQ=255)
                                                                                                                                                                   gagTGTGATACATCAGCCCCTGGTGCTCACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAG‑ACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGta                                                                                   <  1:268149/251‑1 (MQ=255)
                                                                                                                                                                          agacagcaacccctcttGCTCACCACATTCATCGCCATCGATGGTAAAGACATTGCCCGTTGGCGCCGGGAGATAATGACCCCGGAACTCTCTGAATTTACGTTCGCCAATGAATCTAATGCCGGTATAGTCTTTTTTCTTCGCGGTGACCAG‑ACCAAGATCTTCAGCAATCTTTCTCACCTTCTGTTTTTCCAGTACTAAGACAGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGa                                                                            <  2:233279/235‑1 (MQ=255)
                                                                                                                                                                                        ggTGCTCACCAATTTCATCGCCAGCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAG‑ACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCt                                                              <  1:52250/251‑1 (MQ=255)
                                                                                                                                                                                                             cATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAG‑ACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCAc                                           <  1:64898/249‑1 (MQ=255)
                                                                                                                                                                                                             cATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAG‑ACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCAc                                           >  2:64898/1‑249 (MQ=255)
                                                                                                                                                                                                                tggcGGTAATGGTTTTTCCCGGTTGCGCCGTGAGATTATTACCCAGAGACTCGCGGAATTTTCGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGATCAG‑ACAAAGATCTTCAGCTATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGTCTGAGCGTATAAAGGAAGTAGCTCTGGTCATTATTGCGGTctattccactcag                                       <  2:282623/247‑13 (MQ=255)
                                                                                                                                                                                                                                  ccGGTTGCGACGGGAGATATCGACCTAGGAACTCTCTGAAGTCACGTTCGCCAATGAATAAAATGCCGGTAGAGTCTTTCTTATTCTCGGTGACCAG‑ACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTGTTATAGTTCGCCGACCGGGAACAGTCTTTGCGCAATCTGCGCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCTCGCATCAGGCGGCTCTTGCCATCg                    <  2:202269/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    aaCGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAG‑ACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTGCCATCGACATCGGCCCGacgtacg  <  1:268948/251‑1 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       
AGTGCAACTGCTGGGCAATCAACCCGACAGACATCAGCCGCGGGTGTTCATGGCCCTGAGCGACAACCAGAATGTTGTTTTCGACGTCTTTGTCCACCACATACCACGGTTCTTCGGTACCTTCTTTGGTGCCAC‑CGATACCCAGACCTTTACGCTGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGGCCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAG‑ACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTGCCATCGACATCGGCCCGACGTACG  >  NC_000913/1192894‑1193387

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