Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,733,525 T→C 33.7% G40G (GGT→GGC yheL → predicted intracellular sulfur oxidation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,733,5250TC33.7% 114.5 / 87.0 151G40G (GGT→GGCyheLpredicted intracellular sulfur oxidation protein
Reads supporting (aligned to +/- strand):  ref base T (45/55);  new base C (33/18);  total (78/73)
Fisher's exact test for biased strand distribution p-value = 2.57e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 153 total aligned reads displayed.

GCGTCTGCTCAGTGAAGGAGACGAACTGCTATTATTGCAAGATGGCGTAACTGCCGCAGTTGACGGTAACCGCTACCTTGAAAGTCTGCGTAATGCCCCCATTAAGGTCTATGCCCTGAACGAAGACCTTATTG  >  minE/2733459‑2733592
                                                                  |                                                                   
gcGTCTGCTCAGTGAAGGAGACGAACTGCTATTATTGCAAGATGGCGTAACTGCCGCAGTTGACGGt                                                                     <  1:689667/67‑1 (MQ=255)
    cTGCTCAGTGAAGGAGACGAACTGCTATTATTGCAAGATGGCGTAACTGCCGCAGTTGACGGTAACCGc                                                               <  1:906048/69‑1 (MQ=255)
    cTGCTCAGTGAAGGAGACGAACTGCTATTATTGCAAGATGGCGTAACTGCCGCAGTTGACGGTAACCGc                                                               <  1:3255259/69‑1 (MQ=255)
    cTGCTCAGCGAAGGAGACGAACTGCTATTATTGCAAGATGGCGTAACTGCCGCAGTTGACGGTAACCGc                                                               <  1:850104/69‑1 (MQ=255)
     tGCTCAGTGAAGGAGACGAACTGCTATTATTGCAAGATGGCGTAACTGCCGCAGTTGACGGTTACCGc                                                               <  1:904017/68‑1 (MQ=255)
       cTCAGTGAAGGTGACGAACTGCTATTATTGCAAGATGGCGTAACTGCCGCAGTTGACGGCAACTGCt                                                              >  1:2865823/1‑67 (MQ=255)
       cTCAGTGAAGGTGACGAACTGCTATTATTGCAAGATGGCGTAACTGCCGCAGTTGACGGCAACCGCt                                                              >  1:1791844/1‑67 (MQ=255)
       cTCAGTGAAGGTGACGAACTGCTATTATTGCAAGATGGCGTAACTGCCGCAGTTGACGGCAACCGCt                                                              >  1:2325499/1‑67 (MQ=255)
       cTCAGTGAAGGTGACGAACTGCTATTATTGCAAGATGGCGTAACTGCCGCAGTTGACGGCAACCGCt                                                              >  1:3134981/1‑67 (MQ=255)
       cTCAGTGAAGGTGACGAACTGCTATTATTGCAAGATGGCGTAACTGCCGCAGTTGACGGCAACCGCt                                                              >  1:2271421/1‑67 (MQ=255)
       cTCAGTGAAGGTGACGAACTGCTATTATTGCAAGATGGCGTAACTGCCGCAGTTGACGGCAACCGCt                                                              >  1:2221177/1‑67 (MQ=255)
       cTCAGTGAAGGTGACGAACTGCTATTATTGCAAGATGGCGTAACTGCCGCAGTTGACGGCAACCGCt                                                              >  1:2217210/1‑67 (MQ=255)
       cTCAGTGAAGGTGACGAACTGCTATTATTGCAAGATGGCGTAACTGCCGCAGTTGACGGCAACCGCt                                                              >  1:2339602/1‑67 (MQ=255)
       cTCAGTGAAGGTGACGAACTGCTATTATTGCAAGATGGCGTAACTGCCGCAGTTGACGGCAACCGCt                                                              >  1:3238111/1‑67 (MQ=255)
       cTCAGTGAAGGTGACGAACTGCTATTATTGCAAGATGGCGTAACTGCCGCAGTTGACGGCAACCGCt                                                              >  1:2714241/1‑67 (MQ=255)
       cTCAGTGAAGGTGACGAACTGCTATTATTGCAAGATGGCGTAACTGCCGCAGTTGACGGCAACCGCt                                                              >  1:1693370/1‑67 (MQ=255)
       cTCAGTGAAGGTGACGAACTGCTATTATTGCAAGATGGCGTAACTGCCGCAGTTGACGGCAACCGCt                                                              >  1:1675836/1‑67 (MQ=255)
       cTCAGTGAAGGTGACGAACTGCTATTATTGCAAGATGGCGTAACTGCCGCAGTTGACGGCAACCGCt                                                              >  1:2527555/1‑67 (MQ=255)
       cTCAGTGAAGGTGACGAACTGCTATTATTGCAAGATGGCGTAACTGCCGCAGTTGACGGCAACCGCt                                                              >  1:3064548/1‑67 (MQ=255)
       cTCAGTGAAGGTGACGAACTGCTATTATTGCAAGATGGCGTAACTGCCGCAGTTGACGGCAACCGCt                                                              >  1:656200/1‑67 (MQ=255)
       cTCAGTGAAGGTGACGAACTGCTATTATTGCAAGATGGCGTAACTGCCGCAGTTGACGGCAACCGCt                                                              >  1:566472/1‑67 (MQ=255)
       cTCAGTGAAGGTGACGAACTGCTATTATTGCAAGATGGCGTAACTGCCGCAGTTGACGGCAACCGCt                                                              >  1:510258/1‑67 (MQ=255)
             gAAGGAGACGAACTGCTATTATTGCAAGATGGCGTAACTGCCGCAGTTGACGGTAACCGCt                                                              >  1:1946438/1‑61 (MQ=255)
             gAAGGAGACGAACTGCTATTATTGCAAGATGGCGTAACTGCCGCAGTTGACGGTAACCGCt                                                              >  1:2279110/1‑61 (MQ=255)
             gAAGGAGACGAACTGCTATTATTGCAAGATGGCGTAACTGCCGCAGTTGACGGTAACCGCt                                                              >  1:1292544/1‑61 (MQ=255)
             gAAGGAGACGAACTGCTATTATTGCAAGATGGCGTAACTGCCGCAGTTGACGGTAACCGCt                                                              >  1:60673/1‑61 (MQ=255)
             gAAGGAGACGAACTGCTATTATTGCAAGATGGCGTAACTGCCGCAGTTGACGGTAACCGCt                                                              >  1:3260062/1‑61 (MQ=255)
             gAAGGAGACGAACTGCTATTATTGCAAGATGGCGTAACTGCCGCAGTTGACGGTAACCGCt                                                              >  1:425171/1‑61 (MQ=255)
             gAAGGAGACGAACTGCTATTATTGCAAGATGGCGTAACTGCCGCAGTTGACGGTAACCGCt                                                              >  1:1737285/1‑61 (MQ=255)
                  tgaCGAACTGCTATTATTGCAAGATGGCGTAACTGCCGCAGTTGACGGCAACCGCTACCTTGGAAGTCt                                                 <  1:3020482/68‑1 (MQ=255)
                  tgaCGAACTGCTATTATTGCAAGATGGCGTAACTGCCGCAGTTGACGGCAACCGCTACCTTGAAAGTCt                                                 <  1:289886/68‑1 (MQ=255)
                  tgaCGAACTGCTATTATTGCAAGATGGCGTAACTGCCGCAGTTGACGGCAACCGCTACCTTGAAAGTCt                                                 <  1:328964/68‑1 (MQ=255)
                  tgaCGAACTGCTATTATTGCAAGATGGCGTAACTGCCGCAGTTGACGGCAACCGCTACCTTGAAAGTCt                                                 <  1:1825655/68‑1 (MQ=255)
                  tgaCGAACTGCTATTATTGCAAGATGGCGTAACTGCCGCAGTTGACGGCAACCGCTACCTTGAAAGTCt                                                 <  1:1061656/68‑1 (MQ=255)
                  tgaCGAACTGCTATTATTGCAAGATGGCGTAACTGCCGCAGTTGACGGCAACCGCTACCTTGAAAGTCt                                                 <  1:2604934/68‑1 (MQ=255)
                  tgaCGAACTGCTATTATTGCAAGATGGCGTAACTGCCGCAGTTGACGGCAACCGCTACCTTGAAAGTCt                                                 <  1:1616529/68‑1 (MQ=255)
                  tgaCGAACTGCTATTATTGCAAGATGGCGTAACTGCCGCAGTTGACGGCAACCGCTACCTTGAAAGTCt                                                 <  1:880795/68‑1 (MQ=255)
                  tgaCGAACTGCTATTATTGCAAGATGGCGTAACTGCCGCAGTTGACGGCAACCGCTACCTTGAAAGTCt                                                 <  1:1591731/68‑1 (MQ=255)
                  tgaCGAACTGCTATTATTGCAAGATGGCGTAACTGCCGCAGTTGACGGCAACCGCTACCTTGAAAGTCt                                                 <  1:1260893/68‑1 (MQ=255)
                  tgaCGAACTGCTATTATTGCAAGATGGCGTAACTGCCGCAGTTGACGGCAACCGCTACCTTGAAAGTCt                                                 <  1:1537981/68‑1 (MQ=255)
                   gaCGAACTGCTATTATTGCAAGATGGCGTAACTGCCGCAGTTGACGGCAACCGCTACCTTGAAAGTCt                                                 <  1:244802/68‑1 (MQ=255)
                      gAACTGCTATTATTGCAAGATGGCGTAACTGCCGCAGTTGACGGTAACCGCTACCTTGAAAGTCTGCGt                                             <  1:3002348/69‑1 (MQ=255)
                      gAACTGCTATTATTGCAAGATGGCGTAACTGCCGCAGTTGACGGTAACCGCTACCTTGAAAGTCTGCGt                                             <  1:2172298/69‑1 (MQ=255)
                      gAACTGCTATTATTGCAAGATGGCGTAACTGCCGCAGTTGACGGTAACCGCTACCTTGAAAGTCTGCGt                                             <  1:2966130/69‑1 (MQ=255)
                      gAACTGCTATTATTGCAAGATGGCGTAACTGCCGCAGTTGACGGTAACCGCTACCTTGAAAGTCTGCGt                                             <  1:1504752/69‑1 (MQ=255)
                      gAACTGCTATTATTGCAAGATGGCGTAACTGCCGCAGTTGACGGTAACCGCTACCTTGAAAGTCTGCGt                                             <  1:1174318/69‑1 (MQ=255)
                       aaCTGCTATTATTGCAAGATGGCGTAACTGCCGCAGTTGACGGTAACCGCTACCTTGAAAGTCTGGGt                                             <  1:1465208/68‑1 (MQ=255)
                       aaCTGCTATTATTGCAAGATGGCGTAACTGCCGCAGTTGACGGTAACCGCTACCTTGAAAGTCTGCGt                                             <  1:3193013/68‑1 (MQ=255)
                       aaCTGCTATTATTGCAAGATGGCGTAACTGCCGCAGTTGACGGTAACCGCTACCTTGAAAGTCTGCGt                                             <  1:2095934/68‑1 (MQ=255)
                        aCTGCTATTATTGCAAGATGGCGTAACTGCCGCAGTTGACGGTAACCGCTACCTTGGAAGTCTGCGTaa                                           >  1:3277683/1‑69 (MQ=255)
                        aCTGCTATTATTGCAAGATGGCGTAACTGCCGCAGTTGACGGTAACCGCTACCTTGAAAGTCTGCGTaa                                           >  1:3100524/1‑69 (MQ=255)
                        aCTGCTATTATTGCAAGATGGCGTAACTGCCGCAGTTGACGGTAACCGCTACCTTGAAAGTCTGCGTaa                                           >  1:759392/1‑69 (MQ=255)
                        aCTGCTATTATTGCAAGATGGCGTAACTGCCGCAGTTGACGGTAACCGCTACCTTGAAAGTCTGCGTaa                                           >  1:377161/1‑69 (MQ=255)
                        aCTGCTATTATTGCAAGATGGCGTAACTGCCGCAGTTGACGGTAACCGCTACCTTGAAAGTCTGCGTaa                                           >  1:55102/1‑69 (MQ=255)
                        aCTGCTATTATTGCAAGATGGCGTAACTGCCGCAGTTGACGGTAACCGCTACCTTGAAAGTCTGCGTaa                                           >  1:958631/1‑69 (MQ=255)
                        aCTGCTATTATTGCAAGATGGCGTAACTGCCGCAGTTGACGGTAACCGCTACCTTGAAAGTCTGCGTaa                                           >  1:1112252/1‑69 (MQ=255)
                        aCTGCTATTATTGCAAGATGGCGTAACTGCCGCAGTTGACGGTAACCGCTACCTTGAAAGTCTGCGTaa                                           >  1:2694590/1‑69 (MQ=255)
                        aCTGCTATTATTGCAAGATGGCGTAACTGCCGCAGTTGACGGTAACCGCTACCTTGAAAGTCTGCGTaa                                           >  1:253417/1‑69 (MQ=255)
                        aCTGCTATTATTGCAAGATGGCGTAACTGCCGCAGTTGACGGTAACCGCTACCTTGAAAGTCTGCGTa                                            >  1:3436458/1‑68 (MQ=255)
                                 attGCAAGATGGCGTAACTGCCGCAGTTGACGGTAACCGCTACCTTGAAAGTCTGCGTAATGccccc                                    >  1:3099422/1‑67 (MQ=255)
                                 attGCAAGATGGCGTAACTGCCGCAGTTGACGGTAACCGCTACCTTGAAAGTCTGCGTAATGCCCCCa                                   >  1:1361/1‑68 (MQ=255)
                                 attGCAAGATGGCGTAACTGCCGCAGTTGACGGTAACCGCTACCTTGAAAGTCTGCGTAATGCCCCCa                                   >  1:1524139/1‑68 (MQ=255)
                                 attGCAAGATGGCGTAACTGCCGCAGTTGACGGTAACCGCTACCTTGAAAGTCTGCGTAATGCCCCCa                                   >  1:424087/1‑68 (MQ=255)
                                 attGCAAGATGGCGTAACTGACGCAGTTGACGGTAACCGCTACCTTGAAAGTCTGCGTAATGccacca                                   >  1:2132661/1‑68 (MQ=255)
                                        gATGGCGTAACTGCCGCAGTTGACGGTAACCGCTACCTTGAAAGTCTGCGt                                             <  1:1330602/51‑1 (MQ=255)
                                        gATGGCGTAACTGCCGCAGTTGACGGTAACCGCTACCTTGAAAGTCTGCGTAATGCCCCCATTAAGGTc                           <  1:2493997/69‑1 (MQ=255)
                                              gTAACTGCCGCAGTTGACGGTAACCGCTACCTTGAAAGTCTGCGt                                             <  1:2907570/45‑1 (MQ=255)
                                                 aCTGCCGCAGTTGACGGTAACCGCTACCTTGAAAGTCTGCGTAATGCCCCCATTAAGGTCTATGCCCTg                  >  1:2022479/1‑69 (MQ=255)
                                                       gCAGTTGACGGTAACCGCTACCTTGAAAGTCTGCGTAATGCCCCCATTAAGGTCTATGCCCTGAACg              <  1:694323/67‑1 (MQ=255)
                                                            tGACGGTAACCGCTACCTTGAAAGTCTGCGTAATGCCCCCATTAAGGTCTATGCCCTGAACGAAGAcc        >  1:2504534/1‑68 (MQ=255)
                                                            tGACGGTAACCGCTACCTTGAAAGTCTGCGTAATGCCCCCATTAAGGTCTATGCCCTGAACGAAGACCt       <  1:373127/69‑1 (MQ=255)
                                                            tGACGGTAACCGCTACCTTGAAAGTCTGCGTAATGCCCCCATTAAGGTCTATGCCCTGAACGAAGACCt       <  1:1425360/69‑1 (MQ=255)
                                                            tGACGGTAACCGCTACCTTGAAAGTCTGCGTAATGCCCCCATTAAGGTCTATGCCCTGAACGAAGACCt       <  1:1262843/69‑1 (MQ=255)
                                                            tGACGGTAACCGCTACCTTGAAAGTCTGCGTAATGCCCCCATTAAGGTCTATGCCCTGAACGAAGACCt       <  1:592011/69‑1 (MQ=255)
                                                            tGACGGTAACCGCTACCTTGAAAGTCTGCGTAATGCCCCCATTAAGGTCTATGCCCTGAACGAAGACCt       <  1:1360303/69‑1 (MQ=255)
                                                            tGACGGTAACCGCTACCTTGAAAGTCTGCGTAATGCCCCCATTAAGGTCTATGCCCTGAACGAAGACCt       <  1:2870785/69‑1 (MQ=255)
                                                            tGACGGTAACCGCTACCTTGAAAGTCTGCGTAATGCCCCCATTAAGGTCTATGCCCTGAACGAAGACCt       <  1:765109/69‑1 (MQ=255)
                                                            tGACGGTAACCGCTACCTTGAAAGTCTGCGTAATGCCCCCATTAAGGTCTATGCCCTGAACGAAGACCt       >  1:1009097/1‑69 (MQ=255)
                                                            tGACGGTAACCGCTACCTTGAAAGTCTGCGTAATGCCCCCATTAAGGTCTATGCCCTGAACGAAGACCt       <  1:100747/69‑1 (MQ=255)
                                                            tGACGGTAACCGCTACCTTGAAAGTCTGCGTAATGCCCCCATTAAGGTCTATGCCCTGAACGAAGACCt       <  1:3427655/69‑1 (MQ=255)
                                                            tGACGGTAACCGCTACCTTGAAAGTCTGCGTAATGCCCCCATTAAGGTCTATGCCCTGAACGAAGACCt       <  1:3360103/69‑1 (MQ=255)
                                                            tGACGGTAACCGCTACCTTGAAAGTCTGCGTAATGCCCCCATTAAGGTCTATGCCCTGAACGAAGACCt       <  1:1933954/69‑1 (MQ=255)
                                                            tGACGGTAACCGCTACCTTGAAAGTCTGCGTAATGCCCCCATTAAGGTCTATGCCCTGAACGAAGACCt       <  1:2335806/69‑1 (MQ=255)
                                                            tGACGGTAACCGCTACCTTGAAAGTCTGCGTAATGCCCCCATTAAGGTCTATGCCCTGAACGAAGACCt       <  1:2441472/69‑1 (MQ=255)
                                                             gACGGTAACCGCTACCTTGAAAGTCTGCGTAATGCCCCCATTAAGGTCTATGCCCTGAACGaa            >  1:2774138/1‑63 (MQ=255)
                                                             gACGGTAACCGCTACCTTGAAAGTCTGCGTAATGCCCCCATTAAGGTCTATGCCCTGAACGaa            >  1:2037062/1‑63 (MQ=255)
                                                             gACGGTAACCGCTACCTTGAAAGTCTGCGTAATGCCCCCATTAAGGTCTATGCCCTGAACGaa            >  1:3061845/1‑63 (MQ=255)
                                                             gACGGTAACCGCTACCTTGAAAGTCTGCGTAATGCCCCCATTAAGGTCTATGCCCTGAACGaa            >  1:1172264/1‑63 (MQ=255)
                                                             gACGGTAACCGCTACCTTGAAAGTCTGCGTAATGCCCCCATTAAGGTCTATGCCCTGAACGaa            >  1:618326/1‑63 (MQ=255)
                                                             gACGGTAACCGCTACCTTGAAAGTCTGCGTAATGCCCCCATTAAGGTCTATGCCCTGAACGaa            >  1:3118469/1‑63 (MQ=255)
                                                             gACGGTAACCGCTACCTTGAAAGTCTGCGTAATGCCCCCATTAAGGTCTATGCCCTGAACGAAGACCtt      <  1:3366396/69‑1 (MQ=255)
                                                             gACGGTAACCGCTACCTTGAAAGTCTGCGTAATGCCCCCATTAAGGTCTATGCCCTGAACGAAGACCtt      <  1:2963257/69‑1 (MQ=255)
                                                             gACGGTAACCGCTACCTTGAAAGTCTGCGTAATGCCCCCATTAAGGTCTATGCCCTGAACGAAGACCtt      <  1:1168223/69‑1 (MQ=255)
                                                             gACGGTAACCGCTACCTTGAAAGTCTGCGTAATGCCCCCATTAAGGTCTATGCCCTGAACGAAGACCtt      <  1:1953991/69‑1 (MQ=255)
                                                             gACGGTAACCGCTACCTTGAAAGTCTGCGTAATGCCCCCATTAAGGTCTATGCCCTGAACGAAGACCtt      <  1:188448/69‑1 (MQ=255)
                                                             gACGGTAACCGCTACCTTGAAAGTCTGCGTAATGCCCCCATTAAGGTCTATGCCCTGAACGAAGACCt       <  1:1163913/68‑1 (MQ=255)
                                                                 gTAACCGCTACCTTGAAAGTCTGCGTAATGCCCCCATTAAGGTCTATGCCCTGAACGAAGACCTTATTg  <  1:3161780/69‑1 (MQ=255)
                                                                 gTAACCGCTACCTTGAAAGTCTGCGTAATGCCCCCATTAAGGTCTATGCCCTGAACGAAGACCTTATTg  <  1:2298194/69‑1 (MQ=255)
                                                                 gTAACCGCTACCTTGAAAGTCTGCGTAATGCCCCCATTAAGGTCTATGCCCTGAACGAAGACCTTATTg  <  1:1848117/69‑1 (MQ=255)
                                                                 gTAACCGCTACCTTGAAAGTCTGCGTAATGCCCCCATTAAGGTCTATGCCCTGAACGAAGACCTTATTg  <  1:2804048/69‑1 (MQ=255)
                                                                  |                                                                   
GCGTCTGCTCAGTGAAGGAGACGAACTGCTATTATTGCAAGATGGCGTAACTGCCGCAGTTGACGGTAACCGCTACCTTGAAAGTCTGCGTAATGCCCCCATTAAGGTCTATGCCCTGAACGAAGACCTTATTG  >  minE/2733459‑2733592

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: