Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,641,131 T→C L64L (TTG→CTG)  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,1310TC85.0% 39.3 / 1.0 20L64L (TTG→CTG) yjtDputative tRNA/rRNA methyltransferase YjtD
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base C (7/10);  total (8/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.18e-01

CGTAATAAGTTGACAGTGATCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTTAGCGATCTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTGGTGATATTATTGATAATATTAAAGTTTTCCCGACATTGGCTGAATCGTTACACGATGTCGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCCACGCCAGTTGAACTGGTGCCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGTTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCGGATTATCCTTCGCTCAATCTGGGG  >  NC_000913/4640883‑4641376
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      
agTAATAAGTTGAAAGTGATCACCCGTTTCGCGGTTATTTGATCAAGAAGAGGGGAAATAGGTGTATAAGTATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTTAGCGATCTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGTTGGGTCGCACATGGATCTGGTGATATTATTGATAATTTTAAAGTTTTCCCGACATTg                                                                                                                                                                                                                                                      <  2:203577/250‑1 (MQ=255)
                                                        aaTATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGCGAACTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTGGTGATATTATTGATAATATTAAAGTTTTCCCGACACTGGCTGAATCGTTACACGATGTCGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGc                                                                                                                                                                                             <  1:133157/251‑1 (MQ=255)
                                                           atGCGTATAACGATTATTCTGGCCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGCGAACTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGTGTCGCACATGGATCTGGTGAGATTATTGATAATATTAAAGTTTTCCCGACACTGGCTGAATCGTTACACGATGTCGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGCGa                                                                                                                                                                                           <  1:344134/250‑1 (MQ=255)
                                                                                   gCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGCGAACTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTGGTGATATTATTGATAATATTAAAGTTTTCCCGACACTGGCTGAATCGTTACACGATGTCGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGt                                                                                                                                                                  >  1:350921/1‑251 (MQ=255)
                                                                                      cccGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGTTTTAGCGAACTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTGGTGATATTATTGATAATATTAAAGTTTTCCCGACACTGGCTGAATCGTTACACGATGTCGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGtaga                                                                                                                                                               <  1:208743/251‑4 (MQ=255)
                                                                                        cGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTTAGCGATCTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTGGTGATATTATTGATAATATTAAAGTTTTCCCGACATTGGCTGAATCGTTACACGATGTCGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCCACGCCAGTTGAAc                                                                                                                                                             >  1:158857/1‑251 (MQ=255)
                                                                                                   aaaaTATTGGGGCAGCGGCGCGGGCAATAAAACCGATGGGTTTTAGCGTACTGCGGATTGGCGATAGTCAGGCACACCTGTAGCCAGCCACCCGCTGGGTCGCACAGTGATCTGGTGATATTATTGATAATATTAAAGTTTTCCCGACACTGGCTGAATCGTTACACGATGTCGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATAGCATTACTACGCTACTCCAGTAGAACTGGTACCGCTg                                                                                                                                                  <  2:350921/251‑1 (MQ=255)
                                                                                                                         ggCAATGAAAACGATGGGTTTTAGCGAACTGCGGATTGTCGATAGTCAGGCACACCTTGAGCCAGCCACCCGCTGGGTCGCACATGGATCTGGTGATATTATTGATAATATTAAAGTTTTCCCGACACTGGCTGAATCGTTACACGATGTCGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCATTAGAACTGGTTCCGCTGTTAGAGGAAAAATCTTCATGGa                                                                                                                            <  1:384754/251‑1 (MQ=255)
                                                                                                                           aaaTGAAAACGATGGGGTTTAGCGATCTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTGGTGATATTATTGATAATATTAAAGTTTTCCCGACATTGGCTGAATCGTTACACGATGTCGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCCACGCCAGTTGAACTGGTGCCGCTGTTAGAGGAAAAATCTTCATGGATg                                                                                                                          <  2:158857/250‑1 (MQ=255)
                                                                                                                                                                       cAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTGGTGATATTATTGATAATATTAAAGTTTTCCCGACACTGGCTGAATCGTTACACGATGTCGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGAc                                                                              >  1:66593/1‑251 (MQ=255)
                                                                                                                                                                        aGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTGGTGATATTATTGATAATATTAAAGTTTTCCCGACACTGGCTGAATCGTTACACGATGTCGATTTCACTGTCGCCACCACTGCGCGCCGTCGGGCGAAATATCATTACTACGCTACGCCCGTAGAACTGGTACCGCTGTTAGCGGGAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACt                                                                             >  2:47813/1‑251 (MQ=255)
                                                                                                                                                                                      cctgccaCCCTCTGGGTCGCACAGGGATCTGGTGATATTATTGATAATATTAAAGTTTTCCCGACACTGGCTGAATCGTTACACGATGTCGATTTCACTGTCGCCACCACTGCGCGCAGTCGTTCGAAATATCATTTCTACGCTACGCCAGTAGAACTGGTACCGCGGTTAGAGGAAAAATCTTCATGGAGGGGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGc                                                               <  1:236985/248‑1 (MQ=255)
                                                                                                                                                                                            aCCCGCTGGGTCGCACATGGATCTGGTGATATTATTGATAATATTAAAGTTTTCACGACACTGGCTGAATCGTTACACGATGTCGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGACGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGTAGAGTTAGCGTTGGc                                                         >  1:306682/1‑251 (MQ=255)
                                                                                                                                                                                                   gggTCGCCCATGGATCGGGTGATATTATTGATAATATTAAAGTTTTCCCGACACTGGCTGAATCGTTACACGATGTCGATTTCACTGTCGCCACCTCTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGACGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGt                                                   <  1:46100/250‑1 (MQ=255)
                                                                                                                                                                                                           caTGGATCTGGTGATATTATTGATAATATTAAAGTTTTCCCGACACTGGCTGAATCGTTACACGATGTCGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTgg                                          >  1:313326/1‑251 (MQ=255)
                                                                                                                                                                                                                    ggtattatatatgatatatttAAAGTTTTCCCGACACTGGCTGAATCGTTACGCGATGTCGATTTCACTGTCTCCACGACGGCGCGCAGTCGGGCGAAATAGCATTACTACGCTACGCCAGTCGAACGGTTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGAGCGGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGAt                                 <  2:142672/232‑1 (MQ=255)
                                                                                                                                                                                                                          attattTAAAATATTAAAGTTGTCCCGACACGGGTTGAATCGGTACACGATGTCGTTTCTACTTTCGCCACCTCTGCGCGCAGTCGGGCGAAATATCATTACTACGCTTCGCCATTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGTCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGAtggtgg                            <  2:304579/250‑1 (MQ=255)
                                                                                                                                                                                                                          attattGATAATATTAAAGTTTTCCCGACACTGGCTGAATCGTTACACGATGTCGATTTCACTGTCGCAACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGATAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGTAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGc                           >  1:84805/1‑251 (MQ=255)
                                                                                                                                                                                                                                           aGTTTTCCCGACACTGGCTGAATCGTTACACGATGTCGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCa          <  1:206155/251‑1 (MQ=255)
                                                                                                                                                                                                                                                   cGACACTGGCTGAATCGTTACACGATGTCGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTgggg  >  1:207890/1‑251 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      
CGTAATAAGTTGACAGTGATCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTTAGCGATCTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTGGTGATATTATTGATAATATTAAAGTTTTCCCGACATTGGCTGAATCGTTACACGATGTCGATTTCACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCCACGCCAGTTGAACTGGTGCCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGTTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCGGATTATCCTTCGCTCAATCTGGGG  >  NC_000913/4640883‑4641376

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