Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,591,039 G→A 37.8% 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,0390GA37.8% 29.8 / 36.4 37intergenic (+28/+13)yihG/yihFpredicted endonuclease/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (7/16);  new base A (6/8);  total (13/24)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.06e-01

CTCTCATGTCATCACAGCGTCAGAATAAGTAAGTTTTACTCGCCTCTGAATGAGCAGAGGCGAGTGAGTATTTTAATGAACAAAACGTCCGGCACGAGACATAAATTCTTCTTCGCTCATCT  >  W3110S.gb/3590980‑3591101
                                                           |                                                              
ctctcATGTCATCACAGCGTCAGAATAAGTAAGTTTTACTCGCCTCTGAATGAGCAGAGACGAGTGAGTAt                                                     >  1:899026/1‑71 (MQ=255)
  ctcATGTCATCACAGCGTCAGAATAAGTAAGTTTTACTCGCCTCTGAATGAGCAGAGACGAGTGAGTAtt                                                    >  1:477650/1‑70 (MQ=255)
     aTGTCATCACAGCGTCAGAATAAGTAAGTTTTACTCGCCTCTGAATGAGCAGAGGCGAGTGAGTATTTTaa                                                >  1:2104192/1‑71 (MQ=255)
     aTGTCATCACAGCGTCAGAATAAGTAAGTTTTACTCGCCTCTGAATGAGCAGAGACGAGTGAGTAtttt                                                  <  1:480256/69‑1 (MQ=255)
     aTGTCATCACAGCGTCAGAATAAGTAAGTTTTACTCGCCTCTGAATGAGCAGAGACGAGTGAGTAttt                                                   <  1:878017/68‑1 (MQ=255)
     aTGTCATCACAGCGGCAGAATAAGTAAGTTTTACTCGCCTCTGAATGAGCAGAGGCGAGTGAGTAtttt                                                  <  1:4513262/69‑1 (MQ=255)
        tcatcaCAGCGTCAGAATAAGTAAGTTTTACTCGCCTCTGAATGAGCAGAGGc                                                               <  1:1125974/53‑1 (MQ=255)
        tcatcaCAGCGTCAGAATAAGTAAGTTTTACTCGCCTCTGAATGAGCAGAGGc                                                               <  1:3459756/53‑1 (MQ=255)
            cacaGCGTCAGAATAAGTAAGTTTTACTCGCCTCTGAATGAGCAGAGGCGAGTGAGTATTTTAATGAACa                                          >  1:870375/1‑70 (MQ=255)
             acaGCGTCAGAATAAGTAAGTTTTACTCGCCTCTGAATGAGCAGAGGCGAGTGAGTATTTTAATGAACaa                                         <  1:814973/70‑1 (MQ=255)
                  gTCAGAATAAGTAAGTTTTACTCGCCTCTGAATGAGCAGAGGCGAGTGAGTATTTTAATGAACAAAACGTc                                   <  1:1043231/71‑1 (MQ=255)
                   tCAGAATAAGTAAGTTTTACTCGCCTCTGAATGTGCAGAGGCGAGTGAGTATTTTAATGAACAAAACGt                                    <  1:2492807/69‑1 (MQ=255)
                   tCAGAATAAGTAAGTTTTACTCGCCTCTGAATGAGCAGAGACGAGTGAGTATTTTAATGAACaa                                         <  1:4193990/64‑1 (MQ=255)
                      gAATAAGTAAGTTTTACTCGCCTCTGAATGAGCAGAGACGAGTGAGTATTTTAATGAACAAAACg                                     <  1:528039/65‑1 (MQ=255)
                       aaTAAGTAAGTTTTACTCGCCTCTGAATGAGCAGAGGCGAGTGAGTATTTTaa                                                <  1:4002230/53‑1 (MQ=255)
                        aTAAGTAAGTTTTACTCGCCTCTGAATGAGCAGAGACGAGTGAGTATTTTAATGAACAAAACGTCCGGCAc                             >  1:544306/1‑71 (MQ=255)
                        aTAAGTAAGTTTTACTCGCCTCTGAATGAGCAGAGACGAGTGAGTATTTTAATGAACAAAACGTCCGGCAc                             <  1:3013447/71‑1 (MQ=255)
                         taagtaagTTTTACTCGCCTCTGAATGAGCAGAGGCGAGTGAGTATTTTAATGAACAAAAc                                      >  1:2439206/1‑61 (MQ=255)
                         taagtaagTTTTACTCGCCTCTGAATGAGCAGAGGCGAGTGAGTATTTTAATGAACAAAACGTCCGGCAc                             <  1:4364360/70‑1 (MQ=255)
                         taagtaagTTTTACTCGCCTCTGAATGAGCAGAGGCGAGTGAGTATTTTAATGAACAAAACGTCCGGCAc                             >  1:3878501/1‑70 (MQ=255)
                               agtTTTACTCGCCTCTGAATGAGCAGAGGCGAGTGAGTATTTTAATGAACAAAACGTCCGGCACGAGACAt                      <  1:55315/71‑1 (MQ=255)
                                 ttttACTCGCCTCTGAATGAGCAGAGGCGAGTGAGTATTTTAATg                                              <  1:276667/45‑1 (MQ=255)
                                   ttACTCGCCTCTGAATGAGCAGAGGCGAGTGAGTATTTTAATGAACAAAACGTCCGGCACGAGACATaaa                   <  1:2929907/70‑1 (MQ=255)
                                    tACTCGCCTCTGAATGAGCAGAGGCGAGTGAGTATTTTAATGAACAAAACGTCCGGCACGAGACATAAAtt                 >  1:4228644/1‑71 (MQ=255)
                                      cTCGCCTCTGAATGAGCAGAGGCGAGTGAGTATTTTAATGAACAAAACGTCCGGCACGAGACATAAAttct               <  1:4065041/71‑1 (MQ=255)
                                      cTCGCCTCTGAATGAGCAGAGGCGAGTGAGTATTTTAATGAACAAAACGTCCGGCACGAGACATAAAttct               <  1:3566755/71‑1 (MQ=255)
                                      cTCGCCTCTGAATGACCAGAGACGAGTGAGTATTTTAATGAACaa                                         <  1:3636750/45‑1 (MQ=255)
                                         gCCTCTGAATGAGCAGAGACGAGTGAGTATTTTAATGAACAAAACGTCCGGCACGAGACATAAAttctt              <  1:1101047/69‑1 (MQ=255)
                                          ccTCTGAATGAGCAGAGGCGAGTGAGTATTTTAATGAACAAAACGTCCGGCACGAGACATAAAttcttctt           <  1:1983888/71‑1 (MQ=255)
                                            tctGAATGAGCAGAGACGAGTGAGTATTTTAATGAACAAAACGTCCGGCACGAGACATAAAttcttcttc          >  1:921021/1‑70 (MQ=255)
                                              tGAATGAGCAGAGGCGAGTGAGTATTTTAATGAACAAAACGTCCGGCACGAGACATAAATTCTTCTTCGCt       >  1:3975274/1‑71 (MQ=255)
                                              tGAATGAGCAGAGACGAGTGAGTATTTTAATGAACAAAAc                                      >  1:2674149/1‑40 (MQ=255)
                                              tGAATGAGCAGAGACGAGTGAGTATTTTAATGAACAAAACGTCCGGCACGAGACATAAATTCTTCTTCg         <  1:221287/69‑1 (MQ=255)
                                                aaTGAGCAGAGGCGAGTGAGTATTTTAATGAACAAAACGTCCGGCACGAGACATAAATTCTTCTTCGCTc      <  1:26249/70‑1 (MQ=255)
                                                aaTGAGCAGAGACGAGTGAGTATTTTAATGAACAAAACGTCCGGCACGAGACATAAATTCTTCTTCGCTCa     >  1:1412125/1‑71 (MQ=255)
                                                   gAGCAGAGGCGAGTGAGTATTTTAATGAACAAAACGTCCGGCACGAGACATAAATTCTTCTTCGGTCAtct  >  1:277026/1‑71 (MQ=255)
                                                         agGCGAGTGAGTATTTTAATGAACAAAACGTCCGGCAc                             <  1:1190096/38‑1 (MQ=255)
                                                           |                                                              
CTCTCATGTCATCACAGCGTCAGAATAAGTAAGTTTTACTCGCCTCTGAATGAGCAGAGGCGAGTGAGTATTTTAATGAACAAAACGTCCGGCACGAGACATAAATTCTTCTTCGCTCATCT  >  W3110S.gb/3590980‑3591101

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: