Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,591,039 G→A 25.0% intergenic (+28/+13) yihG → / ← yihF predicted endonuclease/conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,591,0390GA25.0% 48.2 / 14.4 28intergenic (+28/+13)yihG/yihFpredicted endonuclease/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (13/8);  new base A (3/4);  total (16/12)
Fisher's exact test for biased strand distribution p-value = 4.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTTCTGACGTCTCTCATGTCATCACAGCGTCAGAATAAGTAAGTTTTACTCGCCTCTGAATGAGCAGAGGCGAGTGAGTATTTTAATGAACAAAACGTCCGGCACGAGACATAAATTCTTCTTCGCTCATCTC  >  W3110S.gb/3590969‑3591102
                                                                      |                                                               
gCTTCTGACGTCTCTCATGTCATCACAGCGTCAGAATAAGTAAGTTTTACTCGCCTCTGAATGAGCAGAgg                                                                 >  1:3239149/1‑71 (MQ=255)
      gACGTCTCTCATGTCATCACAGCGTCAGAATAAGTAAGTTTTACTCGCCTCTGAATGAGCAGAGGCgagt                                                            >  1:1143974/1‑70 (MQ=255)
                aTGTCATCACAGCGTCAGAATAAGTAAGTTTTACTCGCCTCTGAATGAGCAGAGGCGAGTTAGTAtttt                                                   <  1:2539557/69‑1 (MQ=255)
                aTGTCATCACAGCGTCAGAATAAGTAAGTTTTACTCGCCTCTGAATGAGCAGAGGCGAGTGAGTAtttt                                                   <  1:93672/69‑1 (MQ=255)
                aTGGCATCACAGCGTCAGAATAAGTAAGATTTACTCGCCTCTGAATGAGCAGAGACGAGTGAGTATTTTa                                                  <  1:3865248/70‑1 (MQ=255)
                   tcatcaCAGCGTCAGAATAAGTAAGTTTTACTCGCCTCTGAATGAGCAGAGGCGAGTGAGTATTTTAATg                                               >  1:1565027/1‑70 (MQ=255)
                        acaGCGTCAGAATAAGTAAGTTTTACTCGCCTCTGAATGAGCAGAGGCGAGTGAGTATTTTAATGAACaa                                          <  1:921282/70‑1 (MQ=255)
                           gCGTCAGAATAAGTAAGTTTTACTCGCCTCTGAATGAGCAGAGGCGAGTGAGTATTTTAATGAACAAAACg                                      <  1:1627310/71‑1 (MQ=255)
                           gCGTCAGAATAAGTAAGTTTTACTCGCCTCTGAATGAGCAGAGGCGAGTGAGTATTTTAATGAACAAAACg                                      <  1:594165/71‑1 (MQ=255)
                                aGAATAAGTAAGTTTTACTCGCCTCTGAATGAGCAGAGGCGAGTGAGTATTTTAATGAACAAAACGTCCgg                                 >  1:1517847/1‑71 (MQ=255)
                                aGAATAAGTAAGTTTTACTCGCCTCTGAATGAGCAGAGACGAGTGAGTATTTTAATGAACAAAACGTCCgg                                 >  1:1424972/1‑71 (MQ=255)
                                aGAATAAGTAAGTTTTACTCGCCTCTGAATGAGCAGAGACGAGTGAGTATTTTAATGAACAAAACGTCCgg                                 >  1:846133/1‑71 (MQ=255)
                                   aTAAGTAAGTTTTACTCGCCTCTGAATGAGCAGAGACGAGTGAGTATTTTAATGAACAAAACGTCCGGCAc                              <  1:1837964/71‑1 (MQ=255)
                                   aTAAGTAAGTTTTACTCGCCTCTGAATGAGCAGAGACGAGTGAGTATTTTAATGAACAAAACGTCCGGCAc                              <  1:3663870/71‑1 (MQ=255)
                                      agtaagtTTTACTCGCCTCTGAATGAGCAGAGGCGAGTGAGTATTTTAATGAACAAAACGTCCGGCAc                              >  1:3710482/1‑68 (MQ=255)
                                          agtTTTACTCGCCTCTGAATGAGCAGAGGCGAGTGAGTATTTTAATGAACAAAACGTCCGGCACGAGACa                        >  1:3990911/1‑70 (MQ=255)
                                             tttACTCGCCTCTGAATGAGCAGAGGCGAGTGAGTATTTTAATGAACAAAACGTCCGGCACgaga                          >  1:146318/1‑65 (MQ=255)
                                             tttACTCGCCTCTGAATGAGCAGAGGCGAGTGAGTATTTTAATGAACAAAACGTCCGGCACgaga                          >  1:1907774/1‑65 (MQ=255)
                                                 cTCGCCTCTGAATGAGCAGAGGCGAGTGAGTATTTTAATGAACAAAACGTCCGGCACGAGACATAAAttc                 <  1:1784704/70‑1 (MQ=255)
                                                       tctGAATGAGCAGAGGCGAGTGAGTATTTTAATGAACAAAACGTCCGGCACGAGACATAAATTCTTCTTCg          >  1:2620687/1‑71 (MQ=255)
                                                       tctGAATGAGCAGAGGCGAGTGAGTATTTTAATGAACAAAACGTCCGGCACGAGACATAAATTCTTCTTCg          >  1:602400/1‑71 (MQ=255)
                                                         tGAATGAGCAGAGGCGAGTGAGTATTTTAATGAAc                                            >  1:1720386/1‑35 (MQ=255)
                                                         tGAATGAGCAGAGGCGAGTGAGTATTTTAATGAAc                                            >  1:2395685/1‑35 (MQ=255)
                                                         tGAATGAGCAGAGGCGAGTGAGTATTTTAATGAAc                                            >  1:1438776/1‑35 (MQ=255)
                                                         tGAATGAGCAGAGGCGAGTGAGTATTTTAATGAACa                                           <  1:780786/36‑1 (MQ=255)
                                                             tGAGCAGAGACGAGTGAGTATTTTAATGAACAAAACGTCCg                                  >  1:3942244/1‑41 (MQ=255)
                                                               aGCAGAGGCGAGTGAGTATTTTAATGAACAAAACGTCCGGCACGAGACATAAATTCTTCTTCGCTCAtctc  >  1:3491425/1‑71 (MQ=255)
                                                                gCAGAGACGAGTGAGTATTTTAATGAACAAAACGTCCGTCACGAGACATAAATTCTTCTTCGCTCAtctc  <  1:3618243/70‑1 (MQ=255)
                                                                 cAGAGGCGAGTGAGTATTTTAATGAACAAAACGTCCGGCAc                              <  1:66150/41‑1 (MQ=255)
                                                                      |                                                               
GCTTCTGACGTCTCTCATGTCATCACAGCGTCAGAATAAGTAAGTTTTACTCGCCTCTGAATGAGCAGAGGCGAGTGAGTATTTTAATGAACAAAACGTCCGGCACGAGACATAAATTCTTCTTCGCTCATCTC  >  W3110S.gb/3590969‑3591102

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: