Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,122,267 A→C S35A (TCG→GCG)  wcaF ← colanic acid biosynthesis acetyltransferase WcaF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,122,2670AC100.0% 56.0 / NA 18S35A (TCG→GCG) wcaFcolanic acid biosynthesis acetyltransferase WcaF
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (7/11);  total (7/11)

CGGACGAATAACTACGTTTTTTCCTATTTTTGCTCCGAATAAACGTAATAAAAAAGCCCGCCAGCGATACAATACTTGTGGCGACCAGGCAAATATTGTTGCCTGTACTGCCCACCATAATTGCACTTTAATAGCGTTGCCGCCCCGGAACCCTTTCGGCACCGAGAAA  >  NZ_CP009273/2122184‑2122352
                                                                                   |                                                                                     
cggacgAATAACTACGTTTTTTCCTATTTTTGCTCCGAATAAACGTAATAAAAAAGCCCGCCAGCGATACAATACTTGTGGCGCCCAGGc                                                                                 >  1:530764/1‑90 (MQ=255)
        tAACTACGTTTTTTCCTATTTTTGCTCCGAATAAACGTAATAAAAAAGCCCGCCAGCGATACAATACTTGTGGCGCCCAGGCAAATAttg                                                                         <  1:201218/90‑1 (MQ=255)
        tAACTACGTTTTTTCCTATTTTTGCTCCGAATAAACGTAATAAAAAAGCCCGCCAGCGATACAATACTTGTGGCGCCCAGGCAAATAttg                                                                         <  1:378754/90‑1 (MQ=255)
             aCGTTTTTTCCTATTTTTGCTCCGAATAAACGTAATAAAAAAGCCCGCCAGCGATACAATACTTGTGGCGCCCAGGCAAATATTGTTGcc                                                                    <  1:17583/90‑1 (MQ=255)
               gTTTTTTCCTATTTTTGCTCCGAATAAACGTAATAAAAAAGCCCGCCAGCGATACAATACTTGTGGCGCCCAGGCAAATATTGTTGCCTg                                                                  >  1:467773/1‑90 (MQ=255)
                ttttttCCTATTTTTGCTCCGAATAAACGTAATAAAAAAGCCCGCCAGCGATACAATACTTGTGGCGCCCAGGCAAATATTGTTGCCTGt                                                                 >  1:359850/1‑90 (MQ=255)
                       cTATTTTTGCTCCGAATAAACGTAATAAAAAAGCCCGCCAGCGATACAATACTTGTGGCGCCCAGGCAAATATTGTTGCCTGTACTGccc                                                          <  2:157547/90‑1 (MQ=255)
                       cTATTTTTGCTCCGAATAAACGTAATAAAAAAGCCCGCCAGCGATACAATACTTGTGGCGCCCAGGCAAATATTGTTGCCTGTACTGccc                                                          <  2:49612/90‑1 (MQ=255)
                       cTATTTTTGCTCCGAATAAACGTAATAAAAAAGCCCGCCAGCGATACAATACTTGTGGCGCCCAGGCAAATATTGTTGCCTGTACTGccc                                                          <  2:378772/90‑1 (MQ=255)
                                  ccGAATAAACGTAATAAAAAAGCCCGCCAGCGATACAATACTTGTGGCGCCCAGGCAAATATTGTTGCCTGTACTGCCCACCATAATTGc                                               <  1:276534/90‑1 (MQ=255)
                                  ccGAATAAACGTAATAAAAAAGCCCGCCAGCGATACAATACTTGTGGCGCCCAGGCAAATATTGTTGCCTGTACTGCCCACCATAATTGc                                               <  2:81431/90‑1 (MQ=255)
                                              aataaAAAAGCCCGCCAGCGATACAATACTTGTGGCGCCCAGGCAAATATTGTTGcc                                                                    >  2:373829/1‑57 (MQ=255)
                                              aataaAAAAGCCCGCCAGCGATACAATACTTGTGGCGCCCAGGCAAATATTGTTGcc                                                                    <  1:373829/57‑1 (MQ=255)
                                              aataaAAAAGCCCGCCAGCGATACAATACTTGTGGCGCCCAGGCAAATATTGTTGCCTGTACTGCCCACCATAATTGCACTTTAATAGCg                                   >  2:492384/1‑90 (MQ=255)
                                                        cccGCCAGCGATACAATACTTGTGGCGCCCAGGCAAATATTGTTGCCTGTACTGCCCACCATAATTGCACTTTAATAGCGTTGCCGcccc                         <  2:389863/90‑1 (MQ=255)
                                                                               ggCGCCCAGGCAAATATTGTTGCCTGTACTGCCCACCATAATTGCACTTTAATAGCGTTGCCGCCCCGGAAcc                   <  1:556502/73‑1 (MQ=255)
                                                                               ggCGCCCAGGCAAATATTGTTGCCTGTACTGCCCACCATAATTGCACTTTAATAGCGTTGCCGCCCCGGAAcc                   >  2:556502/1‑73 (MQ=255)
                                                                               ggCGCCCAGGCAAATATTGTTGCCTGTACTGCCCACCATAATTGCACTTTAATAGCGTTGCCGCCCCGGAACCCTTTCGGCACCGAGaaa  >  2:3726/1‑90 (MQ=255)
                                                                                   |                                                                                     
CGGACGAATAACTACGTTTTTTCCTATTTTTGCTCCGAATAAACGTAATAAAAAAGCCCGCCAGCGATACAATACTTGTGGCGACCAGGCAAATATTGTTGCCTGTACTGCCCACCATAATTGCACTTTAATAGCGTTGCCGCCCCGGAACCCTTTCGGCACCGAGAAA  >  NZ_CP009273/2122184‑2122352

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