Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,124,493 T→C 47.6% 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,4930TC47.6% 2.2 / 56.0 42intergenic (+95/‑12)damX/damhypothetical protein/DNA adenine methylase
Reads supporting (aligned to +/- strand):  ref base T (8/14);  new base C (3/17);  total (11/31)
Fisher's exact test for biased strand distribution p-value = 1.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCGGGAGTTGCCTGAAGCGCTGGATGCTGTCGGAGCTTTCTCCACAGCCGGAGAAGGTGTAATTAGTTAGTCAGCATGAAGAAAAATCGCGCTTTTTTGAAGTGGGCAGGGGGCAAGTATCCCCTGC  >  W3110S.gb/4124429‑4124556
                                                                |                                                               
tGCGGGAGTTGCCTGAAGCGCTGGATGCTGTCGGAGCTTTCTCCACAGCCGGAGAAGGTGTAAttagt                                                              <  1:3003051/68‑1 (MQ=255)
tGCGGGAGTTGCCTGAAGCGCTGGATGCTGTCGGAGCTTTCTCCACAGCCGGAGAAGGTGTAAttagt                                                              <  1:2283000/68‑1 (MQ=255)
tGCGGGAGTTGCCTGAAGCGCTGGATGCTGTCGGAGCTTTCTCCACAGCCGGAGAAGGTGTAAtcagt                                                              <  1:683055/68‑1 (MQ=255)
 gCGGGAGTTGCCTGAAGCGCTGGATGCTGTCGGAGCTTTCTCCACAGCCGGAGAAGGTGTAAttagt                                                              <  1:2673190/67‑1 (MQ=255)
          tccTGAAGCGCTGGATGCTGTCGGAGCTTTCTCCACAGCCGTAGAAGGTGTAAtcagt                                                              <  1:840385/57‑1 (MQ=255)
          gCCTGAAGCGCTGGATGCTGTCGGAGCTTTCTCCACAGCCTGAGAAGGTGTAAtcagt                                                              <  1:1728201/58‑1 (MQ=255)
          gCCTGAAGCGCTGGATGCTGTCGGAGCTTTCTCCACAGCCGGAGAAGGTGTAAtcagt                                                              <  1:967978/58‑1 (MQ=255)
          gCCTGAAGCGCTGGATGCTGTCGGAGCTTTCTCCACAGCCGGAGAAGGTGTAAtcagt                                                              <  1:3719652/58‑1 (MQ=255)
          gCCTGAAGCGCTGGATGCTGTCGGAGCTTTCTCCACAGCCGGAGAAGGTGTAAtcagt                                                              >  1:124859/1‑58 (MQ=255)
          gCCTGAAGCGCTGGATGCTGTCGGAGCTTTCTCCACAGCCGGAGAAGGTGTAAtcagt                                                              <  1:320461/58‑1 (MQ=255)
          gCCTGAAGCGCTGGATGCTGTCGGAGCTTTCTCCACAGCCGGAGAAGGTGTAAtcagt                                                              <  1:2257564/58‑1 (MQ=255)
          gCCTGAAGCGCTGGATGCTGTCGGAGCTTTCTCCACAGCCGGAGAAGGTGTAATCAGTTAGTCAGCa                                                     >  1:3069635/1‑67 (MQ=255)
                 gcgcTGGATGCTGTCGGAGCTTTCTCCACAGCCGGAGAAGGTGTAATCAGTTAGTCAGCATGAAGaaaa                                            <  1:843572/69‑1 (MQ=255)
                       gATGCTGTCGGAGCTTTCTCCACAGCCGGAGAAGGTGTAAtcagt                                                              <  1:2568513/45‑1 (MQ=255)
                         tGCTGTCGGAGCTTTCTCCACAGCCGGAGAAGGTGTAATCAGTTAGTCAGCATGAAGAAAAATCGCGCt                                    <  1:811527/69‑1 (MQ=255)
                         tGCTGTCGGAGCTTTCTCCACAGCCGGAGAAGGTGTAATCAGTTAGTCAGCATGAAGAAAAATCGCGCt                                    <  1:2327480/69‑1 (MQ=255)
                               cGGAGCTTTCTCCACAGCCGGAGAAGGTGTAATCAGTTAGTCAGCATGAAGAAAAATTGCGCttt                                  <  1:2626417/65‑1 (MQ=255)
                               cGGAGCTTTCTCCACAGCCGGAGAAGGTGTAATCAGTTAGTCAGCATGAAGAAAAATCGCGCttt                                  <  1:908727/65‑1 (MQ=255)
                               cGGAGCTTTCTCCACAGCCGGAGAAGGTGTAATCAGTTAGTCAGCATGAAGAAAAATCGCGCttt                                  <  1:2775663/65‑1 (MQ=255)
                                tgAGCTTTCTCCACAGCCGGAGAAGGTGTAATTAGTTAGTCAGCATGAAGAAAAATCGCGCTTTTTTg                              <  1:685203/67‑1 (MQ=255)
                                ggAGCTTTCTCCACAGCCGGAGACGGTGTAATTAGTTAGTCAGCATGAAGAAAAATCGCGCTTTTTTGa                             <  1:495391/69‑1 (MQ=255)
                                ggAGCTTTCTCCACAGCCGGAGAAGGTGTAATTAGTTAGTCAGCATGAAGAAAAATCGCGCTTTTTTg                              <  1:2025143/68‑1 (MQ=255)
                                ggAGCTTTCTCCACAGCCGGAGAAGGTGTAATTAGTTAGTCAGCATGAAGAAAAATCGCGCTTTTTTg                              <  1:87298/68‑1 (MQ=255)
                                ggAGCTTTCTCCACAGCCGGAGAAGGTGTAATTAGTTAGTCAGCATGAAGAAAAATCGCGCTTTTTTg                              <  1:623568/68‑1 (MQ=255)
                                ggAGCTTTCTCCACAGCCGGAGAAGGTGTAATTAGTTAGTCAGCATGAAGAAAAATCGCGCTTTTTTg                              <  1:3635947/68‑1 (MQ=255)
                                ggAGCTTTCTCCACAGCCGGAGAAGGTGTAATTAGTTAGTCAGCATGAAGAAAAATCGCGCTTTTTTg                              <  1:3405984/68‑1 (MQ=255)
                                ggAGCTTTCTCCACAGCCGGAGAAGGTGTAATTAGTTAGTCAGCATGAAGAAAAATCGCGCTTTTTTg                              <  1:2644966/68‑1 (MQ=255)
                                ggAGCTTTCTCCACAGCCGGAGAAGGTGTAATCAGTTAGTCAGCATGAAGAAAAATCGCTCTTTTTTGa                             <  1:1630233/69‑1 (MQ=255)
                                ggAGCTTTCTCCACAGCCGGAGAAGGTGTAATCAGTTAGTCAGCATGAAGAAAAATCGCGCttttt                                >  1:399727/1‑66 (MQ=255)
                                ggAGCTTTCTCCACAGCCGGAGAAGGTGTAATCAGTTAGTCAGCATGAAGAAAAATCGCGCTTTTTTg                              <  1:3302632/68‑1 (MQ=255)
                                      ttCTCCACAGCCGGAGAAGGTGTAATTAGTTAGTCAGCATGAAGAAAAATCGCGCt                                    >  1:950212/1‑56 (MQ=255)
                                       tctcCACAGCCGGAGAAGGTGTAATCAGTTAGTCAGCATGAAGAAAAATCGCGCTTTTTTGAAGTGGGc                      >  1:3259147/1‑69 (MQ=255)
                                           cacaGCCGGAGAAGGTGTAATCAGTTAGTCAGCAt                                                    <  1:3102101/35‑1 (MQ=255)
                                           cacaGCCGGAGAAGGTGTAATCAGTTAGTCAGCAt                                                    <  1:250196/35‑1 (MQ=255)
                                             caGCCGGAGAAGGTGTAATTAGTTAGTCAGCATGAAGaaaa                                            <  1:1642623/41‑1 (MQ=255)
                                             caGCCGGAGAAGGTGTAATTAGTTAGTCAGCATGAAGaaaa                                            >  1:2204652/1‑41 (MQ=255)
                                              aGCCGGAGAAGGTGTAATTAGTTAGTCAGCATGAAGAAAAATCGCGCTTTTTTg                              <  1:833694/54‑1 (MQ=255)
                                              aGCCGGAGAAGGTGTAATTAGTTAGTCAGCATGAAGAAAAATCGCGCTTTTTTg                              <  1:852894/54‑1 (MQ=255)
                                              aGCCGGAGAAGGTGTAATTAGTTAGTCAGCATGAAGAAAAATCGCGCTTTTTTg                              <  1:442717/54‑1 (MQ=255)
                                              aGCCGGAGAAGGTGTAATTAGTTAGTCAGCATGAAGAAAAATCGCGCTTTTTTGa                             >  1:267253/1‑55 (MQ=255)
                                              aGCCGGAGAAGGTGTAATTAGTTAGTCAGCATGAAGAAAAATCGCGCTTTTTTGa                             >  1:1263537/1‑55 (MQ=255)
                                                 cGGAGAAGGTGTAATTAGTTAGTCAGCATGAAGAAAAATCGCGCTTTTTTGAAGTGGGCAGGGGGCAAg            <  1:2243290/69‑1 (MQ=255)
                                                  ggAGAAGGTGTAATTAGTTAGTCAGCATGAAGAAAAATCGCGCTTTTTTGAAGTGGGCAGGGGGCAAg            <  1:1170551/68‑1 (MQ=255)
                                                  ggAGAAGGTGTAATCAGTTAGTCAGCATGAAGAAAAATCGCGCTTTTTTGAAGTGGGCAGGGGGCAAGt           <  1:2619111/69‑1 (MQ=255)
                                                  ggAGAAGGTGTAATCAGTTAGTCAGCATGAAGAAAAATCGCGCTTTTTTGAAGTGGGCAGGGGGCAAGt           <  1:1051764/69‑1 (MQ=255)
                                                  ggAGAAGGTGTAATCAGTTAGTCAGCATGAAGAAAAATCGCGCTTTTTTGAAGTGGGCAGGGGGCAAGt           <  1:876234/69‑1 (MQ=255)
                                                      aaGGTGTAATTAGTTAGTCAGCATGAAGAAAAATCGCGCTTTTTTGa                             >  1:1595168/1‑47 (MQ=255)
                                                      aaGGTGTAATTAGTTAGTCAGCATGAAGAAAAATCGCGCTTTTTTGa                             >  1:200191/1‑47 (MQ=255)
                                                          tgtAATCAGTTAGTCAGCATGAAGAAAAATCGCGCTTTTTTGAAGTGGGCAGGGGGCAAGTATCCCCTg   <  1:3596797/69‑1 (MQ=255)
                                                           gtAATTAGTTAGTCAGCATGAAGAAAAATCGCGCTTTTTTGAAGTGGGCAGGGGGCAAGTATCCCCTGc  >  1:1196325/1‑69 (MQ=255)
                                                           gtAATTAGTTAGTCAGCATGAAGAAAAATCGCGCTTTTTTGAAGTGGGCAGGGGGCAAGTATCCCCTGc  >  1:1422332/1‑69 (MQ=255)
                                                           gtAATCAGTTAGTCAGCATGAAGAAAAATCGCGCTTTTTTGaa                            <  1:891304/43‑1 (MQ=255)
                                                           gtAATCAGTTAGTCAGCATGAAGAAAAATCGCGCTTTTTTGaa                            <  1:1069115/43‑1 (MQ=255)
                                                           gtAATCAGTTAGTCAGCATGAAGAAAAATCGCGCTTTTTTGaa                            <  1:1148459/43‑1 (MQ=255)
                                                               ttagttagTCAGCATGAAGAAAAATCGCGCTTTTTTGAAGTGGGc                      >  1:397650/1‑45 (MQ=255)
                                                               ttagttagTCAGCATGAAGAAAAATCGCGCTTTTTTGAAGTGGGc                      >  1:3530330/1‑45 (MQ=255)
                                                               ttagttagTCAGCATGAAGAAAAATCGCGCTTTTTTGAAGTGGGc                      >  1:3433964/1‑45 (MQ=255)
                                                               ttagttagTCAGCATGAAGAAAAATCGCGCTTTTTTGAAGTGGGc                      >  1:246412/1‑45 (MQ=255)
                                                                |                                                               
TGCGGGAGTTGCCTGAAGCGCTGGATGCTGTCGGAGCTTTCTCCACAGCCGGAGAAGGTGTAATTAGTTAGTCAGCATGAAGAAAAATCGCGCTTTTTTGAAGTGGGCAGGGGGCAAGTATCCCCTGC  >  W3110S.gb/4124429‑4124556

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