Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,648,720 C→T A19V (GCC→GTC)  ydfB → uncharacterized protein YdfB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,648,7200CT94.7% 51.3 / NA 19A19V (GCC→GTC) ydfBuncharacterized protein YdfB
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (11/7);  minor base G (0/1);  total (11/8)
Fisher's exact test for biased strand distribution p-value = 4.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATGGATTTCGACACAATCATGGAAAAGGCTTACGAAGAATACTTCGAAGGCCTTGCCGAAGGCGAAGAAGCTCTCAGCTTCAGTGAGTTTAAACAGGCGCTTTCCAGTTCGGCAAAATCTAACGGCTGATAAGCGAAACAGCACCGCGAGGAATCAGTATGCAGAAACGAGAACCCGTCATCATCGCGCCAGACTATACCGATGATGAACTTTATGAGTGGATGCGCCAGAAAATTAATGCAGCGCAGGATCTGAAATG  >  NC_000913/1648525‑1648923
                                                                                                                                                                                                   |                                                                                                                                                                                                           
tGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTTAGTGATATGTATTTCGACACAATCATGGAAAAGGCTTACGAAGAATACTGCGAAGGCCTGGTCGAAGGCGAAGAAGCGCTCAGCTTCAGTGAGTTTAAACATGCGCTTTCCAGtcc                                                                                                                                                       <  1:476672/250‑3 (MQ=255)
   caacaaCGAATCTCTGGTGTACGGCGTGTTTCCAAACCAAGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAACCCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATGGATTTCGACACAATCATGGAAAAGGCTTACGAAGAATACTTCGAAGGCCTTGTCGAAGGCga                                                                                                                                                                                                    >  1:425994/1‑202 (MQ=255)
                                           gACGGCACGTTCACCGCAATGACGTATACCAAAAGCGAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATGGATTTCGACACAATCATGGAAAAGGCTTACGAAGAATACTTCGAAGGCCTTGTCGAAGGCGAAGAAGCTCTCAGCTTCAGTGAGTTTAAACAGGCGCTTGCCAGTTCGGCAAAATCTAACGGCTGATAAGc                                                                                                                               <  1:401797/231‑1 (MQ=255)
                                                         cGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATGGATTTCGACACAATCATGGAAAAGGCTTACGAAGAATACTTCGAAGGCCTTGTCGAAGGCGAAGAAGCTCTCAGCTTCAGTGAGTTTAAACAGGCGCTTTCCAGTTCGGCAAAATCTAACGGCTGATAAGCGAAACANc                                                                                                                       <  1:396162/225‑1 (MQ=255)
                                                                    ataCCAACAGCAAACCGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGAGATATGGATTTCGACACAATCATGGAAAAGGCTTACGAAGAATACTTCGAAGGCCTTGTCGACGGCGACGAAGCACTCAGCTTCAGTGAGTTTAAACAGGCGCTTTCCAGTTCGGCAAAATCTCACGGCTGCTAAg                                                                                                                                >  1:296194/1‑205 (MQ=255)
                                                                             gCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATGGATTTCGACACAATCATGGAAAAGGCTTACGAAGAATACTTCGAAGGCCTTGTCGAAGGCGAAGAAGCTCTCAGCTTCAGTGAGTTTAAACAGGCGCTTTCCAGTTCGGCAAAATCTAACGGCTGATAAGCGAAACAGCACCGCGAGGAATCAGTATGCAGAAACGAGAAccc                                                                                     <  1:166157/239‑1 (MQ=255)
                                                                                                                 ggctggAAAGAAACTCAGGTGAGTGANATGGATNNCGNCACAATCATGGAAAAGGCTNACNNNGNNNANTTCGANGGNNTTGTCGAAGNCGAAGAAGCTCTCAGCTTCAGTGAGTTTAAACAGGCGCTTTCCTGTTc                                                                                                                                                       >  2:118981/1‑137 (MQ=255)
                                                                                                                 ggctggAAAGAAACTCAGGTGAGTGAGATGGATTTCGACACAATCATGGAAAAGGCTTACGAAGAATACTTCGAAGGCCTTGGCGAAGGCGAAGAAGCTCTCAGCTTCAGTGAGTTTAAACAGGCGCTTTCCAGTTc                                                                                                                                                       <  1:118981/137‑1 (MQ=255)
                                                                                                                      gaaagaaaCTCAGGTGAGTGATATGGATTTCGACACAATCATGGAAAAGGCTTACGAAGAATACTTCGAAGGCCTTGTCGAAGGCGAAGAAGCTCTCAGCTTCAGTGAGTTTAAACAGGCGCTTTCCAGTTCGGCAAAATCTAACGGCTGATAAGCGAAACAGCACCGCGAGGAATCAGTATGCAGAAACga                                                                                           >  1:248472/1‑192 (MQ=255)
                                                                                                                             aCTCAGGTGAGTGATATGGATTTCGACACAATCATGGAAAAGGCTTACGAAGAATACTTCGAAGGCCTTGTCGAAGGCGAAGAAGCTCTCAGCTTCAGTGAGTTTAAACAGGCGCTTTCCAGTTCGGCAAAATCTAACGGCTGATAAGCGAAACAGCAcc                                                                                                                    >  1:36762/1‑160 (MQ=255)
                                                                                                                                                 tttCGACACAATCATGGAAAACGCTTACGAAGAATACTTCGAAGGCCTTGTCGAAGGCGAAGAAGCTCTCAGCt                                                                                                                                                                                      <  1:296257/74‑1 (MQ=255)
                                                                                                                                                    ngACACAATCATGGAAAAGGCTTACGAAGAATACTTCGAAGGCCTTGTCGAAGGCGACGAAGCTCTCAGCTTCAGTGAGTTTAAACAGGCGCTTTCCAGTTCGGCAAAATCTAACGGCTGATAAGCGAAACAGCACCGCGAGGAATCAGTATGCAGAAACGAGAACCCGTCatcatc                                                                            >  1:57097/2‑177 (MQ=255)
                                                                                                                                                    cGACACAATCATGGAAAAGGCTTACGAAGAATACTTCGAAGGCCTTGTCGAAGGCGAAGAAGCTCTCAGCTTCAGTGAGTTTAAACAGGCGCTTTCCAGTTCGGCAAAATCTAACGGCTGATAAGCGAAACAGCACCGCGAGGAATCAGTATGCAGAAACGAGAACCCGTCATCATCGCGCCAGACTATACCGATGATGAACTTTATGAGTGGATGCGCCAGAAAAt                          >  1:350371/1‑227 (MQ=255)
                                                                                                                                                    cGACACAATAATGGAAAAGGCTTACGAAGAATACTTCGAAGGCCTTGTCGAAGGCGAAGAAGCTCTCAGCTTCAGTGAGTTTAAACAGGCGCTTTCCAGTTCGGCAAAATCTAACGGCTGATAAGCGAAACAGCACCGCGAGGAATCAGTATGCAGAAACGAGAACCCGTCATCATCGCGCCAGACTATACCGATGATGAACTTTATGAGTGGATGCGCCAGAAAATTAATGCAGCGCAGGATCTGAAATg  >  1:109324/1‑251 (MQ=255)
                                                                                                                                                      acacAATCATGGAAAAGGCTTACGAAGAATACTTCGAAGGCATTGTCGAAGGCGAAGAAGCTCTCAGCTTCAGAGAGTTTaaa                                                                                                                                                                        >  1:222706/1‑83 (MQ=255)
                                                                                                                                                             cATGGAAAAGGNTTACGAAGAANACTNNGAAGNCNTNGTCAAAGGCGAAGAAGCTCGCAGCTTCAGNGAGTTTAAACAGGCGCTTTCCAGNNCGNNNATNNCTAAc                                                                                                                                          <  2:374055/106‑1 (MQ=255)
                                                                                                                                                             cATGGAAAAGGCTTACGAAGAATACTTCGAAGGCCTTGTCGAAGGCGAAGAAGCTCTCAGCTTAAGTGAGTTTAAACAGGCGCTTTCCAGTTCGGCAAAATCTAAc                                                                                                                                          >  1:374055/1‑106 (MQ=255)
                                                                                                                                                               tGGAANAGGCTNANNAAGAATACTTCGAAGGCCTNGTCGAAGGCGAAGAANCTCTNAGCTTNAGNg                                                                                                                                                                                >  2:216660/1‑66 (MQ=255)
                                                                                                                                                                         tnANGAAGAATACTTCGAAGGCCTTGTCGAAGGCGAAGAAGCNCTCANCNNCANTNATNTTaaa                                                                                                                                                                        <  2:222706/64‑1 (MQ=255)
                                                                                                                                                                                                   |                                                                                                                                                                                                           
TGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATGGATTTCGACACAATCATGGAAAAGGCTTACGAAGAATACTTCGAAGGCCTTGCCGAAGGCGAAGAAGCTCTCAGCTTCAGTGAGTTTAAACAGGCGCTTTCCAGTTCGGCAAAATCTAACGGCTGATAAGCGAAACAGCACCGCGAGGAATCAGTATGCAGAAACGAGAACCCGTCATCATCGCGCCAGACTATACCGATGATGAACTTTATGAGTGGATGCGCCAGAAAATTAATGCAGCGCAGGATCTGAAATG  >  NC_000913/1648525‑1648923

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: