Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,970,077 G→A *349* (TAG→TAA)  wzzE → Entobacterial Common Antigen (ECA) polysaccharide chain length modulation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,970,0770GA100.0% 65.4 / NA 22*349* (TAG→TAA) wzzEEntobacterial Common Antigen (ECA) polysaccharide chain length modulation protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (11/11);  total (11/11)

TATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTC‑GGGGG‑GCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATGGCGCCGTTGGTGCATGCGTTG  >  NC_000913/3969943‑3970201
                                                                                                                                        |                                                                                                                            
tATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGG‑GCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACaa                                                                                                                           <  1:311350/139‑1 (MQ=255)
  tCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGG‑GCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAacaaca                                                                                                                         <  2:376119/139‑1 (MQ=255)
  tCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGG‑GCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAacaaca                                                                                                                         >  1:50505/1‑139 (MQ=255)
            cGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGG‑GCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGt                                                                                                               >  1:564446/1‑139 (MQ=255)
                   cGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGG‑GCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCa                                                                                                        <  1:55420/139‑1 (MQ=255)
                                                     tGCCTTCCTGATGATTATGCGGGGCATTGTC‑GGGGGAGCTGATCGGGACTGGTGTCGCATTAAGCCGCCGTTGCTCGAAATAACAACGCTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGAtg                                                                      >  1:722650/1‑139 (MQ=255)
                                                                    tATGCGGGGCATTGTCGGGGGG‑GCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCg                                                                                                                 <  1:1024154/81‑1 (MQ=255)
                                                                    tATGCGGGGCATTGTCGGGGGG‑GCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCg                                                                                                                 >  2:1024154/1‑81 (MQ=255)
                                                                      tgcggggCATTGTCGGGGGG‑GCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTa                                                     >  2:1033210/4‑139 (MQ=255)
                                                                           ggCATTGTCGGGGGG‑GCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTgg                                                >  1:833459/1‑139 (MQ=255)
                                                                            gCATTGTCGGGGGG‑GCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGt                                               >  1:263857/1‑139 (MQ=255)
                                                                               ttGTCGGGGGG‑GCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTAcg                                            <  2:263857/139‑1 (MQ=255)
                                                                                tgtcgggggg‑gCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGt                                                               <  1:387240/114‑1 (MQ=255)
                                                                                tgtcgggggg‑gCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGt                                                               >  2:387240/6‑119 (MQ=255)
                                                                                       ggg‑gCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCt                                                                                                <  2:1021451/79‑1 (MQ=255)
                                                                                       ggg‑gCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCt                                                                                                >  1:1021451/1‑79 (MQ=255)
                                                                                                        ggTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATGgcg                    <  2:564446/139‑1 (MQ=255)
                                                                                                         gtgtCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTcgccg                                                                                           >  1:649043/1‑67 (MQ=255)
                                                                                                         gtgtCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTcgccg                                                                                           <  2:649043/67‑1 (MQ=255)
                                                                                                           gtCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactg                                                           <  1:719589/97‑1 (MQ=255)
                                                                                                           gtCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactg                                                           >  2:719589/1‑97 (MQ=255)
                                                                                                            tCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATGGCGCCGt                <  2:833459/139‑1 (MQ=255)
                                                                                                                   aaCCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATGGCGCCGTTGGTGCa         <  1:1033210/139‑1 (MQ=255)
                                                                                                                          cgTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATGGCGCCGTTGGTGCATGCGTTg  >  1:1070102/1‑139 (MQ=255)
                                                                                                                                        |                                                                                                                            
TATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTC‑GGGGG‑GCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATGGCGCCGTTGGTGCATGCGTTG  >  NC_000913/3969943‑3970201

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑