Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,021,281 C→A 43.5% R253R (CGC→CGA fliK → flagellar hook‑length control protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,021,2810CA43.5% 72.3 / 60.7 62R253R (CGC→CGAfliKflagellar hook‑length control protein
Reads supporting (aligned to +/- strand):  ref base C (25/10);  new base A (12/15);  total (37/25)
Fisher's exact test for biased strand distribution p-value = 3.97e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTTGAGTGCACCGCTGGGTTCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAGTGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCA  >  W3110S.gb/2021211‑2021348
                                                                      |                                                                   
ttttGAGTGCACCGCTGGGTTCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGc                                                                     <  1:1661460/71‑1 (MQ=255)
  ttGAGTGCACCGCTGGGTTCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGc                                                                     >  1:561195/1‑69 (MQ=255)
  ttGAGTGCACCGCTGGGTTCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGc                                                                     >  1:4826155/1‑69 (MQ=255)
  ttGAGTGCACCGCTGGGTTCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGc                                                                     >  1:4168907/1‑69 (MQ=255)
    gAGTGCACCGCTGGGTTCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAg                                                                  <  1:256182/70‑1 (MQ=255)
    gAGTGCACCGCTGGGTTCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAg                                                                  <  1:4953119/70‑1 (MQ=255)
           ccGCTGGGTTCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAg                                                                  >  1:3977129/1‑63 (MQ=255)
           ccGCTGGGTTCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAg                                                                  >  1:1804495/1‑63 (MQ=255)
           ccGCTGGGTTCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGACAGGGGcaaca                                                          <  1:3843056/71‑1 (MQ=255)
            cGCTGGGTTCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGcaacaa                                                         <  1:3932515/71‑1 (MQ=255)
            cGCTGGGTTCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGACAGGGGcaaca                                                          >  1:206433/1‑70 (MQ=255)
            cGCTGGGTTCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGACAGGGGcaac                                                           >  1:1403391/1‑69 (MQ=255)
            cGCTGGGTTCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACACGCCAGGGGcaaca                                                          >  1:2981722/1‑70 (MQ=255)
             gCTGGGTTCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGCGGCAACaaa                                                        <  1:762156/71‑1 (MQ=255)
              cTGGGTTCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACaaa                                                        <  1:2319085/70‑1 (MQ=255)
                 ggTTCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGACAGGGGCAACAAAGTg                                                     >  1:3108547/1‑70 (MQ=255)
                 ggTTCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGACAGGGGCAACAAAGTg                                                     >  1:277065/1‑70 (MQ=255)
                 ggTTCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGACAGGGGCAACAAAGTg                                                     >  1:2487668/1‑70 (MQ=255)
                   ttCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAGTg                                                     >  1:3232146/1‑68 (MQ=255)
                   ttCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGACAGGGGCAACAAAGTGCag                                                  <  1:4792362/71‑1 (MQ=255)
                   ttCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGACAGGGGCAACAAAGTGCag                                                  <  1:4025817/71‑1 (MQ=255)
                     ctcACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAGTGCagag                                                >  1:2275644/1‑71 (MQ=255)
                     ctcACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAGTGCagag                                                >  1:520427/1‑71 (MQ=255)
                     ctcACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGACAGGGGCAACAAAGTGCagag                                                >  1:2495781/1‑71 (MQ=255)
                     ctcACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGACAGGGGCAACAAAGTGCagag                                                >  1:3230181/1‑71 (MQ=255)
                      tcACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGACAGGGGCAACAAAGTGCag                                                  <  1:2376532/68‑1 (MQ=255)
                        aCGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGcaac                                                           <  1:2950661/57‑1 (MQ=255)
                        aCGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAGTGCAGAGtt                                              >  1:746307/1‑70 (MQ=255)
                        aCGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAGTGCAGAGt                                               <  1:5126307/69‑1 (MQ=255)
                          gAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGcaacaa                                                         >  1:156304/1‑57 (MQ=255)
                          gAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGcaacaa                                                         >  1:1983287/1‑57 (MQ=255)
                          gAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGACAGGGGCAACAAAGTGCAGAGTTGCg                                           <  1:2846381/71‑1 (MQ=255)
                                  acaaTCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAGTGCAGAGTTGCGTCTGCAcc                                   >  1:1090187/1‑71 (MQ=255)
                                  acaaTCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAGTGCAGAGTTGCGTCTGCAcc                                   >  1:2083001/1‑71 (MQ=255)
                                  acaaTCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAGTGCAGAGTTGCGTCTGCAcc                                   >  1:1101366/1‑71 (MQ=255)
                                   caaTCATTATGCCAGCATATTTCGCTGTTCACCCGCCAGGGGcaaca                                                          <  1:3735534/47‑1 (MQ=255)
                                   caaTCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGcaaca                                                          <  1:2613132/47‑1 (MQ=255)
                                     aTCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAGTGCAGAGTTGCGTCTGCACCCg                                 >  1:53781/1‑70 (MQ=255)
                                     aTCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAGTGCAGAGTTGCGTCTGCACCCg                                 >  1:1448836/1‑70 (MQ=255)
                                     aTCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAGTGCAGAGTTGCGTCTGCACCCg                                 >  1:698443/1‑70 (MQ=255)
                                      tCATTAAGCCAGCATATTTAGCTGTTCACCCGACAGGGGcaaca                                                          <  1:4847527/44‑1 (MQ=255)
                                        aTTAAGCCAGCATATTTCGCTGTTCACCCGCCAgg                                                                 >  1:3459010/1‑35 (MQ=255)
                                        aTTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACaaa                                                        >  1:4831484/1‑44 (MQ=255)
                                        aTTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACaaa                                                        >  1:1111056/1‑44 (MQ=255)
                                        aTTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAGTGCAGAGTTGCGt                                          >  1:4887262/1‑58 (MQ=255)
                                        aTTAAGCCAGCATATTTCGCTGTTCACCCGACAGGGGCAACAAAGTg                                                     <  1:3299722/47‑1 (MQ=255)
                                          tAAGCCAGCATATTTCGCTGTTCACCCGACAGGGGCAACtaa                                                        >  1:2089546/1‑42 (MQ=38)
                                          tAAGCCAGCATATTTCGCTGTTCACCCGACAGGGGCAACaaa                                                        >  1:2606780/1‑42 (MQ=255)
                                                aGCATATTTCGCTGTTCACCCGACAGGGGCAACAAAGTGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGt                     <  1:308384/71‑1 (MQ=255)
                                                 gCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAGTGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGt                     >  1:3637427/1‑70 (MQ=255)
                                                  cATATTTCGCTGTTCACCCGCCAGGGGCAACAAAGTGCaga                                                 >  1:320027/1‑41 (MQ=255)
                                                   atatTTCGCTGTTCACCCGCCAGGGGCAACAAAGTGCAGAGTTGCGt                                          >  1:3042790/1‑47 (MQ=255)
                                                    tatTTCGCTGTTCACCCGACAGGGGCAACAAAGTGCa                                                   <  1:1968632/37‑1 (MQ=255)
                                                    tatTTCGCTGTTCACCCGACAGGGGCAACAAAGTGCa                                                   <  1:3266119/37‑1 (MQ=255)
                                                    tatTTCGCTGTTCACCCGACAGGGGCAACAAAGTGCa                                                   <  1:2842877/37‑1 (MQ=255)
                                                    tatTTCGCTGTTCACCCGACAGGGGCAACAAAGTGCa                                                   <  1:1947326/37‑1 (MQ=255)
                                                    tatTTCGCTGTTCACCCGACAGGGGCAACAAAGTGCAGAGTTGCGt                                          <  1:2062904/46‑1 (MQ=255)
                                                    tatTTCGCTGTTCACCCGACAGGGGCAACAAAGTGCAGAGTTGCGTCTGCACCCGc                                <  1:1034013/56‑1 (MQ=255)
                                                    tatTTCGCTGTTCACCCGACAGGGGCAACAAAGTGCAGAGTTGCGTCTGCACCCGc                                <  1:1956320/56‑1 (MQ=255)
                                                      tttCGCTGTTCACCCGCCAGGGGCAACAAAGTGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTg               >  1:1807731/1‑71 (MQ=255)
                                                      tttCGCTGTTCACCCGACAGGGGCAACAAAGTGCAGAGTTGCGTCt                                        >  1:3433849/1‑46 (MQ=255)
                                                      tttCGCTGTTCACCCGACAGGGGCAACAAAGTGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGt                >  1:202760/1‑70 (MQ=255)
                                                       ttCGCTGTTCACCCGCCAGGGGCAACAAAGTTCAGAGTTGCGTCTGc                                      <  1:3065202/47‑1 (MQ=255)
                                                          gCTGTTCACCCGCCAGGGGCAACAAAGTGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCaaa           >  1:5061820/1‑71 (MQ=255)
                                                                cACCCGACAGGTGCAACAAAGTGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTccc     >  1:4190731/1‑71 (MQ=255)
                                                                    cGCCAGGGGCAACAAAGTGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCa  >  1:2345650/1‑70 (MQ=255)
                                                                      |                                                                   
TTTTGAGTGCACCGCTGGGTTCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAGTGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCA  >  W3110S.gb/2021211‑2021348

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