Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,620,900 C→A 47.4% V124V (GTG→GTT eutQ ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,620,9000CA47.4% 12.4 / 61.7 59V124V (GTG→GTTeutQconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (28/3);  new base A (10/18);  total (38/21)
Fisher's exact test for biased strand distribution p-value = 2.27e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.07e-01

CAGTCAAGCCCACGCAATGCGGCTCTGCGCCGTCGAAACGACCAAACTTGACGCTACTGCCATCAATCACTTTGATGCCGCCTTTGCCGGTCACCGATTTGAAGCTCGGCTGCATTGCGCCCTGCTCCAGCGA  >  minE/1620833‑1620965
                                                                   |                                                                 
cAGTCAAGCCCACGCAATGCGGCTCTGCGCCGTCGAAACGACCAAACTTGACGCTACTGCCatcaatca                                                                  >  1:444957/1‑69 (MQ=255)
         ccACGCAATGCGGCTCTGCGCCGTCGAAACGACCAAACTTGACGCTACTGCCATCAATCACTTTGATgc                                                         <  1:1287546/69‑1 (MQ=255)
         ccACGCAATGCGGCTCTGCGCCGTCGAAACGACCAAACTTGACGCTACTGCCATCAATCACTTTGATgc                                                         <  1:302631/69‑1 (MQ=255)
             gCAATGCGGCTCTGCGCCGTCGAAACGACCAAACTTGACGCTACTGCCATCAATCACTTTGATgccgcc                                                     <  1:1715299/69‑1 (MQ=255)
                  gCGGTTGTGCGCCGTCGAAACGACCAAACTTGACGCTGCTGCCATCAATAACTTTGATGCCGCCTTTGc                                                <  1:2088529/69‑1 (MQ=255)
                  gCGGTTGTGCGCCGTCGAAACGACCAAACTTGACGCTGCTGCCATCAATAACTTTGATGCCGCCTTTGc                                                <  1:344602/69‑1 (MQ=255)
                  gCGGTTGTGCGCCGTCGAAACGACCAAACTTGACGCTGCTGCCATCAATAACTTTGATGCCGCCTTTGc                                                <  1:3397466/69‑1 (MQ=255)
                  gCGGTTGTGCGCCGTCGAAACGACCAAACTTGACGCTGCTGCCATCAATAACTTTGATGCCGCCTTTGc                                                <  1:2154487/69‑1 (MQ=255)
                  gCGGTTGTGCGCCGTCGAAACGACCAAACTTGACGCTGCTGCCATCAATAACTTTGATGCCGCCTTTGc                                                <  1:3293319/69‑1 (MQ=255)
                  gCGGTTGTGCGCCGTCGAAACGACCAAACTTGACGCTGCTGCCATCAATAACTTTGATGCCGCCTTTGc                                                <  1:705601/69‑1 (MQ=255)
                  gCGGTTGTGCGCCGTCGAAACGACCAAACTTGACGCTGCTGCCATCAATAACTTTGATGCCGCCTTTGc                                                <  1:789157/69‑1 (MQ=255)
                  gCGGTTGTGCGCCGTCGAAACGACCAAACTTGACGCTGCTGCCATCAATAACTTTGATGCCGCCTTTGc                                                <  1:160381/69‑1 (MQ=255)
                  gCGGTTGTGCGCCGTCGAAACGACCAAACTTGACGCTGCTGCCATCAATAACTTTGATGCCGCCTTTGc                                                <  1:1543700/69‑1 (MQ=255)
                  gCGGTTGTGCGCCGTCGAAACGACCAAACTTGACGCTGCTGCCATCAATAACTTTGATGCCGCCTTTGc                                                <  1:1535369/69‑1 (MQ=255)
                  gCGGTTGTGCGCCGTCGAAACGACCAAACTTGACGCTGCTGCCATCAATAACTTTGATGCCGCCTTTGc                                                <  1:845436/69‑1 (MQ=255)
                  gCGGTTGTGCGCCGTCGAAACGACCAAACTTGACGCTGCTGCCATCAATAACTTTGATGCCGCCTTTGc                                                <  1:2991575/69‑1 (MQ=255)
                  gCGGTTGTGCGCCGTCGAAACGACCAAACTTGACGCTGCTGCCATCAATAACTTTGATGCCGCCTTTGc                                                <  1:144704/69‑1 (MQ=255)
                  gCGGTTGTGCGCCGTCGAAACGACCAAACTTGACGCTGCTGCCATCAATAACTTTGATGCCGCCTTTGc                                                <  1:2393324/69‑1 (MQ=255)
                  gCGGTTGTGCGCCGTCGAAACGACCAAACTTGACGCTGCTGCCATCAATAACTTTGATGCCGCCTTTGc                                                <  1:2612399/69‑1 (MQ=255)
                  gCGGCTCTGCGCCGTCGAAACGACCAAACTTGACGCTACTGCCATCAATCACTTTGATGCCGCCTTTGc                                                >  1:2261177/1‑69 (MQ=255)
                  gCGGCTCTGCGCCGTCGAAACGACCAAACTTGACGCTACTGCCATCAATCACTTTGATGCCGCCTTTGc                                                >  1:2402135/1‑69 (MQ=255)
                  gCGGCTCTGCGCCGTCGAAACGACCAAACTTGACGCTACTGCCATCAATCACTTTGATGCCGCCTTTGc                                                >  1:1874625/1‑69 (MQ=255)
                  gCGGCTCTGCGCCGTCGAAACGACCAAACTTGACGCTACTGCCATCAATCACTTTGATGCCGCCTTTGc                                                >  1:492798/1‑69 (MQ=255)
                  gCGGCTCTGCGCCGTCGAAACGACCAAACTTGACGCTACTGCCATCAATCACTTTGATGCCGCCTTTGc                                                >  1:1844988/1‑69 (MQ=255)
                  gCGGCTCTGCGCCGTCGAAACGACCAAACTTGACGCTACTGCCATCAATCACTTTGATGCCGCCTTTGc                                                >  1:1466429/1‑69 (MQ=255)
                  gCGGCTCTGCGCCGTCGAAACGACCAAACTTGACGCTACTGCCATCAATCACTTTGATGCCGCCTTTGc                                                >  1:861215/1‑69 (MQ=255)
                  gCGGCTCTGCGCCGTCGAAACGACCAAACTTGACGCTACTGCCATCAATCACTTTGATGCCGCCTTTGc                                                >  1:1294986/1‑69 (MQ=255)
                  gCGGCTCTGCGCAGTCGAAACGACCAAACTTGACGCTACTGCCATCAATCACTTTGATGCCGCCTTTGc                                                >  1:689561/1‑69 (MQ=255)
                                  gAAACGACCAAACTTGACGCTTCTGCCATCAATCACTTTGATGCCGCCTTTGCCGGTCACCGATTTGaa                                >  1:3022885/1‑69 (MQ=255)
                                  gAAACGACCAAACTTGACGCTACTGCCATCAATCACTTTGATGCCGCCt                                                    >  1:2456102/1‑49 (MQ=255)
                                  gAAACGACCAAACTTGACGCTACTGCCATCAATCACTTTGATGCCGCCTTTGCCGGTCACCGATTTGaa                                >  1:290462/1‑69 (MQ=255)
                                  gAAACGACCAAACTTGACGCTACTGCCATCAATCACTTTGATGCCGCCTTTGCCGGTCACCGATTTGaa                                >  1:509079/1‑69 (MQ=255)
                                  gAAACGACCAAACTTGACGCTACTGCCATCAATCACTTTGATGCCGCCTTTGCCGGTCACCGATTTGaa                                >  1:29260/1‑69 (MQ=255)
                                  gAAACGACCAAACTTGACGCTACTGCCATCAATCACTTTGATGCCGCCTTTGCCGGTCACCGATTTGaa                                >  1:907879/1‑69 (MQ=255)
                                  gAAACGACCAAACTTGACGCTACTGCCATCAATCACTTTGATGCCGCCTTTGCCGGTCACCGATTTGaa                                >  1:1013994/1‑69 (MQ=255)
                                  gAAACGACCAAACTTGACGCTACTGCCATCAATCACTTTGATGCCGCCTTTGCCGGTCACCGATTTGaa                                >  1:1519100/1‑69 (MQ=255)
                                  gAAACGACCAAACTTGACGCTACTGCCATCAATCACTTTGATGCCGCCTTTGCCGGTCACCGATTTGaa                                >  1:2514313/1‑69 (MQ=255)
                                  gAAACGACCAAACTTGACGCTACTGCCATCAATCACTTTGATGCCGCCTTTGCCGGTCACCGATTTGaa                                >  1:198648/1‑69 (MQ=255)
                                  gAAACGACCAAACTTGACGCTACTGCCATCAATCACTTTGATGCCGCCTTTGCCGGTCACCGATTTGaa                                >  1:134466/1‑69 (MQ=255)
                                  gAAACGACCAAACTTGACGCTACTGCCATCAATCACTTTGATGCCGCCTTTGCCGGTCACCGATTTGa                                 >  1:2489715/1‑68 (MQ=255)
                                  gAAACGACCAAACTTGACGCTACTGCCATCAATCACTTTGATGCCGCCTTTGCCGGTCACCGATTTGa                                 >  1:1077956/1‑68 (MQ=255)
                                  gAAACGACCAAACTTGACGCTACTGCCATCAATCACTTTGATCCCGCCTTTGCCGGTCACCGATTTGa                                 >  1:1856517/1‑68 (MQ=255)
                                           aaaCTTGACGCTGCTGCCATCAATAACTTTGATGCCGCCTTTGc                                                <  1:3247464/44‑1 (MQ=38)
                                            aaCTTGACGCTGCTGCCATCAATAATTTTGATGCCGCCTTTGCCGGTCACCGATTTAAAGc                              >  1:32838/1‑61 (MQ=255)
                                            aaCTTGACGCTGCTGCCATCAATAACTTTGATGCCGCCTTTGCCGGTCACCGATTTAAAg                               >  1:2285500/1‑60 (MQ=255)
                                            aaCTTGACGCTGCTGCCATCAATAACTTTGATGCCGCCTTTGCCGGTCACCGATTTAAAGc                              >  1:3193420/1‑61 (MQ=255)
                                            aaCTTGACGCTGCTGCCATCAATAACTTTGATGCCGCCTTTGCCGGTCACCGATTTAAAGc                              >  1:3215100/1‑61 (MQ=255)
                                            aaCTTGACGCTGCTGCCATCAATAACTTTGATGCCGCCTTTGCCGGTCACCGATTTAAAGc                              >  1:2281609/1‑61 (MQ=255)
                                            aaCTTGACGCTGCTGCCATCAATAACTTTGATGCCGCCTTTGCCGGTCACCGATTTAAAGc                              >  1:2103939/1‑61 (MQ=255)
                                            aaCTTGACGCTGCTGCCATCAATAACTTTGATGCCGCCTTTGCCGGTCACCGATTTAAAGc                              >  1:2028218/1‑61 (MQ=255)
                                            aaCTTGACGCTGCTGCCATCAATAACTTTGATGCCGCCTTTGCCGGTCACCGATTTAAAGc                              >  1:2369835/1‑61 (MQ=255)
                                            aaCTTGACGCTGCTGCCATCAATAACTTTGATGCCGCCTTTGCCGGTCACCGATTTAAAGc                              >  1:839475/1‑61 (MQ=255)
                                            aaCTTGACGCTGATGCCATCAATAACTTTGATGCCGCCTTTGCCGGTCACCGATTTAAAGc                              >  1:2288352/1‑61 (MQ=255)
                                               ttGACGCTGCTGCCATCAATAACTTTGTTGCCGCCTTTGc                                                <  1:3029402/40‑1 (MQ=37)
                                                  aCGCTGCTGCCATCAATAACTTTGATGCCGCCTTTGc                                                <  1:183039/37‑1 (MQ=38)
                                                          gCCATCAATCACTTTGATGCCGCCTTTGCCGGTCAc                                         >  1:188968/1‑36 (MQ=255)
                                                          gCCATCAATCACTTTGATGCCGCCTTTGCCGGTCACCGa                                      >  1:1818554/1‑39 (MQ=255)
                                                          gCCATCAATCACTTTGATGCCGCCTTTGCCGGTCACCGa                                      >  1:2624385/1‑39 (MQ=255)
                                                           ccATCAATCACTTTGATGCCGCCTTTGCCGGTCACCGATTTGAAGCTCGGCTGCATTGCGCCCTGCTc        >  1:1932577/1‑68 (MQ=255)
                                                           ccATCAATCACTTTGATGCCGCCTTTGCCGGTCACCGATTTGAAGCACGGCTGCATTGCGCCCTGCTcc       >  1:2711168/1‑69 (MQ=255)
                                                               caatcaCTTTGATGCCGCCTTTGCCGGTCACCGATTTGAAGCTCGGCTGCATTGCGCCCTGCTCCAGCg   <  1:3360010/69‑1 (MQ=255)
                                                               caatcaCTTTGATGCCGCCTTTGCCGGTCACCGATTTGAAGCTCGGCTGCATTGCGCCCTGCTCCAGCg   <  1:3118942/69‑1 (MQ=255)
                                                                aatcaCTTTGATGCCGCCTTTGCCGGTCACCGATTTGAAGCTCGGCTGCATTGCGCCCTGCTCCAGCGa  >  1:843683/1‑69 (MQ=255)
                                                                   |                                                                 
CAGTCAAGCCCACGCAATGCGGCTCTGCGCCGTCGAAACGACCAAACTTGACGCTACTGCCATCAATCACTTTGATGCCGCCTTTGCCGGTCACCGATTTGAAGCTCGGCTGCATTGCGCCCTGCTCCAGCGA  >  minE/1620833‑1620965

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: