Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 786,694 Δ1 bp intergenic (+9/+149) aroG → / ← gpmA 3‑deoxy‑7‑phosphoheptulonate synthase, Phe‑sensitive/2,3‑bisphosphoglycerate‑dependent phosphoglycerate mutase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913786,6930T.100.0% 77.6 / NA 19intergenic (+8/+150)aroG/gpmA3‑deoxy‑7‑phosphoheptulonate synthase, Phe‑sensitive/2,3‑bisphosphoglycerate‑dependent phosphoglycerate mutase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (9/10);  total (9/10)

GAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAG  >  NC_000913/786562‑786813
                                                                                                                                   |                                                                                                                        
gagCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATGTCGGATg                                                                                                                 <  2:229737/140‑1 (MQ=255)
gagCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATGTCGGATg                                                                                                                 <  1:359482/140‑1 (MQ=255)
          agCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATACCGATGCCCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATGTCGGATGCGCCGTCaga                                                                                                       <  2:275628/140‑1 (MQ=255)
                      cTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGGAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGGAAGGGTTAATGTCGGATGCGCCGTCAGAGGGGCGGATccc                                                                                           >  2:227478/1‑139 (MQ=255)
                            tACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAAt                                                                                     <  2:9339/140‑1 (MQ=255)
                              cGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATGTCGGATGCGCCGTCAGAGTGGCGTATCCGa                                                                                          <  1:318871/133‑1 (MQ=255)
                              cGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATGTCGGATGCGCCGTCAGAGTGGCGTATCCGa                                                                                          >  2:318871/1‑133 (MQ=255)
                                   agagCATCACCGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACcaca                                                                              >  2:375101/1‑140 (MQ=255)
                                         tCACCGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTg                                                                        >  2:2740/1‑140 (MQ=255)
                                           aCCGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGAt                                                                      >  1:269096/1‑140 (MQ=255)
                                            ccGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATa                                                                     <  1:375101/140‑1 (MQ=255)
                                                      cATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAg                                                           >  2:361631/1‑140 (MQ=255)
                                                             tGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCa                                                    <  1:227478/140‑1 (MQ=255)
                                                                    aTACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCa                                             <  1:2740/140‑1 (MQ=255)
                                                                    aTACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCa                                             <  2:269096/140‑1 (MQ=255)
                                                                           tGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCt                                      >  1:346362/1‑140 (MQ=255)
                                                                                                  aTGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAACAACAAAAAAGcc               <  2:346362/140‑1 (MQ=255)
                                                                                                      aGTAAAAGCGCGTCGCGGGTAAGGTTTAATGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACt           >  1:136579/1‑140 (MQ=255)
                                                                                                              cgcgTCGCGGGTAAGGTTTAATGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCa   >  2:123504/1‑140 (MQ=255)
                                                                                                               gcgTCGCGGGTAAGGTTTAATGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCTCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAg  >  2:322121/1‑140 (MQ=255)
                                                                                                                                   |                                                                                                                        
GAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAG  >  NC_000913/786562‑786813

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