Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_pae_pgi 4,234,667 Δ1 bp intergenic (‑1434/‑499) lysC ← / → yjbE lysine‑sensitive aspartokinase 3/extracellular polysaccharide production threonine‑rich protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_pgi4,234,6660A.100.0% 75.5 / NA 17intergenic (‑1433/‑500)lysC/yjbElysine‑sensitive aspartokinase 3/extracellular polysaccharide production threonine‑rich protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (8/9);  total (8/9)

CGCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATA  >  NC_000913_3_pae_pgi/4234528‑4234801
                                                                                                                                          |                                                                                                                                       
cGCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGat                                                                                                                             >  1:109253/1‑150 (MQ=255)
                       gCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCa                                                                                                      >  3:121128/1‑150 (MQ=255)
                                ccTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATc                                                                                             <  3:8749/150‑1 (MQ=255)
                                                             gTAAAACTTGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGTCGCATCAGGCATCGGTTGCCGGATgcggcgtgaac                                                                <  4:99117/150‑6 (MQ=37)
                                                                aaaCTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTcgcg                                                             <  2:109253/150‑1 (MQ=255)
                                                                      gCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTAt                                                       >  1:185152/1‑150 (MQ=255)
                                                                       cAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATc                                                      <  3:6257/150‑1 (MQ=255)
                                                                         gagaTTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGa                                                                                                                 >  2:60943/1‑89 (MQ=255)
                                                                         gagaTTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGa                                                                                                                 <  1:60943/89‑1 (MQ=255)
                                                                                             ccTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGat                                <  2:185152/150‑1 (MQ=255)
                                                                                             ccTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGat                                <  1:154366/150‑1 (MQ=255)
                                                                                                    aCCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGa                                                                                                                              <  3:60495/49‑1 (MQ=255)
                                                                                                    aCCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGa                                                                                                                              >  4:60495/1‑49 (MQ=255)
                                                                                                                  tGCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAAt           >  1:9717/1‑150 (MQ=255)
                                                                                                                   gCATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAAtt          >  2:6256/1‑150 (MQ=255)
                                                                                                                    cATTGAAGAGGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTa         >  1:85741/1‑150 (MQ=255)
                                                                                                                           gagGCGCTTGCAGTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCAGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATa  <  4:121128/150‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                       
CGCAGTAGAGATGCACCCCAAGGGCGTTTTATCCTTTTTACGAGATGTGGGCATCGCCTCTGTAAAACTGGCAGAGATTAGAGGCGTAAATCCCCTGGCCACCCCGAGAATAGATGCATTGAAGAGGCGCTTGCAGTGAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATA  >  NC_000913_3_pae_pgi/4234528‑4234801

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