Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,613,635 C→G 39.3% L247L (CTG→CTC eutH ← predicted inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,613,6350CG39.3% 56.3 / 83.4 89L247L (CTG→CTCeutHpredicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (34/20);  new base G (20/15);  total (54/35)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCTCACCGGGTTTATCGCCAGGGGCCATAAAGATAGGATCCAGACCGGGGATCAGTTCCCAGCCAAGCAGGAATTTCACTACCGCAGCGGCAAGACCGAGGGTGATCAATGCAACGAGGAATTTGGCGAAGA  >  minE/1613567‑1613699
                                                                    |                                                                
acctcaccGGGTTTATCGCCAGGGGCCATAAAGATAGGATCCAGACCGGGGATCAGTTCCCAGCCAAGc                                                                  >  1:743106/1‑69 (MQ=255)
acctcaccGGGTTTATCGCCAGGGGCCATAAAGATAGGATCCAGACCGGGGATCAGTTCCCAGCCAAGc                                                                  >  1:2669300/1‑69 (MQ=255)
   tcaccGGGTTTATCGCCAGGGGCCATAAAGATAGGATCCAGACCGGGGATCAGTTCCCAGCCAAGCAgg                                                               <  1:2887955/69‑1 (MQ=255)
   tcaccGGGTTTATCGCCAGGGGCCATAAAGATAGGATCCAGACCGGGGATCAGTTCCCAGCCAAGCAgg                                                               <  1:1524875/69‑1 (MQ=255)
   tcaccGGGTTTATCGCCAGGGGCCATAAAGATAGGATCCAGACCGGGGATCAGTTCCCAGCCAAGCAgg                                                               <  1:1521043/69‑1 (MQ=255)
   tcaccGGGTTTATCGCCAGGGGCCATAAAGATAGGATCCAGACCGGGGATCAGTTCCCAGCCAAGCAgg                                                               <  1:2787462/69‑1 (MQ=255)
        gggTTTATCGCCAGGGGCCATAAAGATAGGATCCAGACCGGGGATCAGTTCCCAGCCAAGCAGGAAttt                                                          <  1:2510324/69‑1 (MQ=255)
         ggTTTATCGCCAGGGGCCATAAAGATAGGATCCAGACCGGGGATCAGTTCCCAGCCAAGCAGGAATTTc                                                         <  1:2914808/69‑1 (MQ=255)
         ggTTTATCGCCAGGGGCCATAAAGATAGGATCCAGACCGGGGATCAGTTCCCAGCCAAGCAGGAATTTc                                                         <  1:1833311/69‑1 (MQ=255)
         ggTTTATCGCCAGGGGCCATAAAGATAGGATCCAGACCGGGGATCAGTTCCCAGCCAAGCAGGAATTTc                                                         <  1:497472/69‑1 (MQ=255)
         ggTTTATCGCCAGGGGCCATAAAGATAGGATCCAGACCGGGGATCAGTTCCCAGCCAAGCAGGAATTTc                                                         <  1:2862731/69‑1 (MQ=255)
         ggTTTATCGCCAGGGGCCATAAAGATAGGATCCAGACCGGGGATCAGTTCCCAGCCAAGCAGGAATTTc                                                         <  1:1367232/69‑1 (MQ=255)
           tttATCGCCAGGGGCCATAAAGATAGGATCCAGACCGGGGATCAGTTCCCAGCCAAGCAGGAATTTCAc                                                       >  1:1747647/1‑69 (MQ=255)
           tttATCGCCAGGGGCCATAAAGATAGGATCCAGACCGGGGATCAGTTCCCAGCCAAGCAGGAATTTCAc                                                       >  1:1535757/1‑69 (MQ=255)
                 gCCAGGGGCCATAAAGATAGGATCCAGACCGGGGATCAGTTCCCAGCCAAGCAGGAATTTCACTACCGc                                                 >  1:3368201/1‑69 (MQ=255)
                 gCCAGGGGCCATAAAGATAGGATCCAGACCGGGGATCAGTTCCCAGCCAAGCAGGAATTTCACTACCGc                                                 >  1:3338354/1‑69 (MQ=255)
                 gCCAGGGGCCATAAAGATAGGATCCAGACCGGGGATCAGTTCCCAGCCAAGCAGGAATTTCACTACCGc                                                 >  1:644159/1‑69 (MQ=255)
                       ggCCATAAAGATAGGATCCAGACCGGGGATCAGTTCCCAGCCAAGCAGGAATTTCACTACCGCAGCGGc                                           <  1:156607/69‑1 (MQ=255)
                       ggCCATAAAGATAGGATCCAGACCGGGGATCAGTTCCCAGCCAAGCAGGAATTTCACTACCGCAGCGGc                                           <  1:1570667/69‑1 (MQ=255)
                       ggCCATAAAGATAGGATCCAGACCGGGGATCAGTTCCCAGCCAAGCAGGAATTTCACTACCGCAGCGGc                                           <  1:2239291/69‑1 (MQ=255)
                        gCCATAAAGATAGGATCCAGACCGGGGATCAGTTCCCAGCCAAGCAGGAATTTCACTACCGCAGCg                                             >  1:1815400/1‑66 (MQ=255)
                             aaaGATAGGATCCAGACCGGGGATCAGTTCCCAGCCAAGCAGGAATTTCACTACCGCAGCGGCaa                                         >  1:2429667/1‑65 (MQ=255)
                             aaaGATAGGATCCAGACCGGGGATCAGTTCCCAGCCAAGCAGGAATTTCACTACCGCAGCGGCa                                          >  1:2822939/1‑64 (MQ=255)
                             aaaGATAGGATCCAGACCGGGGATCAGTTCCCAGCCAAGCAGGAATTTCACTACCGCAGCGGCAAg                                        >  1:1576962/1‑66 (MQ=255)
                             aaaGATAGGATCCAGACCGGGGATCAGTTCCCAGCCAAGCAGGAATTTCACTACCGCAGCGGCAAg                                        >  1:2172434/1‑66 (MQ=255)
                             aaaGATAGGATCCAGACCGGGGATCAGTTCCCAGCCAAGCAGGAATTTCACTACCGCAGCGGCAAg                                        >  1:1139332/1‑66 (MQ=255)
                             aaaGATAGGATCCAGACCGGGGATCAGTTCCCAGCCAAGCAGGAATTTCACTACCGCAGCGGCAAg                                        >  1:683212/1‑66 (MQ=255)
                             aaaGATAGGATCCAGACCGGGGATCAGTTCCCAGCCAAGCAGGAATTTCACTACCGCAGCGGCAAg                                        >  1:3235726/1‑66 (MQ=255)
                             aaaGATAGGATCCAGACCGGGGATCAGTTCCCAGCCAAGCAGGAATTTCACTACCGCAGCGGCAAg                                        >  1:2776538/1‑66 (MQ=255)
                             aaaGATAGGATCCAGACCGGGGATCAGTTCCCAGCCAAGCAGGAATTTCACTACCGCAGCGGCAAg                                        >  1:1371430/1‑66 (MQ=255)
                              aaGATAGGATCCGGACCGGGGATCAGTTCCCAGCCAAGCAGGAATTTCACTACCGCAGCGGCAAGACCg                                    >  1:1887195/1‑69 (MQ=255)
                              aaGATAGGATCCAGACCGGGGATCAGTTCCCAGCCAAGCAGGAATTTCACTACCGCAGCGGCAAGACCg                                    >  1:3280911/1‑69 (MQ=255)
                              aaGATAGGATCCAGACCGGGGATCAGTTCCCAGCCAAGCAGGAATTTCACTACCGCAGCGGCAAGACCg                                    >  1:2656518/1‑69 (MQ=255)
                              aaGATAGGATCCAGACCGGGGATCAGTTCCCAGCCAAGCAGGAATTTCACTACCGCAGCGGCAAGACCg                                    >  1:1014896/1‑69 (MQ=255)
                              aaGATAGGATCCAGACCGGGGATCAGTTCCCAGCCAAGCAGGAATTTCACTACCGCAGCGGCAAGACCg                                    >  1:2103817/1‑69 (MQ=255)
                                    ggATCCAGACCGGGGATCAGTTCCCAGCCAAGCAGGAATTTCACTACCGCAGCGGCAAGACCGAGGGTg                              >  1:3093873/1‑69 (MQ=255)
                                    ggATCCAGACCGGGGATCAGTTCCCAGCCAAGCAGGAATTTCACTACCGCAGCGGCAAGACCGAGGGTg                              >  1:2679847/1‑69 (MQ=255)
                                    ggATCCAGACCGGGGATCAGTTCCCAGCCAAGCAGGAATTTCACTACAGCAGCGGCAAGACCGAGGGTg                              >  1:1302316/1‑69 (MQ=255)
                                       tCCAGACCGGGGATCAGTTCCCAGCCAAGCAGGAATTTCACTACCGCAGCGGCAAGACCGAGGGTGAt                            >  1:2880336/1‑68 (MQ=255)
                                            aCCGGTGATCAGTTCCCAGCCAAGCAGGAATTTCACTACCGCAGCGGCAAGACCGAGGGTGATCATTGc                      <  1:77495/69‑1 (MQ=255)
                                            aCCGGGGATCAGTTCCCAGCCAAGCAGGAATTTCACTACCGCAGCGGCAAGACCGAGGGTGATCAATGc                      <  1:1308541/69‑1 (MQ=255)
                                            aCCGGGGATCAGTTCCCAGCCAAGCAGGAATTTCACTACCGCAGCGGCAAGACCGAGGGTGATCAATGc                      <  1:3204734/69‑1 (MQ=255)
                                            aCCGGGGATCAGTTCCCAGCCAAGCAGGAATTTCACTACCGCAGCGGCAAGACCGAGGGTGATCAATGc                      <  1:1731772/69‑1 (MQ=255)
                                            aCCGGGGATCAGTTCCCAGCCAAGCAGGAATTTCACTACCGCAGCGGCAAGACCGAGGGTGATCAATGc                      <  1:1666463/69‑1 (MQ=255)
                                                      aGTTCCCAGCCAAGCAGGAATTTCACTACCGCAGCGGCAAGACCGAGGGTGATCGATGCAACGAGGAAt            >  1:3145556/1‑69 (MQ=255)
                                                      aGTTCCCAGCCAAGCAGGAATTTCACTACCGCAGCGGCAAGACCGAGGGTGATCAATGCAACGAGGAAt            >  1:1504234/1‑69 (MQ=255)
                                                      aGTTCCCAGCCAAGCAGGAATTTCACTACCGCAGCGGCAAGACCGAGGGTGATCAATGCAACGAGGAAt            >  1:1482393/1‑69 (MQ=255)
                                                      aGTTCCCAGCCAAGCAGGAATTTCACTACCGCAGCGGCAAGACCGAGGGTGATCAATGCAACGAGGAAt            >  1:2765141/1‑69 (MQ=255)
                                                      aGTTCCCAGCCAAGCAGGAATTTCACTACCGCAGCGGCAAGACCGAGGGTGATCAATGCAACGAGGAAt            >  1:2251644/1‑69 (MQ=255)
                                                      aGTTCCCAGCCAAGCAGGAATTTCACTACCGCAGCGGCAAGACCGAGGGTGATCAATGCAACGAGGAAt            >  1:1609553/1‑69 (MQ=255)
                                                      aGTTCCCAGCCAAGCAGGAATTTCACTACCGCAGCGGCAAGACCGAGGGTGATCAATGCAACGAGGAAt            >  1:2850349/1‑69 (MQ=255)
                                                      aGTTCCCAGCCAAGCAGGAATTTCACTACCGCAGCGGCAAGACCGAGGGTGATCAATGCAACGAGGAAt            >  1:747268/1‑69 (MQ=255)
                                                      aGTTCCCAGCCAAGCAGGAATTTCACTACCGCAGCGGCAAGACCGAGGGTGATCAATGCAACGAGGAAt            >  1:1760879/1‑69 (MQ=255)
                                                      aGTTCCAAGCCAAGCAGGAATTTCACTACCGCAGCGGCAAGACCGAGGGTGATCAATGCAACGAGGAAt            >  1:1992875/1‑69 (MQ=255)
                                                             aGCCAAGGAGGAATTTCACTACCGCAGCGGCAAGACCGAGGGTGATCAATGCAACGAGGAATTTGGCg      <  1:2201496/68‑1 (MQ=255)
                                                             aGCCAAGGAGGAATTTCACTACCGCAGCGGCAAGACCGAGGGTGATCAATGCAACGAGGAATTTGGCg      <  1:1845437/68‑1 (MQ=255)
                                                             aGCCAAGGAGGAATTTCACTACCGCAGCGGCAAGACCGAGGGTGATCAATGCAACGAGGAATTTGGCg      <  1:3078543/68‑1 (MQ=255)
                                                             aGCCAAGGAGGAATTTCACTACCGCAGCGGCAAGACCGAGGGTGATCAATGCAACGAGGAATTTGGCg      <  1:2948097/68‑1 (MQ=255)
                                                             aGCCAAGGAGGAATTTCACTACCGCAGCGGCAAGACCGAGGGTGATCAATGCAACGAGGAATTTGGCg      <  1:315363/68‑1 (MQ=255)
                                                             aGCCAAGGAGGAATTTCACTACCGCAGCGGCAAGACCGAGGGTGATCAATGCAACGAGGAATTTGGCg      <  1:2915654/68‑1 (MQ=255)
                                                             aGCCAAGGAGGAATTTCACTACCGCAGCGGCAAGACCGAGGGTGATCAATGCAACGAGGAATTTGGCg      <  1:214888/68‑1 (MQ=255)
                                                             aGCCAAGGAGGAATTTCACTACCGCAGCGGCAAGACCGAGGGTGATCAATGCAACGAGGAATTTGGCg      <  1:996202/68‑1 (MQ=255)
                                                             aGCCAAGGAGGAATTTCACTACCGCAGCGGCAAGACCGAGGGTGATCAATGCAACGAGGAATTTGGCg      <  1:883437/68‑1 (MQ=255)
                                                             aGCCAAGGAGGAATTTCACTACCGCAGCGGCAAGACCGAGGGTGATCAATGCAACGAGGAATTTGGCg      <  1:3320123/68‑1 (MQ=255)
                                                             aGCCAAGGAGGAATTTCACTACCGCAGCGGCAAGACCGAGGGTGATCAATGCAACGAGGAATTTGGCg      <  1:2671450/68‑1 (MQ=255)
                                                             aGCCAAGGAGGAATTTCACTACCGCAGCGGCAAGACCGAGGGTGATCAATGCAACGAGGAATTTGGCg      <  1:1438765/68‑1 (MQ=255)
                                                             aGCCAAGGAGGAATTTCACTACCGCAGCGGCAAGACCGAGGGTGATCAATGCAACGAGGAATTTGGCg      <  1:2623792/68‑1 (MQ=255)
                                                             aGCCAAGGAGGAATTTCACTACCGCAGCGGCAAGACCGAGGGTGATCAATGCAACGAGGAATTTGGCg      <  1:2592016/68‑1 (MQ=255)
                                                              gCCAAGGAGGAATTTCACTACCGCAGCGGCAAGACCGAGGGTGATCAATGCAACGAGGAATTTGGCg      <  1:818112/67‑1 (MQ=255)
                                                               ccAAGGAGGAATTTCACTACCGCAGCGGCAAGACCGAGGGTGATCAAt                        >  1:2787324/1‑48 (MQ=255)
                                                               ccAAGGAGGAATTTCACTACCGCAGCGGCAAGACCGAGGGTGATCAATGCAACGAGGAATTTgg        >  1:3403810/1‑64 (MQ=255)
                                                               ccAAGGAGGAATTTCACTACCGCAGCGGCAAGACCGAGGGTGATCAATGCAACGAGGAATTTgg        >  1:1100676/1‑64 (MQ=255)
                                                               ccAAGGAGGAATTTCACTACCGCAGCGGCAAGACCGAGGGTGATCAATGCAACGAGGAATTTgg        >  1:1246628/1‑64 (MQ=255)
                                                               ccAAGGAGGAATTTCACTACCGCAGCGGCAAGACCGAGGGTGATCAATGCAACGAGGAATTTgg        >  1:1258234/1‑64 (MQ=255)
                                                               ccAAGGAGGAATTTCACTACCGCAGCGGCAAGACCGAGGGTGATCAATGCAACGAGGAATTTgg        >  1:1304176/1‑64 (MQ=255)
                                                               ccAAGGAGGAATTTCACTACCGCAGCGGCAAGACCGAGGGTGATCAATGCAACGAGGAATTTgg        >  1:1337091/1‑64 (MQ=255)
                                                               ccAAGGAGGAATTTCACTACCGCAGCGGCAAGACCGAGGGTGATCAATGCAACGAGGAATTTgg        >  1:1393518/1‑64 (MQ=255)
                                                               ccAAGGAGGAATTTCACTACCGCAGCGGCAAGACCGAGGGTGATCAATGCAACGAGGAATTTgg        >  1:1550612/1‑64 (MQ=255)
                                                               ccAAGGAGGAATTTCACTACCGCAGCGGCAAGACCGAGGGTGATCAATGCAACGAGGAATTTgg        >  1:2925615/1‑64 (MQ=255)
                                                               ccAAGGAGGAATTTCACTACCGCAGCGGCAAGACCGAGGGTGATCAATGCAACGAGGAATTTgg        >  1:2036539/1‑64 (MQ=255)
                                                               ccAAGGAGGAATTTCACTACCGCAGCGGCAAGACCGAGGGTGATCAATGCAACGAGGAATTTgg        >  1:2100541/1‑64 (MQ=255)
                                                               ccAAGGAGGAATTTCACTACCGCAGCGGCAAGACCGAGGGTGATCAATGCAACGAGGAATTTgg        >  1:333717/1‑64 (MQ=255)
                                                               ccAAGGAGGAATTTCACTACCGCAGCGGCAAGACCGAGGGTGATCAATGCAACGAGGAATTTgg        >  1:2476617/1‑64 (MQ=255)
                                                               ccAAGGAGGAATTTCACTACCGCAGCGGCAAGACCGAGGGTGATCAATGCAACGAGGAATTTgg        >  1:3313349/1‑64 (MQ=255)
                                                               ccAAGGAGGAATTTCACTACCGCAGCGGCAAGACCGAGGGTGATCAATGCAACGAGGAATTTgg        >  1:2576290/1‑64 (MQ=255)
                                                               ccAAGGAGGAATTTCACTACCGCAGCGGCAAGACCGAGGGTGATCAATGCAACGAGGAATTTgg        >  1:2663391/1‑64 (MQ=255)
                                                               ccAAGGAGGAATTTCACTACCGCAGCGGCAAGACCGAGGGTGATCAATGCAACGAGGAATTTg         >  1:997490/1‑63 (MQ=255)
                                                               ccAAGGAGGAATTTCACTACCGCAGCGGCAAGACCGAGGGTGATCAATGCAACGAGGAATTTg         >  1:1733582/1‑63 (MQ=255)
                                                               caaaGGAGGAATTTCACTACCGCAGCGGCAAGACCGAGGGTGATCa                          >  1:807077/3‑46 (MQ=39)
                                                                cAAGCAGGAATTTCACTACCGCAGCGGCAAGACCGAGGGTGATCAATGCAACGAGGAATTTGGCGAAGa  <  1:642592/69‑1 (MQ=255)
                                                                cAACCAGGAATTTCACTACCGCAGCGGCAAGACCGAGGGTGATCAATGCAACGAGGAATTTGGCGAAGa  <  1:1044140/69‑1 (MQ=255)
                                                                    |                                                                
ACCTCACCGGGTTTATCGCCAGGGGCCATAAAGATAGGATCCAGACCGGGGATCAGTTCCCAGCCAAGCAGGAATTTCACTACCGCAGCGGCAAGACCGAGGGTGATCAATGCAACGAGGAATTTGGCGAAGA  >  minE/1613567‑1613699

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: