Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,251,528 A→G 35.9% G180G (GGT→GGC yebC ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,251,5280AG35.9% 37.9 / 37.1 64G180G (GGT→GGCyebCconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (16/25);  new base G (8/15);  total (24/40)
Fisher's exact test for biased strand distribution p-value = 7.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

GACGGGATCATGGAAACTTCCGCGCTGTCTGCTTTCAGACCTGCCGCTTCCAGAGCGTCGCGCACTTTACCCATTTCTTCCCATGCAGTGTAGACATCAATCGCGCCGTCATCATAGGTCACAACGTCTTCAGCAC  >  minE/1251460‑1251595
                                                                    |                                                                   
gACGGGATCATGGAAACTTCCGCGCTGTCTGCTTTCAGACCTGCCGCTTCCAGAGCGTCGCGCACTTTa                                                                     <  1:338537/69‑1 (MQ=255)
gACGGGATCATGGAAACTTCCGCGCTGTCTGCTTTCAGACCTGCCGCTTCCAGAGCGTCGCGCACTTTa                                                                     <  1:3063806/69‑1 (MQ=255)
    ggATCATGGAAACTTCCGCGCTGTCTGCTTTCAGACCTGCCGCTTCCAGAGCGTCGCGCACTTTACCCa                                                                 <  1:879880/69‑1 (MQ=255)
    ggATCATGGAAACTTCCGCGCTGTCTGCTTTCAGACCTGCCGCTTCCAGAGCGTCGCGCACTTTACCCa                                                                 <  1:3383772/69‑1 (MQ=255)
                      cgcTGTCTGCTTTCAGACCTGCCGCTTCCAGAGCGTCGCGCACTTTACCCATTTCTTCCCATGCAgtg                                                <  1:2534670/68‑1 (MQ=255)
                      cgcTGTCTGCTTTCAGACCTGCCGCTTCCAGAGCGTCGCGCACTTTACCCATTTCTTCCCATGCAgtg                                                <  1:2865285/68‑1 (MQ=255)
                      cgcTGTCTGCTTTCAGACCTGCCGCTTCCAGAGCGTCGCGCACTTTACCCATTTCTTCCCATGCAgtg                                                <  1:3111216/68‑1 (MQ=255)
                       gcTGTCTGCTTTCAGACCTGCCGCTTCCAGAGCGTCGCGCACTTTACCCATTTCTTCCCATGCAGTGTa                                              <  1:2298594/69‑1 (MQ=255)
                       gcTGTCTGCTTTCAGACCTGCCGCTTCCAGAGCGTCGCGCACTTTACCCATTTCTTCCCATGCAGTGTa                                              <  1:2081944/69‑1 (MQ=255)
                       gcTGTCTGCTTTCAGACCTGCCGCTTCCAGAGCGTCGCGCACTTTACCCATTTCTTCCCATGCAGTGTa                                              <  1:87254/69‑1 (MQ=255)
                       gcTGTCTGCTTTCAGACCTGCCGCTTCCAGAGCGTCGCGCACTTTACCCATTTCTTCCCATGCAGTGTa                                              <  1:2736017/69‑1 (MQ=255)
                       gcTGTCTGCTTTCAGACCTGCCGCTTCCAGAGCGTCGCGCACTTTACCCATTTCTTCCCATGCAGTGTa                                              <  1:583437/69‑1 (MQ=255)
                             tGCTTTCAGACCTGCCGCTTCCAGAGCGTCGCGCACTTTACCCATTTCTTCCCATGCAg                                                  >  1:2674528/1‑59 (MQ=255)
                             tGCTTTCAGACCTGCCGCTTCCAGAGCGTCGCGCACTTTACCCATTTCTTCCCATGCAGTGTAGACATc                                        >  1:807979/1‑69 (MQ=255)
                              gCTTTCAGACCTGCCGCTTCCAGAGCGTCGCGCACTTTa                                                                     <  1:2345162/39‑1 (MQ=255)
                                tttCAGACCTGCCGCTTCCAGAGCGTCGCGTACTTTGccc                                                                  >  1:1092780/1‑40 (MQ=38)
                                tttCAGACCTGCCGCTTCCAGAGCGTCGCGTACTTTGCCCATttctt                                                           >  1:1408420/1‑47 (MQ=255)
                                tttCAGACCTGCCGCTTCCAGAGCGTCGCGTACTTTGCCCATTTCTTCCCATGCAGTGTAGAcct                                         >  1:1408833/1‑63 (MQ=255)
                                tttCAGACCTGCCGCTTCCAGAGCGTCGCGTACTTTGCCCATTTCTTCCCATGCAGTGTAGACATCaa                                      >  1:2880214/1‑68 (MQ=255)
                                tttCAGACCTGCCGCTTCCAGAGCGTCGCGTACTTTGCCCATTTCTTCCCATGCAGTGTAGACATCaa                                      >  1:1052656/1‑68 (MQ=255)
                                tttCAGACCTGCCGCTTCCAGAGCGTCGCGTACTTTGCCCATTTCTTCCCATGCAGTGTAGACATCAAt                                     >  1:26185/1‑69 (MQ=255)
                                tttCAGACCTGCCGCTTCCAGAGCGTCGCGTACTTTGCCCATTTCTTCCCATGCAGGGTAGAc                                           >  1:2921879/1‑63 (MQ=255)
                                tttCAGACCTGCCGCTTCCAGAGCGTCGCGTACGTTGCCCATTTCTTCCCATGCAGTGTAGACATCaa                                      >  1:1163359/1‑68 (MQ=255)
                                        cTGCCGCTTCCAGAGCGTCGCGCACTTTACCCATTTCTTCCCATg                                                     >  1:2718526/1‑45 (MQ=255)
                                        cTGCCGCTTCCAGAGCGTCGCGCACTTTACCCATTTCTTCCCATGCAGTGTAGACATCAATCGCGCCGt                             >  1:501927/1‑69 (MQ=255)
                                        cTGCCGCTTCCAGAGCGTCGCGCACTTTACCCATTTCTTCCCATGCAGTGTAGACATCAATCGCGCCGt                             >  1:60638/1‑69 (MQ=255)
                                        cTGCCGCTTCCAGAGCGTCGCGCACTTTACCCATTTCTTCCCATGCAGTGTAGACATCAATCGCGCCGt                             >  1:2610219/1‑69 (MQ=255)
                                        cTGCCGCTTCCAGAGCGTCGCGCACTTTACCCATTTCTTCCCATGCAGTGTAGACATCAATCGCGCCGt                             >  1:241210/1‑69 (MQ=255)
                                        cTGCCGCTTCCAGAGCGTCGCGCACTTTACCCATTTCTTCCCATGCAGTGTAGACATCAATCGCGCCGt                             >  1:956321/1‑69 (MQ=255)
                                        cTGCCGCTTCCAGAGCGTCGCGCACTTTACCCATTTCTTCCCATGCAGTGTAGACATCAATCGCGCCGt                             >  1:2524176/1‑69 (MQ=255)
                                        cTGCCGCTTCCAGAGCGTCGCGCACTTTACCCATTTCTTCCCATGCAGTGTAGACATCAATCGCGCCGt                             >  1:2250937/1‑69 (MQ=255)
                                        cTGCCGCTTCCAGAGCGTCGCGCACTTTACCCATTTCTTCCCATGCAGTGTAGACATCAATCGCGCCGt                             >  1:2132232/1‑69 (MQ=255)
                                        cTGCCGCTTCCAGAGCGTCGCGCACTTTACCCATTTCTTCCCATGCAGTGTAGACATCAATCGCGCCGt                             >  1:2110294/1‑69 (MQ=255)
                                        cTGCCGCTTCCAGAGCGTCGCGCACTTTACCCATTTCTTCCCATGCAGTGTAGACATCAATCGCGCCGt                             >  1:2087811/1‑69 (MQ=255)
                                        cTGCCGCTTCCAGAGCGTCGCGCACTTTACCCATTTCTTCCCATGCAGTGTAGACATCAATCGCGCCGt                             >  1:2081628/1‑69 (MQ=255)
                                        cTGCCGCTTCCAGAGCGTCGCGCACTTTACCCATTTCTTCCCATGCAGTGTAGACATCAATCGCGCAGt                             >  1:1248512/1‑69 (MQ=255)
                                             gCTTCCAGAGCGTCGCGCACTTTACCCATTTCTTCCCATGCAGTGTAGACATCAATCGCGCCGtcat                          <  1:2295911/67‑1 (MQ=255)
                                             gCTTCCAGAGCGTCGCGCACTTTACCCATTTCTTCCCATGCAGTGTAGACATCAATCGCGCCGtcat                          <  1:1443770/67‑1 (MQ=255)
                                             gCTTCCAGAGCGTCGCGCACTTTACCCATTTCTTCCCATGCAGTGTAGACATCAATCGCGCCGtcat                          <  1:1030361/67‑1 (MQ=255)
                                             gCTTCCAGAGCGTCGCGCACTTTACCCATTTCTTCCCATGCAGTGTAGACATCAATCGCGCCGtcat                          <  1:42656/67‑1 (MQ=255)
                                                           gcgtACTTTGCCCATTTCTTCCCATGCAGTGTAGACATCAATCGCGCCGTCATCATAGGTCACAACGTc          <  1:1240391/69‑1 (MQ=255)
                                                           gcgtACTTTGCCCATTTCTTCCCATGCAGTGTAGACATCAATCGCGCCGTCATCATAGGTCACAACGTc          <  1:3269824/69‑1 (MQ=255)
                                                           gcgtACTTTGCCCATTTCTTCCCATGCAGTGTAGACATCAATCGCGCCGTCATCATAGGTCACAACGTc          <  1:3381606/69‑1 (MQ=255)
                                                           gcgtACTTTGCCCATTTCTTCCCATGCAGTGTAGACATCAATCGCGCCGTCATCATAGGTCACAACGTc          <  1:1395195/69‑1 (MQ=255)
                                                           gcgtACTTTGCCCATTTCTTCCCATGCAGTGTAGACATCAATCGCGCCGTCATCATAGGTCACAACGTc          <  1:1372615/69‑1 (MQ=255)
                                                           gcgtACTTTGCCCATTTCTTCCCATGCAGTGTAGACATCAATCGCGCCGTCATCATAGGTCACAACGTc          <  1:3092190/69‑1 (MQ=255)
                                                           gcgtACTTTGCCCATTTCTTCCCATGCAGTGTAGACATCAATCGCGCCGTCATCATAGGTCACAACGTc          <  1:1231909/69‑1 (MQ=255)
                                                           gcgtACTTTGCCCATTTCTTCCCATGCAGTGTAGACATCAATCGCGCCGTCATCATAGGTCACAACGTc          <  1:10712/69‑1 (MQ=255)
                                                           gcgtACTTTGCCCATTTCTTCCCATGCAGTGTAGACATCAATCGCGCCGTCATCATAGGTCACAACGTc          <  1:1038149/69‑1 (MQ=255)
                                                           gcgtACTTTGCCCATTTCTTCCCATGCAGTGTAGACATCAATCGCGCCGTCATCATAGGTCACAACGTc          <  1:2929102/69‑1 (MQ=255)
                                                           gcgtACTTTGCCCATTTCTTCCCATGCAGTGTAGACATCAATCGCGCCGTCATCATAGGTCACAACGTc          <  1:19721/69‑1 (MQ=255)
                                                           gcgtACTTTGCCCATTTCTTCCCATGCAGTGTAGACATCAATCGCGCCGTCATCATAGGTCACAACGTc          <  1:2002638/69‑1 (MQ=255)
                                                           gcgtACTTTGCCCATTTCTTCCCATGCAGTGTAGACATCAATCGCGCCGTCATCATAGGTCACAACGTc          <  1:2005791/69‑1 (MQ=255)
                                                           gcgtACTTTGCCCATTTCTTCCCATGCAGTGTAGACATCAATCGCGCCGTCATCATAGGTCACAACGTc          <  1:1569999/69‑1 (MQ=255)
                                                           gcgtACTTTGCCCATTTCTTCCCATGCAGTGTAGACATCAATCGCGCCGTCATCATAGGTCACAACGTc          <  1:2661246/69‑1 (MQ=255)
                                                                cTTTACCCATTTCTTCCCATGCAGTGTAGACATCAATCGCGCCGTCATCATAGGTCACAACGTCTTCAg     <  1:619686/69‑1 (MQ=255)
                                                                 tttACCCATTTCTTCCCATGCAGTGTAGACATCAATCGCGCCGTCATCATAGGTCACAACGTCTTCAGc    <  1:2445113/69‑1 (MQ=255)
                                                                 tttACCCATTTCTTCCCATGCAGTGTAGACATCAATCGCGCCGTCATCATAGGTCACAACGTCTTCAGc    <  1:2501895/69‑1 (MQ=255)
                                                                 tttACCCATTTCTTCCCATGCAGTGTAGACATCAATCGCGCCGTCATCATAGGTCACAACGTCTTCAGc    <  1:1226468/69‑1 (MQ=255)
                                                                 tttACCCATTTCTTCCCATGCAGTGTAGACATCAATCGCGCCGTCATCATAGGTCACAACGTCTTCAGc    <  1:633933/69‑1 (MQ=255)
                                                                 tttACCCATTTCTTCCCATGCAGTGTAGACATCAATCGCGCCGTCATCATAGGTCACAACGTCTTCAGc    <  1:2686211/69‑1 (MQ=255)
                                                                 tttACCCATTTCTTCCCATGCAGTGTAGACATCAATCGCGCCGTCATCATAGGTCACAACGTCTTCAGc    <  1:1307350/69‑1 (MQ=255)
                                                                 tttACCCATTTCTTCCCATGCAGTGTAGACATCAATCGCGCCGTCATCATAGGTCACAACGTCTTCAGc    <  1:2743985/69‑1 (MQ=255)
                                                                 tttACCCATTTCTTCCCATGCAGTGTAGACATCAATCGCGCCGTCATCATAGGTCACAACGTCTTCAGc    <  1:355078/69‑1 (MQ=255)
                                                                 tttACCCATTTCTTCCCATGCAGTGTAGACATCAATCGAGCCGTCATCATAGGTCACAACGTCTTCAGc    <  1:533275/69‑1 (MQ=255)
                                                                  ttACCCATTTCTTCCCATGCAGTGTAGACATCAATCGCGCCGTCATCATAGGTCACAACGCCTTCAGc    <  1:934182/68‑1 (MQ=255)
                                                                   tACCCATTTCTTCCCATGCAGTGTAGACATCAATCGCGCCGTCATCATAGGTCACAACGTCTTCAGCAc  >  1:2048217/1‑69 (MQ=255)
                                                                    |                                                                   
GACGGGATCATGGAAACTTCCGCGCTGTCTGCTTTCAGACCTGCCGCTTCCAGAGCGTCGCGCACTTTACCCATTTCTTCCCATGCAGTGTAGACATCAATCGCGCCGTCATCATAGGTCACAACGTCTTCAGCAC  >  minE/1251460‑1251595

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: