Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,539,733 C→G 16.7% P13A (CCA→GCA)  yjhR → predicted frameshift suppressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,539,7330CG16.7% 83.3 / 9.5 36P13A (CCA→GCA) yjhRpredicted frameshift suppressor
Reads supporting (aligned to +/- strand):  ref base C (18/12);  new base G (4/2);  total (22/14)
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

AGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAACCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAG  >  W3110S.gb/4539665‑4539801
                                                                    |                                                                    
agagGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCa                                                                    <  1:178244/71‑1 (MQ=255)
         gAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATggagg                                                            >  1:294536/1‑70 (MQ=255)
          aaaaaGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATggaggga                                                          <  1:1022805/71‑1 (MQ=255)
           aaaaGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATggaggga                                                          <  1:894386/70‑1 (MQ=255)
            aaaGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGGCAAATg                                                                <  1:1712895/63‑1 (MQ=255)
              aGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGGCAAATgga                                                              >  1:1018565/1‑63 (MQ=255)
              aGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGGCAAATgga                                                              >  1:3676/1‑63 (MQ=255)
              aGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGGCAAATgga                                                              >  1:1778058/1‑63 (MQ=255)
              aGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGGCAAATgga                                                              >  1:1656954/1‑63 (MQ=255)
                       ttCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAAccc                                             >  1:633351/1‑71 (MQ=255)
                           ccGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGGCAAATg                                                                <  1:812258/48‑1 (MQ=255)
                             gCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAACCCTCTTGa                                       <  1:990210/71‑1 (MQ=255)
                                   ggATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAACCCTCTTGAAGCTg                                  >  1:1284690/1‑70 (MQ=255)
                                   ggATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAACCCTCTTGAAGCTg                                  >  1:424840/1‑70 (MQ=255)
                                   ggATATTTACATATTGATGCTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAAcc                                              >  1:1584870/1‑58 (MQ=255)
                                        tttACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAACCCTCTTGAAGCTGAAACaa                            >  1:638324/1‑71 (MQ=255)
                                        tttACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAACCCTCTTGAAGCTGAAACaa                            >  1:1318058/1‑71 (MQ=255)
                                          tACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAACCCTCTTGAAGCTGAAACaa                            >  1:6560/1‑69 (MQ=255)
                                                ttGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCa                                                   <  1:1067672/40‑1 (MQ=255)
                                                ttGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCa                                                   <  1:1439413/40‑1 (MQ=255)
                                                ttGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCa                                                   <  1:1803396/40‑1 (MQ=255)
                                                ttGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAACCCTCTTGAAGCTGAAACAATAgcagca                    <  1:1589302/71‑1 (MQ=255)
                                                       tAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAACCCTCTTGAAGCTGAAAc                              <  1:1471499/54‑1 (MQ=255)
                                                        agagGCATGAAGCCAAATGGAGGGAGCCGTCATAACCctct                                          >  1:1436372/1‑41 (MQ=255)
                                                        agagGCATGAAGCCAAATGGAGGGAGCCGTCATAACCCTCTTGAAGCTGAAACAATa                          >  1:1036516/1‑57 (MQ=255)
                                                        agagGCATGAAGCCAAATGGAGGGAGCCGTCATAACCCTCTTGAAGCTGAAACAATAGCAGCATGGCtggt            >  1:145074/1‑71 (MQ=255)
                                                        agagGCATGAAGCCAAATGGAGGGAGCCGTCATAACCCTCTTGAAGCTGAAACAATAGCAGCATGGCtgg             >  1:1237386/1‑70 (MQ=255)
                                                              aTGAAGCCAAATGGAGGGAGCCGTCATAACCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGcaca      >  1:1752220/1‑71 (MQ=255)
                                                                 aaGCCAAATGGAGGGAGCCGTCATAACCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACAt     >  1:571662/1‑69 (MQ=255)
                                                                 aaGCCAAATGGAGGGAGCCGTCATAACCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACAt     >  1:1723342/1‑69 (MQ=255)
                                                                 aaGCCAAATGGAGGGAGCCGTCATAACCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACAt     >  1:1182620/1‑69 (MQ=255)
                                                                 aaGCCAAATGGAGGGAGCCGTCATAACCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACAt     >  1:856649/1‑69 (MQ=255)
                                                                 aaGCCAAATGGAGGGAGCCGTCATAACCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACAt     >  1:1689611/1‑69 (MQ=255)
                                                                   gCCAAATGGAGGGAGCCGTCATAACCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAg  <  1:371682/70‑1 (MQ=255)
                                                                   gCCAAATGGAGGGAGCCGTCATAACCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAg  <  1:369534/70‑1 (MQ=255)
                                                                   gCCAAATGGAGGGAGCCGTCATAACCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAg  <  1:603276/70‑1 (MQ=255)
                                                                    |                                                                    
AGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAACCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAG  >  W3110S.gb/4539665‑4539801

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: