Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,641,544 A→G Q201Q (CAA→CAG yjtD → putative tRNA/rRNA methyltransferase YjtD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,641,5440AG90.0% 48.0 / ‑1.4 20Q201Q (CAA→CAGyjtDputative tRNA/rRNA methyltransferase YjtD
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base G (7/11);  total (7/13)
Fisher's exact test for biased strand distribution p-value = 5.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.61e-01

GTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCGGATTATCCTTCGCTCAATCTGGGGCAGGCGGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGGCGAAAAGTGATGCAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGCCATGACATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAACAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGAAAAAAATATCACCAAATAAAAAACGCCTTAGTAAGTATTTTTC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4641316‑4641652
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                             
gTTGGCTGACGTTCTTACTGGTGTGCCGAGGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCATCTGCAAGCTTTACGCGAACGAGTCATGGCATTGCGGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGc                                                                                                                                                                                                                         <  1:244754/251‑1 (MQ=255)
          gttatgaCTGGTGCTCCGATGGTGGCAGATTATCCGTCGCTCAATCGGGGGCAGGCAGTGATTGACTATTGCTACCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCAGTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACGGGTCGACTGGTTACAGCAACGCCTGGGTCTTTTAGATCAACGAGacac                                                                                                                                                                                                               <  2:19368/245‑1 (MQ=255)
                                 ggcagATTATCCGTCGCGCAATCTGGGGCAGGCAGGGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACCCAACGGCAGACCAACATCAACTGCAAGCTGTACGCGAACGAGTCATGGCATTGCTGACGGCTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTtgctg                                                                                                                                                                                         <  2:363576/246‑1 (MQ=255)
                                       ttATCCTTCGCTCAATCTGGGGCAGGCGGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGGCGAAAAGTGATGCAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGCCATGACATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAACAACGCCTGGGGCTTTTagag                                                                                                                                                                                                                          <  1:97775/211‑1 (MQ=255)
                                                                        ggTCTATTGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGAAAAAAATATCACCAAATAAAAAACGCCtt                                                                                                                                                 >  1:275031/1‑251 (MQ=255)
                                                                               tGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGacac                                                                                                                                                                                                               <  1:297894/182‑1 (MQ=255)
                                                                               tGCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGacac                                                                                                                                                                                                               >  2:297894/1‑182 (MQ=255)
                                                                                gCTATCAATTAGCAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATTGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGAAAAAAATATCACCAAATAAAAAACGCCTTagtaagt                                                                                                                                          <  1:343633/250‑1 (MQ=255)
                                                                                  tATCAATTAGCAACATTAATACAACAACCGCCGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTtgct                                                                                                                                                                                          >  1:258650/1‑200 (MQ=255)
                                                                                  tATCAATGAGCAACAGTAAGACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTGTACGCGAACGAGTCATGTCATTGCGGACGACTCTGGCAGTGTCAGATGACATAAAACGGGTCGTCTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCTTTtgct                                                                                                                                                                                          <  2:258650/200‑1 (MQ=255)
                                                                                    acaaatAGCAACATTAATACAACAACCGACGAAAATGTACACAACGGCAGACCAACAGCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGAAAAAAATATCACCAAATAAAAAACGCCTTAGTAAGTAtttt                                                                                                                                     <  2:275031/246‑1 (MQ=255)
                                                                                            cAACATTAATACAACAACCGACGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTTCGCGAACGAGTCATGGCATTGCTGACGACGCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACTAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGAAAAAAATATCACCAAATAAAAAACGCCTTAGTAAGTATTTTTCagcttt                                                                                                                             >  1:15895‑M2/1‑245 (MQ=255)
                                                                                                                 cGAAAAGTGACACAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGAAAAAAATATCACCAAATAAAAAACGCCTTAGTAAGTATTTTTCagcttttcattctaactgaaacgggc                                                                                                         <  2:351288‑M2/250‑27 (MQ=255)
                                                                                                                         tgacacaacggCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTCGCTGCATGATATTGAAAAAAATATCACCAAATAAAAAACGCCTTAGTAAGTATTTTTCagattttcattctgactgcaacgggcaatatgtc                                                                                                 <  1:35024‑M2/246‑35 (MQ=255)
                                                                                                                               aaCGGCAGACCAACATCAACTGCAAGCTGTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGaaaa                                                                                                                                                                           <  1:36126/171‑1 (MQ=255)
                                                                                                                                             aTCAACTTCAAGCTTTACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTTGCAGATGACATAAAACTGGTCTACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGAAAAAAATATCACCAAATAAAAAACGCCTTAGTAAGTATTTTTCagcttttcattctgactgcaacgggcaatatgtctctgtgtggattaaaaaaag                                                                             >  1:370817‑M2/1‑197 (MQ=255)
                                                                                                                                                           ttACGCGAACGAGTCATGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGAAAAAAATATCACCAAATAAAAAACGCCTTAGTAAGTATTTTTCagcttttcattctgactgcaacgggcaatatgtctctgtgtggattaaaaaaagagtgtctgatag                                                                 >  1:351288‑M2/1‑183 (MQ=255)
                                                                                                                                                                          aTGGCATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGAAAAAAATATCACCAAATAAAAAACGCCTTAGTAAGTATTTTTCagcttttcattctgactgcaacgggcaatatgtctctgtgtggattaaaaaaagagtgtctgatagcagcttctgaactggtt                                                <  1:326251‑M2/251‑84 (MQ=255)
                                                                                                                                                                            gACATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAACAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGAAAAAAATATCACCAAATAAAAAACGCCTTAGTAAGTATTTTTCagcttttcattctgactgcaacgggcaatatgtctctgtgtggattaaaaaaagagtgtctgatagcagcttctgaactggttac                                              <  1:384032‑M2/251‑86 (MQ=255)
                                                                                                                                                                                                                         cGACTGGTTACAGCAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGAAAAAAATATCACCAAATAAAAAACGCCTTAGTAAGTATTTTTCagcttttcagtctgactgcaacgggcaatatgtctctgtgtggattaaaaaaagagtgtctgatagcagcttctgaactggttacctgccgtgagtaaattaaaattttattgacttaggtcactaaat  >  1:272196‑M2/1‑121 (MQ=255)
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                             
GTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCGGATTATCCTTCGCTCAATCTGGGGCAGGCGGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGGCGAAAAGTGATGCAACGGCAGACCAACATCAACTGCAAGCTTTACGCGAACGAGCCATGACATTGCTGACGACTCTGGCAGTGGCAGATGACATAAAACTGGTCGACTGGTTACAACAACGCCTGGGGCTTTTAGAGCAACGAGACACGGCAATGTTGCACCGTTTGCTGCATGATATTGAAAAAAATATCACCAAATAAAAAACGCCTTAGTAAGTATTTTTC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4641316‑4641652

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