Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,969,748 T→C S240P (TCG→CCG)  wzzE → enterobacterial common antigen polysaccharide co‑polymerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,969,7480TC96.4% 82.2 / NA 28S240P (TCG→CCG) wzzEenterobacterial common antigen polysaccharide co‑polymerase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (13/14);  total (13/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTATGTTGCTTTTGCCAGCCAGCGTGCAGCCAGCCATCTGAATGATGAGCTGAAAGGCGCA‑TGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGTACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTGGGCGTCCAATGCTTCAGGCTCGACTGGAAAATTTACAGGCCGTCGGTCCGGCCTTTGATCTCGACTATGATCAGAATCGGGCCATGTTAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTAC‑GCCGGAAGAACCGGTAAAACGCGATAGCCCACG  >  NC_000913/3969511‑3969992
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                     
gTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTATGTTGCTTTTGCCAGCCAGCGTGCAGCCAGCCATCTGAATGATGAGCTGAAAGGCGCA‑TGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTCCGCGCAGTGCGAc                                                                                                                                                                                                                                          >  1:68989/1‑251 (MQ=255)
                    cgcgCCTGACGCTAATAACCTGTTACGTCAGTATGTTGCTTTTGCCAGCCAGCGTGCAGCCAGCCATCTGAATGATGAGCTGAAAGGCGCA‑TGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTCCGCGCAGTGCGACAGATGTACCTGCCGAGGAAt                                                                                                                                                                                                                      >  1:112394/1‑251 (MQ=255)
                                                                         gtgccgcccactcgcTGAATGAGGAGCTGAAAGG‑GCAATGGGCGGCGATTTCAATCCAGTGGAAAGCTCAGGTGAAGCGTCAGTAGGAGGGGGCGAAACCCATCTACGACCGTCGGATAAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTCCGCGCAGTGCGACAGATGTACCTGCCGAGGAATTACCTTATTCAGAAATGTTCCTGCTTGGGCGTCCAATGCTTCAGGCTCGACTg                                                                                                                                                                 <  2:58760/237‑1 (MQ=255)
                                                                         gTGCAGCCAGCCATCTGAATGATGAGCTGAAAGGCGCA‑TGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTCCGCGCAGTGCGACAGATGTACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTGGGCGTCCAATGCTTCAGGCTCGACTg                                                                                                                                                                 >  1:270912/1‑251 (MQ=255)
                                                                         gTGCAGCCAGCCATCTGAATGATGAACTGACAGGCGCA‑TGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGGGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGATGAGCAGCATAATATTCCGCGCAGTGCGACAGATGGACCTGCCGAGGAATTACCTGCTTCAGAAATGTTCCTGCTTGGGCGTCCAATGCTTCAGGCTCGACTg                                                                                                                                                                 >  2:149129/1‑251 (MQ=255)
                                                                           gCAGCCAGCCATCTGAATGATGAGCTGAAAGGCGCA‑TGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTCCGCGCAGTGCGACAGATGTACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTGGGCGTCCAATGCTTCAGGCTCGACTGGa                                                                                                                                                               <  1:11852/251‑1 (MQ=255)
                                                                                        tGAATGATGAGCTGAAAGGCGCA‑TGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTCCGCGCAGTGCGACAGATGTACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTGGGCGTCCAATGCTTCAGGCTCGACTGGAAAATTTACAGGcc                                                                                                                                                  >  1:323805/1‑251 (MQ=255)
                                                                                         gAATGATGCGCTGAAAGGCGCT‑GGGGCGGCGCGTCCCTTCCAGATGAAAGTGCAGGGGAAGCGACAGGAATAGGTGGCGAAAGCCATCTAGGACCGCCGGATGAACCGAATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATCCCGCGCAGTGCGACAGATGTACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCGTGGGCGTCCAATGCTTCAGGCTCGACTGGAAAATTTACAGGcc                                                                                                                                                  <  2:275633/250‑1 (MQ=255)
                                                                                            tgatgaGCTGAAAGGCGCA‑TGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTCCGCGCAGTGCGACAGATGTACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTggg                                                                                                                                                                                         >  1:22755/1‑208 (MQ=255)
                                                                                            tgatgaGCGAAAAGGCGCA‑TGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCAATCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTCCGCGCAGTGCGACAGATGTACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTggg                                                                                                                                                                                         <  2:22755/208‑1 (MQ=255)
                                                                                                           cgcA‑TGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTCCGCGCAGTGCGACAGATGTACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTGGGCGTcc                                                                                                                                                                                    >  2:332497/1‑198 (MQ=255)
                                                                                                           cgcA‑TGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTCCGCGCAGTGCGACAGATGTACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTGGGCGTcc                                                                                                                                                                                    <  1:332497/198‑1 (MQ=255)
                                                                                                           cgcA‑TGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTCCGCGCAGTGCGACAGATGTACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTGGGCGTCCAATGCTTCAGGCTCGACTGGAAAATTTACAGGCCGTCGGTCCGGCCTTTGAtc                                                                                                                               <  1:305331/251‑1 (MQ=255)
                                                                                                             cA‑TGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTCCGCGCAGTGCGACAGATGTACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTGGGCGTCCAATGCTTCAGGCTCGACTGGAAAATTTACAGGCCGTCGGTCCGGCCTTTGAtctc                                                                                                                             <  1:56374/251‑1 (MQ=255)
                                                                                                                       cgcgTTCCAGCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTCCGCGCAGTGCGACAGATGTACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTGGGCGTCCAATGCTTCAGGCTCGACTGGAAAATTTACAGGCCGTCGGTCCGGCCTTTGATCTCGACTATGAt                                                                                                                    <  1:44612/251‑1 (MQ=255)
                                                                                                                         cgTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTCCGCGCAGTGCGACAGATGTACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTGGGCGTCCAATGCTTCAGGCTCGACTGGAAAATTTACAGGCCGTCGGTCCGGCCTTTGATCTCGACTATGATCa                                                                                                                  <  1:26003/251‑1 (MQ=255)
                                                                                                                            ccctaccAGAGGAAAACTCAGGTGATGCGTAAGGAATAGGTGGCGAAAGCTACCTTCGACCGCCGTATGAACAGCCTTGAGCAGGCGCTGAAAATTGCTGAGCAGTATAATATTTCGCGAAGTGCGACAGATGTACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTGGGCGTCCAATGCTTCAGGCTCGACTGGAAAATTTACAGGCCGTCGGTCCGGCCTTTGATCTCGACTATGATCAGaa                                                                                                               <  2:293782/246‑1 (MQ=255)
                                                                                                                             ccacccAGAGGAAAGCTCAGGGTAAGCGAGAGGAAGAGGGGGCGAAAGCCATCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGAATAATATTCCGCGCAGTTCGACAGAGGTACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTGGGCTTCCAATGCTTCAGGCTCGACTGGAAAATTTACAGGCCGTCGGTCCGGCCTTTGATCTCGACTAGGATCAGAAt                                                                                                              <  2:164803/247‑1 (MQ=255)
                                                                                                                             ccaccaaGATGAAAGCACAGGTGAGGCGTCAGGAAGAGGGGGGGAAAGTCATCTACAACCGGCGGATGAACAGCATGGAGCAGGGGCAGAAAATTGCTGAGCAGCATAATATTCCGCGCAGTGCGACAGATGTACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTGGGCGTCCAATGCTTCAGGCTCGACTGGAAAATTTACAGGCCGTCGGTCCGGCCTTTGATCTCTACTATGATCAGAAt                                                                                                              <  2:141961/245‑1 (MQ=255)
                                                                                                                                     aTGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTCCGCGCAGTGCGACAGATGTACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTGGGCGTCCAATGCTTCAGGCGCGACTGGAAAATTTACAGGCCGTCGGTCCGGCCTTTGATCTCGACTATGATCAGCATCGGGCCAt                                                                                                      >  2:10010/1‑251 (MQ=255)
                                                                                                                                               aGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTCCGCGCAGTGCGACAGATGTACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTGGGCGTCCAATGCTTCAGGCTCGACTGGAAAATTTACAGGCCGTCGGTCCGGCCTTTGATCTCGACTATGATCAGAATCGGGCCATGTTAAACAcc                                                                                            >  1:128092/1‑251 (MQ=255)
                                                                                                                                                              aaGAGGTGGCGAAAGCCTGCTCCGACCGCCGGATGAACAGCATTTAGCAGGCGCTGAAAAATGCTGAGCAGCATAATATTCCGCGCAGTGCGACAGATGTACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTGGTCGTCCAATGCTTCAGGCTCGACTGGAAAATTTACAGGCCGTCGGTCCGGCCTTTGATCTCGACTATGATCAGAATCGGGCCATGTTAAACACCCTGAATGTTGGTCCa                                                                             <  1:11779/251‑1 (MQ=255)
                                                                                                                                                                      gCGAAAGCCATCTACGACCGCGGGATGAACAGCATTGAGCAGGCGCTGAAAATTTCTGAGCAGCATAATATTCCGCGCAGTGCGACAGATGTACCTGCCGAGTAATTACCTGATTCAGAAATGTTCCTGCTTGGGCGTCCAATGCTTCAGGCTCGACTGGAAAAGTTACAGGCCGTCGGTCCGGCCTTTGATCTCGACTATGATCAGAATCGGGCCATGTTAAACACCCTGAATGTTGGTCCAACCCTGGa                                                                     <  1:307206/251‑1 (MQ=255)
                                                                                                                                                                                                     gCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTCCGCGCAGTGCGACAGATGTACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTGGGCGTCCAATGCTTCAGGCTCGACTGGAAAATTTACAGGCCGTCGGTCCGGCCTTTGATCTCGACTATGATCAGAATCGGGCCATGTTAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGt                                      >  1:189714/1‑251 (MQ=255)
                                                                                                                                                                                                                  gcTGAAAATTGCTGGGGAGAATAATATCCCGCGCAGTGCGACAGATGTACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTGGGCGTCCAATGCTTCAGGCTCGACTGGAAAATTTACAGGCCGTCGGTCCGGCCTTTGATCTCGACTATGATCAGAATCGGGCCATGTTAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTAC‑GCCGGAAGAAc                        <  2:128092/251‑1 (MQ=255)
                                                                                                                                                                                                                  gcTGAAAATTGCTGAGCAGAATAATATTCCGCGCAGTGCGACAGATGTACATGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTGGGCGTCCAATGCTTCAGGCTCGACTGGAAAATTTACAGGCCGTCGGTCCGGCCTTTGATCTCGACTATGATCAGAATCGGGCCATGTTAAACACCCTGCATTTTGGTACAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTGC‑GCAGgaagacc                        >  2:112902/1‑249 (MQ=255)
                                                                                                                                                                                                                                        aataTTCCGCGCAGTGCGACAGATGTACATGCCGCGGCATTACCTGATTCAGAAATATTACTGCTTGCGCGTCCAATGCTTCAGGCTCGACTGGAAAATTTACAGGCCGTCGGTCCGGCCTTTGATCTCGACTATGATCAGAATCGGGCCATGTTAAACACCCTGAATGTTGGTCTAACCCTGGATCCGCGTTTTCAGACCAATCGCTATTTG‑GGACTGCCGAAAGACCCGGTAAACAGCGATAGCCCAcg  >  2:291872/1‑251 (MQ=255)
                                                                                                                                                                                                                                        aataTTCCGCGCAGAGCGACAGATGTACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTGGGCGTCCAATGCTTCAGGCTCGACTGGAAAATTTACAGGCCGTCGGTCCGGCCTTTGATCTCGACTATGATCAGAATCGGGCCATGTTAAACACCCTGCATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTTTCGCTATTTGCGTAC‑GCCGGAAGAACCGGTAAAACGCGATAGCCCAcg  >  2:313505/1‑251 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                     
GTGAAGCTTATTGCCGAAACCGCGCCTGACGCTAATAACCTGTTACGTCAGTATGTTGCTTTTGCCAGCCAGCGTGCAGCCAGCCATCTGAATGATGAGCTGAAAGGCGCA‑TGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGTACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTGGGCGTCCAATGCTTCAGGCTCGACTGGAAAATTTACAGGCCGTCGGTCCGGCCTTTGATCTCGACTATGATCAGAATCGGGCCATGTTAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTAC‑GCCGGAAGAACCGGTAAAACGCGATAGCCCACG  >  NC_000913/3969511‑3969992

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