Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 612,948 C→T 31.1% A128A (GCG→GCA ycaO ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE612,9480CT31.1% 57.2 / 29.3 64A128A (GCG→GCAycaOconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (33/11);  new base T (12/8);  total (45/19)
Fisher's exact test for biased strand distribution p-value = 2.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.97e-01

ACCGGATTGTAGGTCAATCAGCATGCTGCCGGTCAGTTCATTCTCCGGATCATAAAACGCGCGCAGACGGTCATCGAGCAGCCCTTCTGGCACATCGTCATTTTCGGTCAGTGGGAACCATTTTTC  >  minE/612891‑613016
                                                         |                                                                    
aCCGGATTGCAGGTCAATCAGCATACTACCGGTCAGTTCATTCTCCGGATCGTAAAATGCGCGCAGCCg                                                           >  1:1498035/1‑69 (MQ=255)
aCCGGATTGCAGGTCAATCAGCATACTACCGGTCAGTTCATTCTCCGGATCGTAAAATGCGCGCAGACg                                                           >  1:3255686/1‑69 (MQ=255)
aCCGGATTGCAGGTCAATCAGCATACTACCGGTCAGTTCATTCTCCGGATCGTAAAATGCGCGCAGACg                                                           >  1:277342/1‑69 (MQ=255)
aCCGGATTGCAGGTCAATCAGCATACTACCGGTCAGTTCATTCTCCGGATCGTAAAATGCGCGCAGACg                                                           >  1:3053805/1‑69 (MQ=255)
aCCGGATTGCAGGTCAATCAGCATACTACCGGTCAGTTCATTCTCCGGATCGTAAAATGCGCGCAGACg                                                           >  1:155391/1‑69 (MQ=255)
aCCGGATTGCAGGTCAATCAGCATACTACCGGTCAGTTCATTCTCCGGATCGTAAAATGCGCGCAGACg                                                           >  1:1446605/1‑69 (MQ=255)
aCCGGATTGCAGGTCAATCAGCATACTACCGGTCAGTTCATTCTCCGGATCGTAAAATGCGCGCAGACg                                                           >  1:3200841/1‑69 (MQ=255)
aCCGGATTGCAGGTCAATCAGCATACTACCGGTCAGTTCATTCTCCGGATCGTAAAATGCGCGCAGACg                                                           >  1:3207657/1‑69 (MQ=255)
aCCGGATTGCAGGTCAATCAGCATACTACCGGTCAGTTCATTCTCCGGATCGTAAAATGCGCGCAGACg                                                           >  1:1813672/1‑69 (MQ=255)
aCCGGATTGCAGGTCAATCAGCATACTACCGGTCAGTTCATTCTCCGGATCGTAAAATGCGCGCAGACg                                                           >  1:2356485/1‑69 (MQ=255)
aCCGGATTGCAGGTCAATCAGCATACTACCGGTCAGTTCATTCTCCGGATCGTAAAATGCGCGCAGACg                                                           >  1:2226707/1‑69 (MQ=255)
aCCGGATTGCAGGTCAATCAGCATACTACCGGTCAGTTCATTCTCCGGATCGTAAAATGCGCGCAGACg                                                           >  1:425304/1‑69 (MQ=255)
           ggTCAATCAGCATGCTGCCGGTCAGTTCATTCTCCGGATCATAAAACGCGCGCAGCCGGTCATCGag                                                  >  1:3277102/1‑67 (MQ=255)
           ggTCAATCAGCATGCTGCCGGTCAGTTCATTCTCCGGATCATAAAACGCGCGCAGACGGTCATCGagca                                                >  1:1277258/1‑69 (MQ=255)
           ggTCAATCAGCATGCTGCCGGTCAGTTCATTATCCGGATCATAAAACGCGCGCAGACGGTCATCGagca                                                >  1:1231903/1‑69 (MQ=255)
                     cATGCTGCCGGTCAGTTCATTCTCCGGATCATAAAACGCGCGCAg                                                              >  1:626100/1‑45 (MQ=255)
                     cATGCTGCCGGTCAGTTCATTCTCCGGATCATAAAACGCGCGCAGACGGTCATCGAGCAGCCCTTCTgg                                      >  1:3168312/1‑69 (MQ=255)
                     cATGCTGCCGGTCAGTTCATTCTCCGGATCATAAAACGCGCGCAGACGGTCATCGAGCAGCCCTTCTgg                                      >  1:180593/1‑69 (MQ=255)
                     cATGCTGCCGGTCAGTTCATTCTCCGGATCATAAAACGCGCGCAGACGGTCATCGAGCAGCCCTTCTgg                                      >  1:2471259/1‑69 (MQ=255)
                     cATGCTGCCGGTCAGTTCATTCTCCGGATCATAAAACGCGCGCAGACGGTCATCGAGCAGCCCTTCTgg                                      >  1:2221094/1‑69 (MQ=255)
                     cATGCTGCCGGTCAGTTCATTCTCCGGATCATAAAACGCGCGCAGACGGTCATCGAGCAGCCCTTCTgg                                      <  1:2061171/69‑1 (MQ=255)
                     cATGCTGCCGGTCAGTTCATTCTCCGGATCATAAAACGCGCGCAGACGGTCATCGAGCAGCCCTTCTgg                                      >  1:1531783/1‑69 (MQ=255)
                     cATGCTGCCGGTCAGTTCATTCTCCGGATCATAAAACGCGCGCAGACGGTCATCGAGCAGCCCTTCTgg                                      <  1:143878/69‑1 (MQ=255)
                     cATGCTGCCGGTCAGTTCATTCTCCGGATCATAAAACGCGCGCAGACGGTCATCGAGCAGCCCTTCTgg                                      >  1:1174477/1‑69 (MQ=255)
                     cATGCTGCCGGTCAGTTCATTCTCCGGATCATAAAACGCGCGCAGACGGTCATCGAGCAGCCCTTCTgg                                      >  1:350302/1‑69 (MQ=255)
                     cATGCTGCCGGTCAGTTCATTCTCCGGATCATAAAACGCGCGCAGACGGTCATCGAGCAGCCCTTCTgg                                      >  1:560765/1‑69 (MQ=255)
                        gctgcCGGTCAGTTCATTCTCCGGATCATAAAACGCGCGCAGACGGTCATCGag                                                  <  1:2929663/54‑1 (MQ=255)
                        gctgcCGGTCAGTTCATTCTCCGGATCATAAAACGCGCGCAGACGGTCATCGag                                                  <  1:1767226/54‑1 (MQ=255)
                        gctgcCGGTCAGTTCATTCTCCGGATCATAAAACGCGCGCAGACGGTCATCGAGCAGCCCTTCTGGcac                                   <  1:821672/69‑1 (MQ=255)
                                     tCATTCTCCGGATCATAAAACGCGCGCAGACGGTCATc                                                     >  1:103485/1‑38 (MQ=255)
                                     tCATTCTCCGGATCATAAAACGCGCGCAGACGGTCATCGAGCAGCCCTTCTGGCACATCGTCATTTTCg                      >  1:902449/1‑69 (MQ=255)
                                     tCATTCTCCGGATCATAAAACGCGCGCAGACGGTCATCGAGCAGCCCTTCTGGCACATCGTCATTTTCg                      >  1:1603213/1‑69 (MQ=255)
                                     tCATTCTCCGGATCATAAAACGCGCGCAGACGGTCATCGAGCAGCCCTTCTGGCACATCGTCATTTTCg                      >  1:2793241/1‑69 (MQ=255)
                                     tCATTCTCCGGATCATAAAACGCGCGCAGACGGTCATCGAGCAGCCCTTCTGGCACATCGTCATTTTCg                      >  1:2655686/1‑69 (MQ=255)
                                     tCATTCTCCGGATCATAAAACGCGCGCAGACGGTCATCGAGCAGCCCTTCTGGCACATCGTCATTTTCg                      >  1:2444260/1‑69 (MQ=255)
                                     tCATTCTCCGGATCATAAAACGCGCGCAGACGGTCATCGAGCAGCCCTTCTGGCACATCGTCATTTTCg                      >  1:3332061/1‑69 (MQ=255)
                                     tCATTCTCCGGATCATAAAACGCGCGCAGACGGTCATCGAACAGCCCTTCTGGCACATCGTCATTTTc                       >  1:2716672/1‑68 (MQ=255)
                                                   aTAAAACGCGCGCAGACTGTCATCGATCAGCCCTTCTGGCACATCGTCATTTTCGGTCAGTGGTAACCa        <  1:1574725/69‑1 (MQ=255)
                                                   aTAAAACGCGCGCAGACGGTCATCGAGCAGCCCTTCTGGCACATCGTCATTTTCGGTCAGTGGGAACCa        <  1:1157026/69‑1 (MQ=255)
                                                   aTAAAACGCGCGCAGACGGTCATCGAGCAGCCCTTCTGGCACATCGTCATTTTCGGTCAGTGGGAACCa        <  1:2372235/69‑1 (MQ=255)
                                                   aTAAAACGCGCGCAGACGGTCATCGAGCAGCCCTTCTGGCACATCGTCATTTTCGGTCAGTGGGAACCa        <  1:1444329/69‑1 (MQ=255)
                                                   aTAAAACGCGCGCAGACGGTCATCGAGCAGCCCTTCTGGCACATCGTCATTTTCGGTCAGTGGGAACCa        <  1:156683/69‑1 (MQ=255)
                                                   aTAAAACGCGCGCAGACGGTCATCGAGCAGCCCTTCTGGCACATCGTCATTTTCGGTCAGTGGGAACCa        <  1:2970305/69‑1 (MQ=255)
                                                      aaaTGCGCGCAGACGGTCATCGAGCAGCCCTTCTGGtaca                                  <  1:3387609/40‑1 (MQ=38)
                                                      aaaTGCGCGCAGACGGTCATCGAGCAGCCCTTCTGGtaca                                  <  1:1143180/40‑1 (MQ=38)
                                                      aaaTGCGCGCAGACGGTCATCGAGCAGCCCTTCTGGtaca                                  <  1:2795263/40‑1 (MQ=38)
                                                      aaaTGCGCGCAGACGGTCATCGAGCAGCCCTTCTGGtaca                                  <  1:457605/40‑1 (MQ=38)
                                                      aaaTGCGCGCAGACGGTCATCGAGCAGCCCTTCTGGtaca                                  <  1:1084269/40‑1 (MQ=38)
                                                      aaaTGCGCGCAGACGGTCATCGAGCAGCCCTTCTGGtaca                                  <  1:2969469/40‑1 (MQ=38)
                                                      aaaTGCGCGCAGACGGTCATCGAGCAGCCCTTCTGGtaca                                  <  1:1707741/40‑1 (MQ=38)
                                                      aaaTGCGCGCAGACGGTCATCGAGCAGCCCTTCTGGtaca                                  <  1:2314570/40‑1 (MQ=38)
                                                      aaaCGCGCGCAGACGGTCATCGAGCATCCCTTCTGGCACATCGTCATTTTCGGTc                   >  1:2149787/1‑55 (MQ=255)
                                                      aaaCGCGCGCAGACGGTCATCGAGCAGCCCTTCTGGCACATCGTCATTTTCGGTc                   >  1:2742902/1‑55 (MQ=255)
                                                      aaaCGCGCGCAGACGGTCATCGAGCAGCCCTTCTGGCACATCGTCATTTTCGGTc                   >  1:2192035/1‑55 (MQ=255)
                                                      aaaCGCGCGCAGACGGTCATCGAGCAGCCCTTCTGGCACATCGTCATTTTCGGTc                   >  1:1125550/1‑55 (MQ=255)
                                                      aaaCGCGCGCAGACGGTCATCGAGCAGCCCTTCTGGCACATCGTCATTTTCGGTc                   >  1:2007734/1‑55 (MQ=255)
                                                      aaaCGCGCGCAGACGGTCATCGAGCAGCCCTTCTGGCACATCGTCATTTTCGGTc                   >  1:1959644/1‑55 (MQ=255)
                                                      aaaCGCGCGCAGACGGTCATCGAGCAGCCCTTCTGGCACATCGTCATTTTCGGTc                   >  1:1952924/1‑55 (MQ=255)
                                                      aaaCGCGCGCAGACGGTCATCGAGCAGCCCTTCTGGCACATCGTCATTTTCGGTc                   >  1:144299/1‑55 (MQ=255)
                                                      aaaCGCGCGCAGACGGTCATCGAGCAGCCCTTCTGGCACATCGTCATTTTCGGTc                   >  1:2374168/1‑55 (MQ=255)
                                                      aaaCGCGCGCAGACGGTCATCGAGCAGCCCTTCTGGCACATCGTCATTTTCGGTc                   >  1:104742/1‑55 (MQ=255)
                                                      aaaCGCGCGCAGACGGTCATCGAGCAGCCCTTCTGGCACATCGTCATTTTCGGTc                   >  1:2484304/1‑55 (MQ=255)
                                                      aaaCGCGCGCAGACGGTCATCGAGCAGCCCTTCTGGCACATCGTCATTTTCGGTc                   >  1:1023915/1‑55 (MQ=255)
                                                      aaaCGCACGCAGACGGTCATCGAGCAGCCCTTCTGGCACATCGTCATTTTCGCTc                   >  1:436419/1‑55 (MQ=255)
                                                         cgcgcgCAGACGGTCATCGAGCAGCCCTTCTGTCACATCGTCATTTTCGGTCAGTGGGAACCATTTTTc  <  1:1380163/69‑1 (MQ=255)
                                                         cgcgcgCAGACGGTCATCGAGCAGCCCTTCTGGCACATCGTCATTTTCGGTCAGTGGGAACCTTTTTTc  <  1:435708/69‑1 (MQ=255)
                                                         cgcgcgCAGACGGTCATCGAGCAGCCCTTCTGGCACATCGTCATTTTCGGTCAGTGGGAACCATTTTTc  <  1:619375/69‑1 (MQ=255)
                                                         cgcgcgCAGACGGTCATCGAGCAGCCCTTCTGGCACATCGTCATTTTCGGTCAGTGGGAACCATTTTTc  <  1:2188618/69‑1 (MQ=255)
                                                         cgcgcgCAGACGGTCATCGAGCAGCCCTTCTGGCACATCGTCATTTTCGGTCAGTGGGAACCATTTTTc  <  1:2190063/69‑1 (MQ=255)
                                                         cgcgcgCAGACGGTCATCGAGCAGCCCTTCTGGCACATCGTCATTTTCGGTCAGTGGGAACCATTTTTc  <  1:3305878/69‑1 (MQ=255)
                                                         cgcgcgCAGACGGTCATCGAGCAGCCCTTCTGGCACATCGTCATTTTCGGTCAGTGGGAACCATTTTTc  <  1:146636/69‑1 (MQ=255)
                                                         cgcgcgCAGACGGTCATCGAGCAGCCCTTCTGGCACATCGTCATTTTCGGTCAGTGGGAACCATTTTTc  <  1:2506568/69‑1 (MQ=255)
                                                         cgcgcgCAGACGGTCATCGAGCAGCCCTTCTGGCACATCGTCATTTTCGGTCAGTGGGAACCATTTTTc  <  1:259014/69‑1 (MQ=255)
                                                         cgcgcgCAGACGGTCATCGAGCAGCCCTTCTGGCACATCGTCATATTCGGTCAGTGGGAACCATTTTTc  <  1:851501/69‑1 (MQ=255)
                                                         |                                                                    
ACCGGATTGTAGGTCAATCAGCATGCTGCCGGTCAGTTCATTCTCCGGATCATAAAACGCGCGCAGACGGTCATCGAGCAGCCCTTCTGGCACATCGTCATTTTCGGTCAGTGGGAACCATTTTTC  >  minE/612891‑613016

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: