Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,068,568 T→A 24.0% 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,5680TA24.0% 44.7 / 12.6 25G22G (GGT→GGAymdFconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (10/9);  new base A (2/4);  total (12/13)
Fisher's exact test for biased strand distribution p-value = 6.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

AATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGC  >  W3110S.gb/1068502‑1068630
                                                                  |                                                              
aaTGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTa                                                               >  1:790954/1‑68 (MQ=255)
   ggCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTaaaa                                                            <  1:419809/68‑1 (MQ=255)
   ggCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTaaaa                                                            <  1:321050/68‑1 (MQ=255)
   ggCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTaaaa                                                            <  1:1856713/68‑1 (MQ=255)
   ggCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTaaaa                                                            <  1:1075111/68‑1 (MQ=255)
   ggCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCCGAAAGAGCATCAGAAGCAGGTaaaaaa                                                          <  1:3025575/69‑1 (MQ=255)
    gCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTaaaaaa                                                          <  1:2048657/69‑1 (MQ=255)
    gCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTaaaaaa                                                          >  1:3197140/1‑69 (MQ=255)
          cATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGAaaaaaa                                                          <  1:2307882/63‑1 (MQ=255)
          cATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGAaaaaaa                                                          <  1:305534/63‑1 (MQ=255)
           aTCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGAAAAAAAGGTGGac                                                   <  1:3627179/69‑1 (MQ=255)
                   ggTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGacagcacagc                                           >  1:154971/1‑69 (MQ=255)
                      tCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGAAAAAAAGGTGGACAGCACAGCggg                                        >  1:3084143/1‑69 (MQ=255)
                      tCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGAAAAAAAGGTGGACAGCACAGCggg                                        >  1:327650/1‑69 (MQ=255)
                       ccGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCggg                                        >  1:1907659/1‑68 (MQ=255)
                                ttGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTaaaa                                                            >  1:3354395/1‑39 (MQ=255)
                                  gCAGAAGACCGCGAAAGAGCATCAGAAGCAGTTaaaaaa                                                          <  1:1185250/39‑1 (MQ=255)
                                   cAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGTTAATTTCAAg                            <  1:206647/68‑1 (MQ=255)
                                        gaCCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTg                                                      >  1:2219936/1‑37 (MQ=255)
                                                         aGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGa     >  1:1754300/1‑69 (MQ=255)
                                                         aGAAGCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAg                            <  1:2587521/46‑1 (MQ=255)
                                                            aGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGc  >  1:3325954/1‑69 (MQ=255)
                                                            aGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGc  >  1:392995/1‑69 (MQ=255)
                                                            aGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGc  >  1:575828/1‑69 (MQ=255)
                                                                ggTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAg                            <  1:2767686/39‑1 (MQ=255)
                                                                  |                                                              
AATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGC  >  W3110S.gb/1068502‑1068630

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: