Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,169,960:1 +T 100% intergenic (‑3/+166) yegZ ← / ← yegR ECK2080:JW5886:b2083; hypothetical protein fragment/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,169,9601.T94.7% 50.9 / ‑2.6 19intergenic (‑3/+166)yegZ/yegRECK2080:JW5886:b2083; hypothetical protein fragment/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base T (11/7);  total (11/8)
Fisher's exact test for biased strand distribution p-value = 4.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

CGGTGTTTCCGTGTAAATATCTGCCCGACCGGTAGCCAGGCTGATGGAGAACTCTGCTACACGGGCTTATTTTTTTATGCATAAGCCCTATCTCTGGTAACCGTCTTCCATTGACCACATCGATAGAA  >  W3110S.gb/2169892‑2170019
                                                                     |                                                           
cGGTGTTTCCGTGTAAATATCTGCCCGACCGGTAGCCAGGCTGATGGAGAACTCTGCTACACGGGCTTAtt                                                           <  1:140773/71‑1 (MQ=255)
     tttCCGTGTAAATATCTGCCCGACCGGTAGCCAGGCTGATGGAGAACTCTGCTACACGGGCTTAttttttt                                                      <  1:3028398/71‑1 (MQ=255)
             tAAATATCTGCCCGACCGGTAGCCAGGCTGATGGAGAACTCTGCTACACGGGCTTATTTTTTTTATGCAt                                                >  1:1774842/1‑70 (MQ=255)
               aaTATCTGCCCGACCGGTAGCCAGGCTGATGGAGATCTCTGCTACACGGGCTTATTTTTTTTATGCATAAg                                             >  1:1748486/1‑71 (MQ=255)
                 tatCTGCCCGACCGGTAGCCAGGCTGATGGAGAACTCTGCTACACGGGCTTATTTTTTTTATGCAt                                                <  1:1681462/66‑1 (MQ=255)
                 tatCTGCCCGACCGGTAGCCAGGCTGATGGAGAACTCTGCTACACGGGCTTATTTTTTTTATGCAt                                                <  1:2769493/66‑1 (MQ=255)
                  atCTGCCCGACCGGTAGCCAGGCTGATGGAGAACTCTGCTACACGGGCTTATTTTTTTTATGCATAAGccc                                          >  1:892890/1‑71 (MQ=255)
                   tCTGCCCGACCGGTAGCCAGGCTGATGGAGAACTCTGCTACACGGGCTTAttttttt                                                      <  1:530505/57‑1 (MQ=255)
                   tCTGCCCGACCGGTAGCCAGGCTGATGGAGAACTCTGCTACACGGGCTTAttttttt                                                      <  1:442489/57‑1 (MQ=255)
                                      ggCTGATGGAGAACTCTGCTACACGGGCTTATTTTTTTTATGCATAAGCCCTATCTCTGGTAACCGt                          >  1:1080574/1‑67 (MQ=255)
                                      ggCTGATGGAGAACTCTGCTACACGGGCTTATTTTTTTTATGCATAAGCCCTATCTCTGGTAACCGt                          >  1:2544499/1‑67 (MQ=255)
                                      ggCTGATGGAGAACTCTGCTACACGGGCTTATTTTTTTTATGCATAAGCCCTATCTCTGGTAACCGTCTTc                      <  1:3046729/71‑1 (MQ=255)
                                      ggCTGATGGAGAACTCTGCTACACGGGCTTATTTTTTTTATGCATAAGCCCTATCTCTGGTAACCGTCTTc                      <  1:471145/71‑1 (MQ=255)
                                       gCTGATGGAGAACTCTGCTACACGGGCTTATTTTTTTTATGCa                                                 >  1:98573/1‑43 (MQ=37)
                                       gCTGATGGAGAACTCTGCTACACGGGCTTATTTTTTTTATGCATAAGCCCTATCTCTGGTAACCGTCtt                       <  1:3536891/69‑1 (MQ=255)
                                              gagaACTCTGCTACACGGGCTTATTTTTTTTATGCATAAGCCCTATCTCTGGTAACCGTCTTCCATTGa                >  1:414659/1‑69 (MQ=255)
                                                gaACTCTGCTACACGGGCTTATTTTTTTTATGCATAAGCCCTATCTCTGGTAACCGTCTTCCATTGACca             >  1:3258556/1‑70 (MQ=255)
                                                   ctctGCTACACGGGCTTATTTTTTTTATGCATAAGCCCTATCTCTGGTAACAGTCTTCCATTGACCACATc         <  1:2607741/71‑1 (MQ=255)
                                                      tGCTACACGGGCTTATTTTTTTTATGCATAAGCCCTATCTCTGGTAACCGTCTTCCATTGACCACATCGa       >  1:3477443/1‑70 (MQ=255)
                                                      tGCTACACGGGCTTATTTTTTTTATGCATAAGCCCTATCTCTGGTAACCGTCTTCCATTGACCACATCGa       >  1:1080523/1‑70 (MQ=255)
                                                      tGCTACACGGGCTTATTTTTTTTATGCATAAGCCCTATCTCTGGTAACCGTCTTCCATTGACCACATCGa       >  1:1750695/1‑70 (MQ=255)
                                                           cacGGGCTTATTTTTTTATGCATAAGCCCTATCTCTGGTAACCGTCTTCCATTGACCACATCGATAGaa  <  1:3497366/69‑1 (MQ=255)
                                                             cGGGCTTATTTTTTTTATGCATAAGCCCTATCTCTGGTAACCGTCTTCCAtt                  <  1:482169/52‑1 (MQ=255)
                                                                     |                                                           
CGGTGTTTCCGTGTAAATATCTGCCCGACCGGTAGCCAGGCTGATGGAGAACTCTGCTACACGGGCTTATTTTTTTATGCATAAGCCCTATCTCTGGTAACCGTCTTCCATTGACCACATCGATAGAA  >  W3110S.gb/2169892‑2170019

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: