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,0640GA83.3% 69.8 / 8.3 30I567I (ATC→ATTflhApredicted flagellar export pore protein
Reads supporting (aligned to +/- strand):  ref base G (5/0);  new base A (16/9);  total (21/9)
Fisher's exact test for biased strand distribution p-value = 2.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GTTAATTCATGTGGATCGCTTTGGATGGGCGCATGTTCCGCCAGCGTTTCGAGAATGGTGCGCATATCGCGAATCGGTACTT  >  W3110S.gb/1965041‑1965122
                       |                                                          
gTTAATTCATGTGGATCGCTTTGGATGGGCGCATGTTCCGCCAGCGTTTCGAGAATGGTGCGCATATcgc              >  1:953446/1‑70 (MQ=255)
 ttAATTCATGTGGATCGCTTTGAATGGGCGCATGTTCCGCCAGCGTTTCGAGAATGGTgcg                      <  1:1090011/61‑1 (MQ=255)
   aaTTCATGTGGATCGCTTTGAATGGGCGCATGTTCCGCCAGCGTTTCGAGAATGGTGCGCa                    >  1:372609/1‑61 (MQ=255)
   aaTTCATGTGGATCGCTTTGAATGGGCGCATGTTCCGCCAGCGTTTCGAGAATGGTGCGCa                    >  1:888346/1‑61 (MQ=255)
       cATGTGGATCGCTTTGAATGGGCGCATGTTCCGCCAGCGTTTCgag                               >  1:297555/1‑46 (MQ=255)
       cATGTGGATCGCTTTGAATGGGCGCATGTTCCGCCAGCGTTTCgag                               >  1:1683325/1‑46 (MQ=255)
       cATGTGGATCGCTTTGAATGGGCGCATGTTCCGCCAGCGTTTCgag                               >  1:629492/1‑46 (MQ=255)
       cATGTGGATCGCTTTGAATGGGCGCATGTTCCGCCAGCGTTTCgag                               >  1:1621184/1‑46 (MQ=255)
       cATGTGGATCGCTTTGAATGGGCGCATGTTCCGCCAGCGTTTCgag                               >  1:158529/1‑46 (MQ=255)
       cATGTGGATCGCTTTGAATGGGCGCATGTTCCGCCAGCGTTTCgag                               >  1:1322088/1‑46 (MQ=255)
         tgtgGATCGCTTTGAATGGGCGCATGTTCCGCCAGCGTTTCGAGAATGGTGCGCata                  <  1:713746/57‑1 (MQ=255)
         tgtgGATCGCTTTGAATGGGCGCATGTTCCGCCAGCGTTTCGAGAATGGTGCGCata                  <  1:757806/57‑1 (MQ=255)
         tgtgGATCGCTTTGAATGGGCGCATGTTCCGCCAGCGTTTCGAGAATGGTGCGCata                  <  1:1466632/57‑1 (MQ=255)
         tgtgGATCGCTTTGAATGGGCGCATGTTCCGCCAGCGTTTCGAGAATGGTGCGCata                  <  1:1763770/57‑1 (MQ=255)
         tgtgGATCGCTTTGAATGGGCGCATGTTCCGCCAGCGTTTCGAGAATGGTGCGCata                  <  1:624570/57‑1 (MQ=255)
         tgtgGATCGCTTTGAATGGGCGCATGTTCCGCCAGCGTTTCGAGAATGGTGCGCata                  <  1:1801440/57‑1 (MQ=255)
         tgtgGATCGCTTTGAATGGGCGCATGTTCCGCCAGCGTTTCGAGAATGGTGCGCata                  <  1:546653/57‑1 (MQ=255)
         tgtgGATCGCTTTGAATGGGCGCATGTTCCGCCAGCGTTTCGAGAATGGTGCGCata                  <  1:1952315/57‑1 (MQ=255)
           tgGATCGCTTTGGATGGGCGCATGTTCCGCCAGCGTTTCGAGAATGGTGCGCATATCGCGAATCGGTACtt  >  1:1160716/1‑71 (MQ=255)
           tgGATCGCTTTGGATGGGCGCATGTTCCGCCAGCGTTTCGAGAATGGTGCGCATATCGCGAATCGGTACtt  >  1:1675230/1‑71 (MQ=255)
           tgGATCGCTTTGGATGGGCGCATGTTCCGCCAGCGTTTCGAGAATGGTGCGCATATCGCGAATCGGTACtt  >  1:1568759/1‑71 (MQ=255)
           tgGATCGCTTTGGATGGGCGCATGTTCCGCCAGCGTTTCGAGAATGGTGCGCATATCGCGAATCGGTACtt  >  1:1422149/1‑71 (MQ=255)
           tgGATCGCTTTGAATGGGCGCATGTTCCGCCAGCGTTTCGAGAATGGTGCGCATATc                >  1:259383/1‑57 (MQ=255)
           tgGATCGCTTTGAATGGGCGCATGTTCCGCCAGCGTTTCGAGAATGGTGCGCATATc                >  1:533435/1‑57 (MQ=255)
           tgGATCGCTTTGAATGGGCGCATGTTCCGCCAGCGTTTCGAGAATGGTGCGCATATc                >  1:1935412/1‑57 (MQ=255)
           tgGATCGCTTTGAATGGGCGCATGTTCCGCCAGCGTTTCGAGAATGGTGCGCATATc                >  1:1782383/1‑57 (MQ=255)
           tgGATCGCTTTGAATGGGCGCATGTTCCGCCAGCGTTTCGAGAATGGTGCGCATATc                >  1:1499129/1‑57 (MQ=255)
           tgGATCGCTTTGAATGGGCGCATGTTCCGCCAGCGTTTCGAGAATGGTGCGCATATc                >  1:924393/1‑57 (MQ=255)
           tgGATCGCTTTGAATGGGCGCATGTTCCGCCAGCGTTTCGAGAATGGTGCGCATATc                >  1:1227620/1‑57 (MQ=255)
           tgGATCACTTTGAATGGGCGCATGTTCCGCCAGCGATTCGAGAATGGTGCGCATATc                >  1:252108/1‑57 (MQ=255)
                       |                                                          
GTTAATTCATGTGGATCGCTTTGGATGGGCGCATGTTCCGCCAGCGTTTCGAGAATGGTGCGCATATCGCGAATCGGTACTT  >  W3110S.gb/1965041‑1965122

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: