Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,993,568 C→A 16.6% intergenic (‑9/+25) ygeK ← / ← ygeL predicted DNA‑binding transcriptional regulator/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,993,5680CA16.6% 83.8 / 4.2 30intergenic (‑9/+25)ygeK/ygeLpredicted DNA‑binding transcriptional regulator/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (9/16);  new base A (2/3);  total (11/19)
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.48e-01

ATATGAGCATCTATCTTATGCGATTTAACCCATTTTACGAGTTCGGCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGA  >  W3110S.gb/2993507‑2993635
                                                             |                                                                   
atatGAGCATCTATCTTATGCGATTTAACCCATTTTACGAGTTCGGCCCCCATCATTCCACCTGCAAATtc                                                            >  1:890863/1‑71 (MQ=255)
         tctatctTATGCGATTTAACCCATTTTACGAGTGCGGCCCCCATCATTCCACCTGCaaa                                                               <  1:1847048/59‑1 (MQ=255)
           tatctTATGCGATTTAACCCATTTTACGAGTTCGGCCCCCATCATTCCACCTGCAAATTCTCCTAaaata                                                  >  1:2720015/1‑70 (MQ=255)
             tctTATGCGATTTAACCCATTTTACGAGTTCGGCCCCCATCATTCCACCTGCAAATTCTCCTAaaata                                                  <  1:702829/68‑1 (MQ=255)
             tctTATGCGATTTAACCCATTTTACGAGTTCGGCCCCCATCATTCCACCTGCAAATTCTCCTAaaata                                                  <  1:333520/68‑1 (MQ=255)
             tctTATGCGATTTAACCCATTTTACGAGTTCGGCCCCCATCATTCCACCTGCAAATTCTCCTAaaata                                                  <  1:2738795/68‑1 (MQ=255)
               ttATGCGATTTAACCCATTTTACGAGTTCGGCCCCCATCATTCCACCTGCAAATTCTCCTa                                                       >  1:448684/1‑61 (MQ=255)
                 aTGCGATTTAACCCATTTTACGAGTTCGGCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATAtt                                            <  1:2270332/70‑1 (MQ=255)
                 aTGCGATTTAACCCATTTTACGAGTTCCGCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATAtt                                            <  1:1166619/70‑1 (MQ=255)
                      aTTTAACCCATTTTACGAGTTCGGCCCCCATCATTCCACATGCaaa                                                               <  1:510130/46‑1 (MQ=255)
                        ttAACCCATTTTACGAGTTCGGCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTg                                     >  1:2120566/1‑70 (MQ=255)
                           aCCCATTTTACGAGTTCGGCCCCCATCATTCCACATGCAAATTCTCCTAAAATAAATATTTAGGCTGTAg                                  <  1:1271878/70‑1 (MQ=255)
                                ttttACGAGTTCGGCCCCCATCATTCCACCTGCAAATTCTCCTAaaata                                                  <  1:838473/49‑1 (MQ=255)
                                ttttACGAGTTCGGCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGt                            <  1:1631846/71‑1 (MQ=255)
                                ttttACGAGTTCGGCCCCCATCATTCCACATGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGt                            <  1:11707/71‑1 (MQ=255)
                                  ttACGAGTTCGGCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTa                          >  1:138100/1‑71 (MQ=255)
                                   tACGAGTTCGGCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTa                          <  1:1835288/70‑1 (MQ=255)
                                    aCGAGTTCGGCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTac                         <  1:556017/70‑1 (MQ=255)
                                       aGTTCGGCCCCCATCATTCCACCTGCAAATTCTCCTaa                                                      <  1:573675/38‑1 (MQ=255)
                                                 ccATCATTCCACCTGCAAATTCTCCTAaaataaat                                               <  1:1928233/35‑1 (MQ=255)
                                                 ccATCATTCCACCTGCAAATTCTCCTAaaataaat                                               <  1:1837663/35‑1 (MQ=255)
                                                   atcatTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGAt                                >  1:2477255/1‑48 (MQ=255)
                                                   atcatTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGt                            <  1:202981/52‑1 (MQ=255)
                                                    tcatTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTa                          <  1:580051/53‑1 (MQ=255)
                                                      atTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTg                             >  1:2220003/1‑48 (MQ=255)
                                                      atTCCACATGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTttctt      >  1:450004/1‑71 (MQ=255)
                                                      atTCCACATGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTttctt      >  1:675493/1‑71 (MQ=255)
                                                         ccACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTttctt      >  1:2383659/1‑68 (MQ=255)
                                                         ccACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTttctt      >  1:137628/1‑68 (MQ=255)
                                                          cACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGa  <  1:1995714/71‑1 (MQ=255)
                                                             |                                                                   
ATATGAGCATCTATCTTATGCGATTTAACCCATTTTACGAGTTCGGCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGA  >  W3110S.gb/2993507‑2993635

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: