Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,441,198 G→A 15.2% T434I (ACC→ATC)  ydbK ← fused predicted Fe‑S subunit of pyruvate‑flavodoxin oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,441,1980GA15.2% 152.9 / 17.1 59T434I (ACC→ATC) ydbKfused predicted Fe‑S subunit of pyruvate‑flavodoxin oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (33/17);  new base A (4/5);  total (37/22)
Fisher's exact test for biased strand distribution p-value = 2.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTAAACAAAATAGCCCTGTGCGTACCACGGCGTGGAATTACCGATAATCTTGATATTGTTTTTGGTCGCGGAAACGCTGCCATCACTACCAAGGCCATAAAACAAGGCTTCCAGTTTCGCCGAGTT  >  W3110S.gb/1441134‑1441259
                                                                |                                                             
gTAAACAAAATAGCCCTGTGCGTACCACGGCGTGGAATTACCGATAATCTTGATATTGTTTTTGGTcg                                                            >  1:105373/1‑68 (MQ=255)
gTAAACAAAATAGCCCTGTGCGTACCACGGCGTGGAATTACCGATAATCTTGATATTGTTTTTGGTcg                                                            >  1:3519877/1‑68 (MQ=255)
   aaCAAAATAGCCCTGTGCGTACCACGGCGTGGAATTACCGATAATCTTGATATTGTTTTTGGTCGCGGa                                                        >  1:837584/1‑69 (MQ=255)
     cAAAATAGCCCTGTGCGTACCACGGCGTGGAATTACCGATAATCTTGATATTGTTTTTGGTCGCGGaaa                                                      <  1:2554859/69‑1 (MQ=255)
     cAAAATAGCCCTGTGCGTACCACGGCGTGGAATTACCGATAATCTTGATATTGTTTTTGGTCGCGGaaa                                                      >  1:2836271/1‑69 (MQ=255)
     cAAAATAGCCCTGTGCGTACCACGGCGTGGAATTACCGATAATCTTGATATTGTTTTTGGTCGCGGaaa                                                      <  1:3341276/69‑1 (MQ=255)
     cAAAATAGCCCTGTGCGTACCACGGCGTGGAATTACCGATAATCTTGATATTGTTTTTGGTCGCGGaaa                                                      <  1:1396168/69‑1 (MQ=255)
      aaaaTAGCCCTGTGCGTACCACGGCGTGGAATTACCGATAATCTTGATATTGTTTTTGGTCGCGGAAAc                                                     <  1:1490316/69‑1 (MQ=255)
      aaaaTAGCCCTGTGCGTACCACGGCGTGGAATTACCGATAATCTTGATATTGTTTTTGGTCGCGGAAAc                                                     <  1:2637613/69‑1 (MQ=255)
      aaaaTAGCCCTGTGCGTACCACGGCGTGGAATTACCGATAATCTTGATATTGTTTTTGGTCGCGGAAAc                                                     <  1:3019951/69‑1 (MQ=255)
              ccTCTGCGTACCACGGCGTGGAATTACCGATAATCTTGATATTGTTTTTGGTCGCGGAAAAGCTGCCAt                                             <  1:2187757/69‑1 (MQ=255)
                  tgCGTACCACGGCGTGGAATTACTGATAATCTTGATATTGTTTTTGGTCGCGGAAACGCTGCCATCACt                                         >  1:3683372/1‑69 (MQ=255)
                  tgCGTACCACGGCGTGGAATTACCGATAATCTTGATATTGTTTTTGGTCGCGGAAACGCTGc                                                >  1:516765/1‑62 (MQ=255)
                  tgCGTACCACGGCGTGGAATTACCGATAATCTTGATATTGTTTTTGGTCGCGGAAACGCTGc                                                <  1:2556450/62‑1 (MQ=255)
                  tgCGTACCACGGCGTGGAATTACCGATAATCTTGATATTGTTTTTGGTCGCGGAAACGCTGc                                                <  1:3060997/62‑1 (MQ=255)
                  tgCGTACCACGGCGTGGAATTACCGATAATCTTGATATTGTTTTTGGTCGCGGAAACGCTGc                                                >  1:327603/1‑62 (MQ=255)
                  tgCGTACCACGGCGTGGAATTACCGATAATCTTGATATTGTTTTTGGTCGCGGAAACGCTGc                                                <  1:699090/62‑1 (MQ=255)
                  tgCGTACCACGGCGTGGAATTACCGATAATCTTGATATTGTTTTTGGTCGCGGAAACGCTGc                                                >  1:1999912/1‑62 (MQ=255)
                  tgCGTACCACGGCGTGGAATTACCGATAATCTTGATATTGTTTTTGGTCGCGGAAACGCTGc                                                <  1:692830/62‑1 (MQ=255)
                  tgCGTACCACGGCGTGGAATTACCGATAATCTTGATATTGTTTTTGGTCGCGGAAACGCTGc                                                >  1:1917897/1‑62 (MQ=255)
                  tgCGTACCACGGCGTGGAATTACCGATAATCTTGATATTGTTTTTGGTCGCGGAAACGCTGc                                                <  1:3605155/62‑1 (MQ=255)
                  tgCGTACCACGGCGTGGAATTACCGATAATCTTGATATTGTTTTTGGTCGCGGAAACGCTGc                                                >  1:514940/1‑62 (MQ=255)
                  tgCGTACCACGGCGTGGAATTACCGATAATCTTGATATTGTTTTTGGTCGCGGAAACGCTGc                                                >  1:1506772/1‑62 (MQ=255)
                  tgCGTACCACGGCGTGGAATTACCGATAATCTTGATATTGTTTTTGGTCGCGGAAACGCTGc                                                >  1:574745/1‑62 (MQ=255)
                  tgCGTACCACGGCGTGGAATTACCGATAATCTTGATATTGTTTTTGGTCGCGGAAACGCTGCCATCACt                                         >  1:3295100/1‑69 (MQ=255)
                  tgCGTACCACGGCGTGGAATTACCGATAATCTTGATATTGTTTTTGGTCGCGGAAACGCTGCCATCACt                                         >  1:1772071/1‑69 (MQ=255)
                  tgCGTACCACGGCGTGGAATTACCGATAATCTTGATATTGTTTTTGGTCGCGGAAACGCTGCCATCACt                                         >  1:1642043/1‑69 (MQ=255)
                  tgCGTACCACGGCGTGGAATTACCGATAATCTTGATATTGTTTTTGGTCGCGGAAACGCTGCCATCACt                                         >  1:699216/1‑69 (MQ=255)
                    cGTACCACGGCGTGGAATTACCGATAATCTTGATATTGTTTTTGGTCGCGGAAACGCTGCCAt                                             <  1:1029606/63‑1 (MQ=255)
                    cGTACCACGGCGTGGAATTACCGATAATCTTGATATTGTTTTTGGTCGCGGAAACGCTGCCAt                                             <  1:1606324/63‑1 (MQ=255)
                    cGTACCACGGCGTGGAATTACCGATAATCTTGATATTGTTTTTGATCGCGGa                                                        <  1:1952481/52‑1 (MQ=255)
                          aCGGCGTGGAATTACCGATAATCTTGATATTGTTTTTGGt                                                              >  1:1459205/1‑40 (MQ=255)
                          aCGGCGTGGAATTACCGATAATCTTGATATTGTTTTTGGTCGCGGa                                                        >  1:624384/1‑46 (MQ=255)
                          aCGGCGTGGAATTACCGATAATCTTGACATTGTTTTTGGTCGCGGa                                                        >  1:556252/1‑46 (MQ=255)
                               gTGGAATTACCGATAATCTTGATATTGTTTTTGGTCGCGG‑‑ACGCTGCCATCACTACCAAGGCCATaaaa                          <  1:2997521/69‑1 (MQ=255)
                                 ggAATTACCGATAATCTTGATATTGTTTTTGGTCGCGGAAACGCTGCCATCACTACCAAGGCCATaaaa                          <  1:10221/69‑1 (MQ=255)
                                  gAATTACCGATAATCTTGATATTGTTTTTGGTCGCGGAAACGCTGCCATc                                            >  1:2305654/1‑50 (MQ=255)
                                  gAATTACCGATAATCTTGATATTGTTTTTGGTCGCGGAAACGCTGCCATc                                            >  1:1240471/1‑50 (MQ=255)
                                  gAATTACCGATAATCTTGATATTGTTTTTGATCGCGGAAACGCTGCCATCACTACCAAGGCCATaaa                           >  1:2216129/1‑67 (MQ=255)
                                  gAATTACCGATAATCTTGATATTGTTTTTGATCGCGGAAACGCTGCCATCACTACCAAGGCCATaaa                           >  1:2337802/1‑67 (MQ=255)
                                   aaTTACCGATAATCTTGATATTGTTTTTGGTCGCGGAAACGCTGCCATCACTACCAAGGCCATAAAACa                        >  1:1039517/1‑69 (MQ=255)
                                   aaTTACCGATAATCTTGATATTGTTTTTGGTCGCGGAAACGCTGCCATCACTACCAAGGCCATAAAACa                        >  1:2527279/1‑69 (MQ=255)
                                    aTTACCGATAATCTTGATATTGTTTTTGATCGCGGAAACGCTGCCATCACTACCAAGGCCATAAAACaa                       <  1:634260/69‑1 (MQ=255)
                                    aTTACCGATAATCTTATATTGTTTCTGATCGCGGAAACGCTGCCATCACTACCAAGGCCATAAAACaa                       <  1:3278533/68‑1 (MQ=255)
                                     ttACCGATAATCTTGATATTGTTTTTGGTCGCGGAAACGCTGCCATCACTAcc                                      <  1:2300401/53‑1 (MQ=255)
                                     ttACCGATAATCTTGATATTGTTTTTGATCGCGGAAACGCTGCCATCACTACCAAGGCCATAAAACaa                       <  1:894463/68‑1 (MQ=255)
                                       aCCGATAATCTTGATATTGTTTTTGGTCGCGGAAACGCTGCCATCACTACCAAGGCCATAAAACAAGGc                    >  1:1031374/1‑69 (MQ=255)
                                             aaTCTTGATATTGTTTTTGGTCGCGGAAACGCTGc                                                >  1:1828837/1‑35 (MQ=255)
                                             aaTCTTGATATTGTTTTTGGTCGCGGAAACGCTGc                                                >  1:2359476/1‑35 (MQ=255)
                                              aTCTTGATATTGTTTTTGGTCGCGGAAACGCTGCCATCACTACCAAGGCCATAAAACAAGGCTTCCAg              >  1:1088659/1‑68 (MQ=255)
                                              aTCTTGATATTGTTTTTGGTCGCGGAAACGCTGCCATCACTACCAAGGCCATAAAACAAGGCTTCCAg              >  1:736749/1‑68 (MQ=255)
                                                 ttGATATTGTTTTTGGTCGCGGAACCGCTGCCATCACTACCAAGGCCATaa                            >  1:2859344/1‑51 (MQ=255)
                                                 ttGATATTGTTTTTGGTCGCGGAAACGCTGCCATCACTACCAAGGCCATaa                            >  1:1810754/1‑51 (MQ=255)
                                                 ttGATATTGTTTTTGGTCGCGGAAACGCTGCCATCACTACCAAGGCCATaa                            >  1:2266899/1‑51 (MQ=255)
                                                 ttGATATTGTTTTTGGTCGCGGAAACGCTGCCATCACTACCAAGGCCATaa                            >  1:3175238/1‑51 (MQ=255)
                                                       ttGTTTTTGATCGCGGAAACGCTGCCATCACTACCAAGGCCATAAAACAAGGCTTCCAGTTTCGCCGAg    <  1:2337249/69‑1 (MQ=255)
                                                         gTTTTTGGTCGCGGAAACGCTGCCATCACTACCAAGGCCATAAAACAAGGCTTCCAGTTTCGCCGAGtt  >  1:3171873/1‑69 (MQ=255)
                                                         gTTTTTGATCGCGGAAACGCTGGCATCACTACCAAGGCCATAAAACAAGGCTTCCAGTTTCGCCGAGtt  >  1:2701122/1‑69 (MQ=255)
                                                         gTTTTTGATCGCGGAAACGCTGCCATCACTACCAAGGCCATAAAACAAGGCTTCCAGTTTCGCCGAGtt  >  1:2841777/1‑69 (MQ=255)
                                                                |                                                             
GTAAACAAAATAGCCCTGTGCGTACCACGGCGTGGAATTACCGATAATCTTGATATTGTTTTTGGTCGCGGAAACGCTGCCATCACTACCAAGGCCATAAAACAAGGCTTCCAGTTTCGCCGAGTT  >  W3110S.gb/1441134‑1441259

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: