Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 960,227 A→G intergenic (+132/‑37) aroA → / → ycaL 3‑phosphoshikimate 1‑carboxyvinyltransferase/periplasmic protease YcaL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913960,2270AG100.0% 47.6 / NA 14intergenic (+132/‑37)aroA/ycaL3‑phosphoshikimate 1‑carboxyvinyltransferase/periplasmic protease YcaL
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (8/6);  total (8/6)

GTCGGCGTTTTTTTTCGGACCTTGTGAGTCATTTTGATTAATGGTAGCGTCGCTTGTCAATGTAAGTTGTTGATACATAATATTTATATATGATTAATCAACGGATGATTCACATGAAGAATACTAAATTACTGCTGGCGATTGCGACCTC  >  NC_000913/960151‑960301
                                                                            |                                                                          
gTCGGCGTTTTTTTTCGGACCTTGTGAGTCATTTTGATTAATGGTAGCGTCGCTTGTCAATGTAAGTTGTTGATACgtaa                                                                         >  1:193183/1‑80 (MQ=255)
 tCGGCGTTTTTTTTCGGACCTTGTGAGTCATTTTGATTAATGGTAGCGTCGCTTGTCAATGTAAGTTGTTGATACgtaata                                                                       >  2:193186/1‑81 (MQ=255)
 tCGGCGTTTTTTTTCGGACCTTGTGAGTCATTTTGATTAATGGTAGCGTCGCTTGTCAATGTAAGTTGTTGATACgtaat                                                                        <  1:193159/80‑1 (MQ=255)
  cGGCGTTTTTTTTCGGACCTTGTGAGTCATTTTGATTAATGGTAGCGTCGCTTGTCAATGTAAGTTGTTGATACGTAatat                                                                      >  1:193184/1‑81 (MQ=255)
      gTTTTTTTTCGGACCTTGTGAGTCATTTTGATTAATGGTAGCGTCGCTTGTCAATGTAAGTTGTTGATACGTAATATTta                                                                   >  1:193179/1‑80 (MQ=255)
                       gtgAGTCATTTTGATTAATGGTAGCGTCGCTTGTCAATGTAAGTTGTTGATACGTAATATTTATATATGATTAATCAACg                                                  >  2:193182/1‑80 (MQ=255)
                         gAGTCATTTTGATTAATGGTAGCGTCGCTTGTCAATGTAAGTTGTTGATACGTAATATTTATATATGATTAATCAACGgat                                               >  2:193185/1‑81 (MQ=255)
                          aGTCATTTTGATTAATGGTAGCGTCGCTTGTCAATGTAAGTTGTTGATACGTAATATTTATATATGATTAATCAACGgatg                                              <  1:193160/81‑1 (MQ=255)
                                 ttGATTAATGGTAGCGTCGCTTGTCAATGTAAGTTGTTGATACGTAATATTTATATATGATTAATCAACGGATGATTcac                                        <  1:193161/80‑1 (MQ=255)
                                           gTAGCGTCGCTTGTCAATGTAAGTTGTTGATACGTAATATTTATATATGATTAATCAACGGATGATTCACATGAAGAATAc                             >  1:193188/1‑81 (MQ=255)
                                             aGCGTCGCTTGTCAATGTAAGTTGTTGATACGTAATATTTATATATGATTAATCAACGGATGATTCACATGAAGAATACt                            >  1:193189/1‑80 (MQ=255)
                                                   gCTTGTCAATGTAAGTTGTTGATACGTAATATTTATATATGATTAATCAACGGATGATTCACATGAAGAATACTAAATTAc                     >  1:193178/1‑81 (MQ=255)
                                                    cTTGTCAATGTAAGTTGTTGATACGTAATATTTATATATGATTAATCAACGGATGATTCACATGAAGAATACTAAATTAct                    <  1:193162/81‑1 (MQ=255)
                                                        tCAATGTAAGTTGTTGATACGTAATATTTATATATGATTAATCAACGGATGATTCACATGAAGAATACTAAATTActgct                 <  2:193163/80‑1 (MQ=255)
                                                              tAAGTTGTTGATACGTAATATTTATATATGATTAATCAACGGATGATTCACATGAAGAATACTAAATTACTGCTGGCGAtt          <  2:193164/81‑1 (MQ=255)
                                                                 gttgttGATACGTAATATTTATATATGATTAATCAACGGATGATTCACATGAAGAATACTAAATTACTGCTGGCGATTGc        >  2:193193/1‑80 (MQ=255)
                                                                       gATACGTAATATTTATATATGATTAATCAACGGATGATTCACATGAAGAATACTAAATTACTGCTGGCGATTGCGACctc  <  1:193165/80‑1 (MQ=255)
                                                                            |                                                                          
GTCGGCGTTTTTTTTCGGACCTTGTGAGTCATTTTGATTAATGGTAGCGTCGCTTGTCAATGTAAGTTGTTGATACATAATATTTATATATGATTAATCAACGGATGATTCACATGAAGAATACTAAATTACTGCTGGCGATTGCGACCTC  >  NC_000913/960151‑960301

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