Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,675,194 G→C 45.6% A93A (GCC→GCG gntY ← predicted gluconate transport associated protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,675,1940GC45.6% 17.2 / 89.5 79A93A (GCC→GCGgntYpredicted gluconate transport associated protein
Reads supporting (aligned to +/- strand):  ref base G (8/35);  new base C (19/17);  total (27/52)
Fisher's exact test for biased strand distribution p-value = 1.94e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

TACTCCACGCGCTCCATCAGCGGTGCATCGTCTGCCACTTTACGCATTTTGGCGTTCGGGGCTTTCAGCGTTAACTGGGAACCCAACTGGTCGGTAACAAAATCGATCTCTGCATCTTCCAGGTA  >  minE/2675135‑2675259
                                                           |                                                                 
taCTCCACGCGCTCCATCAGCGGTGCATCGTCTGCCACTTTACGCATTTTGGCGTTCggg                                                                   <  1:1950608/60‑1 (MQ=255)
taCTCCACGCGCTCCATCAGCGGTGCATCGTCTGCCACTTTACGCATTTTGGCGTTCggg                                                                   <  1:1115675/60‑1 (MQ=255)
taCTCCACGCGCTCCATCAGCGGTGCATCGTCTGCCACTTTACGCATTTTGGCGTTCggg                                                                   <  1:1185640/60‑1 (MQ=255)
taCTCCACGCGCTCCATCAGCGGTGCATCGTCTGCCACTTTACGCATTTTGGCGTTCGGGGCTTTCAGc                                                          >  1:68112/1‑69 (MQ=255)
taCTCCACGCGCTCCATCAGCGGTGCATCGTCTGCCACTTTACGCATTTTGGCGTTCGGGGCTTTCAGc                                                          >  1:1394442/1‑69 (MQ=255)
  cTCCACGCGCTCCATCAGCGGTGCATCGTCTGCCACTTTACGCATTTTGGCGTTCGGGGCTTTCAGCg                                                         >  1:1569446/1‑68 (MQ=255)
  cTCCACGCGCTCCATCAGCGGTGCATCGTCTGCCACTTTACGCATTTTGGCGTTCGGGGCTTTCAGCg                                                         >  1:2225107/1‑68 (MQ=255)
    ccACGCGCTCCATCAGCGGTGCATCGTCTGCCACTTTACGCATTTTGGCGTTCGGGGCTTTCAGCGtt                                                       <  1:66910/68‑1 (MQ=255)
    ccACGCGCTCCATCAGCGGTGCATCGTCTGCCACTTTACGCATTTTGGCGTTCGGGGCTTTCAGCGtt                                                       <  1:1200639/68‑1 (MQ=255)
    ccACGCGCTCCATCAGCGGTGCATCGTCTGCCACTTTACGCATTTTGGCGTTCGGGGCTTTCAGCGtt                                                       <  1:480733/68‑1 (MQ=255)
    ccACGCGCTCCATCAGCGGTGCATCGTCTGCCACTTTACGCATTTTGGCGTTCGGGGCTTTCAGCGtt                                                       <  1:697695/68‑1 (MQ=255)
    ccACGCGCTCCATCAGCGGTGCATCGTCTGCCACTTTACGCATTTTGGCGTTCGGGGCTTTCAGCGtt                                                       <  1:1414476/68‑1 (MQ=255)
    ccACGCGCTCCATCAGCGGTGCATCGTCTGCCACTTTACGCATTTTGGCGTTCGGGGCTTTCAGCGTTa                                                      <  1:227000/69‑1 (MQ=255)
    ccACGCGCTCCATCAGCGGTGCATCGTCTGCCACTTTACGCATTTTGGCGTTCGGGGCTTTCAGCGTTa                                                      <  1:652204/69‑1 (MQ=255)
    ccACGCGCTCCATCAGCGGTGCATCGTCTGCCACTTTACGCATTTTGGCGTTCGGGGCTTTCAGCGTTa                                                      <  1:125482/69‑1 (MQ=255)
    ccACGCGCTCCATCAGCGGTGCATCGTCTGCCACTTTACGCATTTTGGCGTTCGGGGCTTTCAGCGTTa                                                      <  1:3174121/69‑1 (MQ=255)
          gcTCCATCAGCGGTGCATCGTCTGCCACTTTACGCATTTTGGCGTTCGGGGCTTTCAGCGTTa                                                      <  1:794560/63‑1 (MQ=255)
                tCAGCGGTGCATCGTCTGCCACTTTACGCATTTTGGCGTTCGGGGCTTTCAGCGTTAACTGGGAAccc                                           <  1:1319068/68‑1 (MQ=255)
                tCAGCGGTGCATCGTCTGCCACTTTACGCATTTTGGCGTTCGGGGCTTTCAGCGTTAACTGGGAAccc                                           <  1:509645/68‑1 (MQ=255)
                tCAGCGGTGCATCGTCTGCCACTTTACGCATTTTGGCGTTCGGGGCTTTCAGCGTTAACTGGGAAccc                                           <  1:2642276/68‑1 (MQ=255)
                tCAGCGGTGCATCGTCTGCCACTTTACGCATTTTGGCGTTCGGGGCTTTCAGCGTTAACTGGGAAccc                                           <  1:538016/68‑1 (MQ=255)
                tCAGCGGTGCATCGTCTGCCACTTTACGCATTTTGGCGTTCGGGGCTTTCAGCGTTAACTGGGAAccc                                           <  1:1930203/68‑1 (MQ=255)
                tCAGCGGTGCATCGTCTGCCACTTTACGCATTTTGGCGTTCGGGGCTTTCAGCGTTAACTGGGAACCCa                                          <  1:2296864/69‑1 (MQ=255)
                tCAGCGGTGCATCGTCTGCCACTTTACGCATTTTGGCGTTCGGGGCTTTCAGCGTTAACTGGGAACCCa                                          <  1:1956173/69‑1 (MQ=255)
                tCAGCGGTGCATCGTCTGCCACTTTACGCATTTTGGCGTTCGGGGCTTTCAGCGTTAACTGGGAACCCa                                          <  1:1542603/69‑1 (MQ=255)
                tCAGCGGTGCATCGTCTGCCACTTTACGCATTTTGGCGTTCGGGGCTTTCAGCGTTAACTGGGAACCCa                                          <  1:2556972/69‑1 (MQ=255)
                tCAGCGGTGCATCGTCTGCCACTTTACGCATTTTGGCGTTCGGGGCTTTCAGCGTTAACTGGGAACCCa                                          <  1:265218/69‑1 (MQ=255)
                tCAGCGGTGCATCGGCTGCCACTTTACGCATTTTGGCGTTCGGGGCTTTCAGCGTTAACTGGGAACCCa                                          <  1:87586/69‑1 (MQ=255)
                               cagCCACTTTACGCATTTTGGCGTTAGGCGCTTTCAGCGTCAGCTGCGAACCCAACTGGTCGGt                                <  1:2424146/62‑1 (MQ=255)
                               cagCCACTTTACGCATTTTGGCGTTAGGCGCTTTCAGCGTCAGCTGCGAACCCAACTGGTCGGt                                <  1:1024157/62‑1 (MQ=255)
                               cagCCACTTTACGCATTTTGGCGTTAGGCGCTTTCAGCGTCAGCTGCGAACCCAACTGGTCGGt                                <  1:3029231/62‑1 (MQ=255)
                               cagCCACTTTACGCATTTTGGCGTTAGGCGCTTTCAGCGTCAGCTGCGAACCCAACTGGTCGGt                                <  1:1071363/62‑1 (MQ=255)
                               cagCCACTTTACGCATTTTGGCGTTAGGCGCTTTCAGCGTCAGCTGCGAACCCAACTGGTCGGt                                <  1:2400887/62‑1 (MQ=255)
                               cagCCACTTTACGCATTTTGGCGTTAGGCGCTTTCAGCGTCAGCTGCGAACCCAACTGGTCGGt                                <  1:2926157/62‑1 (MQ=255)
                               cagCCACTTTACGCATTTTGGCGTTAGGCGCTTTCAGCGTCAGCTGCGAACCCAACTGGTCGGt                                <  1:2835651/62‑1 (MQ=255)
                               cagCCACTTTACGCATTTTGGCGTTAGGCGCTTTCAGCGTCAGCTGCGAACCCAACTGGTCGGt                                <  1:2730767/62‑1 (MQ=255)
                               cagCCACTTTACGCATTTTGGCGTTAGGCGCTTTCAGCGTCAGCTGCGAACCCAACTGGTCGGt                                <  1:2243344/62‑1 (MQ=255)
                               cagCCACTTTACGCATTTTGGCGTTAGGCGCTTTCAGCGTCAGCTGCGAACCCAACTGGTCGGt                                <  1:469617/62‑1 (MQ=255)
                               cagCCACTTTACGCATTTTGGCGTTAGGCGCTTTCAGCGTCAGCTGCGAACCCAACTGGTCGGt                                <  1:1916666/62‑1 (MQ=255)
                               cagCCACTTTACGCATTTTGGCGTTAGGCGCTTTCAGCGTCAGCTGCGAACCCAACTGGTCGGt                                <  1:550283/62‑1 (MQ=255)
                               cagCCACTTTACGCATTTTGGCGTTAGGCGCTTTCAGCGTCAGCTGCGAACCCAACTGGTCGGt                                <  1:617873/62‑1 (MQ=255)
                               cagCCACTTTACGCATTTTGGCGTTAGGCGCTTTCAGCGTCAGCTGCGAACCCAACTGGTCGGt                                <  1:1676440/62‑1 (MQ=255)
                               cagCCACTTTACGCATTTTGGCGTTAGGCGCTTTCAGCGTCAGCTGCGAACCCAACTGGTCGGt                                <  1:1547720/62‑1 (MQ=255)
                               cagCCACTTTACGCATTTTGGCGTTAGGCGCTTTCAGCGTCAGCTGCGAACCCAACTGGTCGGt                                <  1:1421565/62‑1 (MQ=255)
                               cagCCACTTTACGCATTTTGGCGTTAGGCGCTTTCAGCGTCAGCTGCGAACCCAACTGGTCGGt                                <  1:927629/62‑1 (MQ=255)
                                        tACGCATTTTGGCGTTCGGGGCTTTCAGCGTTAACTGGGAACCCAACTGGTCGGTAACAAAATCGAtct                  >  1:2828389/1‑69 (MQ=255)
                                            cATTTTGGCGTTCGGGGCTTTCAGCGTTAACTGGGAACCCAACTGGTCGGTAACAAAATCGATCTCTg               <  1:3407570/68‑1 (MQ=255)
                                            cATTTTGGCGTTCGGGGCTTTCAGCGTTAACTGGGAACCCAACTGGTCGGTAACAAAATCGATCTCTg               <  1:734159/68‑1 (MQ=255)
                                            cATTTTGGCGTTCGGGGCTTTCAGCGTTAACTGGGAACCCAACTGGTCGGTAACAAAATCGATCTCTg               <  1:2564706/68‑1 (MQ=255)
                                            cATTTTGGCGTTCGGGGCTTTCAGCGTTAACTGGGAACCCAACTGGTCGGTAACAAAATCGATCTCTg               <  1:1116449/68‑1 (MQ=255)
                                            cATTTTGGCGTTCGGGGCTTTCAGCGTTAACTGGGAACCCAACTGGTCGGTAACAAAATCGATCTCTg               <  1:1458141/68‑1 (MQ=255)
                                            cATTTTGGCGTTCGGGGCTTTCAGCGTTAACTGGGAACCCAACTGGTCGGTAACAAAATCGATCTCTg               <  1:2133060/68‑1 (MQ=255)
                                             aTTTTGGCGTTCGGGGCTTTCAGCGTTAACTGGGAACCCAACTGGTCGGTATCAAAATCGATCTCTGCa             <  1:888360/69‑1 (MQ=255)
                                             aTTTTGGCGTTCGGGGCTTTCAGCGTTAACTGGGAACCCAACTGGTCGGTAACAAAATCGATCTCTg               <  1:309279/67‑1 (MQ=255)
                                             aTTTTGGCGTTCGGGGCTTTCAGCGTTAACTGGGAACCCAACTGGTCGGTAACAAAATCGATCTCTGCa             <  1:539326/69‑1 (MQ=255)
                                               tttGGCGTTAGGCGCTTTCAGCGTCAGCTGCGAACCCAACTGGTCGGTAACAAAATCGATCTCTGCAt            >  1:1318861/1‑68 (MQ=255)
                                               tttGGCGTTAGGCGCTTTCAGCGTCAGCTGCGAACCCAACTGGTCGGTAACAAAATCGATCTCTGCAt            >  1:1511469/1‑68 (MQ=255)
                                               tttGGCGTTAGGCGCTTTCAGCGTCAGCTGCGAACCCAACTGGTCGGTAACAAAATCGATCTCTGCAt            >  1:1909103/1‑68 (MQ=255)
                                               tttGGCGTTAGGCGCTTTCAGCGTCAGCTGCGAACCCAACTGGTCGGTAACAAAATCGATCTCTGCATc           >  1:1914720/1‑69 (MQ=255)
                                               tttGGCGTTAGGCGCTTTCAGCGTCAGCTGCGAACCCAACTGGTCGGTAACAAAATCGATCTCTGCATc           >  1:966685/1‑69 (MQ=255)
                                               tttGGCGTTAGGCGCTTTCAGCGTCAGCTGCGAACCCAACTGGTCGGTAACAAAATCGATCTCTGCATc           >  1:955663/1‑69 (MQ=255)
                                               tttGGCGTTAGGCGCTTTCAGCGTCAGCTGCGAACCCAACTGGTCGGTAACAAAATCGATCTCTGCATc           >  1:1073158/1‑69 (MQ=255)
                                               tttGGCGTTAGGCGCTTTCAGCGTCAGCTGCGAACCCAACTGGTCGGTAACAAAATCGATCTCTGCATc           >  1:1194917/1‑69 (MQ=255)
                                               tttGGCGTTAGGCGCTTTCAGCGTCAGCTGCGAACCCAACTGGTCGGTAACAAAATCGATCTCTGCATc           >  1:558374/1‑69 (MQ=255)
                                               tttGGCGTTAGGCGCTTTCAGCGTCAGCTGCGAACCCAACTGGTCGGTAACAAAATCGATCTCTGCATc           >  1:552767/1‑69 (MQ=255)
                                               tttGGCGTTAGGCGCTTTCAGCGTCAGCTGCGAACCCAACTGGTCGGTAACAAAATCGATCTCTGCATc           >  1:2325995/1‑69 (MQ=255)
                                               tttGGCGTTAGGCGCTTTCAGCGTCAGCTGCGAACCCAACTGGTCGGTAACAAAATCGATCTCTGCATc           >  1:2725638/1‑69 (MQ=255)
                                               tttGGCGTTAGGCGCTTTCAGCGTCAGCTGCGAACCCAACTGGTCGGTAACAAAATCGATCTCTGCATc           >  1:2889616/1‑69 (MQ=255)
                                               tttGGCGTTAGGCGCTTTCAGCGTCAGCTGCGAACCCAACTGGTCGGTAACAAAATCGATCTCTGCATc           >  1:2655202/1‑69 (MQ=255)
                                               tttGGCGTTAGGCGCTTTCAGCGTCAGCTGCGAACCCAACTGGTCGGTAACAAAATCGATCTCTGCATc           >  1:2612127/1‑69 (MQ=255)
                                               tttGGCGTTAGGCGCTTTCAGCGTCAGCTGCGAACCCAACTGGTCGGTAACAAAATCGATCTCTGCATc           >  1:3162576/1‑69 (MQ=255)
                                               tttGGCGTTAGGCGCTTTCAGCGTCAGCTGCGAACCCAACTGGTCGGTAACAAAATCGATCTCTGCATc           >  1:2248276/1‑69 (MQ=255)
                                               tttGGCGTTAGGCGCTTTCAGCGTCAGCTGCGAACCCAACTGGTCGGTAACAAAATCGATCTCTGCATc           >  1:2061941/1‑69 (MQ=255)
                                               tttGGCGTTAGGCGCTTTCAGCGTCAGCTGCGAACCCAACTGGTCGGTAACAAAATCGATCTCTGCATc           >  1:2036987/1‑69 (MQ=255)
                                                      ttCGGGGCTTTCAGCGTTAACTGGGAACCCAACTGGTCGGTAACAAAATCGATCTCTGCATCTTCCAgg    <  1:692650/69‑1 (MQ=255)
                                                      ttCGGGGCTTTCAGCGTTAACTGGGAACCCAACTGGTCGGTAACAAAATCGATCTCTGCATCTTCCAgg    <  1:2786161/69‑1 (MQ=255)
                                                      ttCGGGGCTTTCAGCGTTAACTGGGAACCCAACTGGTCGGTAACAAAATCGATCTCTGCATCTTCCAgg    <  1:1185204/69‑1 (MQ=255)
                                                      ttCGGGGCTTTCAGCGTTAACTGGGAACCCAACTGGTCGGTAACAAAATCGATCTCTGCATCTTCCAgg    <  1:83444/69‑1 (MQ=255)
                                                      ttCGGGGCTTTCAGCGTTAACTGGGAACCCAACTGGTCGGTAACAAAATCGATCTCTGCATCTTCCAgg    <  1:1955076/69‑1 (MQ=255)
                                                        cGGGGCTTTCAGCGTTAACTGGGAACCCAACTGGTCGGTAACAAAATCGATCTCTGCATCTTCCAGGTa  >  1:1360174/1‑69 (MQ=255)
                                                        cGGGGCTTTCAGCGTTAACTGGGAACCCAACTGGTCGGTAACAAAATCGATCTCTGCATCTTCCAGGTa  >  1:3318156/1‑69 (MQ=255)
                                                        cGGGGCTTTCAGCGTTAACTGGGAACCCAACTGGTCGGTAACAAAATCGATCTCTGCATCTTCCAGGTa  >  1:3116670/1‑69 (MQ=255)
                                                           |                                                                 
TACTCCACGCGCTCCATCAGCGGTGCATCGTCTGCCACTTTACGCATTTTGGCGTTCGGGGCTTTCAGCGTTAACTGGGAACCCAACTGGTCGGTAACAAAATCGATCTCTGCATCTTCCAGGTA  >  minE/2675135‑2675259

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: