Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,124,493 T→C 57.1% intergenic (+95/‑12) damX → / → dam hypothetical protein/DNA adenine methylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,124,4930TC57.1% 5.7 / 21.8 21intergenic (+95/‑12)damX/damhypothetical protein/DNA adenine methylase
Reads supporting (aligned to +/- strand):  ref base T (6/3);  new base C (5/7);  total (11/10)
Fisher's exact test for biased strand distribution p-value = 3.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGTTGCGGGAGTTGCCTGAAGCGCTGGATGCTGTCGGAGCTTTCTCCACAGCCGGAGAAGGTGTAATTAGTTAGTCAGCATGAAGAAAAATCGCGCTTTTTTGAAGTGGGCAGGGGGCAAGTATCCCCTGCTTGATGA  >  W3110S.gb/4124425‑4124563
                                                                    |                                                                      
tCGTTGCGGGAGTTGCCTGAAGCGCTGGATGCTGTCGGAGCTTTCTCCACAGCCGGAGAAGGTGTAAtcag                                                                      <  1:1695378/71‑1 (MQ=255)
     gCGGGAGTTGCCTGAAGCGCTGGATGCTGTCGGAGCTTTCTCCACAGCCGGAGAAGGTGTAATcagttagt                                                                 >  1:3348119/1‑71 (MQ=255)
       gggAGTTGCCTGAAGCGCTGGATGCTGTCGGAGCTTTCTCCACAGCCGGAGAAGGTGTAAtcag                                                                      <  1:336732/64‑1 (MQ=255)
              gCCTGAAGCGCTGGATGCTGTCGGAGCTTTCTCCACAGCCGGAGAAGGTGTAATTAGTTAGTCAGCATg                                                          >  1:589956/1‑69 (MQ=255)
              gCCTGAAGCGCTGGATGCTGTCGGAGCTTTCTCCACAGCCGGAGAAGGTGTAATTAGTTAGTCAGCATg                                                          >  1:2078155/1‑69 (MQ=255)
              gCCTGAAGCGCTGGATGCTGTCGGAGCTTTCTCCACAGCCGGAGAAGGTGTAATTAGTTAGTCAGCATg                                                          >  1:251922/1‑69 (MQ=255)
                    aGCGCTGGATGCTGTCGGAGCTTTCTCCACAGCCGGAGAAGGTGTAATCAGTTAGTCAGCATGAAGaaa                                                    <  1:347440/69‑1 (MQ=255)
                        cTGGATGCTGTCGGAGCTTTCTCCACAGCCGGAGAAGGTGTAATTAGTTAGTCAGCATGAAGAAAAATcgc                                              >  1:381775/1‑71 (MQ=255)
                           gATGCTGTCGGAGCTTTCTCCACAGCCGGAGAAGGTGTAATTAGTTAGTCAGCATGAAGAAAAATCGCGCt                                           >  1:1183451/1‑71 (MQ=255)
                           gATGCTGTCGGAGCTTTCTCCACAGCCGGAGAAGGTGTAATTAGTTAGTCAGCATGAAGAAAAATCGCGCt                                           >  1:904627/1‑71 (MQ=255)
                                  tCGGAGCTTTCTCCACAGCCGGAGAAGGTGTAATTAGTTAGTCAGCATGAAGAAAAATCGCGCTTTTTTg                                     <  1:2504138/70‑1 (MQ=255)
                                  tCGGAGCTTTCTCCACAGCCGGAGAAGGTGTAATCAGTTAGTCAGCATGAAGAAAAATCGCGCTTTTTTGa                                    <  1:2943158/71‑1 (MQ=255)
                                    ggAGCTTTCTCCACAGCCGGAGAAGGTGTAATTAGTTAGTCAGCATGAAGAAAAATCGCGCTTTTTTGaa                                   <  1:2073009/70‑1 (MQ=255)
                                        cTTTCTCCACAGCCGGAGAAGGTGTAATTAGTTAGTCAGCATGAAGAAAAATCGCGCTTTTTTGAAGTggg                              <  1:2104974/71‑1 (MQ=255)
                                        cTTTCTCCACAGCCGGAGAAGGTGTAATCAGTTAGTCAg                                                              <  1:762684/39‑1 (MQ=255)
                                        cTTTCTCCACAGCCGGAGAAGGTGTAATCAGTTAGTCAg                                                              <  1:1749505/39‑1 (MQ=255)
                                        cTTTCTCCACAGCCGGAGAAGGTGTAATCAGTTAGTCAg                                                              <  1:1778263/39‑1 (MQ=255)
                                           tctcCACAGCCGGAGAAGGTGTAATCAGTTAGTCAGCATGAAGAAAAATCGCGCTTTTTTGAAGTGGGCAg                           >  1:2468937/1‑71 (MQ=255)
                                           tctcCACAGCCGGAGAAGGTGTAATCAGTTAGTCAGCATGAAGAAAAATCGCGCTTTTTTGAAGTGGGCAg                           >  1:2583006/1‑71 (MQ=255)
                                           tctcCACAGCCGGAGAAGGTGTAATCAGTTAGTCAGCATGAAGAAAAATCGCGCTTTTTTGAAGTGGGCAg                           >  1:3595027/1‑71 (MQ=255)
                                           tctcCACAGCCGGAGAAGGTGTAATCAGTTAGTCAGCATGAAGAAAAATCGCGCTTTTTTGAAGTGGGCAg                           >  1:2106133/1‑71 (MQ=255)
                                                        agaAGGTGTAATCAGTTAGTCAGCATGAAGAAAAATCGCGCTTTTTTGAAGTGGGCAGGGGGCAAGTATcc              <  1:1104139/71‑1 (MQ=255)
                                                                tAATCAGTTAGTCAGCATGAAGAAAAATCGCGCttttt                                       >  1:2380679/1‑38 (MQ=255)
                                                                tAATCAGTTAGTCAGCATGAAGAAAAATCGCGCTTTTTTGAAGTg                                <  1:2143592/45‑1 (MQ=255)
                                                                   ttagttagTCAGCATGAAGAAAAATCGCGCTTTTTTGa                                    <  1:2956028/38‑1 (MQ=255)
                                                                   ttagttagTCAGCATGAAGAAAAATCGCGCTTTTTTGa                                    <  1:2891856/38‑1 (MQ=255)
                                                                   ttagttagTCAGCATGAAGAAAAATCGCGCTTTTTTGa                                    <  1:446754/38‑1 (MQ=255)
                                                                   ttagttagTCAGCATGAAGAAAAATCGCGCTTTTTTGa                                    <  1:2501129/38‑1 (MQ=255)
                                                                    tagttagtCAGCATGAAGAAAAATCGCGCTTTTTTGAAGTGGGCAGGGGGCAAGTATCCCCTGCTtgatga  >  1:3202488/1‑71 (MQ=255)
                                                                    tagttagtCAGCATGAAGAAAAATCGCGCTTTTTTGAAGTGGGCAGGGGGCAAGTATCCCCTGCTtgatga  >  1:2730475/1‑71 (MQ=255)
                                                                    tagttagtCAGCATGAAGAAAAATCGCGCTTTTTTGAAGTGGGCAGGGGGCAAGTATCCCCTGCTtgatga  >  1:2726429/1‑71 (MQ=255)
                                                                    |                                                                      
TCGTTGCGGGAGTTGCCTGAAGCGCTGGATGCTGTCGGAGCTTTCTCCACAGCCGGAGAAGGTGTAATTAGTTAGTCAGCATGAAGAAAAATCGCGCTTTTTTGAAGTGGGCAGGGGGCAAGTATCCCCTGCTTGATGA  >  W3110S.gb/4124425‑4124563

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: