Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,913,761 G→A L414L (CTG→TTG)  rlmD ← 23S rRNA m(5)U1939 methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,913,7610GA92.9% 80.1 / ‑5.3 28L414L (CTG→TTG) rlmD23S rRNA m(5)U1939 methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (10/16);  total (11/17)
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

GTTGCAGACAATACGCCCAGGTTTCGGCTAAGCACTCACACGACTTCTGGCTGGTAATACCCAGACTTGCGATCCATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACGCTTCGCTATCCCGAGCCAGCGTTGCAGGGTTACAGGATACATAAACTATACGAATAGGTTCCAGTTTTATAATTTGCTGCATAACACCTGCGGCACCTGCTCGCGCCGGGTCCAGCAACA  >  NC_000913/2913519‑2913922
                                                                                                                                                                                                                                                   |                                                                                                                                                                 
gTTGCAGACAATACGCCCAGGTTTCGGCTAAGCACTCACACGACTTCTGGCTGGTAATACCCAGACTTGCGATCCATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCTCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAACATCGcc                                                                                                                                                            <  1:165823/250‑1 (MQ=255)
      gACAATACGCCCGGGTTCCGGCTAAGCACCCCCACGACTTCTGGCTGGTAATACCCAGACTTGCGATCCATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAACATCGCCAGTcgcg                                                                                                                                                     <  2:252131/251‑1 (MQ=255)
      gACAATACGCCCAGGTTTCGGCTAAGCACTCACACGACTTCTGGCTGGTAATACCCAGACTTGCGATCCATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAACATCGCCAGTcgcg                                                                                                                                                     >  1:6553/1‑251 (MQ=255)
                       tCGGCTAAGCACTCACACGACTTCTGGTTGGGAATACCCAGACTTGGGATCCATTTTCCCGGATCAAATTCACCAGCCTTAGTGATATGGGCTCTTCTTACCGCAACCATCGTCCTCTCCTGTAGGAACCAGCCCAGCCGAAATCGGCAAATCGCATCTATTTAACGCTCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTTTGTGAACATAGCTAACATCGCCAGTCGCGCAATGGTATATCCTGc                                                                                                                                     <  1:296152/250‑1 (MQ=255)
                                      cacGACTTCTGGATGGTAATACCCAGACTTTCTATCCATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAAGCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAACAGCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACGCTTCg                                                                                                                     <  1:39930/251‑1 (MQ=255)
                                      cacGACTTCGGGTTGGTAATACCAAGACTTGCGATCAATTTTTCCGGAACAATTTCACCAGCCTTATTGATAGGTGCACTTCTTCCCCCAACCATCGTCCCCCCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATGCCAGATGTCCCGTGTGTGGGAACATATCCAACAGCGCCAGTCGCGCAAGGGTATATCCTGCTTTTAATAACGCTTCg                                                                                                                     <  2:129782/251‑1 (MQ=255)
                                      cacGACTTATGGCTGGTAATACCCAGACGTGCGATCCATTTTTCCGGAGCAAATCCACCAGCCTTATGGATAGTGGCGCTTCTTACCGCAACCATCTTCCTCTCCTTTAGGAACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGAGTCCAGACGTCCCGTGTGTGGGAACATAGCCAACATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACTCTTCg                                                                                                                     <  2:281316/251‑1 (MQ=255)
                                              gtggctggAAATACACAGACTTGTGAGCCATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGAACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGAACCATCGATTCCAGATGTCCCGTGAGTGGGAACATATCAAACATCCCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACGCTTCGCTATCCCg                                                                                                             >  2:217704/2‑251 (MQ=255)
                                                         accccatccTTGCGATCCATTTTTCGGTACCAAATTCGCAAGCTTTATTGCAATGTGCACTTCTTACCGCACCTATCTTCCCCCCCTGTAGGGACCTGCCCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGTGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAACATCGCCAGGCGCGCAATGGTATATCCTGCTTTTAATACCGATTCGTTATCCCGAGCCAGCGTTg                                                                                                  <  2:302188/243‑1 (MQ=255)
                                                          aCCCAGACTTGCGATACATTTTTACGGAGCAAATTCAACAGCCTTATCGCTATGTGCACTTCTTACCGCAACCATCGTCCTCTACTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCAAACATCCACATTCGCGCAATGGAATATCCTGCTTTTAATAACGCTTCGCTATcccactccagcgtttc                                                                                                 >  2:175105/1‑238 (MQ=255)
                                                                actggCGATCCATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGAAACCAGCGGCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAACATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACGCTTCGCTATCCCGAGCCAGCGTTGCAGGGtt                                                                                           <  1:40989/247‑1 (MQ=255)
                                                                      gctCCATTCTGCCGGATCAAATCCACCAGCCTTTGTGACAGGTCCACTTCTTACCGCACCCATCGTCCTCACCTTTAGGGACCAGGCTTGCCGAAATCGGCTAAGCGCAACGATTTAACGCGCGAGAAAAGTTCCATCGATTCCATATGTCCCGTGTGTTGGAACATAGCCAGCAACGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACGCTTCGCGATTCCGAGCTAACGGTTCAGGGTTTCAGGa                                                                                     <  2:226811/249‑1 (MQ=255)
                                                                                tCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCCACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAACATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACGCTTCGCTATCCCGAGCCAGCGTTGCAGGGTTACAGGATACATAAACt                                                                           >  2:316243/1‑251 (MQ=255)
                                                                                   ggATCAAATTCCCCAGCCTTATTTATATGTGCACTTCTTCCCGCAACCATCGTCCTCCCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAACCTCGCCAGTCGCGCAATGGTATAACCTGCTTTTAATAACGCTTCGCTATCCCGAGCCAGCGTTGCAGGGTTACAGGATACATAAACtata                                                                        <  2:270959/251‑1 (MQ=255)
                                                                                               ccAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAACATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACGCTTCGCTATCCCGAGCCAGCGTTGCAGGGTTACAGGATACATAAACTATACGAATAGGTTcc                                                            >  2:76072/1‑251 (MQ=255)
                                                                                                   ccTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAACATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACGCTTCGCTATCCCGAGCCAGc                                                                                                      >  1:10372/1‑205 (MQ=255)
                                                                                                   ccTTATTGATAGGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATAGCCAACATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACGCTTCGCTATCCCGAGCCAGc                                                                                                      <  2:10372/205‑1 (MQ=255)
                                                                                                         tGATATGTGCACTTCTTACCGCAACAATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATAGGCAAATCGCAACTATATAACGCGCAAGCAAAGTACCATCGATTACAGATGTCCAGTGTGTGGGAACATATCCAACATCGCCAGTCGCGCACTGGTATATCCTGCTTTTAATAACGCTTCGCTATCCCGAGCCAGCGTTGCAGGGTTACAGGATACATAATCTACACGAATAGGTTCCAGGTTtata                                                   >  1:85428/1‑250 (MQ=255)
                                                                                                          gATATGTGCACTTCTTACCGCAACCATCGTCCTCTACTTTAGGGACCAGACATGCCGAACTCGGCCAATCGCCACTATTTAACGCGCGCGAAAAGTCCCATAGACTCCCGATGTCCCGTGTGTGGGAACATATCCAACATCGCCAGTCGCGCAATGGAATATCCTGCTTTTAATAACGCTTCGCTATCCCGCGCCAGCGTTGCAGGGTTACAGGATACATAAACTATACGAATAGGTTCCAGTTTTATAAt                                                 >  1:24053/1‑251 (MQ=255)
                                                                                                           atatGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACGCTTCGCTATCCCGAGCCAGCGTTGCAGGGTTACAGTATACATAAACTATACGAATAGGTTCCAGTTTTATAAtt                                                >  2:101619/1‑251 (MQ=255)
                                                                                                                tgCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGACAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAACATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACGCTTCGCTATCCCGAGCCAGCGTTGCAGGGTTACAGGAt                                                                                    >  2:14556/1‑210 (MQ=255)
                                                                                                                tgCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAACATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACGCTTCGCTATCCCGAGCCAGCGTTGCAGGGTTACAGGAt                                                                                    <  1:14556/210‑1 (MQ=255)
                                                                                                                        cttACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCACATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAACATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACGCTTCGCTATCCCGAGCCa                                                                                                        >  2:329139/1‑183 (MQ=255)
                                                                                                                        cttACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAACATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACGCTTCGCTATCCCGAGCCa                                                                                                        <  1:329139/183‑1 (MQ=255)
                                                                                                                          tACCGCAACCATCGTCCGCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAACATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACGCTTCGCTATCCCGAGCCAGCGTTGCAGGGTTACAGGATACATAAACTATACGAATAGGTTCCAGTTTTATAATTTGCTGCATAac                                     <  1:316243/248‑1 (MQ=255)
                                                                                                                                              cgttgaGGGCCCAGCCCTGCCGAAATCGGCAAAGCGCAACTATTTAACGCGCGAGAAAAGGACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAACAGCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACGCTTCGCTATCCCGAGCCAGCGTTGCAGGGTTACAGGATACATAAACTATACGAATAGGTTCCAGTTTTATAATTTGCTGCATAACACCTGCGGCACCTGCTCGCGCCg              <  2:307523/246‑1 (MQ=255)
                                                                                                                                                    gggACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTGCCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAACATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACGCTTCGCTATCCCGAGCCAGCGTTGCAGGGTTACAGGATACATAAACTATACGAATAGGTTCCAGTTTTATAATTTGCTGCATAACACCTGCGGCACCTGCTCGCGCCGGGTCCa        <  2:133881/251‑1 (MQ=255)
                                                                                                                                                          aGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAACATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACGCTTCGCTATCCCGAGCCAGCGTTGCAGGGTTACAGGATACATAAACTATACGAATAGGTTCCAGTTTTATAATTTGCTGCATAACACCTGCGGCACCTGCTCGCGCCGGGTCCAGCAaca  <  1:76072/251‑1 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                 
GTTGCAGACAATACGCCCAGGTTTCGGCTAAGCACTCACACGACTTCTGGCTGGTAATACCCAGACTTGCGATCCATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACGCTTCGCTATCCCGAGCCAGCGTTGCAGGGTTACAGGATACATAAACTATACGAATAGGTTCCAGTTTTATAATTTGCTGCATAACACCTGCGGCACCTGCTCGCGCCGGGTCCAGCAACA  >  NC_000913/2913519‑2913922

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