Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,224,919 G→A 64.3% intergenic (‑423/‑123) ymgG ← / → ymgH hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,224,9190GA64.3% 22.2 / 26.1 28intergenic (‑423/‑123)ymgG/ymgHhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (6/4);  new base A (11/7);  total (17/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATATTCATCAGCAACGATTACATTAGTCATTTTATTTTGCCGACGGCCTCATTGTCGAAAGATAAGCGTACGACAGTATTATCAGAAAAGAGTGATTTTTTATCCAACTACACTTCAGCGCACTGCGTGTAAA  >  W3110S.gb/1224851‑1224984
                                                                    |                                                                 
gaTATTCATCAGCAACGATTACATTAGTCATTTTATTTTGCCGACGGCCTCATTGTCGAAAGATAAGCGTa                                                                 >  1:1532184/1‑71 (MQ=255)
gaTATTCATCAGCAACGATTACATTAGTCATTTTATTTTGCCGACGGCCTCATTGTCGAAAGATAAGCGTa                                                                 >  1:598656/1‑71 (MQ=255)
      catcaGCAACGATTACATTAGTCATTTTATTTTGCCGACGGCCTCATTGTCGAAAGATAAGCGTACGACa                                                            <  1:29583/70‑1 (MQ=255)
             aaCGATTACATTAGTCATTTTATTTTGCCGACGGCCTCATTGTCGAAAGATAAGCGTACGACAGTATTATc                                                    <  1:1576192/71‑1 (MQ=255)
             aaCGATTACATTAGTCATTTTATTTTGCCGACGGCCTCATTGTCGAAAGATAAGCATACGACAGtattat                                                     >  1:588093/1‑70 (MQ=255)
             aaCGATTACATTAGTCATTTTATTTTGCCGACGGCCTCATTGTCGAAAGATAAGCATACGACAGtattat                                                     >  1:1548666/1‑70 (MQ=255)
                 aTTACATTAGTCATTTTATTTTGCCGACGGCCTCATTGTCGAAAGATAAGCGTACGACAGTATTATCAGaa                                                <  1:1076439/71‑1 (MQ=255)
                       ttAGTCATTTTATTTTGCCGACGGCCTCATTGTCGAAAGATAAGCATACGACAGTATTATCAGAAAAGAGt                                          >  1:1607423/1‑71 (MQ=255)
                             attttattttGCCGACGGCCTCATTGTCGAAAGATAAGCATACGACa                                                            >  1:1169449/1‑47 (MQ=255)
                             attttattttGCCGACGGCCTCATTGTCGAAAGATAAGCATACGACAGTATTATCAGAAAAGAGTGAtttt                                    <  1:1448325/71‑1 (MQ=255)
                              ttttattttGCCGACGGCCTCATTGTCGAAAGATAAGCATACGa                                                              <  1:695421/44‑1 (MQ=255)
                              ttttattttGCCGACGGCCTCATTGTCGAAAGATAAGCATACGa                                                              <  1:506148/44‑1 (MQ=255)
                               tttattctGCCGACGGCCTCATTGTCGAAAGATAAGCGTACGACAGTATTATCAGAAAAGAGTGAttttt                                   >  1:1569423/1‑70 (MQ=255)
                                ttattttGCCGACGGCCTCATTGTCGAAAGATAAGCGTACGACAGTATTATCAGAAAagag                                           >  1:849435/1‑61 (MQ=255)
                                  attttGCCGACGGCCTCATTGTCGAAAGATAAGCGTACGACAGTATTATCAGAAAAGAGTGATTTTTTATc                               >  1:1365691/1‑71 (MQ=255)
                                  attttGCCGACGGCCTCATTGTCGAAAGATAAGCATACGACAGTATTATCAGAAAAGAGTGATTTTTTATc                               >  1:789421/1‑71 (MQ=255)
                                               ccTCATTGTCGAAAGATAAGCGTACGACAGTATTATCAGAAAAGAGTGATTTTTTATCCAACTACACtt                    >  1:493948/1‑69 (MQ=255)
                                               ccTCATTGTCGAAAGATAAGCGTACGACAGTATTATCAGAAAAGAGTGATTTTTTATCCAACTACACTTCa                  <  1:1222084/71‑1 (MQ=255)
                                               ccTCATTGTCGAAAGATAAGCATACGACAGTATTATCAGAAAAGAGTGAtt                                      <  1:1135274/51‑1 (MQ=255)
                                                 tCATTGTCGAAAGATAAGCATACGACAGTATTATCAGAAAAGAGTGATTTTTTATCCAACTACACTTCAgc                >  1:1764514/1‑71 (MQ=255)
                                                   aTTGTCGAAAGATAAGCATACGACAGTATTATCAGAAAAGAGTGATTTTTTATCCaa                            >  1:875366/1‑57 (MQ=255)
                                                    ttGTCGAAAGATAAGCATACGACAGTATTATCAGAAAag                                             <  1:580461/39‑1 (MQ=255)
                                                    ttGTCGAAAGATAAGCATACGACAGTATTATCAGAAAag                                             <  1:303341/39‑1 (MQ=255)
                                                      gTCGAAAGATAAGCATACGACAGTATTATCAGAAAAGAGTGATTTTTTATCCAACTACACTTCAGCGCAc            >  1:955348/1‑70 (MQ=255)
                                                      gTCGAAAGATAAGCATACGACAGTATTATCAGAAAAGAGTGATTTTTTATCCAACTACACTTCAGCGCAc            >  1:1516764/1‑70 (MQ=255)
                                                           aaGATAAGCATACGACAGTATTATCAGAAAAGAGTGATTTTTTATCCAACTacac                      >  1:1734197/1‑55 (MQ=255)
                                                            aGATAAGCATACGACAGTATTATCAGAAAAGAGTGATTTTTTATCCAACTACACTTCAGCGCACTGCgtgt     <  1:689327/71‑1 (MQ=255)
                                                                aaGCATACGACAGTATTATCAGAAAAGAGTGATTTTATATCCAACTACACTTCAGCGCACTGCGTGTaaa  >  1:180456/1‑70 (MQ=255)
                                                                    |                                                                 
GATATTCATCAGCAACGATTACATTAGTCATTTTATTTTGCCGACGGCCTCATTGTCGAAAGATAAGCGTACGACAGTATTATCAGAAAAGAGTGATTTTTTATCCAACTACACTTCAGCGCACTGCGTGTAAA  >  W3110S.gb/1224851‑1224984

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