Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,512,090 +GC coding (334/522 nt) aroK ← shikimate kinase AroK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,512,0881.G100.0% 43.1 / NA 13A112A (GCA→GCCaroKshikimate kinase AroK
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (5/8);  total (5/8)
*NZ_CP0092733,512,0882.C100.0% 43.9 / NA 13A112A (GCA→GCGaroKshikimate kinase AroK
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (5/8);  total (5/8)

TTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGT‑‑GCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTT  >  NZ_CP009273/3512002‑3512168
                                                                                       ||                                                                                
ttGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCgtgcgt‑‑gcg                                                                               >  2:437212/1‑89 (MQ=255)
               aaCGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCGCAAGTTGCTTTTCGa                                                                  <  1:271016/90‑1 (MQ=255)
                    ttCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCGCAAGTTGCTTTTCGATGGTc                                                             <  1:352375/90‑1 (MQ=255)
                    ttCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCGCAAGTTGCTTTTCGATGGTc                                                             <  1:482810/90‑1 (MQ=255)
                         gAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCGCAAGTTGCTTTTCGATGGTCGTTTc                                                        <  1:166742/90‑1 (MQ=255)
                             ttCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCGCAAGTTGCTTTTCGATGGTCGTTTCAAGa                                                    >  2:357863/1‑90 (MQ=255)
                              tCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCGCAAGTTGCTTTTCGATGGTCGTTTCAAGAt                                                   <  2:120080/90‑1 (MQ=255)
                                        gtgtTTCAACGTGCAGCAACGGGCGTTGTTTATCACGCTTCGTGCGTGCGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAAc                                         >  1:400414/1‑90 (MQ=255)
                                                 cGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGc                                <  1:84384/90‑1 (MQ=255)
                                                  gTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGTGCGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCg                               >  2:34553/1‑90 (MQ=255)
                                                         aaCGGGCGTTTTTTATCACGCTGCGTGCGTGCGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGAc                        <  1:138441/90‑1 (MQ=255)
                                                            gggCGTTTTTTATCACGCTGCGTGCGTGCGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGt                     >  1:109953/1‑90 (MQ=255)
                                                            gggCGTTTTTTATCACGCTGCGTGCGTGCGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGt                     >  1:565356/1‑90 (MQ=255)
                                                                          aCGCTGCGTGCGTGCGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACgg       <  2:385507/90‑1 (MQ=255)
                                                                               gcgtgcgtgcGCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGAttt  >  1:27508/1‑90 (MQ=255)
                                                                                       ||                                                                                
TTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGCGTTTTTTATCACGCTGCGTGCGT‑‑GCAAGTTGCTTTTCGATGGTCGTTTCAAGATAAACGACAACGCCACGAGCGGAAAGACGGTTACGCGTTTCACGGGATTT  >  NZ_CP009273/3512002‑3512168

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