Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,675,212 G→C 52.9% S87S (TCC→TCG gntY ← predicted gluconate transport associated protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,675,2120GC52.9% 6.9 / 78.4 68S87S (TCC→TCGgntYpredicted gluconate transport associated protein
Reads supporting (aligned to +/- strand):  ref base G (7/25);  new base C (19/17);  total (26/42)
Fisher's exact test for biased strand distribution p-value = 1.24e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCAGCGGTGCATCGTCTGCCACTTTACGCATTTTGGCGTTCGGGGCTTTCAGCGTTAACTGGGAACCCAACTGGTCGGTAACAAAATCGATCTCTGCATCTTCCAGGTATGGTGCGCTTAACTC  >  minE/2675151‑2675274
                                                             |                                                              
tCAGCGGTGCATCGTCTGCCACTTTACGCATTTTGGCGTTCGGGGCTTTCAGCGTTAACTGGGAAccc                                                          <  1:538016/68‑1 (MQ=255)
tCAGCGGTGCATCGTCTGCCACTTTACGCATTTTGGCGTTCGGGGCTTTCAGCGTTAACTGGGAAccc                                                          <  1:509645/68‑1 (MQ=255)
tCAGCGGTGCATCGTCTGCCACTTTACGCATTTTGGCGTTCGGGGCTTTCAGCGTTAACTGGGAAccc                                                          <  1:1930203/68‑1 (MQ=255)
tCAGCGGTGCATCGTCTGCCACTTTACGCATTTTGGCGTTCGGGGCTTTCAGCGTTAACTGGGAAccc                                                          <  1:1319068/68‑1 (MQ=255)
tCAGCGGTGCATCGTCTGCCACTTTACGCATTTTGGCGTTCGGGGCTTTCAGCGTTAACTGGGAAccc                                                          <  1:2642276/68‑1 (MQ=255)
tCAGCGGTGCATCGTCTGCCACTTTACGCATTTTGGCGTTCGGGGCTTTCAGCGTTAACTGGGAACCCa                                                         <  1:1956173/69‑1 (MQ=255)
tCAGCGGTGCATCGTCTGCCACTTTACGCATTTTGGCGTTCGGGGCTTTCAGCGTTAACTGGGAACCCa                                                         <  1:2296864/69‑1 (MQ=255)
tCAGCGGTGCATCGTCTGCCACTTTACGCATTTTGGCGTTCGGGGCTTTCAGCGTTAACTGGGAACCCa                                                         <  1:2556972/69‑1 (MQ=255)
tCAGCGGTGCATCGTCTGCCACTTTACGCATTTTGGCGTTCGGGGCTTTCAGCGTTAACTGGGAACCCa                                                         <  1:1542603/69‑1 (MQ=255)
tCAGCGGTGCATCGTCTGCCACTTTACGCATTTTGGCGTTCGGGGCTTTCAGCGTTAACTGGGAACCCa                                                         <  1:265218/69‑1 (MQ=255)
tCAGCGGTGCATCGGCTGCCACTTTACGCATTTTGGCGTTCGGGGCTTTCAGCGTTAACTGGGAACCCa                                                         <  1:87586/69‑1 (MQ=255)
               cagCCACTTTACGCATTTTGGCGTTAGGCGCTTTCAGCGTCAGCTGCGAACCCAACTGGTCGGt                                               <  1:2835651/62‑1 (MQ=255)
               cagCCACTTTACGCATTTTGGCGTTAGGCGCTTTCAGCGTCAGCTGCGAACCCAACTGGTCGGt                                               <  1:2243344/62‑1 (MQ=255)
               cagCCACTTTACGCATTTTGGCGTTAGGCGCTTTCAGCGTCAGCTGCGAACCCAACTGGTCGGt                                               <  1:2400887/62‑1 (MQ=255)
               cagCCACTTTACGCATTTTGGCGTTAGGCGCTTTCAGCGTCAGCTGCGAACCCAACTGGTCGGt                                               <  1:2424146/62‑1 (MQ=255)
               cagCCACTTTACGCATTTTGGCGTTAGGCGCTTTCAGCGTCAGCTGCGAACCCAACTGGTCGGt                                               <  1:2730767/62‑1 (MQ=255)
               cagCCACTTTACGCATTTTGGCGTTAGGCGCTTTCAGCGTCAGCTGCGAACCCAACTGGTCGGt                                               <  1:1071363/62‑1 (MQ=255)
               cagCCACTTTACGCATTTTGGCGTTAGGCGCTTTCAGCGTCAGCTGCGAACCCAACTGGTCGGt                                               <  1:2926157/62‑1 (MQ=255)
               cagCCACTTTACGCATTTTGGCGTTAGGCGCTTTCAGCGTCAGCTGCGAACCCAACTGGTCGGt                                               <  1:3029231/62‑1 (MQ=255)
               cagCCACTTTACGCATTTTGGCGTTAGGCGCTTTCAGCGTCAGCTGCGAACCCAACTGGTCGGt                                               <  1:469617/62‑1 (MQ=255)
               cagCCACTTTACGCATTTTGGCGTTAGGCGCTTTCAGCGTCAGCTGCGAACCCAACTGGTCGGt                                               <  1:550283/62‑1 (MQ=255)
               cagCCACTTTACGCATTTTGGCGTTAGGCGCTTTCAGCGTCAGCTGCGAACCCAACTGGTCGGt                                               <  1:617873/62‑1 (MQ=255)
               cagCCACTTTACGCATTTTGGCGTTAGGCGCTTTCAGCGTCAGCTGCGAACCCAACTGGTCGGt                                               <  1:927629/62‑1 (MQ=255)
               cagCCACTTTACGCATTTTGGCGTTAGGCGCTTTCAGCGTCAGCTGCGAACCCAACTGGTCGGt                                               <  1:1024157/62‑1 (MQ=255)
               cagCCACTTTACGCATTTTGGCGTTAGGCGCTTTCAGCGTCAGCTGCGAACCCAACTGGTCGGt                                               <  1:1916666/62‑1 (MQ=255)
               cagCCACTTTACGCATTTTGGCGTTAGGCGCTTTCAGCGTCAGCTGCGAACCCAACTGGTCGGt                                               <  1:1421565/62‑1 (MQ=255)
               cagCCACTTTACGCATTTTGGCGTTAGGCGCTTTCAGCGTCAGCTGCGAACCCAACTGGTCGGt                                               <  1:1676440/62‑1 (MQ=255)
               cagCCACTTTACGCATTTTGGCGTTAGGCGCTTTCAGCGTCAGCTGCGAACCCAACTGGTCGGt                                               <  1:1547720/62‑1 (MQ=255)
                        tACGCATTTTGGCGTTCGGGGCTTTCAGCGTTAACTGGGAACCCAACTGGTCGGTAACAAAATCGAtct                                 >  1:2828389/1‑69 (MQ=255)
                            cATTTTGGCGTTCGGGGCTTTCAGCGTTAACTGGGAACCCAACTGGTCGGTAACAAAATCGATCTCTg                              <  1:2133060/68‑1 (MQ=255)
                            cATTTTGGCGTTCGGGGCTTTCAGCGTTAACTGGGAACCCAACTGGTCGGTAACAAAATCGATCTCTg                              <  1:734159/68‑1 (MQ=255)
                            cATTTTGGCGTTCGGGGCTTTCAGCGTTAACTGGGAACCCAACTGGTCGGTAACAAAATCGATCTCTg                              <  1:1116449/68‑1 (MQ=255)
                            cATTTTGGCGTTCGGGGCTTTCAGCGTTAACTGGGAACCCAACTGGTCGGTAACAAAATCGATCTCTg                              <  1:3407570/68‑1 (MQ=255)
                            cATTTTGGCGTTCGGGGCTTTCAGCGTTAACTGGGAACCCAACTGGTCGGTAACAAAATCGATCTCTg                              <  1:2564706/68‑1 (MQ=255)
                            cATTTTGGCGTTCGGGGCTTTCAGCGTTAACTGGGAACCCAACTGGTCGGTAACAAAATCGATCTCTg                              <  1:1458141/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:1511469/1‑68 (MQ=255)
                               tttGGCGTTAGGCGCTTTCAGCGTCAGCTGCGAACCCAACTGGTCGGTAACAAAATCGATCTCTGCAt                           >  1:1318861/1‑68 (MQ=255)
                               tttGGCGTTAGGCGCTTTCAGCGTCAGCTGCGAACCCAACTGGTCGGTAACAAAATCGATCTCTGCAt                           >  1:1909103/1‑68 (MQ=255)
                               tttGGCGTTAGGCGCTTTCAGCGTCAGCTGCGAACCCAACTGGTCGGTAACAAAATCGATCTCTGCATc                          >  1:2889616/1‑69 (MQ=255)
                               tttGGCGTTAGGCGCTTTCAGCGTCAGCTGCGAACCCAACTGGTCGGTAACAAAATCGATCTCTGCATc                          >  1:966685/1‑69 (MQ=255)
                               tttGGCGTTAGGCGCTTTCAGCGTCAGCTGCGAACCCAACTGGTCGGTAACAAAATCGATCTCTGCATc                          >  1:955663/1‑69 (MQ=255)
                               tttGGCGTTAGGCGCTTTCAGCGTCAGCTGCGAACCCAACTGGTCGGTAACAAAATCGATCTCTGCATc                          >  1:2036987/1‑69 (MQ=255)
                               tttGGCGTTAGGCGCTTTCAGCGTCAGCTGCGAACCCAACTGGTCGGTAACAAAATCGATCTCTGCATc                          >  1:2248276/1‑69 (MQ=255)
                               tttGGCGTTAGGCGCTTTCAGCGTCAGCTGCGAACCCAACTGGTCGGTAACAAAATCGATCTCTGCATc                          >  1:1073158/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:1194917/1‑69 (MQ=255)
                               tttGGCGTTAGGCGCTTTCAGCGTCAGCTGCGAACCCAACTGGTCGGTAACAAAATCGATCTCTGCATc                          >  1:1914720/1‑69 (MQ=255)
                               tttGGCGTTAGGCGCTTTCAGCGTCAGCTGCGAACCCAACTGGTCGGTAACAAAATCGATCTCTGCATc                          >  1:3162576/1‑69 (MQ=255)
                               tttGGCGTTAGGCGCTTTCAGCGTCAGCTGCGAACCCAACTGGTCGGTAACAAAATCGATCTCTGCATc                          >  1:2061941/1‑69 (MQ=255)
                               tttGGCGTTAGGCGCTTTCAGCGTCAGCTGCGAACCCAACTGGTCGGTAACAAAATCGATCTCTGCATc                          >  1:2655202/1‑69 (MQ=255)
                               tttGGCGTTAGGCGCTTTCAGCGTCAGCTGCGAACCCAACTGGTCGGTAACAAAATCGATCTCTGCATc                          >  1:2612127/1‑69 (MQ=255)
                                      ttCGGGGCTTTCAGCGTTAACTGGGAACCCAACTGGTCGGTAACAAAATCGATCTCTGCATCTTCCAgg                   <  1:2786161/69‑1 (MQ=255)
                                      ttCGGGGCTTTCAGCGTTAACTGGGAACCCAACTGGTCGGTAACAAAATCGATCTCTGCATCTTCCAgg                   <  1:1185204/69‑1 (MQ=255)
                                      ttCGGGGCTTTCAGCGTTAACTGGGAACCCAACTGGTCGGTAACAAAATCGATCTCTGCATCTTCCAgg                   <  1:692650/69‑1 (MQ=255)
                                      ttCGGGGCTTTCAGCGTTAACTGGGAACCCAACTGGTCGGTAACAAAATCGATCTCTGCATCTTCCAgg                   <  1:1955076/69‑1 (MQ=255)
                                      ttCGGGGCTTTCAGCGTTAACTGGGAACCCAACTGGTCGGTAACAAAATCGATCTCTGCATCTTCCAgg                   <  1:83444/69‑1 (MQ=255)
                                        cGGGGCTTTCAGCGTTAACTGGGAACCCAACTGGTCGGTAACAAAATCGATCTCTGCATCTTCCAGGTa                 >  1:3318156/1‑69 (MQ=255)
                                        cGGGGCTTTCAGCGTTAACTGGGAACCCAACTGGTCGGTAACAAAATCGATCTCTGCATCTTCCAGGTa                 >  1:3116670/1‑69 (MQ=255)
                                        cGGGGCTTTCAGCGTTAACTGGGAACCCAACTGGTCGGTAACAAAATCGATCTCTGCATCTTCCAGGTa                 >  1:1360174/1‑69 (MQ=255)
                                                        aaCTGGGAACCCAACTGGTCGGTAACAAAATCGATCTCTGCATCTTCCAGGTATGGTGCGCTTAACtc  >  1:3240870/1‑68 (MQ=255)
                                                        aaCTGGGAACCCAACTGGTCGGTAACAAAATCGATCTCTGCATCTTCCAGGTATGGTGCGCTTAACtc  >  1:3133174/1‑68 (MQ=255)
                                                        aaCTGGGAACCAAACTGGTCGGTAACAAAATCGAtctc                                >  1:1105884/1‑38 (MQ=255)
                                                             |                                                              
TCAGCGGTGCATCGTCTGCCACTTTACGCATTTTGGCGTTCGGGGCTTTCAGCGTTAACTGGGAACCCAACTGGTCGGTAACAAAATCGATCTCTGCATCTTCCAGGTATGGTGCGCTTAACTC  >  minE/2675151‑2675274

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: