Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 108,284 T→C 10.7% S176S (TCT→TCC aceE → pyruvate dehydrogenase, decarboxylase component E1, thiamin‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE108,2840TC10.7% 151.6 / 8.8 56S176S (TCT→TCCaceEpyruvate dehydrogenase, decarboxylase component E1, thiamin‑binding
Reads supporting (aligned to +/- strand):  ref base T (26/24);  new base C (3/3);  total (29/27)
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.83e-01

CTGGAAGGTCGTCTGACTCAGGAGCAGCTGGATAACTTCCGTCAGGAAGTTCACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGTTCCCGACCGTATCTATGGGTCTGGGTCCGA  >  minE/108216‑108354
                                                                    |                                                                      
cTGGAAGGTCGTCTGACTCAGGAGCAGCTGGATAACTTCCGTCAGGAAGTTCACGGCAATGGCCTCTCTTc                                                                      >  1:3975400/1‑71 (MQ=255)
             tGACTCAGGAGCAGCTGGATAACTTCCGTCAGGAAGTTCACGGCAATGGCCTCTCTTCCTATCCGCAccc                                                          >  1:2777565/1‑70 (MQ=255)
                   aagAGCAGCTGCTGAACTTCCGCCGCGAAGTCGACGGCAAGGGCCTGTCCTCGTACCCGCACCCGAAACt                                                    >  1:1829558/3‑70 (MQ=255)
                   aGGAGCAGCTGGATAACTTCCGTCAGGAAGTTCACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACt                                                    <  1:4074345/70‑1 (MQ=255)
                    ggAGCAGCTGGATAACTTCCGTCAGGAAGTTCACGGCAATGGCCTCTCTTCCTa                                                                   <  1:2315943/54‑1 (MQ=255)
                       gcagcTGGATAACTTCCGTCAGGAAGTTCACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGc                                               >  1:912929/1‑71 (MQ=255)
                             ggATAACTTCCGTCAGGAAGTTCACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAAt                                         >  1:4479318/1‑71 (MQ=255)
                              gATAACTTCCGTCAGGAAGTTCACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAAtt                                        <  1:5081467/71‑1 (MQ=255)
                              gATAACTTCCGTCAGGAAGTTCACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAAtt                                        <  1:2364702/71‑1 (MQ=255)
                              gATAACTTCCGTCAGGAAGTTCACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAAt                                         <  1:1249502/70‑1 (MQ=255)
                              gATAACTTCCGTCAGGAAGTTCACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAAt                                         <  1:1688861/70‑1 (MQ=255)
                              gATAACTTCCGTCAGGAAGTTCACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAAt                                         <  1:2058886/70‑1 (MQ=255)
                              gATAACTTCCGTCAGGAAGTTCACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAAt                                         <  1:3102709/70‑1 (MQ=255)
                               aTAACTTCCGTCAGGAAGTTCACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTc                                       >  1:4948167/1‑71 (MQ=255)
                                tAACTTCCGTCAGGAAGTTCACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCt                                      >  1:1417434/1‑71 (MQ=255)
                                tAACTTCCGTCAGGAAGTTCACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCt                                      >  1:1650229/1‑71 (MQ=255)
                                tAACTTCCGTCAGGAAGTTCACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCt                                      >  1:3270176/1‑71 (MQ=255)
                                gaaCTTCCGCCGCGAAGTCGACGGCAAGGGCCTGTCCTCGTACCCGCACCCGAAACTGATGCCGGACTTCt                                      >  1:1308823/2‑71 (MQ=255)
                                gaaCTTCCGCCGCGAAGTCGACGGCAAGGGCCTGTCCTCGTACCCGCACCCGAAACTGATGCCGGACTTCt                                      >  1:1142796/2‑71 (MQ=255)
                                         tCAGGAAGTTCACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGt                                <  1:2377000/68‑1 (MQ=255)
                                         tCAGGAAGTTCACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGt                                <  1:2446882/68‑1 (MQ=255)
                                           aGGAAGTTCACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGTTc                              >  1:3796696/1‑68 (MQ=255)
                                             gAAGTTCACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAAt                                         <  1:5028764/55‑1 (MQ=255)
                                             gAAGTTCACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAAt                                         <  1:96439/55‑1 (MQ=255)
                                                gTTCTCGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGtt                               <  1:803461/62‑1 (MQ=255)
                                                gTTCACGGCAATGGCCTCTCTTCCTATCCGCCCCCGAAACTGATGCCGGAATTCTGGCAGTTCCCGAcc                        >  1:3341770/1‑69 (MQ=255)
                                                gTTCACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGc                                               >  1:226187/1‑46 (MQ=255)
                                                gTTCACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGtt                               <  1:4356506/62‑1 (MQ=255)
                                                gTTCACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGtt                               <  1:1681075/62‑1 (MQ=255)
                                                gTTCACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGtt                               <  1:4049966/62‑1 (MQ=255)
                                                gTTCACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGtt                               <  1:3827363/62‑1 (MQ=255)
                                                  tCACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGTTccc                            >  1:544506/1‑63 (MQ=255)
                                                  tCACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGTTccc                            >  1:2118353/1‑63 (MQ=255)
                                                   cACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGTTCCCGACCg                       >  1:4633758/1‑67 (MQ=255)
                                                   cACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGTTCCCGACCGt                      >  1:4345289/1‑68 (MQ=255)
                                                   cACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGTTCCCGACCGt                      >  1:944058/1‑68 (MQ=255)
                                                   cACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGTTCCCGACCGTATc                   >  1:3405465/1‑71 (MQ=255)
                                                   cACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGTTCCCGACCGTATc                   >  1:1188865/1‑71 (MQ=255)
                                                   cACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGTTCCCGACCGTATc                   >  1:1731088/1‑71 (MQ=255)
                                                   cACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGTTCCCGACCGTATc                   >  1:1973453/1‑71 (MQ=255)
                                                   cACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGTTCCCGACCGTATc                   >  1:2543948/1‑71 (MQ=255)
                                                    aCGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGTTCCCGACCGTATc                   <  1:4580337/70‑1 (MQ=255)
                                                    aCGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGTTCCCGACCGTATc                   <  1:4015569/70‑1 (MQ=255)
                                                      ggCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGTTCCCGACCGTATCTAt                >  1:901876/1‑71 (MQ=255)
                                                      ggCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGTTCCCGACCGTATCTAt                >  1:2389253/1‑71 (MQ=255)
                                                         aaTGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGTTCCCGACCGTATCTAtggg             <  1:4713268/71‑1 (MQ=255)
                                                           tGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGt                                <  1:2142976/50‑1 (MQ=255)
                                                             gCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGTTcc                             >  1:4759201/1‑51 (MQ=255)
                                                             gCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGTTCCCGACCGTATCTAtgggtctg         >  1:3543206/1‑71 (MQ=255)
                                                             gCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGTTCCCGACCGTATCTAtgggtctg         >  1:2655520/1‑71 (MQ=255)
                                                             gCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGTTCCCGACCGTATCTAtgggtct          <  1:2425248/70‑1 (MQ=255)
                                                             gCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGTTCCCGACCGTATCTAtgggtct          <  1:4802160/70‑1 (MQ=255)
                                                              ccTGTCCTCGTACCCGCACCCGAAACTGATGCCGGACTTCTGGCAGTTCCCGACCGTGTCGAtgggcctgg        <  1:5077737/71‑1 (MQ=255)
                                                              ccTGTCCTCGTACCCGCACCCGAAACTGATGCCGGACTTCTGGCAGTTCCCGACCGTGTCGAtgggcctgg        <  1:5996/71‑1 (MQ=255)
                                                              ccTGTCCTCGTACCCGCACCCGAAACTGATGCCGGACTTCTGGCAGTTCCCGACCGTGTCGAtgggcctgg        <  1:928690/71‑1 (MQ=255)
                                                              ccTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGTTCCCGAcc                        <  1:1216028/55‑1 (MQ=255)
                                                                tctctTACTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGTTCCCGACCGTATCTAtgggtctgggt      >  1:3214/1‑71 (MQ=255)
                                                                 ctctTCCTATCCGCACCCGAAACTGATGCCGGAAt                                         <  1:3545755/35‑1 (MQ=255)
                                                                    ttCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGTTCCCGACCGTATCTAtgggtctgggt      >  1:274237/1‑67 (MQ=255)
                                                                    ttCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGTTCCCGACCGTATCTATGGGTCTGGGTcc    <  1:684846/69‑1 (MQ=255)
                                                                    ttCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGTTCCCGACCGTATCTATGGGTCTGGGTcc    <  1:2824243/69‑1 (MQ=255)
                                                                    ttCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGTTCCCGACCGTATCTATGGGTCTGGGTCCGa  >  1:542396/1‑71 (MQ=255)
                                                                    |                                                                      
CTGGAAGGTCGTCTGACTCAGGAGCAGCTGGATAACTTCCGTCAGGAAGTTCACGGCAATGGCCTCTCTTCCTATCCGCACCCGAAACTGATGCCGGAATTCTGGCAGTTCCCGACCGTATCTATGGGTCTGGGTCCGA  >  minE/108216‑108354

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: