Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 2,169,967 (T)7→8 intergenic (‑10/+159) 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.T87.0% 53.8 / 5.0 23intergenic (‑3/+166)yegZ/yegRECK2080:JW5886:b2083; hypothetical protein fragment/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/3);  new base T (3/17);  total (3/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

CGGTGTTTCCGTGTAAATATCTGCCCGACCGGTAGCCAGGCTGATGGAGAACTCTGCTACACGGGCTTATTTTTTTATGCATAAGCCCTATCTCTGGTAACCGTCTTCCATTGACCACATCGATAGAATCCTC  >  W3110S.gb/2169892‑2170024
                                                                     |                                                                
cGGTGTTTCCGTGTAAATATCTGCCCGACCGGTAGCCAGGCTGATGGAGAACTCTGCTACACGGGCTTAtt                                                                <  1:2913509/71‑1 (MQ=255)
cGGTGTTTCCGTGTAAATATCTGCCCGACCGGTAGCCAGGCTGATGGAGAACTCTGCTACACGGGCTTAtt                                                                <  1:2761784/71‑1 (MQ=255)
        ccGTGTAAATATCTGCCCGACCGGTAGCCAGGCTGATGGAGAACTCTGCTACACGGGCTTAtttttttta                                                         <  1:3252477/70‑3 (MQ=255)
        ccGTGTAAATATCTGCCCGACCGGTAGCCAGGCTGATGGAGAACTCTGCTACACGGGCTTAtttttttta                                                         <  1:3935762/70‑3 (MQ=255)
        ccGTGTAAATATCTGCCCGACCGGTAGCCAGGCTGATGGAGAACTCTGCTACACGGGCTTAtttttttta                                                         <  1:3635649/70‑3 (MQ=255)
        ccGTGTAAATATCTGCCCGACCGGTAGCCAGGCTGATGGAGAACTCTGCTACACGGGCTTAtttttttta                                                         <  1:2286137/70‑3 (MQ=255)
        ccGTGTAAATATCTGCCCGACCGGTAGCCAGGCTGATGGAGAACTCTGCTACACGGGCTTAtttttttta                                                         <  1:3279851/70‑3 (MQ=255)
        ccGTGTAAATATCTGCCCGACCGGTAGCCAGGCTGATGGAGAACTCTGCTACACGGGCTTAtttttttta                                                         <  1:2574222/70‑3 (MQ=255)
        ccGTGTAAATATCTGCCCGACCGGTAGCCAGGCTGATGGAGAACTCTGCTACACGGGCTTAtttttttta                                                         <  1:3392163/70‑3 (MQ=255)
                atatCTGCCCGACCGGTAGCCAGGCTGATGGAGAACTCTGCTACACGGGCTTATTTTTTTTATGCATAAGc                                                 <  1:2465469/71‑1 (MQ=255)
                atacCTGCCCGACCGGTAGCCAGGCTGATGGAGAACTCTGCTACACGGGCTTATTTTTTTTATGCATAAGc                                                 <  1:348346/71‑1 (MQ=255)
                         cGACCGGTAGCCAGGCTGATGGAGAACTCTGCTACACGGGCTTATTTTTTTTATGCATAAGCCCTATctct                                        <  1:3222734/71‑1 (MQ=255)
                         cGACCGGTAGCCAGGCTGATGGAGAACTCTGCTACACGGGCTTATTTTTTTTATGCATAAGACCTATctct                                        <  1:2349417/71‑1 (MQ=255)
                               gTAGCCAGGCTGATGGAGAACTCTGCTACACGGGCTTAtt                                                                <  1:4072324/40‑1 (MQ=255)
                                    cAGGCTGATGGAGAACTCTGCTACACGGGCTTAtttttttta                                                         <  1:2422052/42‑3 (MQ=255)
                                      ggCTGATGGAGAACTCTGCTACACGGGCTTATTTTTTTTATGCATAAGCCCTATCTCTGGTAACCGTCTTc                           <  1:2278314/71‑1 (MQ=255)
                                       gCTGATGGAGAACTCTGCTACATGGGCTTATTTTTTTTATGCATAAGCCCTATCTCTGGTAACCGTCTTc                           <  1:2103438/70‑1 (MQ=255)
                                       gCTGATGGAGAACTCTGCTACACGGGCTTATTTTTTTTATGCATAAGCCCTAtctc                                         <  1:1089822/56‑1 (MQ=39)
                                       gCTGATGGAGAACTCTGCTACACGGGCTTATTTTTTTTATGCATAAGCCCTATCTCTGGTAACCGTCTTc                           <  1:2301973/70‑1 (MQ=255)
                                       gCTGATGGAGAACTCTGCTACACGGGCTTATTTTTTTTATGCATAAGCCCTATCTCTGGTAACCGTCTTc                           <  1:2297405/70‑1 (MQ=255)
                                         tGATGGAGAACTCTGCTACACGGGCTTATTTTTTTTATGCATAAGCCCTATCTCTGGTAACCGTCTTc                           <  1:2371890/68‑1 (MQ=255)
                                                gaACTCTGCTACACGGGCTTATTTTTTTTATGCATAAGCCCTATCTCTGGTAACCGTCTTCCATTGACca                  >  1:1130059/1‑70 (MQ=255)
                                                gaACTCTGCTACACGGGCTTATTTTTTTTATGCATAAGCCCTATCTCTGGTAACCGTCTTCCATTGACca                  >  1:2932385/1‑70 (MQ=255)
                                                gaACTCTGCTACACGGGCTTATTTTTTTTATGCATAAGCCCTATCTCTGGTAAACGTCTTCCATTGACca                  >  1:1223750/1‑70 (MQ=255)
                                                 aaCTCTGCTACACGGGCTTATTTTTTTTATGCATAAGCCCTATCTCTGGTAACCGTCTTCCATTGACcaca                <  1:3157136/71‑1 (MQ=255)
                                                 aaCTCTGCTACACGGGCTTATTTTTTTTATGCATAAGCCCTATCTCTGGTAACCGTCTTCCATTGACcaca                <  1:3284782/71‑1 (MQ=255)
                                                 aaCTCTGCTACACGGGCTTATTTTTTTTATGCATAAGCCCTATCTCTGGTAACCGTCTTCCATTGACcaca                <  1:1903071/71‑1 (MQ=255)
                                                 aaCTCTGCTACACGGGCTTATTTTTTTTATGCATAAGCCCTATCTCTGGTAACCGTCTTCCATTGACcaca                <  1:1897885/71‑1 (MQ=255)
                                                 aaCTCTGCTACACGGGCTTATTTTTTTTATGCATAAGCCCTATCTCTGGTAACCGTCTTCCATTGACcaca                <  1:4012355/71‑1 (MQ=255)
                                                   ctctGCTACACGGGCTTATTTTTTTTATGCATAAGCCCTATCTCTGGTAACCGTCTTCCATTGACcaca                <  1:2558095/69‑1 (MQ=255)
                                                      tGCTACACGGGCTTATTTTTTTATGCATAAGCCCTATCTCTGGTAACCGTCTTCCATTGACCACATCGa            <  1:2727927/69‑1 (MQ=255)
                                                      tGCTACACGGGCTTATTTTTTTATGCATAAGCCCTATCTCTGGTAACCGTCTTCCATTGACCACATCGa            <  1:2047789/69‑1 (MQ=255)
                                                      tGCTACACGGGCTTATTTTTTTATGCATAAGCCCTATCTCTGGTAACCGTCTTCCATTGACCACATCGa            <  1:3609492/69‑1 (MQ=255)
                                                               ggCTTATTTTTTTTATGCATAAGCCCTATCTCTGGTAACCGTCTTCCATTGACCACATCGATAGAAtcctc  <  1:3815806/71‑1 (MQ=255)
                                                                     |                                                                
CGGTGTTTCCGTGTAAATATCTGCCCGACCGGTAGCCAGGCTGATGGAGAACTCTGCTACACGGGCTTATTTTTTTATGCATAAGCCCTATCTCTGGTAACCGTCTTCCATTGACCACATCGATAGAATCCTC  >  W3110S.gb/2169892‑2170024

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