Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 2,860,110 G→A V9I (GTC→ATC)  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,1100GA76.4% 102.6 / 32.2 55V9I (GTC→ATC) ygbJpredicted dehydrogenase, with NAD(P)‑binding Rossmann‑fold domain
Reads supporting (aligned to +/- strand):  ref base G (5/8);  new base A (3/39);  total (8/47)
Fisher's exact test for biased strand distribution p-value = 1.34e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACAATATCAAACAAAATATCACTTAAATTAACAAGGAGAGCAGATGAAAACGGGATCTGAGTTTCATGTCGGTATCGTTGGCTTAGGGTCAATGGGAATGGGAGCAGCACTGTCATATGTCCGCGCAGGT  >  W3110S.gb/2860043‑2860172
                                                                   |                                                              
acaATATCAAACAAAATATCACTTAAATTAACAAGGAGAGCAGATGAAAACGGGATCTGAGTTTCATGTc                                                              <  1:958810/70‑1 (MQ=255)
acaATATCAAACAAAATATCACTTAAATTAACAAGGAGAGCAGATGAAAACGGGATCTGAGTTTCATATc                                                              <  1:984958/70‑1 (MQ=255)
acaATATCAAACAAAATATCACTTAAATTAACAAGGAGAGCAGATGAAAACGGGATCTGAGTTTCATATc                                                              <  1:2768981/70‑1 (MQ=255)
acaATATCAAACAAAATATCACTTAAATTAACAAGGAGAGCAGATGAAAACGGGATCTGAGTTTCATATc                                                              <  1:367891/70‑1 (MQ=255)
   atatCAAACAAAATATCACTTAAATTAACAAGGAGAGCAGATGAAAACGGGATCTGAGTTTCATATCGGTa                                                          <  1:2424547/71‑1 (MQ=255)
      tCAAACAAAATATCACTTAAATTAACAAGGAGAGCAGATGAAAACGGGATCTGAGTTTCATATCGGTATc                                                        <  1:98378/70‑1 (MQ=255)
      tCAAACAAAATATCACTTAAATTAACAAGGAGAGCAGATGAAAACGGGATCTGAGTTTCATATCGGTATc                                                        <  1:1756478/70‑1 (MQ=255)
      tCAAACAAAATATCACTTAAATTAACAAGGAGAGCAGATGAAAACGGGATCTGAGTTTCATATCGGTATc                                                        <  1:373583/70‑1 (MQ=255)
        atacaaaATATCACTTAAATTAACAAGGAGAGCAGATGAAAACGGGATCTGAGTTTCATATCGGTATc                                                        <  1:3961127/66‑1 (MQ=255)
        aaacaaaATATCACTTAAATTAACAAGGAGAGCAGATGAAAACGGGATCTGAGTTTCATATCGGTATCGtt                                                     <  1:1080857/71‑1 (MQ=255)
              aaTATCACTTAAATTAACAAGGAGAGCAGATGAAAACGGGATCTGAGTTTCATATCGGTATCGTTGGCTTa                                               <  1:1330526/71‑1 (MQ=255)
              aaTATCACTTAAATTAACAAGGAGAGCAGATGAAAACGGGATCTGAGTTTCATATCGGTATCGTTGGCTTa                                               <  1:337701/71‑1 (MQ=255)
              aaTATCACTTAAATTAACAAGGAGAGCAGATGAAAACGGGATCTGAGTTTCATATCGGTATCGTTGGCTTa                                               <  1:2161782/71‑1 (MQ=255)
              aaTATCACTTAAATTAACAAGGAGAGCAGATGAAAACGGGATCTGAGTTTCATATCGGTATCGTTGGCTTa                                               <  1:2069597/71‑1 (MQ=255)
              aaTATCACTTAAATTAACAAGGAGAGCAGATGAAAACGGGATCTGAGTTTCATATCGGTATCGTTGGCTTa                                               <  1:2783176/71‑1 (MQ=255)
              aaTATCACTTAAATTAACAAGGAGAGCAGATGAAAACGGGATCTGAGTTTCATATCGGTATCGTTGGCTTa                                               <  1:4119261/71‑1 (MQ=255)
              aaTATCACTTAAATTAACAAGGAGAGCAGATGAAAACGGGATCTGAGTTTCATATCGGTATCGTTGGCTTa                                               <  1:767486/71‑1 (MQ=255)
              aaTATCACTTAAATTAACAAGGAGAGCAGATGAAAACGGGAGCTGAGTTTCATATCGGTATCGTTGGCTTa                                               <  1:1867761/71‑1 (MQ=255)
               atatCACTTAAATTAACAAGGAGAGCAGATGAAAACGGGATCTGAGTTTCATGTCGGTATCGTTGGCTTAg                                              <  1:179944/71‑1 (MQ=255)
               atatCACTTAAATTAACAAGGAGAGCAGATGAAAACGGGATCTGAGTTTCATGTCGGTATCGTTGGCTTAg                                              <  1:977983/71‑1 (MQ=255)
               atatCACTTAAATTAACAAGGAGAGCAGATGAAAACGGGATCTGAGTTTCATGTCGGTATCGTTGGCTTAg                                              <  1:1773244/71‑1 (MQ=255)
               atatCACTTAAATTAACAAGGAGAGCAGATGAAAACGGGATCTGAGTTTCATGTCGGTATCGTTGGCTTAg                                              <  1:872419/71‑1 (MQ=255)
               atatCACTTAAATTAACAAGGAGAGCAGATGAAAACGGGATCTGAGTTTCATGTCGGTATCGTTGGCTTAg                                              <  1:901232/71‑1 (MQ=255)
                tatCACTTAAATTAACAAGGAGAGCAGATGAAAACGGGATCTGAGTTTCATGTCGGTATCGTTGGCTTAg                                              <  1:3865540/70‑1 (MQ=255)
                tatCACTTAAATTAACAAGGAGAGCAGATGAAAACGGGATCTGAGTTTCATATCGGTATCGTTGGCTTAgg                                             <  1:1869845/71‑1 (MQ=255)
                    aCTTAAATTAACAAGGAGAGCAGATGAAAACGGGATCTGAGTTTCATATCGGTATCGTTGGCTTAGGGTc                                          <  1:242936/70‑1 (MQ=255)
                    aCTTAAATTAACAAGGAGAGCAGATGAAAACGGGATCTGAGTTTCATATCGGTATCGTTGGCTTAGGGTc                                          <  1:3166286/70‑1 (MQ=255)
                     cTTAAATTAACAAGGAGAGCAGATGAAAACGGGATCTGAGTTTCATATCGGTATCGTTGGCTTAGGGTc                                          <  1:2835235/69‑1 (MQ=255)
                         aaTTAACAAGGAGAGCAGATGAAAACGGGATCTGAGTTTCATATCGGTATCGTTGGCTTAGGGTCaatggg                                    <  1:4108481/71‑1 (MQ=255)
                         aaTTAACAAGGAGAGCAGATGAAAACGGGATCTGAGTTTCATATCGGTATCGTTGGCTTAGGGTCaatggg                                    <  1:2612730/71‑1 (MQ=255)
                         aaTTAACAAGGAGAGCAGATGAAAACGGGATCTGAGTTTCATATCGGTATCGTTGGCTTAGGGTCaatggg                                    <  1:1231815/71‑1 (MQ=255)
                         aaTTAACAAGGAGAGCAGATGAAAACGGGATCTGAGTTTCATATCGGTATCGTTGGCTTAGGGTCaatggg                                    <  1:1450730/71‑1 (MQ=255)
                         aaTTAACAAGGAGAGCAGATGAAAACGGGATCTGAGTTTCATATCGGTATCGTTGGCTTAGGGTCaatggg                                    <  1:1653206/71‑1 (MQ=255)
                          aTTAACAAGGAGAGCAGATGAAAACGGGATCTGAGTTTCATATCGGTATCGTTGGCTTAGGGTCgatggg                                    <  1:1997070/70‑1 (MQ=255)
                               cAAGGAGAGCAGATGAAAACGGGATCTGAGTTTCATATCGGTATCGTTGGCTTAGGGTCaatgggaatggg                              <  1:1135599/71‑1 (MQ=255)
                                 aGGAGAGCAGATGAATACGGGATCTGAGTTTCATGTCGGTATCGTTGGCTTAg                                              <  1:2304226/53‑1 (MQ=255)
                                 aGGAGAGCAGATGAAAACGGGATCTGAGTTTCATATCGGTATCGTTGGCTTAGGGTCaatgggaatggg                              <  1:2706725/69‑1 (MQ=255)
                                       gCAGATGAAAACGGGATCTGAGTTTCATGTCGGTATCGTTGGCTTAGGGTCaatgggaa                                  >  1:2976509/1‑59 (MQ=255)
                                       gCAGATGAAAACGGGATCTGAGTTTCATGTCGGTATCGTTGGCTTAGGGTCaatgggaa                                  >  1:2256513/1‑59 (MQ=255)
                                       gCAGATGAAAACGGGATCTGAGTTTCATGTCGGTATCGTTGGCTTAGGGTCaatgggaa                                  >  1:2086265/1‑59 (MQ=255)
                                       gCAGATGAAAACGGGATCTGAGTTTCATGTCGGTATCGTTGGCTTAGGGTCaatgggaa                                  >  1:2574684/1‑59 (MQ=255)
                                       gCAGATGAAAACGGGATCTGAGATTCATGTCGGTATCGTTgg                                                   >  1:3293763/1‑42 (MQ=255)
                                         aGATGAAAACGGGATCTGAGTTTCATATCGGTATCGtt                                                     <  1:739902/38‑1 (MQ=255)
                                                   gggATCTGAGTTTCATATCGGTATCGTTGGCTTAGGGTCaatggg                                    <  1:2599528/45‑1 (MQ=255)
                                                           aGTTTCATATCGGTATCGTTGGCTTAGGGTCAATGGGAATGGGAGCAGCACTGTCATATGTCCGCGCAgg   <  1:3807864/70‑1 (MQ=255)
                                                           aGTTTCATATCGGTATCGTTGGCTTAGGGTCAATGGGAATGGGAGCAGCACTGTCATATGTCCGCGCAgg   <  1:1736618/70‑1 (MQ=255)
                                                           aGTTTCATATCGGTATCGTTGGCTTAGGGTCAATGGGAATGGGAGCAGCACTGTCATATGTCCGCGCAgg   <  1:3346834/70‑1 (MQ=255)
                                                           aGTTTCATATCGGTATCGTTGGCTTAGGGTCAATGGGAATGGGAGCAGCACTGTCATATGTCCGCGCAgg   <  1:2942029/70‑1 (MQ=255)
                                                           aGTTTCATATCGGTATCGTTGGCTTAGGGTCAATGGGAATGGGAGCAGCACTGTCATATGTCCGCGCAgg   <  1:1436061/70‑1 (MQ=255)
                                                           aGTTTCATATCGGTATCGTTGGCTTAGGGTCAATGGGAATGGGAGCAGCACTGTCATATGTCCGCGCAGGt  <  1:1918762/71‑1 (MQ=255)
                                                           aGTTTCATATCGGTATCGTTGGCTTAGGGTCAATGGGAATGGGAGCAGCACTGTCATATGTCCGCGCAGGt  <  1:3961377/71‑1 (MQ=255)
                                                           aGTTTCATATCGGTATCGTTGGCTTAGGGTCAATGGGAATGGGAGCAGCACTGTCATATGTCCGCGCAGGt  <  1:1764156/71‑1 (MQ=255)
                                                                cATATCGGTATCGTTGGCTTAGGGTCAATGGGAATGGGagcagc                        >  1:3191608/1‑44 (MQ=255)
                                                                cATATCGGTATCGTTGGCTTAGGGTCAATGGGAATGGGAGCAGCAc                      >  1:186213/1‑46 (MQ=255)
                                                                cATATCGGTATCGTTGGCTTAGGGTCAATGGGAATGGGAGCAGCAc                      >  1:1529374/1‑46 (MQ=255)
                                                                   |                                                              
ACAATATCAAACAAAATATCACTTAAATTAACAAGGAGAGCAGATGAAAACGGGATCTGAGTTTCATGTCGGTATCGTTGGCTTAGGGTCAATGGGAATGGGAGCAGCACTGTCATATGTCCGCGCAGGT  >  W3110S.gb/2860043‑2860172

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: