Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 1,068,568 T→A G22G (GGT→GGA ymdF → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,068,5680TA98.3% 200.8 / ‑3.9 60G22G (GGT→GGAymdFconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base A (21/38);  total (22/38)
Fisher's exact test for biased strand distribution p-value = 3.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

CAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCAT  >  W3110S.gb/1068507‑1068623
                                                             |                                                       
cAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGaaa                                                       >  1:1599311/1‑64 (MQ=255)
  aaCCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGaaaaaa                                                    >  1:965190/1‑65 (MQ=255)
  aaCCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGaaaaaa                                                    >  1:2067517/1‑65 (MQ=255)
  aaCCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGaaaaaa                                                    >  1:2167428/1‑65 (MQ=255)
  aaCCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGaaaaaa                                                    >  1:1615419/1‑65 (MQ=255)
  aaCCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGaaaaaa                                                    >  1:2909688/1‑65 (MQ=255)
  aaCCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGaaaaaa                                                    >  1:1189929/1‑65 (MQ=255)
   aCCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGaaaaaa                                                    >  1:2732636/1‑64 (MQ=255)
   aCCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGaaaaaa                                                    >  1:989421/1‑64 (MQ=255)
   aCCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGAaaaaaa                                                   >  1:3007761/1‑65 (MQ=255)
   aCCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGAaaaaaa                                                   >  1:2955536/1‑65 (MQ=255)
   aCCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGAaaaaaa                                                   >  1:3059638/1‑65 (MQ=255)
   aCCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGAaaaaaa                                                   >  1:717965/1‑65 (MQ=255)
   aCCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGAaaaaaa                                                   >  1:82197/1‑65 (MQ=255)
   aCCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGAaaaaaa                                                   >  1:2715966/1‑65 (MQ=255)
   aCCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGAaaaaaa                                                   >  1:2675068/1‑65 (MQ=255)
   aCCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGAaaaaaa                                                   >  1:517857/1‑65 (MQ=255)
   aCCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGAaaaaaa                                                   >  1:651913/1‑65 (MQ=255)
   aCCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGAaaaaaa                                                   >  1:2643753/1‑65 (MQ=255)
   aCCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGAaaaaaa                                                   >  1:264031/1‑65 (MQ=255)
   aCCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGAaaaaaa                                                   >  1:2288200/1‑65 (MQ=255)
   aCCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGAaaaaaa                                                   >  1:702719/1‑65 (MQ=255)
   aCCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGAaaaaaa                                                   >  1:2087057/1‑65 (MQ=255)
   aCCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGAaaaaaa                                                   >  1:1954610/1‑65 (MQ=255)
   aCCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGAaaaaaa                                                   >  1:962330/1‑65 (MQ=255)
   aCCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGAaaaaaa                                                   >  1:1935276/1‑65 (MQ=255)
   aCCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGAaaaaaa                                                   >  1:1076453/1‑65 (MQ=255)
   aCCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGAaaaaaa                                                   >  1:1142661/1‑65 (MQ=255)
   aCCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGAaaaaaa                                                   >  1:97930/1‑65 (MQ=255)
   aCCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGAaaaaaa                                                   >  1:1528169/1‑65 (MQ=255)
   aCCATCGAGGCGCTTCCCGCAATATTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGCaaaaaa                                                    >  1:1662653/1‑64 (MQ=255)
                    ggCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGAAAAAAAGGTGGAAGCACAGCgg                                  <  1:1386517/64‑1 (MQ=255)
                     gCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGAAAAAAAGGTGGACAGCACAGCggg                                 <  1:575567/65‑1 (MQ=255)
                     gCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGAAAAAAAGGTGGACAGCACAGCggg                                 <  1:784856/65‑1 (MQ=255)
                     gCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGAAAAAAAGGTGGACAGCACAGCggg                                 <  1:817118/65‑1 (MQ=255)
                     gCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGAAAAAAAGGTGGACAGCACAGCgg                                  <  1:1744709/64‑1 (MQ=255)
                     gCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGAAAAAAAGGTGGACAGCACAGCgg                                  <  1:1085436/64‑1 (MQ=255)
                     gCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGAAAAAAAGGTGGACAGCACAGCgg                                  <  1:1097019/64‑1 (MQ=255)
                     gCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGAAAAAAAGGTGGACAGCACAGCgg                                  <  1:222330/64‑1 (MQ=255)
                     gCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGAAAAAAAGGTGGACAGCACAGCgg                                  <  1:3005422/64‑1 (MQ=255)
                     gCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGAAAAAAAGGTGGACAGCACAGCgg                                  <  1:1218988/64‑1 (MQ=255)
                     gCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGAAAAAAAGGTGGACAGCACAGCgg                                  <  1:2835163/64‑1 (MQ=255)
                     gCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGAAAAAAAGGTGGACAGCACAGCgg                                  <  1:2493647/64‑1 (MQ=255)
                     gCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGAAAAAAAGGTGGACAGCACAGCgg                                  <  1:2656027/64‑1 (MQ=255)
                     gCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGAAAAAAAGGTGGACAGCACAGCgg                                  <  1:878996/64‑1 (MQ=255)
                     gCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGAAAAAAAGGTGGACAGCACAGCgg                                  <  1:269864/64‑1 (MQ=255)
                      cAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGAAAAAAAGGTGGACAGCACAGCggg                                 <  1:2923667/64‑1 (MQ=255)
                      cAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGAAAAAAAGGTGGACAGCACAGCggg                                 <  1:3109072/64‑1 (MQ=255)
                      cAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGAAAAAAAGGTGGACAGCACAGCgg                                  <  1:793783/63‑1 (MQ=255)
                          tttGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTaa                            >  1:955224/1‑65 (MQ=255)
                                          aaaGAGCATCAGAAGCAGGAAAAAAAGGTGGACAGCACAGCgg                                  <  1:1267423/43‑1 (MQ=255)
                                              agCATCAGAAGCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAgc        <  1:2208104/65‑1 (MQ=255)
                                              agCATCAGAAGCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAgc        <  1:1582944/65‑1 (MQ=255)
                                              agCATCAGAAGCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAgc        <  1:479749/65‑1 (MQ=255)
                                              agCATCAGAAGCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAgc        <  1:963730/65‑1 (MQ=255)
                                              agCATCAGAAGCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAgc        <  1:1792341/65‑1 (MQ=255)
                                              agCATCAGAAGCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAgc        <  1:2173899/65‑1 (MQ=255)
                                                   cAGAAGCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCa   <  1:1075101/65‑1 (MQ=255)
                                                   cAGAAGCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCa   <  1:2723499/65‑1 (MQ=255)
                                                   cAGAAGCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCa   <  1:1424451/65‑1 (MQ=255)
                                                    aGAAGCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCAt  <  1:336255/65‑1 (MQ=255)
                                                    aGAAGCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCAt  <  1:3077844/65‑1 (MQ=255)
                                                    aGAAGCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCAt  <  1:2846281/65‑1 (MQ=255)
                                                    aGAAGCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCAt  <  1:2242086/65‑1 (MQ=255)
                                                    aGAAGCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCAt  <  1:1941799/65‑1 (MQ=255)
                                                    aGAAGCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCAt  <  1:1172589/65‑1 (MQ=255)
                                                    aGAAGCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCAt  <  1:1003930/65‑1 (MQ=255)
                                                    aGAAGCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCAt  <  1:994108/65‑1 (MQ=255)
                                                    aGAAGCAGGAAAAAAACGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCAt  <  1:1404101/65‑1 (MQ=255)
                                                     gAAGCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCAt  <  1:1996349/64‑1 (MQ=255)
                                                             |                                                       
CAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCAT  >  W3110S.gb/1068507‑1068623

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: