Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 338,985 G→A 19.6% intergenic (+18/‑32) yahG → / → yahH conserved hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb338,9850GA19.6% 102.9 / 18.5 46intergenic (+18/‑32)yahG/yahHconserved hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (16/21);  new base A (6/3);  total (22/24)
Fisher's exact test for biased strand distribution p-value = 2.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.68e-01

TGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGT  >  W3110S.gb/338916‑339048
                                                                     |                                                               
tGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGa                                                                >  1:3178700/1‑71 (MQ=255)
 gCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGAt                                                               >  1:3319722/1‑71 (MQ=255)
  ccTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGa                                                                <  1:1094914/69‑1 (MQ=255)
  ccTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGAt                                                               <  1:3571295/70‑1 (MQ=255)
     gCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGa                                                                >  1:343105/1‑66 (MQ=255)
         cGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATgcggcg                                                         <  1:355908/69‑1 (MQ=255)
          gagaAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTg                                                       >  1:1606689/1‑70 (MQ=255)
          gagaAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTg                                                       <  1:4196226/70‑1 (MQ=255)
             aaaGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATgcggc                                                          <  1:2885639/64‑1 (MQ=255)
             aaaGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACg                                                   >  1:1046795/1‑71 (MQ=255)
              aaGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGc                                                  <  1:2715722/71‑1 (MQ=255)
              aaGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGc                                                  <  1:966934/71‑1 (MQ=255)
               aGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGcc                                                 >  1:3086115/1‑71 (MQ=255)
               aGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGAATGCGGCGTGAACGcc                                                 >  1:2802019/1‑71 (MQ=255)
                 ccATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCtt                                               <  1:1377214/71‑1 (MQ=255)
                 ccATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCtt                                               <  1:3978488/71‑1 (MQ=255)
                 ccATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCtt                                               <  1:1894502/71‑1 (MQ=255)
                    ttCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATgcggcg                                                         <  1:4543342/58‑1 (MQ=255)
                    ttCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAAc                                                    <  1:4497057/63‑1 (MQ=255)
                    ttCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATc                                            <  1:2368279/71‑1 (MQ=255)
                    ttCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGAATGCGGCGTGa                                                      <  1:1691404/61‑1 (MQ=255)
                      cTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCg                                          >  1:2339213/1‑71 (MQ=255)
                      cTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCg                                          >  1:1030980/1‑71 (MQ=255)
                              gTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGaa                                                     >  1:3407648/1‑52 (MQ=255)
                                    aCGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGaa                                                     >  1:3610407/1‑46 (MQ=255)
                                     cGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTAt                                             >  1:52706/1‑53 (MQ=255)
                                       ttACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGgcgc                           >  1:4296519/1‑69 (MQ=255)
                                         aCGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGaa                                                     >  1:548129/1‑41 (MQ=255)
                                         aCGGCGTCTGACCTACATGTTCATGCCGAATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAg                        >  1:3294173/1‑70 (MQ=255)
                                          cGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATc                                            <  1:4610852/49‑1 (MQ=255)
                                          cGGCGTCTGACCTACATGTTCATGCCGAATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGa                       <  1:967651/70‑1 (MQ=255)
                                           ggCGTCTGACCTACATGTTCATGCCGAATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAAt                     <  1:1071273/71‑1 (MQ=255)
                                              gTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTg                  <  1:488100/71‑1 (MQ=255)
                                               tCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGa                       <  1:751418/65‑1 (MQ=255)
                                                 tGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAg               >  1:2116043/1‑71 (MQ=255)
                                                   aCCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGc             <  1:3271811/71‑1 (MQ=255)
                                                     cTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAg               <  1:3783959/67‑1 (MQ=255)
                                                       aCATGTTCATGCCGAATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGc             >  1:262360/1‑67 (MQ=255)
                                                       aCATGTTCATGCCGAATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGa         >  1:3214458/1‑71 (MQ=255)
                                                       aCATGTTCATGCCGAATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGa         >  1:208749/1‑71 (MQ=255)
                                                       aCATGTTCATGCCGAATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGAAGAATCTGTAGGCCTGa         >  1:4183707/1‑70 (MQ=14)
                                                        cATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGc                                        >  1:2292852/1‑39 (MQ=37)
                                                        cATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGa         <  1:2571254/70‑1 (MQ=255)
                                                          tGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATaa      <  1:2536850/71‑1 (MQ=39)
                                                            ttCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGc    >  1:249879/1‑71 (MQ=35)
                                                               aTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGt  <  1:2939614/70‑1 (MQ=31)
                                                                     |                                                               
TGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGT  >  W3110S.gb/338916‑339048

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: