Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,748,841 C→T R240K (AGA→AAA)  ydhT ← uncharacterized protein YdhT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,748,8410CT85.7% 45.4 / ‑1.0 21R240K (AGA→AAA) ydhTuncharacterized protein YdhT
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base T (10/8);  total (11/10)
Fisher's exact test for biased strand distribution p-value = 5.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.38e-01

TTCGTTCAGGGACTTACGGCATGTGCAGAAATTGGTGAGGCGATGGCTCGGGCGGTCGTAAAGCCTGCATCCCGTGCTCGTCGGCGTCTTTCGTTTGATTAAAGGTAAAGCTATCAGCACTCAACAAACTGATGCTGTTCATTAAGACGATAGCGGACTCCGGCTCGAATATTGTTTTCACCTTCAATGGCGACGGCGAAGCGAACTCTTTCGCCGTCATGATATGCCAGTGCGGCGCTTCCTCCTGGCCCGAGAATAATTGAGTCGACCACGCCCGTGCTCACTACCGTTGAATTTTCACCGCTGGCAATGATGTGTACATTGTCGCCGCTGTTGGCAATTCTGGCGTTTGCGCCAAAGCTGGCGATTTGCACCAGGTCACCATTACTGGCGACATGACACCGTTCACCCAATGTACAAACGCGTACTCGCATTCCGGTATTCGCAATCCGAC  >  NC_000913/1748634‑1749087
                                                                                                                                                                                                               |                                                                                                                                                                                                                                                      
ttCGTTCAGGGACTTACGGCATGTGCAGAAATTGGTGAGGCGATGGCTCGGGCGGTCGTAAAGCCTGCATCCCGTGCTCGTCGGCGTCTTTCGTTTGATTAAAGGTAAAGCTATCAGCACTCAACAAACTGATGCTGTTCATTAAGACGATAGCGGACTCCGGCTCGAATATTGTTTTCACCTTCAATGGCGACGGCGAAGCGAACTTTTTCGCCGTCATGATATGCCAGTGCGGCGCTTCCTCCTGGcc                                                                                                                                                                                                              <  1:117748/250‑1 (MQ=255)
                                   tGAGGCGATGGCTCGGGCGGTCGTAAAGCCCGCATCCCGTGCTCGTCGGCGTCTTTCGTTTGATTAAAGCAAAAGCTCTCAGCACTCAACAACCTGATGCTGGTCATTAAAACGGTAGCGGACTACGGCTCGTATATTGATTTTGCCTTCAATGACCACGGCGACGCGAACTTTTTCCCCGTCATGAGAGACCAGTGCGGCGCTTCGTCCTGGTCCGATCATAATTGCGTCGGCCACGACCGTGCTCACt                                                                                                                                                                           >  2:233834/1‑250 (MQ=255)
                                                      gTCGTAAAGCCTGCATCCCGTGCTCGTCGGCGTCTTTCGTTTGATTAAAGGTAAAGCTATCAGCACTCAACAAACTGATGCTGTTCATTAAGACGATAGCGGACTCCGGCTCGAATATTGTTTTCACCTTCAATGGCGACGGCGAAGCGAACTTTTTCGCCGTCATGATATGCCAGTGCGGCGCTTCCTCCTGGCCCGAGAATAATTGAGTCGACCACGCCCGTGCTCAATACCGTTGAATTTTCACCGCt                                                                                                                                                       >  1:210083/1‑251 (MQ=255)
                                                      gTCGTAAAGCCTGCATCCCGTGCTCGTAGGCGTCTTTCGTTTGATTAAAGGTAAAGCTATCAGCACTCAACAAACTGATGCTGTTAATTAAGACGATAGCGGACTCCGGCTCGAATATTGTTTTCACCTTCAATGGCGACGGCGAAGCGAACTTTTTCGCCGTCATGATATGCCAGTGCGGCGCTTCCTCCTGGCCCGAGAATAATTGAGTCGACCACGCCCGTGCTCACTACCGTTGAATTTTCACCGCt                                                                                                                                                       >  1:107115/1‑251 (MQ=255)
                                                        cGTAAAGCCTGCATCCCGTGCTCGTCGGCGTCTTTCGTTTGATTAAAGGTAAAGCTATCAGCACTCAACAAACTGATGCTGTTCATTAAGACGATAGCGGACTCCGGCTCGAATATTGTTTTCACCTTCAATTGCGACGGCGAAGCGAACTTTTTCGCCGTCATGATATGCCAGTGCGGCGCTTCCTCCTGGCCCGAGAATAATTGAGTCGACCACGCCCGTGCTCACTACCGTTGAATTTTCACCGCTg                                                                                                                                                      <  1:124480/250‑1 (MQ=255)
                                                                cTGCATCCCGTGCTCGTCGGCGTCTTTCGTTTGATTAAAGGTAAAGCTATCAGCACTCAACAAACTGATGCTGTTCATTAAGACGATAGCGGACTCCGGCTCGAATATTGTTTTCACCTTCAATGGCGACGGCGAAGCGAACTTTTTCGCCGTCATGATATGCCAGTGCGGCGCTTCCTCCTGGCCCGAGAATAATTGAGTCGACCACGCCCGAGCTCACTACCGTTGAATTTTCACCGCTGGCAatgatg                                                                                                                                             >  1:92891/1‑251 (MQ=255)
                                                                     atcctgcctttttggcttCTTTCGTTTGAATTAAGGTAAAGTTATCTGTACTCAATAAACTGATTTGGTTCATTATGCCGATAGCCTACTCCTGCTCTAATATTGTTTTCTCCTTCAATGGCGACGGCTCAGCGAACGCTTTCGCCGTCATGCTATGCCAGTGCTGCGCTTCCTCCTGGCCCGAGAATAATTGAGTCGACCACGCCCTTGCTCACTACCGTTGAATTTTCCCCTCTGGCAatgaggtttac                                                                                                                                        <  2:264589/234‑8 (MQ=255)
                                                                                  gtcGTCTTGCGTTGGATTAAATGTAAAACTATCATCACTCAACAATCTGTTGCGTTTCATTAAGACGAGAGCGGACTCCGGCTCGAATAATTTTTTCACCGGCAATATCGACGGCGAAGCGAACTTTTTCGCCGTCATGATATGCCAGGGCGGCGCTTCCTCCTTGCCCGAGAATAATTGAGTCGACCACGCCCGTGCTCACTACCGTTGAATTTTCACCGCTGGTAATGCTGTGTACATTGTCGCCTCTg                                                                                                                           <  2:207717/249‑1 (MQ=255)
                                                                                  ggCGTCTTTCGTTTGATTAAAGGTAAAGCTATCAGCACTCAACAAACTGATGCTGTTCATTAAGACGATAGCGGACTCCGGCTCGAATATTGTTTTCACCTTCAATGGCGACGGCGAAGCGAACTTTTTCGCCGTCATGATATGCCAGTGCGGCGCTTCCTCCTGGCCCGAGAATAATTGAGTCGACCACGCCCGTGCTCACTACCGTTGAATTTTCACCGCTGGCAATGATGTGTACATTTTCGCCGCTg                                                                                                                           >  2:6748/1‑251 (MQ=255)
                                                                                                         tAAAGCTATCAGCAATCAACAAACTGAAGCTGTTCATTAAGACGATAGCGGACTCCGGCTCGAATATTGTTTTCACATTCAATGTCGACGGCGAAGCGAACTTTTTCGCCGTCATGAAATGCCAGCGCGGCGCTTACTCCGGGCCCGAGCATAATTGAGTCAACCAAGCCCGTGCTCACTACCGTTGAATTTACACCGATGGCAATGCTGTGTACCTTGTCGCCGCTGTTGTCAATTCTGGCGTTCGCGcc                                                                                                    >  2:285196/1‑251 (MQ=255)
                                                                                                                              aaCTGATGCTGTTCATTAAGACGATAGCGGACTCCGGCTCGAATATTGTTTTCACCTTCAATGGCGACGGCGAAGCGAACTTTTTCGCCGTCATGATATGCCAGTGCGGCGCTTCCTCCTGGCCCGAGAATAATTGAGTCGACCACGCCCGTGCTCACTACCGTTGAATTTTCACCGCTGGCAATGATGTGTACATTGTCGCCGCTGTTGGCAATTCTGGCGTTTGCGCCAAAGCTGGCGATTTGCAcc                                                                                 >  1:42901/1‑249 (MQ=255)
                                                                                                                               aCTGATGCTGTTCATTAAGACGATAGCGGACTCCGGCTCGAATATTGTTTTCACCTTCAATGGCGACGGCGAAGCGAACTCTTTCGCCGTCATGATATGCCAGTGCGGCACTTCCTCCTGGCCCGAGAATAATTGAGTCGACCACGCCCGTGCTCACTACCGTTGAATTTTCACCGCTGGCAATGATGTGTACATTGTCGCCGCTGTTGGCAATTCTGGCGTTTGCGCCAAAGCTGGCGATTCGCACCAgg                                                                              >  1:152700/1‑251 (MQ=255)
                                                                                                                                cTGATGCTGTTCATTAAGACGATAGCGGACTCCGGCTCGCATATTGTTTTCACCTTCAATGGCGACGGCGAAGCGAACTTTTTCGCCGTCATGATATGCCAGTGCGGCGCTTCCTCCTGGCCCGAGAATAATTGAGTCGACCACGCCCGTGCTCACTACCGTTGAATTTTCACCGCTGGCAATGATGTGTACATTGTCGCCGCTGTTGGCAATTCTGGCGTTTGCGCCAAAGCTGGCGATTTGCACCAGGt                                                                             >  1:260007/1‑251 (MQ=255)
                                                                                                                                          tcctttaGACTATAGCGGACTCCGGATCGAATAGTGTTTTCACCTTCAATGGCGACGGCGACTCGATCTTTTTCGCCGTCATGATATGCCAGTTCGGCGCTTCCTCCTGGCCCGAGAATAAGTGAGTCGACCACGCCCGTGCTCACTACCGTTGAAGTTTCACCGCTGGCAATGATGTGTACATTGTCGCCGCTGTTGGCAATTCTGGCGTTTGCGCCAAAGCTGGCGATTTGCACCAGGTCTCCATTAc                                                                    <  2:263356/244‑1 (MQ=255)
                                                                                                                                                 gacgacAGCGGACTCCGGCTCGAATATGTTTTTCACCTTCAATGGCGACTGCGAAGCGAACTTTTTCGCCGTCATGATATGCCAGTGCGGCGCTTCCTCCTGGCCCTAGAATAATTGAGTCGTCCACGCCCGTGTTCACTACCGTTGAATTTTCACCGCTGGCAATGATGTGTACATTGTCGCCGCTGTTGGCAATTCTGGCGTGTGCGCCAATGCTGGCGATTTGCACCAGGTCACCATTACTGGCGACa                                                            <  2:207115/251‑1 (MQ=255)
                                                                                                                                                       aGCGGACTCCGGCTCGAATATTGTTTTCACCTTCAATGGCGACGGCGAAGCGCACTTTTTCGCCGTCATAATATGCCAGTGCGGCGCTTCCTCCTGGCCCGAGAATAATTGAGGCGACCACGCCCGTGCTCACTACCGTTGAATTTTCACCGCTGGCAATGATGTGTACATTGTCGCCGCTGt                                                                                                                          >  1:287061/1‑183 (MQ=255)
                                                                                                                                                                          atTGTTTTCACCTTCAATGGCGACGGCGAAGCGAACTTTTTCGCCGTCATGATATGCCAGTGCGGCGCTTCCTCCTGGCCCGAGAATAATTGAGTCGACCACGCCCGTGCTCACTACCGTTGAATTTTCACCGCTGGCAATGATGTGTACATTGTCGCCGCTGTTGGCAATTCTGGCGTTTGCGCCAAAGCTGGCGATTTGCACCAGGTCACCATTACCGGCGCCATAACACCGTTCACCCAATGTACac                                    >  2:289757/1‑249 (MQ=255)
                                                                                                                                                                                                 cGGCGAAGCGAACTTTTTCGCCGACATGATATGCCAGTGCGGCGCTTCCTTCTGGCCCGAGAATAATTGAGTCGACCACGCCCGTGCTCACTACCGTTGAATTTTCACCGCTGGCAATGATGTGTACATTGTCGCCGCTGTTGGCAATTCTGGCGTTTGCGCCAAAGCTGGCGATTTGCACCAGGTCACCATTACTGGCGACATGACACCGTTCACCCAATGTACAAACGCGTACTCGCATTCCGGTAttt            <  1:228213/251‑2 (MQ=255)
                                                                                                                                                                                                  ggCGAAGCGAACTCTTTCGCCGTCATGATATGCCAGTGCGGCGCTTCCTCCTGGCCCGAGAATAATTGAGTCGACCACGCCCGTGCTCACTACCGTTGAATTTTCACCGCTGGCAATGATGTGTACATTGTCGCCGCTGTTGGCAATTCTGGCGTTTGCGCCAAAGCTGGCGATTTGCACCAGGTCACCATTACTGGCGACATGACACCGTTCACCCAATGTACAAACGCGTACTCGCATTCCGGTATTCg           <  1:51861/251‑1 (MQ=255)
                                                                                                                                                                                                   gcgaagcgaaCTTTTGCGCCGGCATGATATGCCAGTGCGGCGCTTACTCCTGGCCCGAGAATAATTGAGTCGACCACGCCCGTTCTCACTTCCGTTGAATTTTCACCGCTGGCATGGATGTGTACATTGTCGCCGCTGTTGGCAATTCTGGCGTTTGCGCCAAAGCTGGCGATTTGCACCAGGTCACCATTACTGGCTACAGGACACCGTTCACCCAATGTACAAACGCGTACTCGCATTCCTGTATTCGc          <  2:259607/251‑1 (MQ=255)
                                                                                                                                                                                                           aaCTTTTTCGCCGTCAGTATATGACAGTGCGGCGTCTCCACCTGGCCCGAGAATAATTGAGTCGACCATGCCCGTGCTCACTACCGTGGAATTTTCACCGCTGTCAATGATGTGTACATTGTCGCCGCTGTTGGCAATTCTGGCGTTGGCGCCCAAGCTGGCGATTTGCACCAGGTCACCATTACTGGCGACATGACACCGTTCGCCCAATTGACAAACGCGTCCTCGCTTTCCGGGATTCGCAATCCGAc  <  2:264365/251‑1 (MQ=255)
                                                                                                                                                                                                               |                                                                                                                                                                                                                                                      
TTCGTTCAGGGACTTACGGCATGTGCAGAAATTGGTGAGGCGATGGCTCGGGCGGTCGTAAAGCCTGCATCCCGTGCTCGTCGGCGTCTTTCGTTTGATTAAAGGTAAAGCTATCAGCACTCAACAAACTGATGCTGTTCATTAAGACGATAGCGGACTCCGGCTCGAATATTGTTTTCACCTTCAATGGCGACGGCGAAGCGAACTCTTTCGCCGTCATGATATGCCAGTGCGGCGCTTCCTCCTGGCCCGAGAATAATTGAGTCGACCACGCCCGTGCTCACTACCGTTGAATTTTCACCGCTGGCAATGATGTGTACATTGTCGCCGCTGTTGGCAATTCTGGCGTTTGCGCCAAAGCTGGCGATTTGCACCAGGTCACCATTACTGGCGACATGACACCGTTCACCCAATGTACAAACGCGTACTCGCATTCCGGTATTCGCAATCCGAC  >  NC_000913/1748634‑1749087

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