Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 338,000 A→G 12.0% H151R (CAC→CGC)  yahG → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb338,0000AG12.0% 124.8 / 9.8 50H151R (CAC→CGC) yahGconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (25/19);  new base G (3/3);  total (28/22)
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.78e-01

AAGATCTCGACGAGGCGGCTGAACTGGCGGCTTCCGGGGAGATCACCTTCTCGCCGTGTCACGAGCACGACTGCGTGGGATCGATGGCGGGTGTTACCTCGGCCTCGATGTTTATGCACATCGTAAAAAACAAAAC  >  W3110S.gb/337934‑338069
                                                                  |                                                                     
aaGATCTCGACGAGGCGGCTGAACTGGCGGCTTCCGGGGAGATCACCTTCTCGCCGTGTCACGAGCAc                                                                      <  1:2520892/68‑1 (MQ=255)
                ggcTGAACTGGCGGCTTCCGGGGAGATCACCTTCTCGCCGTGTCACGAGCACGACTGCGTGGGATCGATgg                                                   <  1:3598845/71‑1 (MQ=255)
                ggcTGAACTGGCGGCTTCCGGGGAGATCACCTTCTCGCCGTGTCACGAGCACGACTGCGTGGGATCGATgg                                                   <  1:3131117/71‑1 (MQ=255)
                   tGAACTGGCGGCTTCCGGGGAGATCACCTTCTCGCCGTGTCACGAGCACGACTGCGTGGGATCGATGGCgg                                                >  1:2593254/1‑71 (MQ=255)
                   tGAACTGGCGGCTTCCGGGGAGATCACCTTCTCGCCGTGTCACGAGCACGACTGCGTGGGATCGATGGCgg                                                >  1:749936/1‑71 (MQ=255)
                    gAACTGGCGGCTTCCGGGGAGATCACCTTCTCGCCGTGTCACGAGCACGACTGCGTGGGa                                                          >  1:411710/1‑60 (MQ=255)
                    gAACTGGCGGCTTCCGGGGAGATCACCTTCTCGCCGTGTCACGAGCACGACTGCGTGGGa                                                          >  1:1640087/1‑60 (MQ=255)
                    gAACTGGCGGCTTCCGGGGAGATCACCTTCTCGCCGTGTCACGAGCACGACTGCGTGGGa                                                          >  1:2492152/1‑60 (MQ=255)
                     aaCTGGCGGCTTCCGGGGAGATCACCTTCTCGCCGTGTCACGAGCACGa                                                                    >  1:1036764/1‑49 (MQ=255)
                      aCTGGCGGCTTCCGGGGAGATCACCTTCTCGCCGTGTCACGAGCACGACTGCGTGGGATCGATGGCGGgt                                              <  1:3863250/70‑1 (MQ=255)
                         ggcggcTTCCGGGGAGATCACCTTCTCGCCGTGTCACGAGCACGACTGCGTGGGATCGAt                                                     >  1:3201730/1‑60 (MQ=255)
                         ggcggcTTCCGGGGAGATCACCTTCTCGCCGTGTCACGAGCACGACTGCGTGGGATCGATGGCGGGTGTTa                                          >  1:1359204/1‑71 (MQ=255)
                          gcggcTTCCGGGGAGATCACCTTCTCGCCGTGTCACGAGCACGACTGCGTGGGATCGATGGCGGGTGtt                                           <  1:1492580/69‑1 (MQ=255)
                           cggcTTCCGGGGAGATCACCTTCTCGCCGTGTCACGAGCACGACTGCGTGGGATCGATGGCGGGTGTTAcc                                        <  1:4243662/71‑1 (MQ=255)
                               ttCCGGGGAGATCACCTTCTCGCCGTGTCACGAGCACGACTGCGTGGGATCGATGGCGGGTGTTACCTCgg                                    >  1:131469/1‑71 (MQ=255)
                                tCCGGGGAGATCACCTTCTCGCCGTGTCACGAGCAc                                                                      >  1:4189240/1‑36 (MQ=255)
                                  cGGGGAGATCACCTTCTCGCCGTGTCAAGAGCACGACTGCGTGGGATCGATGGCGGGTGTTACCTCGGCCt                                 >  1:4605858/1‑71 (MQ=255)
                                    gggAGATCACCTTCTCGCCGTGTCACGAGCACGACTGCGTGGGATCGATGGCGGGTGTTACCTCGGCCTCg                               <  1:197680/71‑1 (MQ=255)
                                      gagaTCACCTTCTCGCCGTGTCACGAGCGCGACTGCGTg                                                             <  1:1408755/39‑1 (MQ=255)
                                      gagaTCACCTTCTCGCCGTGTCACGAGCACGACTGCGTGGGATCTATGGCgg                                                >  1:4138647/1‑52 (MQ=255)
                                        gaTCACCTTCTCGCCGTGTCACGAGCGCGACTGCGTggg                                                           >  1:3045210/1‑39 (MQ=255)
                                        gaTCACCTTCTCGCCGTGTCACGAGCACGACTGCGTGGGATCGATGGCGGGTGTTACCTCGGCCTCGATGt                           >  1:519990/1‑71 (MQ=255)
                                         aTCACCTTCTCGCCGTGTCACGAGCACGACTGCGTGGGATCGa                                                      <  1:3217196/43‑1 (MQ=255)
                                         aTCACCTTCTCGCCGTGTCACGAGCACGACTGCGTGGGATCGATg                                                    <  1:365605/45‑1 (MQ=255)
                                         aTCACCTTCTCGCCGTGTCACGAGCACGACTGCGTGGGATCGATGGCGGGTGTTACCTCGGCCTCGATGtt                          <  1:2745427/71‑1 (MQ=255)
                                          tCACCTTCTCGCCGTGTCACGAGCACGACTGCGTGGGATCGATGGCGGGTGtt                                           <  1:2681560/53‑1 (MQ=255)
                                          tCACCTTCTCGCCGTGTACGAGCACGACTGCGTGGGATCGATGGCGGGTGTTACCTCGGCCTCGATGTTTa                        <  1:3265012/71‑1 (MQ=255)
                                           cACCTTCTCGCCGTGTCACGAGCACGACTGCGTGGGATCGATGGCGGGTGTTACCTc                                      >  1:2983871/1‑57 (MQ=255)
                                            aCCTTCTCGCCGTGTCACGAGCACGACTGCGTGGGATc                                                        >  1:4402399/1‑38 (MQ=255)
                                            aCCTTCTCGCCGTGTCACGAGCACGACTGCGTGGGATCGa                                                      >  1:2387896/1‑40 (MQ=255)
                                              cTTCTCGCCGTGTCACGAGCGCGACTGCGTGGGATCGATGGCGGGTGtt                                           <  1:2141817/49‑1 (MQ=255)
                                                tctcGCCGTGTCACGAGCACGACTGCGTGGGATCGATGGCGGGTGTTACCTCGGCCTCGATGTTTATGca                    >  1:1774186/1‑70 (MQ=255)
                                                tctcGCCGTGTCACGAGCACGACTGCGTGGGATCGATGGCGGGTGTTACCTCGGCCTCGATGTTTATGca                    >  1:1651782/1‑70 (MQ=255)
                                                tctcGCCGTGTCACGAGCACGACTGCGTGGGATCGATGGCGGGTGTTACCTCGGCCTCGATGTTTATGca                    >  1:3939465/1‑70 (MQ=255)
                                                 ctcGCCGTGTCACGAGCACGACTGCGTGGGATCGa                                                      >  1:310447/1‑35 (MQ=255)
                                                   cGCCGTGTCACGAGCACGACTGCGTGGGATCGATGGCGGGTGTTACCTCGGCCTCGATGTTTATGCACATc                >  1:4028555/1‑71 (MQ=255)
                                                    gCCGTGTAACGAGCGCGACTGCGTGGGATCGATGGCGGGTGTTACCt                                       >  1:904832/1‑47 (MQ=39)
                                                     ccGTGTCACGAGCACGACTGCGTGGGATCGATGGCGGGTGTTACCTCGGCCTCg                               <  1:124167/54‑1 (MQ=255)
                                                          tCACGAGCACGACTGCGTGGGATCGATGGCGGGTGTTACCTCGGCCTc                                <  1:1504624/48‑1 (MQ=255)
                                                          tCACGAGCACGACTGCGTGGGATCGATGGCGGGTGTTACCTCGGCCTCGATGTTTATGCACATCGTaaaaa         <  1:1980459/71‑1 (MQ=255)
                                                           cACGAGCACGACTGCGTGGGATCGATGGCGGGTGTTACCTCGGCCTCGATGTTTATGCACATCGTaaaaa         <  1:1785599/70‑1 (MQ=255)
                                                            aCGAGCGCGACTGCGTGGGATCGATGGCGGGTGTTAc                                         <  1:634255/37‑1 (MQ=255)
                                                            aCGAGCGCGACTGCGTGGGATCGATGGCGGGTGTTACCTCGGCCTCGATGTTTATGCACATCGTaaaaaa        >  1:846763/1‑70 (MQ=255)
                                                            aCGAGCACGACTGCGTGGGATCGATGGCGGGTGTTACCTCgg                                    >  1:248004/1‑42 (MQ=255)
                                                             cGAGCACTACTGCGTGGGATCGATGGCGGGTGTTACCTCGGCCTCGATGTTTAt                       <  1:275903/54‑1 (MQ=255)
                                                             cGAGCACGACTGCGTGGGATCGATGGCGGGTGTTAc                                         <  1:2729985/36‑1 (MQ=255)
                                                              gAGCACGACTGCGTGGGATCGATGGCGGGTGTTACCTCGGCCTCGATGTTTATGCACATCGTAAaaaaca      >  1:3872857/1‑70 (MQ=255)
                                                              gAGCACGACTGCGTGGGATCGATGGCGGGTGTTACATCGGCCTCGATGTTTATGCACATCGTAAaaaaca      >  1:2130335/1‑70 (MQ=255)
                                                               aGCACGACTGCGTGGGATCGATGGCGGGTGTTACCTCGGCCTCGATGTTTATGCACATCGTAAaaaacaaa    <  1:4619149/71‑1 (MQ=255)
                                                                 cACGACTGCGTGGGATCGATGGCGGGTGTTACCTCGGCCTCGATGTTTATGCACATCGTAAaaaacaaaac  >  1:2677839/1‑71 (MQ=255)
                                                                  |                                                                     
AAGATCTCGACGAGGCGGCTGAACTGGCGGCTTCCGGGGAGATCACCTTCTCGCCGTGTCACGAGCACGACTGCGTGGGATCGATGGCGGGTGTTACCTCGGCCTCGATGTTTATGCACATCGTAAAAAACAAAAC  >  W3110S.gb/337934‑338069

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: