Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,089,667 A→C 11.2% R243R (CGT→CGG yeeY ← predicted DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,089,6670AC11.2% 70.5 / 6.0 36R243R (CGT→CGGyeeYpredicted DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (18/13);  new base C (2/2);  total (20/16)
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.36e-01

AAGAGAGTTTGGGCGACCACTTTCTTCTCGATTAACGCGCCACTATCGAGAAAACGTTGCGCCAGATAACGGGGTAAATAGCCGCAGCCGAGGCCGCTAATTTGCAGCTCCAGTTTGGTTTTAAAATCGAAAACGGT  >  W3110S.gb/2089599‑2089735
                                                                    |                                                                    
aaGAGAGTTTGGGCGACCACTTTCTTCTCGATTAACGCGCCACTATCGAGAAAACGTTGCGCCAGATAACg                                                                    <  1:418307/71‑1 (MQ=255)
   agagTTTGGGCGACCACTTTCTTCTCGATTAACGCGCCACTATCGAGAAAACGTTGCGCCAGATAACggg                                                                  >  1:1969613/1‑70 (MQ=255)
              gACCACTTTCTTCTCGATTAACGCGCCACTATCGAGAAAACGTTGCGCCAGATACCGTGGTAAATAg                                                          >  1:154607/1‑67 (MQ=255)
              gACCACTTTCTTCTCGATTAACGCGCCACTATCGAGAAAACGTTGCGCCAGATACCGTGGTAAATAg                                                          >  1:2526716/1‑67 (MQ=255)
                   cTTTCTTCTCGATTAACGCGCCACTATCGAGAAAACGTTGCGCCAGATAACGGGGTAAATAGCCGCAGcc                                                  >  1:3686977/1‑70 (MQ=255)
                   cTTTCTTCTCGATTAACGCGCCACTATCGAGAAAACGTTGCGCCAGATAACGGGGTAAATAGCCGCAGcc                                                  >  1:3256130/1‑70 (MQ=255)
                    tttCTTCTCGATTAACGCGCCACTATCGAGAAAACGTTGCGCCAGATAACg                                                                    <  1:399266/51‑1 (MQ=255)
                     ttcttcTCGATTAACGCGCCACTATCGAGAAAACGTTGCGCCAGATAACGGGGTAAATAGCCGCAGCCGAg                                               >  1:2624169/1‑71 (MQ=255)
                      tcttctCGATTAACGCGCCACTATCGAGAAAACGTTGCGCCAGATATCGGGGTaa                                                              <  1:1471626/55‑1 (MQ=255)
                      tcttctCGATTAACGCGCCACTATCGAGAAAACGTTGCGCCAGATAACGGGGTAAATAGCCGCAGCCGAgg                                              >  1:2726946/1‑71 (MQ=255)
                      tcttctCGATTAACGCGCCACTATCGAGAAAACGTTGCGCCAGATAACGGGGTAAATAGCCGCAGCCGAgg                                              >  1:1286918/1‑71 (MQ=255)
                          ctcGATTAACGCGCCACTATCGAGAAAACGTTGCGCCAGATAACGGGGTAAAt                                                            >  1:1068286/1‑53 (MQ=255)
                             gATTAACGCGCCACTATCGAGAAAACGTTGCGCCAGATAAc                                                                     >  1:1972807/1‑41 (MQ=255)
                             gATTAACGCGCCACTATCGAGAAAACGTTGCGCCAGATAAc                                                                     >  1:1968490/1‑41 (MQ=255)
                                          cTATCGAGAAAACGTTGCGCCAGATAACGGGGTAAATAGCCGCAGCCGAGGCCGCTAATTTGCAGCTCCAg                          >  1:3248311/1‑71 (MQ=255)
                                           tATCGAGAAAACGTTGCGCCAGATACCGTGGTAAATAGCCGCAGCCAAGGCCGCTAATTTGc                                  <  1:3430695/62‑1 (MQ=255)
                                           tATCGAGAAAACGTTGCGCCAGATACCGTGGTAAATAGCCGCAGCCAAGGCCGCTAATTTGc                                  <  1:2762008/62‑1 (MQ=255)
                                             tCGAGAAAACGTTGCGCCAGATAACGGGGTAAATAGCCGCAGCCGAGGCCGCTAATTTGCa                                 <  1:2555864/61‑1 (MQ=255)
                                               gagaAAACGTTGCGCCAGATAACGGGGTAAATAGCCGCAGCCGAGGCCGCTAATTTGCAGc                               >  1:780322/1‑61 (MQ=255)
                                               gagaAAACGTTGCGCCAGATAACGGGGTAAATAGCCGCAGCCGAGGCCGCTAATTTGCAGc                               >  1:2937682/1‑61 (MQ=255)
                                                   aaaCGTTGCGCCAGATAACGGGGTAAATAGCCGCAGCCGAGGCCGCTAATTTGCAGCTCCAGTTTGGttt                  >  1:1636012/1‑70 (MQ=255)
                                                    aaCGTTGCGCCAGATAACGGGGTAAATAGCCGCAGCCGAGGCCGCTAATTTGCa                                 >  1:2662180/1‑54 (MQ=255)
                                                    aaCGTTGCGCCAGATAACGGGGTAAATAGCCGCAGCCGAGGCCGCTAATTTGCAGCTCCAGTTTGGtttt                 >  1:3617999/1‑70 (MQ=255)
                                                      cGTTGCGCCAGATAACGGGGTAAATAGCCGCAGCCGAGGCCGCTAATTTGCAGCTCCAGTTTGGTTTTaaa              <  1:890281/71‑1 (MQ=255)
                                                      cGTTGCGCCAGATAACGGGGTAAATAGCCGCAGCCGAGGCCGCTAATTTGCAGCTCCAGTTTGGTTTTaa               >  1:1317192/1‑70 (MQ=255)
                                                       gTTGCGCCAGATAACGGGGTAAATAGCCGCAGCCGAGGCCGCTAATTTGCa                                 >  1:3170185/1‑51 (MQ=255)
                                                       gTTGCGCCAGATAACGGGGTAAATAGCCGCAGCCGAGGCCGCTAATTTGCAGCTCCAGTTTGGTTTTaaaa             <  1:1889344/71‑1 (MQ=255)
                                                       gTTGCGCCAGATAACGGGGTAAATAGCCGCAGCCGAGGCCGCTAATTTGCAGCTCCAGTTTGGTTTTaaaa             <  1:636945/71‑1 (MQ=255)
                                                            gcCAGATAACGGGGTAAATAGCCGCAGCCGAGGCCGCTAATTTGCAGCTCCAGTTTGGTTTTAAAATCGaa        <  1:201513/71‑1 (MQ=255)
                                                            gcCAGATAACGGGGTAAATAGCCGCAGCCGAGGCCGCTAATTTGCAGCTCCAGTTTGGTTTTAAAATCGaa        <  1:720464/71‑1 (MQ=255)
                                                             ccAGATAACGGGGTAAATAGCCGCAGCCGAGGCCGCTAATTTGCAGCTCCAGTTTGGTTTTaaa              <  1:117599/64‑1 (MQ=255)
                                                                  tAACGGGGTAAATAGCCGCAGCCGAGGCCGCTAATTTGCAg                                >  1:1144257/1‑41 (MQ=255)
                                                                  tAACGGGGTAAATAGCCGCAGCCGAGGCCGCTAATTTGCAGCTCCAGTTTGGTTTTAAAATCGAAAACGGt  <  1:2958176/71‑1 (MQ=255)
                                                                  tAACGGGGTAAATAGCCGCAGCCGAGGCCGCTAATTTGCAGCTCCAGTTTGGTTTTAAAATCGAAAACGGt  <  1:3190305/71‑1 (MQ=255)
                                                                  tAACGGGGTAAATAGCCGCAGCCGAGGCCGCTAATTTGCAGCTCCAGTTTGGTTTTAAAATCGAAAACGGt  <  1:1057355/71‑1 (MQ=255)
                                                                  tAACGGGGTAAATAGCCGCAGCCGAGGCCGCTAATTTGCAGCTCCAGTTTGGTTTTAAAATCGAAAACGGt  <  1:452377/71‑1 (MQ=255)
                                                                   aaCGGGGTAAATAGCCGCAGCCGAGGCCGCTAATTTGCAGCTCCAGTTTGGTTTTAAAATCGAAAACGGt  <  1:2268684/70‑1 (MQ=255)
                                                                    |                                                                    
AAGAGAGTTTGGGCGACCACTTTCTTCTCGATTAACGCGCCACTATCGAGAAAACGTTGCGCCAGATAACGGGGTAAATAGCCGCAGCCGAGGCCGCTAATTTGCAGCTCCAGTTTGGTTTTAAAATCGAAAACGGT  >  W3110S.gb/2089599‑2089735

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: