Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 693,599 T→C T453A (ACA→GCA)  miaB ← tRNA‑i(6)A37 methylthiotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913693,5990TC94.4% 116.7 / ‑1.2 36T453A (ACA→GCA) miaBtRNA‑i(6)A37 methylthiotransferase
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base C (17/17);  total (17/19)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.11e-01

AAAAAGGCGCAAGCATTGCAATAAAGATGGGGATAAAGAGAGAAAAAACAAGGCCCACCGGAACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCGGCGTGCCTTCGAAGTTGACCACGCGGTTATTTTCCG  >  NC_000913/693453‑693745
                                                                                                                                                  |                                                                                                                                                  
aaaaaGGCGCAAGCATTGCAATAAAGATGGGGATAAAGAGAGAAAAAACAAGGCCCACCGGAACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGGCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGCTTCt                                                                                                                                                >  1:592371/1‑151 (MQ=255)
   aaGGCGCAAGCATTGCAATAAAGATGGGGATAAAGAGAGAAAAAACAAGGCCCACCGGAACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGCTTCTGcc                                                                                                                                             >  2:359253/1‑151 (MQ=255)
          aaGCATTGCAATAAAGATGGGGATAAAGAGAGAAAAAACAAGGCCCACCGGAACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGCTTCTGCCACGCGCa                                                                                                                                      <  1:533419/151‑1 (MQ=255)
           aGCATTGCAATAAAGATGGGGATAAAGAGAGAAAAAACAAGGCCCACCGGAACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGCTTCTGCCACGCGCAg                                                                                                                                     <  2:331709/151‑1 (MQ=255)
              aTTGCAATAAAGATGGGGATAAAGAGAGAAAAAACAAGGCCCACCGGAACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGCTTCTGCCACGCGCAGAcc                                                                                                                                  >  2:363941/1‑151 (MQ=255)
                tGCAATAAAGATGGGGATAAAGAGAGAAAAAACAAGGCCCACCGGAACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGCTTCTGCCACGCGCAGAccc                                                                                                                                 >  1:65086/1‑150 (MQ=255)
                      aaaGATGGGGATAAAGAGAGAAAAAACAAGGCCCACCGGAACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGCTTCTGCCACGCGCAGACCCATTTCa                                                                                                                           <  2:40819/150‑1 (MQ=255)
                      aaaGATGGGGATAAAGAGAGAAAAAACAAGGCCCACCGGAACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGCTTCTGCCACGCGCAGACCCATTTCAt                                                                                                                          <  1:411284/151‑1 (MQ=255)
                           tGGGGATAAAGAGAGAAAAAACAAGGCCCACCGGAACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGCTTCTGCCACGCGCAGACCCATTTCATCTTc                                                                                                                      >  1:464344/1‑150 (MQ=255)
                               gATAAAGAGAGAAAAAACAAGGCCCACCGGAACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGCTTCTGCCACGCGCAGACCCATTTCATCTTCAGTa                                                                                                                  >  2:447643/1‑150 (MQ=255)
                               gATAAAGAGAGAAAAAACAAGGCCCACCGGAACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGCTTCTGCCACGCGCAGACCCATTTCATCTTCAGTAc                                                                                                                 <  1:447638/151‑1 (MQ=255)
                                       gagaAAAAACAAGGCCCACCGGAACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGCTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACt                                                                                                         >  2:67908/1‑151 (MQ=255)
                                          aaaaaaCAAGGCCCACCGGAACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGCTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTa                                                                                                      <  2:485076/151‑1 (MQ=255)
                                           aaaaaCAAGGCCCACCGGAACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTAc                                                                                                     <  1:205638/151‑1 (MQ=255)
                                             aaaCAAGGCCCACCGGAACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGCTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACcg                                                                                                   <  1:495431/151‑1 (MQ=255)
                                              aaCAAGGCCCACCGGAACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGCTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACcgc                                                                                                  >  2:613185/1‑151 (MQ=255)
                                                     cccACCGGAACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGCTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCagaga                                                                                           <  2:436037/151‑1 (MQ=255)
                                                     cccACCGGAACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGCTTCTGCCACGCGCAGACCCATGTCATCTTCAGTACGCACCACTTTACCGCGCagaga                                                                                           <  1:512577/151‑1 (MQ=255)
                                                           ggAACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGCTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCgg                                                                                     >  2:385409/1‑151 (MQ=255)
                                                              aCGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGCTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTa                                                                                  <  1:169853/151‑1 (MQ=255)
                                                                   aGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGCTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGt                                                                             <  1:145009/151‑1 (MQ=255)
                                                                    ggCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTc                                                                            <  2:609315/151‑1 (MQ=255)
                                                                      ccTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGCTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCgg                                                                          >  1:622737/1‑151 (MQ=255)
                                                                         gagaATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGCTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCGCTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGt                                                                         <  2:434971/149‑1 (MQ=255)
                                                                          agaATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGCTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAAt                                                                      <  2:12130/151‑1 (MQ=255)
                                                                                    cTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGCTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATcg                                                            >  1:600840/1‑151 (MQ=255)
                                                                                      gATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGCTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCATCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATcg                                                            <  1:74619/149‑1 (MQ=255)
                                                                                                 aCGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGCTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCg                                                >  1:260900/1‑150 (MQ=255)
                                                                                                      aaGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGCTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCat                                           >  2:444979/1‑150 (MQ=255)
                                                                                                       aGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGCTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCatat                                         <  2:332147/151‑1 (MQ=255)
                                                                                                          tCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGCTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCg                                      >  2:661279/1‑151 (MQ=255)
                                                                                                            gTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGCTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCGGc                                    <  1:599990/151‑1 (MQ=255)
                                                                                                                 cTTTGCGGGTACGGGCAATCACTGATTCCGGTGCTTCTGCCCCCCGCAGCCCAATTTCACCTTCAGAACGAACACCTTTACCGCGCAGAGGTTTCGGGTAGCCGCCGGTATTTTCTCCATCAACGATTTTACCGACCATATccgccggcc                                >  1:212468/1‑147 (MQ=255)
                                                                                                                         gTACGGGCAATCACTGATTCCGGTGCTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCGGCGTGCCTTCGAAGt                       <  1:291175/151‑1 (MQ=255)
                                                                                                                            cGGGCAATCACTGATTCCGGTGCTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCGGCGTGCCTTCGAAGTTg                     >  1:621638/1‑150 (MQ=255)
                                                                                                                                              ggTGCTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCGGCGTGCCTTCGAAGTTGACCACGCGGTTATTTTCCg  >  2:495167/1‑151 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
AAAAAGGCGCAAGCATTGCAATAAAGATGGGGATAAAGAGAGAAAAAACAAGGCCCACCGGAACGGCAGGCCTGAGAATTACGGCTGATAATAACCCACGCCAAGGTCGTTTTCTTTGCGGGTACGGGCAATCACTGATTCCGGTGTTTCTGCCACGCGCAGACCCATTTCATCTTCAGTACGCACCACTTTACCGCGCAGAGAGTTCGGGTAGACGTCGGTAATTTCTACATCGACGAATTTACCGATCATATCCGGCGTGCCTTCGAAGTTGACCACGCGGTTATTTTCCG  >  NC_000913/693453‑693745

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