Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 437,421 G→A L239F (CTC→TTC)  yajO ← 1‑deoxyxylulose‑5‑phosphate synthase YajO

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913437,4210GA96.0% 55.4 / NA 25L239F (CTC→TTC) yajO1‑deoxyxylulose‑5‑phosphate synthase YajO
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (16/8);  minor base C (1/0);  total (17/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

CAGGATGCGGTTTATACGGCGTTTCCAGTTCGGCAATCTGTTCCGGCTTCAAAGTGATATCCACCGCGTTCAATAGCTCATCAAGCTGTTCTTCGCGCGAAGTTCCGATAATCGGTGCGGCAATGCCCGGTTTACTCAACAACCAGGCCAGCGCAACTTGTGCTCGTGTCGCCCCCAGTTCTTCACTGACGCCTGTTAACCGCTCTGCGATCTGCGCGTCATTTTCATCGCTTTCTTTATAGAGATTTTTCCCCACCTCATCAGACACCAGTCGTGCGGTAGTTTCTCCCCACGGACGCGTCAGACGGCCCCTTGCCAGCGGGCTCCATGGAATTACCGCCACGCCCTCCTGATAACACAGTGGTAGCATCTCGCGCTCTTCTTCACGATAAATCAGATTGTAGTGATCCTGCATACTGACAAACTGCGCCCAGCCGTGCTGTTTTTG  >  NC_000913/437178‑437625
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                            
cAGGATGCGGTTTATACGGCGTTTCCAGTTCGGCAATATGTTCCGGCTTCAAAGTGATATCCCACGCGTTCAATAGCTCCTAAAGCTGTTCTTCGCGCGAAGTTCCGATAATCGGTGCGGCAATGCCCGGTTTACTCAACAACCAGGCCAGCGCCACTTATGCTCGTGTCGCCCCCAGTTCTTCACTGACGCCTGTTAACCGCTCTGCGATCTGCGCGTCATTTTCATCGCTTTCTTTATAGAAATTTTTc                                                                                                                                                                                                       >  1:444674/1‑251 (MQ=255)
                         cAGTTCGGCAATCTGTTCCGGCTTCAAAGTGATATCCACCGCGTTCAATAGCTCATCCAGCTGTTCTTCGCGCGAAGTTCCGATAATCGGTGCGGCAATGCCCGGTTTACTCAACAACCAGGCCAGCGCAACTTGTGCTCGTGTCGCCTCCAGTTCTTCACTGACGCCTGTTAACCGCTCTGCGATCTGCGCGTCATTTTCATCGCTTTCTTTATAGAAATTTTTCCCCACCttatc                                                                                                                                                                                            >  1:506293/1‑237 (MQ=255)
                                       gTTCCGGCTTCAAAGTGATATCCACCGCGTTCAATAGCTCATCAAGCTGTTCTTCGCGCGAAGTTCCGATAATCGGTGCGGCAATGCCCGGTTTACTCAACGACCAGGCCAGCGCAACTTGTGCTCGTGTCGCCCCCAGTTCTTCACTGACGCCTGGTAACCGCTCTGCGATCTGCGCGTCATTTTCATCGCTTTCTTTATAGAAATTTTTCCCCACCTTATCAGACACCAGTCGTGCGGTAGTTTCTcc                                                                                                                                                                 <  1:486513/250‑1 (MQ=255)
                                                                                                 cgAAGTTCCGATAATCGGTGAGGGAATGTACGGTTTACTCAACAACCAGGCCAGCGCAACTCGTGCTCGTGTCGCCCCCAGTTCTTCACTGAAGCCTGTTAACCGCTCTGCGATCTGCGCGTCATTTTCATCGCTTTCTTTATAGCAATTTTTCCCCACCTTATCAGACACCAGTCGTGCGGTCGTTTCTCCCCACGGACGCGACAGACGGCCCCTTGCCAGCGGGCTCCATGGACTTACCGCCACGCcct                                                                                                      >  1:504113/1‑251 (MQ=255)
                                                                                                  gAAGTTCCGATAATCGGTGCGGCAATGCCCGGTTTACTCAACAACCAGGCCAGCGCAACTTGTGCTCGTGTCGCCCCCAGTTCTTCACTGACGCCTGTTAACCGCTCTGCGATCTGCGCGTCATTTTCATCGCTTTCTTTATAGAAATTTTTCCCCACCTTATCAGACACCAGTCGTGCGGTAGTTTCTccc                                                                                                                                                                >  1:235141/1‑192 (MQ=255)
                                                                                                     gTTCCGATACTCGGTGCGGCAATGCCCGGTTTACTCAACAACCAGGCCAGCGCAACTTGTGCTCGTGTCGCCCCCAGTTCTTCACTGACGCCTGTTCACCGCTCTGCGATCTGCGCGTCATTTTCATCGCTTTCTTTATAGAAATTTTTCCCCACCTTATCAGACACCAGTCGTGCGGTAGTTTCTCCCCACGGCCGCGTCAGACGGcccc                                                                                                                                          >  1:108623/1‑211 (MQ=255)
                                                                                                                       nccntgcncnntnncccnaacaacCATGCCAGCGCAACTNNTGCTCGGGTNCCCCCCAGTTCTTCACTGACGCCTGTTAACCGCTCTGCGATCTGCGCGTCATTNTCATCGCTTNNTTNATAGANANNNTNNNCCNCCTTATCAGACACCAGGCGNGCNNTAGTTTNTCCCCACGGACGCGTCAGACGGcccc                                                                                                                                          <  2:108623/175‑1 (MQ=255)
                                                                                                                                   ttACTCAACAACCAGGACAGCGCAACTTGTGCTCGTGTCGCCCCCAGTTCTTCACTGACGCCTGTTAACCGCTCTGCGATCTGCGCGTCATTTTCATCGCTTTCTTTATAGAAATTTTTCCCCACCTTATCAGACACCAGTCGTGCGGTAGTTTCTCCCCACGGACGCGTCAGACGGCCCCTTGCCAGCGGGCTCCATGGAATTACCGCCACGCCCTCCTGATAACACAGTGGTAGCATCTCGCGCtcttc                                                                    >  1:421808/1‑251 (MQ=255)
                                                                                                                                         aacaacCAGGCCAGCGCAACTTGTGCTCGTGTCGACCCCAGTTCTTCACTGACGCCTGTTAACCGCTCTGCGCTCTGCGCGTCATTTTCATCGCTTTCTTTATAGAAATTTTTcccc                                                                                                                                                                                                    >  1:424984/1‑117 (MQ=255)
                                                                                                                                            aacCAGGCNANNNCNNATTGTGCTCGTGTCGCCCCCAGTTCTTCACTGACGCCTGTTAACNGATCTGCGATCTGCGNGTNATTTNCATNGCTTTCTTTATAGAAATTTTTCCCCACCTTATCAGACACCAGTCGTGCGGTCGTTTCTCCCCACGGAc                                                                                                                                                         >  2:358003/1‑157 (MQ=255)
                                                                                                                                            aacCAGGCCAGCGCAACTGGGGCTCGTGGCGCCCCCAGTTCTTCACTGACGCCGGTGAACCGCTCTGCGATCGGCGCGTCAGTTGCATCGCTTTCTTAATAGAAATGTTTCCCCCCCATATCAGACACCAGTCGGGCGGTAGTTTCTCCCCACGGAc                                                                                                                                                         <  1:358003/157‑1 (MQ=255)
                                                                                                                                                 ggCCAGCGCAACTTGTGCTCGTGTCGCCCCCAGTTCTTCACTGACGCCTGTTAACCGCTCTGCGATCTGCGCGTCATTTTCATCGCTTTCTTTATAGAAATTTTTCCCCACCTTATCAGACACCAGTCGTGCGGTAGTTTCTCCCCACGGACGCGTCAGACGGCCCCTTGCCAGCGGGCTCCATGGAATTACCGCCACGCCCTCCTGATAACACAGTGGTAGCATCTCGCGCTCTTCTTCACGATAAATc                                                       >  1:497752/1‑250 (MQ=255)
                                                                                                                                                                 gCTCGTGTNGNNNCNNGTTCTTCACTGACGCCTGTTAACCGCTCTGCGATCTGCGCGTCANTTTCATCGCTTTCTTNATNGAAANTTTNCCCCACCTTATCAGACACCAGTCGTGCGGTAGTTTCTCCCCACGGACGCGTCAGACGGCCCCTTGCCAGCg                                                                                                                                 >  2:345564/1‑160 (MQ=255)
                                                                                                                                                                 gCTCGTGTCGCCCCCATTTCTTCACTGACGCCTGTTAACCGCTCTGCGATCTGCGCGTCATTTTCATCGCTTTCTTTATAGAAATTTTGCCCCACCTTATCAGACACCAGTCGTGCGGTAGTTTCTCCCCACGGACGCGTCAGACGGCCCCTTGCCAGCg                                                                                                                                 <  1:345564/160‑1 (MQ=255)
                                                                                                                                                                       ntCGCCCCCAGTTCTTCACTGACGCCTGTTAACCGCTCTGCGATCTGCGCGTCATTTTCATCGCTTTCTTTATAGAAATTTTTCCCCACCTTATCAGACACCAGTCGTGCGGTAGTTTCTCCCCACGGACGCGTCAGACGGCCCCTTGCCAGCGGGCTCCATGGAATTACCGCCACGCCCTCCTGATAACACAGTGGTAGCATCTCGCGCTCTTCTTCACGATAAATCAGATTGTAGTGATCCTGCATACt                                >  1:126500/2‑251 (MQ=255)
                                                                                                                                                                                         cTGACGCCTGTTAACCGCTCTGCGAACTGCGCGTCATTTTCATCGCTTTCTTTATAGAACTTTTTCCCCACCTTATCAGACACCAGTCGTGCGGTAGTTTCTCCCCACGGACGCGTCAGACGGCCCCTTGc                                                                                                                                      >  1:182850/1‑131 (MQ=255)
                                                                                                                                                                                           gACGCCTGTTAACCGCTCTGCGATCTGCGCGTCATTTTCATCGCTTTCTTTATAGAAATTTTTCCCCACCTTATCAGACACCAGTCGGGCGGTAGTTTCTCCCCACGGACGCGTCAGACGGCCCCTTGCCAGCGGGCTCCATGGAATTACCGCCACGCCCTCCTGATAACACAGTGGTAGCATCTCGCGCTCGTCTTCACGATAAATCAGATTGTAGTGATCCTGCATACTGACAAACTGCGCCCAGcc              <  1:456883/249‑1 (MQ=255)
                                                                                                                                                                                             ngCCTGTTAACCGCTCTGCGATCTGCGCGTCATTTTCATCGCTTTCTTTATAGAAATTTTTCCCCACCTTATCAGACACCAGTCGTGCGGTAGTTTCTCCCCACGGACGCGTCAGACTGCCCCTTGCCAGCGGGCTCCATGGCATTACCGCCACGCCCTCCTGATAACACATTGGTAGCATATCGCGCTCTTCTTCACGATAAATCAGATTGTAGTGATCCTGCATACTGACAAACTGCGCCCAGCCGttc          >  1:302609/2‑249 (MQ=255)
                                                                                                                                                                                             cGCCTGTTAACCGCTCTGCGATCTGCGAGTCATTTTCATCGCTTTCTTTATAGAAATTTTTCCCCACCTTATCAGACACCAGTCGTGCGGTAGTTTCTCCCCACGGACGCGTCAGACGGCCCCTTGCCAGCGGGCTCCATGGAATTACCGCCACGCCCTCCTGATAACACAGTGGTAGCATCTCGCGCTCTTCTTCACGATAAATCAGATTGTAGTGAt                                          >  1:136962/1‑219 (MQ=255)
                                                                                                                                                                                              gCCTGTTAACCGCTCTGCGATCTGCGCGTCATTTTAATCGCTTTCTTTATAGAAATTTTTCCCCACCTTATCAGACACAAGTCGTGCGGTAGTTTCTCCCCACGGACGCGTCAGACGGCCCCTTGCCAGCGGGc                                                                                                                              >  1:175755/1‑134 (MQ=255)
                                                                                                                                                                                               ncTGTTAACCGCTCTGCGATCTGCGCGTCATTTTCATCGCTTTCTTTATAGAAATTTTTCCCCACCTTATCAGACACCAGTCGTGCGGTAGTTtc                                                                                                                                                                    >  1:105708/2‑95 (MQ=255)
                                                                                                                                                                                               ccTGTTAACCGCTCTGCGATCTGCGCGTCATTTTCATCGCTTTCTTTATAGAAATTTTTCCCCACCTTATCAGACACCAGTCGTGCGGTAGTTTCTCCCCACGGACGCGTCAGACGGCCCCTTGCCAGCGGGCTCCATGGAATTACCGCCACGCCCTCCTGATAACACAGTGGTAGCATCTCGCGCTCTTCTTCACGATAAATCAGATTGTAGTGAACCTGCATACTGACAAACGGCTCCCAGCCGTGCTg        <  1:152537/251‑1 (MQ=255)
                                                                                                                                                                                                cTGTTAACCGCTCTGCGATCTGCGCGTCATTTTCATCGCTTTCTTTATAGAAATTTTTCCCCACCTTATCAGACACCAGTCGTGCGGTAGTTTCTCCCCACGGACGCGTCAGACGGCCCCTTGCCAGCGGGc                                                                                                                              <  1:389512/132‑1 (MQ=255)
                                                                                                                                                                                                     aaCCGCTCTGCGATCTGCGCGTCATTTTCATCGCTTTCTTTATAGAAATTTTTCCCCACCTTATCAGACACCAGTCGTGCGGTAGTTTCTCCCCACGGACGCGTCAGACGGCCCCTGGCCAGCGGGCTCCATGGAATTACCGCCACGCCCTCCTGATAACACAGTGGTAGCATCTCGCGCTCTTCTTCACGATAAATCAGGTTGTAGTGATCCTGCATACTGACAAACTGCGCCCAGCCGTGCTGTTTTTg  <  1:449800/251‑1 (MQ=255)
                                                                                                                                                                                                              gCGATCTGCGCGTCATTTTCCTCGCTTTCTTTATAGACATTTTTCCCCACCTTATCCGAAACCAGTCGTGCGGTAGTTTCTCCACACGGACGCGTc                                                                                                                                                    >  1:416237/1‑96 (MQ=255)
                                                                                                                                                                                                                     tcgcgTCATTTTCATCGCTTTCTTTATAGAAATTTTGCCCCTCCTTATCAGACACCAGTCGTGCGGTAGTTTCTCCCCACGGACGAGTCAGACGGCCCCTTGCCCGCGGGCTCCAGGTAATTACCGCCACGCCCTCCTGATAACACAGTGGTAGCATCTCGCGCTCTTCTTCACTATAAATCAGATTGGAGTGATCCTGCATACTGACAAACTgcgc                    <  1:162221/216‑1 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                            
CAGGATGCGGTTTATACGGCGTTTCCAGTTCGGCAATCTGTTCCGGCTTCAAAGTGATATCCACCGCGTTCAATAGCTCATCAAGCTGTTCTTCGCGCGAAGTTCCGATAATCGGTGCGGCAATGCCCGGTTTACTCAACAACCAGGCCAGCGCAACTTGTGCTCGTGTCGCCCCCAGTTCTTCACTGACGCCTGTTAACCGCTCTGCGATCTGCGCGTCATTTTCATCGCTTTCTTTATAGAGATTTTTCCCCACCTCATCAGACACCAGTCGTGCGGTAGTTTCTCCCCACGGACGCGTCAGACGGCCCCTTGCCAGCGGGCTCCATGGAATTACCGCCACGCCCTCCTGATAACACAGTGGTAGCATCTCGCGCTCTTCTTCACGATAAATCAGATTGTAGTGATCCTGCATACTGACAAACTGCGCCCAGCCGTGCTGTTTTTG  >  NC_000913/437178‑437625

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