Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,965,064 G→A 100% I567I (ATC→ATT flhA ← predicted flagellar export pore protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,965,0640GA88.2% 88.4 / 5.0 34I567I (ATC→ATTflhApredicted flagellar export pore protein
Reads supporting (aligned to +/- strand):  ref base G (0/4);  new base A (18/12);  total (18/16)
Fisher's exact test for biased strand distribution p-value = 3.92e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

ACTGCTGGGTAATCGCCCGTCCCAACGCCACGCGCACGACGGCGGTTAATTCATGTGGATCGCTTTGGATGGGCGCATGTTCCGCCAGCGTTTCGAGAATGGTGCGCATATCGCGAATCGGTACTTTTTCA  >  W3110S.gb/1964997‑1965127
                                                                   |                                                               
aCTGCTGGGTAATCGCCCGTCCCAACGCCACGCGCACGACGGCGGTTAATTCATGTGGATCGCTTTGAAt                                                               >  1:1277181/1‑70 (MQ=255)
aCTGCTGGGTAATCGCCCGTCCCAACGCCACGCGCACGACGGCGGTTAATTCATGTGGATCGCTTTGAAt                                                               >  1:1402967/1‑70 (MQ=255)
aCTGCTGGGTAATCGCCCGTCCCAACGCCACGCGCACGACGGCGGTTAATTCATGTGGATCGCTTTGAAt                                                               >  1:445676/1‑70 (MQ=255)
    ctgGGTAATCGCCCGTCCCAACGCCACGCGCACGACGGCGGTTAATTCATGTGGATCGCTTTGAATGGgcg                                                          <  1:3071264/71‑1 (MQ=255)
      gggTAATCGCCCGTCCCAACGCCACGCGCACGACGGCGGTTAATTCATGTGGATCGCTTTGAATGGGCGCa                                                        >  1:1897039/1‑71 (MQ=255)
           aTCGCCCGTCCCAACGCCACGCGCACGACGGCGGTTAATTCATGTGGATCGCTTTGAATGGGCGCATGTTc                                                   >  1:629924/1‑71 (MQ=255)
           aTCGCCCGTCCCAACGCCACGCGCACGACGGCGGTTAATTCATGTGGATCGCTTTGAATGGGCGCATGTTc                                                   >  1:2102935/1‑71 (MQ=255)
            tCGCCCGTCCCAACGCCACGCGCACGACGGCGGTTAATTCATGTGGATCGCTTTGAATGGGCGCATGTTcc                                                  <  1:764882/71‑1 (MQ=255)
            tCGCCCGTCCCAACGCCACGCGCACGACGGCGGTTAATTCATGTGGATCGCTTTGAATGGGCGCATGTTcc                                                  <  1:1818594/71‑1 (MQ=255)
            tCGCCCGTCCCAACGCCACGCGCACGACGGCGGTTAATTCATGTGGATCGCTTTGAATGGGCGCATGTTcc                                                  <  1:3337125/71‑1 (MQ=255)
                     ccAACGCCACGCGCACGACGGCGGTTAATTCATGTGGATCGCTTTGAATGGGCGCATGTTCCGCCAGCGtt                                         >  1:856012/1‑71 (MQ=255)
                     ccAACGCCACGCGCACGACGGCGGTTAATTCATGTGGATCGCTTTGAATGGGCGCATGTTCCGCCAGCGtt                                         >  1:2961131/1‑71 (MQ=255)
                     ccAACGCCACGCGCACGACGGCGGTTAATTCATGTGGATCGCTTTGAATGGGCGCATGTTCCGCCAGCGtt                                         >  1:2241232/1‑71 (MQ=255)
                       aaCGCCACGCGCACGACGGCGGTTAATTCATGTGGATCGCTTTGAATGGGCGCATGTTCCGCCAGCGTTTc                                       >  1:2307447/1‑71 (MQ=255)
                       aaCGCCACGCGCACGACGGCGGTTAATTCATGTGGATCGCTTTGAATGGGCGCATGTTCCGCCAGCGTTTc                                       >  1:1482625/1‑71 (MQ=255)
                       aaCGCCACGCGCACGACGGCGGTTAATTCATGTGGATCGCTTTGAATGGGCGCATGTTCCGCCAGCGTTTc                                       >  1:59930/1‑71 (MQ=255)
                       aaCGCCACGCGCACGACGGCGGTTAATTCATGTGGATCGCTTTGAATGGGCGCATGTTCCGCCAGCGTTTc                                       >  1:1192800/1‑71 (MQ=255)
                          gCCACGCGCACGACGGCGGTTAATTCATGTGGATCGCTTTGGATGGGCGCATGTTCCGCCAGCGTTTCgag                                    <  1:2352381/71‑1 (MQ=255)
                          gCCACGCGCACGACGGCGGTTAATTCATGTGGATCGCTTTGGATGGGCGCATGTTCCGCCAGCGTTTCgag                                    <  1:2098097/71‑1 (MQ=255)
                          gCCACGCGCACGACGGCGGTTAATTCATGTGGATCGCTTTGAATGGGCGCATGTTCCGCCAGCGTTTc                                       >  1:1814681/1‑68 (MQ=255)
                          gCCACGCGCACGACGGCGGTTAATTCATGTGGATCGCTTTGAATGGGCGCATGTTCCGCCAGCGTTTCgag                                    <  1:2629704/71‑1 (MQ=255)
                            gaCGCGCACGACGGCGGTTAATTCATGTGGATCGCTTTGAATGGGCGCATGTTCCGCCAGCGTTTCGAGaa                                  <  1:745355/70‑1 (MQ=255)
                            cACGCGCACGACGGCGGTTAATTCATGTGGATCGCTTTGAATGGGCGCATGTTCCGCCAGCGTTTCGAGaa                                  <  1:26980/71‑1 (MQ=255)
                            cACGCGCACGACGGCGGTTAATTCATGTGGATCGCTTTGAATGGGCGCATGTTCCGCCAGCGTTTCGAGaa                                  <  1:72860/71‑1 (MQ=255)
                            cACGCGCACGACGGCGGTTAATTCATGTGGATCGCTTTGAATGGGCGCATGTTCCGCCAGCGTTTCGAGaa                                  <  1:627010/71‑1 (MQ=255)
                                cgcACGACGGCGGTTAATTCATGTGGATCGCTTTGAATGGGCGCATGTTCCGCCAGCGttt                                        >  1:3408901/1‑61 (MQ=255)
                                cgcACGACGGCGGTTAATTCATGTGGATCGCTTTGAATGGGCGCATGTTCCGCCAGCGttt                                        >  1:2078794/1‑61 (MQ=255)
                                            gTTAATTCATGTGGATCGCTTTGGATGGGCGCATGTTCCGCCAGCGTTTCGAGAATGGTGCGCatat                      <  1:1450685/67‑1 (MQ=255)
                                            gTTAATTCATGTGGATCGCTTTGGATGGGCGCATGTTCCGCCAGCGTTTCGAGAATGGTGCGCATATcgcg                  <  1:1222826/71‑1 (MQ=255)
                                             ttAATTCATGTGGATCGCTTTGAATGGGCGCATGTTCCGCCAGCGTTTCGAGAATGGTgcg                           <  1:3506630/61‑1 (MQ=255)
                                                   cATGTGGATCGCTTTGAATGGGCGCATGTTCCGCCAGCGTTTCGAGAATGGTGCGCATATCGCGAATCgg            >  1:3442114/1‑70 (MQ=255)
                                                   cATGTGGATCGCTTTGAATGGGCGCATGTTCCGCCAGCGTTTCGAGAATGGTGCGCATATCGCGAATCgg            <  1:2744486/70‑1 (MQ=255)
                                                   cATGTGGATCGCTTTGAATGGGCGCATGTTCCGCCAGCGTTTCGAGAATGGTGCGCATATCGCGAATCgg            <  1:2740067/70‑1 (MQ=255)
                                                            cGCTTTGAATGGGCGCATGTTCCGCCAGCGTTTCGAGAATGGTGCGCATATCGCGAATCGGTACTTTTTCa  >  1:327671/1‑71 (MQ=255)
                                                                   |                                                               
ACTGCTGGGTAATCGCCCGTCCCAACGCCACGCGCACGACGGCGGTTAATTCATGTGGATCGCTTTGGATGGGCGCATGTTCCGCCAGCGTTTCGAGAATGGTGCGCATATCGCGAATCGGTACTTTTTCA  >  W3110S.gb/1964997‑1965127

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