Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,860,149 G→T 33.3% A22S (GCA→TCA)  ygbJ → predicted dehydrogenase, with NAD(P)‑binding Rossmann‑fold domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,860,1490GT33.3% 83.1 / 29.4 45A22S (GCA→TCA) ygbJpredicted dehydrogenase, with NAD(P)‑binding Rossmann‑fold domain
Reads supporting (aligned to +/- strand):  ref base G (14/16);  new base T (8/7);  total (22/23)
Fisher's exact test for biased strand distribution p-value = 7.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGAAAACGGGATCTGAGTTTCATGTCGGTATCGTTGGCTTAGGGTCAATGGGAATGGGAGCAGCACTGTCATATGTCCGCGCAGGTCTTTCTACCTGGGGCGCAGACCTGAACAGCAATGC  >  W3110S.gb/2860087‑2860207
                                                              |                                                          
tGAAAACGGGATCTGAGTTTCATGTCGGTATCGTTGGCTTAGGGTCAATGGGAATGGGAGCATCACTGTCa                                                    <  1:2870942/71‑1 (MQ=255)
tGAAAACGGGATCTGAGTTTCATGTCGGTATCGTTGGCTTAGGGTCAATGGGAATGGGAGCAGCACTGTCa                                                    <  1:941181/71‑1 (MQ=255)
tGAAAACGGGATCTGAGTTTCATGTCGGTATCGTTGGCTTAGGGTCAATGGGAATGGGAGCAGCACTGTCa                                                    <  1:1335042/71‑1 (MQ=255)
 gAAAACGGGATCTGAGTTTCATGTCGGTATCGTTGGCTTAGGGTCAATGGGAATGGGAgcagca                                                          >  1:3146235/1‑64 (MQ=255)
 gAAAACGGGATCTGAGTTTCATGTCGGTATCGTTGGCTTAGGGTCAATGGGAATGGGAgcagca                                                          >  1:2093760/1‑64 (MQ=255)
 gAAAACGGGATCTGAGTTTCATGTCGGTATCGTTGGCTTAGGGTCAATGGGAATGGGAgcagca                                                          >  1:1587331/1‑64 (MQ=255)
 gAAAACGGGATCTGAGTTTCATGTCGGTATCGTTGGCTTAGGGTCAATGGGAATGGGAGCAGCACTGTCa                                                    >  1:1999870/1‑70 (MQ=255)
 gAAAACGGGATATGAGTTTCATGTCGGTATCGTTGGCTTAGGGTCAATGGGAATGGGAgcagca                                                          >  1:1385880/1‑64 (MQ=255)
  aaaaCGGGATCTGAGTTTCATGTCGGTATCGTTGGCTTAGGGTCAATGGGAATGGGAGCATCACTGTCata                                                  >  1:3554122/1‑71 (MQ=255)
  aaaaCGGGATCTGAGTTTCATGTCGGTATCGTTGGCTTAGGGTCAATGGGAATGGGAGCATCACTGTCata                                                  >  1:2305871/1‑71 (MQ=255)
  aaaaCGGGATCTGAGTTTCATGTCGGTATCGTTGGCTTAGGGTCAATGGGAATGGGAGCATCACTGTCat                                                   <  1:2356533/70‑1 (MQ=255)
      cGGGATCTGAGTTTCATGTCGGTATCGTTGGCTTAGGGTCAATGGGAATGGGAGCATCACTGTCATATGt                                               >  1:1787432/1‑70 (MQ=255)
      cGGGATCTGAGTTTCATGTCGGTATCGTTGGCTTAGGGTCAATGGGAATGGGAGCATCACTGTCATATGTc                                              >  1:3593994/1‑71 (MQ=255)
          aTCTGAGTTTCATGTCGGTATCGTTGGCTTAGGGTCAATGGGAATGGGAgcatca                                                          >  1:997437/1‑55 (MQ=255)
              gAGTTTCATGTCGGTATCGTTGGCTTAGGGTCAATGGGAATGGGAGCATCACTGTCATATg                                                <  1:2531469/61‑1 (MQ=255)
                 tttCATGTCGGTATCGTTGGCTTAGGGTCAATGGGAATGGGAGCAGCACTGTCATATGTCCGCGCAGGTCt                                   <  1:956028/71‑1 (MQ=255)
                 tttCATGTCGGTATCGTTGGCTTAGGGTCAATGGGAATGGGAGCAGCACTGTCATATGTCCGCGCAGGTCt                                   <  1:1531602/71‑1 (MQ=255)
                 tttCATGTCGGTATCGTTGGCTTAGGGTCAATGGGAATGGGAGCAGCACTGTCATATGTCCGCGCAGGTCt                                   <  1:748801/71‑1 (MQ=255)
                   tCATGTCGGTATCGTTGGCTTAGGGTCAATGGGAATGGGAgcagca                                                          >  1:1056570/1‑46 (MQ=255)
                   tCATGTCGGTATCGTTGGCTTAGGGTCAATGGGAATGGGAgcagca                                                          >  1:1608859/1‑46 (MQ=255)
                   tCATGTCGGTATCGTTGGCTTAGGGTCAATGGGAATGGGAGCAGCAc                                                         >  1:1188421/1‑47 (MQ=255)
                   tCATGTCGGTATCGTTGGCTTAGGGTCAATGGGAATGGGAGCAGCAc                                                         >  1:3055827/1‑47 (MQ=255)
                   tCATGTCGGTATCGTTGGCTTAGGGTCAATGGGAATGGGAGCAGCAc                                                         >  1:2571609/1‑47 (MQ=255)
                    cATGTCGGTATCGTTGCCTTAGGGTCAATGGGAATGGGAGCAGCACTGTCATATGTCCGCGCAGGTCttt                                 <  1:3209543/70‑1 (MQ=255)
                     aTGTCGGTATCGTTGGCTTAGGGTCAATGGGAATGGGAGCAGCACTGTCATATGTCCGCGCAGGTCTTTc                                <  1:3834676/70‑1 (MQ=255)
                          ggTATCGTTGGCTTAGGGTCAATGGGAATGGGTGCATCACTGTCat                                                   <  1:1465252/46‑1 (MQ=255)
                          ggTATCGTTGGCTTAGGGTCAATGGGAATGGGAGCATCACTGTCat                                                   <  1:878256/46‑1 (MQ=255)
                           gTATCGTTGGCTTAGGGTCAATGGGAATGGGAGCATCACTGTCat                                                   <  1:185889/45‑1 (MQ=255)
                            tATCGTTGGCTTAGGGTCAATGGGAATGGGAGCAGCACt                                                        <  1:3065511/39‑1 (MQ=255)
                            tATCGTTGGCTTAGGGTCAATGGGAATGGGAGCAGCACt                                                        <  1:2280368/39‑1 (MQ=255)
                                gTTGGCTTAGGGTCAATGGGAATGGGAGCATCACTGTCATATGTCCGCGCAg                                       >  1:2286914/1‑52 (MQ=255)
                                gTTGGCTTAGGGTCAATGGGAATGGGAGCATCACTGTCATATGTCCGCGCAGGTCTTTCTACCTggg                        >  1:3116342/1‑67 (MQ=255)
                                gTTGGCTTAGGGTCAATGGGAATGGGAGCAGCACTGTCata                                                  >  1:1088765/1‑41 (MQ=255)
                                gTTGGCTTAGGGTCAATGGGAATGGGAGCAGCACTGTCATATGTCCGCGCAGGTCt                                   <  1:674641/56‑1 (MQ=255)
                                gTTGGCTTAGGGTCAATGGGAATGGGAGCAGCACTGTCATATGTCCGCGCAGGTCTTTCTAc                             >  1:192534/1‑62 (MQ=255)
                                gTTGGCTTAGGGTCAATGGGAATGGGAGCAGCACTGTCATATGTCCGCGCAGGTCTTTCTAc                             >  1:3987345/1‑62 (MQ=255)
                                gTTGGCTTAGGGTCAATGGGAATGGGAGCAGCACTGTCATATGTCCGCGCAGGTCTTTCTACCTGGGgcg                     <  1:992234/70‑1 (MQ=255)
                                gTTGGCTTAGGGTCAATGGGAATGGGAGCAGCACTGTCATATGTCCGCGCAGGTCTTTCTACCTGGGgcg                     >  1:680623/1‑70 (MQ=255)
                                gTTGGCTTAGGGTCAATGGGAATGGGAGCAGCACTGTCATATGTCCGCGCAGGTCTTTCTACCTGGGgcg                     >  1:2434662/1‑70 (MQ=255)
                                 ttGGCTTAGGGTCAATGGGAATGGGAGCAGCACTGTCATATGTCCGCGCAGGTCTTTCTACCTGGGCGCAg                  <  1:2161246/71‑1 (MQ=255)
                                   ggCTTAGGGTCAATGGGAATGGGAGCAGCACTGTCATATGTCCGCGCAGGTCTTTCTACCTGGGGCGCAg                  <  1:1209010/70‑1 (MQ=255)
                                   ggCTTAGGGTCAATGGGAATGGGAGCAGCACTGTCATATGTCCGCGCAGGTCTTTCTACCTGGGGCGCAg                  <  1:3973728/70‑1 (MQ=255)
                                   ggCTTAGGGTCAATGGGAATGGGAGCAGCACTGTCATATGTCCGCGCAGGTCTTTCTACCTGGGGCGCAGa                 >  1:3918666/1‑71 (MQ=255)
                                        aGGGTCAATGGGAATGGGAGCAGCACTGTCATATGTCCGCGCAGGTCTTTCTACCTGGGGCGCAGACCTg             <  1:432919/70‑1 (MQ=255)
                                         gggTCAATGGGAATGGGAGCAGCACTGTCATATGTCCGCGCAGGt                                     >  1:1757822/1‑45 (MQ=255)
                                         gggTCAATGGGAATGGGAGCAGCACTGTCATATGTCCGCGCAGGt                                     >  1:3692568/1‑45 (MQ=255)
                                              aatgggaatgggAGCAGCACTGTCATATGTCCGCGCAGGTCTTTCTACCTGGGGCGCAGACCTGAACAGc       <  1:1045389/70‑1 (MQ=255)
                                               atgggaatgggaGCATCACTGTCATATGTCCGCGCAGGTCTTTCTACCTGGGGCGCAGACCTGAACAGCaa     <  1:130203/71‑1 (MQ=255)
                                                 gggaatgggaGCAGCACTGTCATATGTCCGCGCAGGTCTTTCTACCTGGGGCGCAGACCTGAACAGCAATg   >  1:1430763/1‑71 (MQ=255)
                                                  ggaatgggaGCATCACTGTCATATGTCCGCGCAGGTCTTTCTACCTGGGGCGCAGACCTGAACAGCAATGc  >  1:1255845/1‑71 (MQ=255)
                                                  ggaatgggaGCAGCACTGTCATATGTCCGCGCAGGTCTTTCTACCTGGGGCGCAGACCTGAACAGCAATGc  >  1:354604/1‑71 (MQ=255)
                                                     atgggaGCATCACTGTCATATGTCCGCGCAGGTCTTTCTAc                             >  1:2649214/1‑41 (MQ=255)
                                                              |                                                          
TGAAAACGGGATCTGAGTTTCATGTCGGTATCGTTGGCTTAGGGTCAATGGGAATGGGAGCAGCACTGTCATATGTCCGCGCAGGTCTTTCTACCTGGGGCGCAGACCTGAACAGCAATGC  >  W3110S.gb/2860087‑2860207

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: