Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 786,791 C→T intergenic (+106/+52) aroG → / ← gpmA 3‑deoxy‑D‑arabino‑heptulosonate‑7‑phosphate synthase, phenylalanine repressible/phosphoglyceromutase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913786,7910CT100.0% 39.0 / NA 13intergenic (+106/+52)aroG/gpmA3‑deoxy‑D‑arabino‑heptulosonate‑7‑phosphate synthase, phenylalanine repressible/phosphoglyceromutase 1
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (7/6);  total (7/6)

TCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGT  >  NC_000913/786676‑786895
                                                                                                                   |                                                                                                        
tCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAATAAAAAAGCCGACTCACTTGCAGt                                                                                   >  2:44078/1‑139 (MQ=255)
          ggTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAATAAAAAAGCCGACTCACTTGCAGTCGGCTTtctc                                                                         <  1:333770/139‑1 (MQ=255)
                                     gTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAATAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGc                                              <  2:103125/139‑1 (MQ=255)
                                                   gAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAATAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTg                                >  2:98482/1‑139 (MQ=255)
                                                     aTCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAGTGTGCTCCATTGTTAGCAATAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCa                              >  2:294953/1‑139 (MQ=255)
                                                       caccacAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAATAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAAc                            <  2:543800/139‑1 (MQ=255)
                                                        accacAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAATAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAAcc                           >  1:531258/1‑139 (MQ=255)
                                                         ccacAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAATAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTaaa                                                                 <  1:62056/100‑1 (MQ=255)
                                                         ccacAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAATAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTaaa                                                                 >  2:62056/1‑100 (MQ=255)
                                                                         tCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAATAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGc          >  1:489491/1‑139 (MQ=255)
                                                                                 cGTCGCATCAGGCAATGTGCTCCATTGTTAGCAATAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGt  <  2:489491/139‑1 (MQ=255)
                                                                                                    cTCCATTGTTAGCAATAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTAccc                                       <  1:201218/83‑1 (MQ=255)
                                                                                                    cTCCATTGTTAGCAATAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTAccc                                       >  2:201218/1‑83 (MQ=255)
                                                                                                                   |                                                                                                        
TCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGT  >  NC_000913/786676‑786895

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