Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,182,950 T→C 42.8% A217A (GCA→GCG mipA ← scaffolding protein for murein synthesizing machinery

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,182,9500TC42.8% 18.8 / 58.2 77A217A (GCA→GCGmipAscaffolding protein for murein synthesizing machinery
Reads supporting (aligned to +/- strand):  ref base T (24/20);  new base C (18/15);  total (42/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

CCAGGATTTATCCACCATCGGGCTGTCAGTAACTTCATCAGACAGACGGGTGTAGCGTGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAG  >  minE/1182893‑1183011
                                                         |                                                             
ccAGGATTTATCCACCATCGGGCTGTCAGTAACTTCATCAGACAGACGGGTGTAGCGTGCGGTACCGTa                                                    <  1:2629533/69‑1 (MQ=255)
ccAGGATTTATCCACCATCGGGCTGTCAGTAACTTCATCAGACAGACGGGTGTAGCGTGCGGTACCGTa                                                    <  1:2451595/69‑1 (MQ=255)
ccAGGATTTATCCACCATCGGGCTGTCAGTAACTTCATCAGACAGACGGGTGTAGCGTGCGGTACCGTa                                                    <  1:2479262/69‑1 (MQ=255)
ccAGGATTTATCCACCATCGGGCTGTCAGTAACTTCATCAGACAGACGGGTGTAGCGTGCGGTACCGTa                                                    <  1:1494909/69‑1 (MQ=255)
ccAGGATTTATCCACCATCGGGCTGTCAGTAACTTCATCAGACAGACGGGTGTAGCGTGCGGTACCGTa                                                    <  1:237952/69‑1 (MQ=255)
ccAGGATTTATCCACCATCGGGCTGTCAGTAACTTCATCAGACAGACGGGTGTAGCGTGCGGTACCGTa                                                    <  1:2248798/69‑1 (MQ=255)
ccAGGATTTATCCACCATCGGGCTGTCAGTAACTTCATCAGACAGACGGGTGTAGCGTGCGGTACCGTa                                                    <  1:1839243/69‑1 (MQ=255)
ccAGGATTTATCCACCATCGGGCTGTCAGTAACTTCATCAGACAGACGGGTGTAGCGTGCGGTACCGTa                                                    <  1:1868423/69‑1 (MQ=255)
ccAGGATTTATCCACCATCGGGCTGTCAGTAACTTCATCAGACAGACGGGTGTAGCGTGCGGTACCGTa                                                    <  1:48203/69‑1 (MQ=255)
ccAGGATTTATCCACCATCGGGCTGTCAGTAACTTCATCAGACAGACGGGTGTAGCGTGCGGTACCGTa                                                    <  1:1334206/69‑1 (MQ=255)
ccAGGATTTATCCACCATCGGGCTGTCAGTAACTTCATCAGACAGACGGGTGTAGCGTGCGGTACCGTa                                                    <  1:2953873/69‑1 (MQ=255)
ccAGGATTTATCCACCATCGGGCTGTCAGTAACTTCATCAGACAGACGGGTGTAGCGTGCGGTACCGTa                                                    <  1:2972648/69‑1 (MQ=255)
ccAGGATTTATCCACCATCGGGCTGTCAGTAACTTCATCAGACAGACGGGTGTAGCGTGCGGTACCGTa                                                    <  1:2957984/69‑1 (MQ=255)
 cAGGATTTATCCACCATCGGGCTGTCAGTAACTTCATCAGACAGACGGGTGTAGCGTGCGGTACCGTa                                                    >  1:2926686/1‑68 (MQ=255)
 cAGGATTTATCCACCATCGGGCTGTCAGTAACTTCATCAGACAGACGGGTGTAGCGTGCGGTACCGTa                                                    >  1:1153542/1‑68 (MQ=255)
 cAGGATTTATCCACCATCGGGCTGTCAGTAACTTCATCAGACAGACGGGTGTAGCGTGCGGTACCGTa                                                    >  1:1517474/1‑68 (MQ=255)
     aTTTATCCACCATTGGGCTGTCAGTAACTTCATCAGACAGACGGGTGTAGCGCGCGATACCGTAAacac                                               <  1:2794879/69‑1 (MQ=255)
     aTTTATCCACCATTGGGCTGTCAGTAACTTCATCAGACAGACGGGTGTAGCGCGCGATACCGTAAacac                                               <  1:22443/69‑1 (MQ=255)
     aTTTATCCACCATTGGGCTGTCAGTAACTTCATCAGACAGACGGGTGTAGCGCGCGATACCGTAAacac                                               <  1:2613047/69‑1 (MQ=255)
     aTTTATCCACCATTGGGCTGTCAGTAACTTCATCAGACAGACGGGTGTAGCGCGCGATACCGTAAacac                                               <  1:1475850/69‑1 (MQ=255)
     aTTTATCCACCATTGGGCTGTCAGTAACTTCATCAGACAGACGGGTGTAGCGCGCGATACCGTAAacac                                               <  1:1311182/69‑1 (MQ=255)
      tttATCCACCATCGGGCTGTCAGTAACTTCATCAGACAGACGGGTGTAGCGTGCGGTACCGTAAacac                                               >  1:1478176/1‑68 (MQ=255)
      tttATCCACCATCGGGCTGTCAGTAACTTCATCAGACAGACGGGTGTAGCGTGCGGTACCGTAAacac                                               >  1:1450061/1‑68 (MQ=255)
      tttATCCACCATCGGGCTGTCAGTAACTTCATCAGACAGACGGGTGTAGCGTGCGGTACCGTAAacac                                               >  1:2129714/1‑68 (MQ=255)
        tATCCCCCATCGGGCTGTCAGTAACTTCATCAGACAGACGGGTGTAGCGTGCGGTACCGTAAACACTcc                                            >  1:1664693/1‑69 (MQ=255)
        tATCCACCATCGGGCTGTCAGTAACTTCATCAGACAGACGGGTGTAGCGTGCGGTACCGTAAACACTcc                                            >  1:2373182/1‑69 (MQ=255)
        tATCCACCATCGGGCTGTCAGTAACTTCATCAGACAGACGGGTGTAGCGTGCGGTACCGTAAACACTcc                                            >  1:2013308/1‑69 (MQ=255)
        tATCCACCATCGGGCTGTCAGTAACTTCATCAGACAGACGGGTGTAGCGTGCGGTACCGTAAACACTc                                             >  1:2468741/1‑68 (MQ=255)
            caccaTTGGGCTGTCAGTAACTTCATCAGACAGACGGGTGTAGCGCGCGATACCGTAAacac                                               <  1:1548284/62‑1 (MQ=255)
             accaTCGGGCTGTCAGTAACTTCATCAGACAGACGGGTGTAGCGTGCGGTACCGTa                                                    <  1:2064951/56‑1 (MQ=255)
                aTCGGGCTGTCAGTAACTTCATCAGACAGACGGGTGTAGCGTGCGGTACCGTAAACACt                                              >  1:2667441/1‑59 (MQ=255)
                aTCGGGCTGTCAGTAACTTCATCAGACAGACGGGTGTAGCGTGCGGTACCGTAAACACTCCAGTCGCCg                                    >  1:279155/1‑69 (MQ=255)
                aTCGGGCTGTCAGTAACTTCATCAGACAGACGGGTGTAGCGTGCGGTACCGTAAACACTCCAGTCGCCg                                    >  1:104366/1‑69 (MQ=255)
                    ggCTGTCAGTAACTTCATCAGACAGACGGGTGTAGCGTGCGGTACCGTa                                                    <  1:1398077/49‑1 (MQ=255)
                               aCTTCATCAGACAGACGGGTGTAGCGTGCGGTACCGTaaa                                                  <  1:2897101/40‑1 (MQ=255)
                               aCTTCATCAGACAGACGGGTGTAGCGTGCGGTACCGTaaa                                                  <  1:3109197/40‑1 (MQ=255)
                               aCTTCATCAGACAGACGGGTGTAGCGTGCGGTACCGTaaa                                                  <  1:2423763/40‑1 (MQ=255)
                               aCTTCATCAGACAGACGGGTGTAGCGTGCGGTACCGTaaa                                                  <  1:2916078/40‑1 (MQ=255)
                               aCTTCATCAGACAGACGGGTGTAGCGTGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCt                      >  1:669579/1‑68 (MQ=255)
                               aCTTCATCAGACAGACGGGTGTAGCGTGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTg                     >  1:3115814/1‑69 (MQ=255)
                               aCTTCATCAGACAGACGGGTGTAGCGTGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTg                     >  1:3264806/1‑69 (MQ=255)
                               aCTTCATCAGACAGACGGGTGTAGCGTGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTg                     >  1:447427/1‑69 (MQ=255)
                               aCTTCATCAGACAGACGGGTGTAGCGTGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTg                     >  1:1316745/1‑69 (MQ=255)
                               aCTTCATCAGACAGACGGGTGTAGCGTGCGGTACCGTAAACACTCAAGTCGCCGAGGAAGTTGTAGCTg                     >  1:2010595/1‑69 (MQ=255)
                                   catcaGACAGACGGGTGTAGCGCGCGATACCGTAAACACTCCAg                                          >  1:3149525/1‑44 (MQ=255)
                                   catcaGACAGACGGGTGTAGCGCGCGATACCGTAAACACTCCAGTCACCGAGGAAg                              >  1:3324622/1‑56 (MQ=255)
                                   catcaGACAGACGGGTGTAGCGCGCGATACCGTAAACACTCCAGTCACCGAGGAAGTTGTAGCTGgcg                  >  1:543487/1‑68 (MQ=255)
                                   catcaGACAGACGGGTGTAGCGCGCGATACCGTAAACACTCCAGTCACCGAGGAAGTTGTAGCTGgcg                  >  1:860211/1‑68 (MQ=255)
                                   catcaGACAGACGGGTGTAGCGCGCGATACCGTAAACACTCCAGTCACCGAGGAAGTTGTAGCTGgcg                  >  1:322780/1‑68 (MQ=255)
                                   catcaGACAGACGGGTGTAGCGCGCGATACCGTAAACACTCCAGTCACCGAGGAAGTTGTAGCTGgcg                  >  1:3095314/1‑68 (MQ=255)
                                   catcaGACAGACGGGTGTAGCGCGCGATACCGTAAACACTCCAGTCACCGAGGAAGTTGTAGCTGgcg                  >  1:2688537/1‑68 (MQ=255)
                                   catcaGACAGACGGGTGTAGCGCGCGATACCGTAAACACTCCAGTCACCGAGGAAGTTGTAGCTGgcg                  >  1:1115731/1‑68 (MQ=255)
                                   catcaGACAGACGGGTGTAGCGCGCGATACCGTAAACACTCCAGTCACCGAGGAAGTTGTAGCTGgcg                  >  1:1271906/1‑68 (MQ=255)
                                   catcaGACAGACGGGTGTAGCGCGCGATACCGTAAACACTCCAGTCACCGAGGAAGTTGTAGCTGgcg                  >  1:1476180/1‑68 (MQ=255)
                                   catcaGACAGACGGGTGTAGCGCGCGATACCGTAAACACTCCAGTCACCGAGGAAGTTGTAGCTGgcg                  >  1:1501490/1‑68 (MQ=255)
                                   catcaGACAGACGGGTGTAGCGCGCGATACCGTAAACACTCCAGTCACCGAGGAAGTTGTAGCTGgcg                  >  1:1621398/1‑68 (MQ=255)
                                   catcaGACAGACGGGTGTAGCGCGCGATACCGTAAACACTCCAGTCACCGAGGAAGTTGTAGCTGgcg                  >  1:181820/1‑68 (MQ=255)
                                   catcaGACAGACGGGTGTAGCGCGCGATACCGTAAACACTCCAGTCACCGAGGAAGTTGTAGCTGgcg                  >  1:1821449/1‑68 (MQ=255)
                                   catcaGACAGACGGGTGTAGCGCGCGATACCGTAAACACTCCAGTCACCGAGGAAGTTGTAGCTGgcg                  >  1:2247874/1‑68 (MQ=255)
                                   catcaGACAGACGGGTGTAGCGCGCGATACCGTAAACACTCCAGTCACCGAGGAAGTTGTAGCTGgcg                  >  1:1821699/1‑68 (MQ=255)
                                   catcaGACAGACGGGTGTAGCGCGCGATACCGTAAACACTCCAGTCACAGAGGAAGTTGTAGCTGgcg                  >  1:1989129/1‑68 (MQ=255)
                                   catcaGACAGACGGGTGTAGCGCGCGATACAGTAAACACTCCAGTCACCGAGGAAGTTGTAGCTGgcg                  >  1:3047824/1‑68 (MQ=255)
                                      cagacagaCGGGTGTAGCGTGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCa              >  1:33833/1‑69 (MQ=255)
                                      cagacagaCGGGTGTAGCGTGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCa              >  1:2844313/1‑69 (MQ=255)
                                      cagacagaCGGGTGTAGCGTGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCa              >  1:638009/1‑69 (MQ=255)
                                      cagacagaCGGGTGTAGCGTGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCa              >  1:129237/1‑69 (MQ=255)
                                      cagacagaCGGGTGTAGCGTGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCa              >  1:3360034/1‑69 (MQ=255)
                                                 gtgtAGTGTGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTaa   <  1:1987286/69‑1 (MQ=255)
                                                  tgtAGCGCGCGATACCGTAAACACTCCAGTCACCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAg  <  1:3232300/69‑1 (MQ=255)
                                                  tgtAGCGCGCGATACCGTAAACACTCCAGTCACCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAg  <  1:1356811/69‑1 (MQ=255)
                                                  tgtAGCGCGCGATACCGTAAACACTCCAGTCACCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAg  <  1:1653402/69‑1 (MQ=255)
                                                  tgtAGCGCGCGATACCGTAAACACTCCAGTCACCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAg  <  1:627958/69‑1 (MQ=255)
                                                  tgtAGCGCGCGATACCGTAAACACTCCAGTCACCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAg  <  1:3208510/69‑1 (MQ=255)
                                                  tgtAGCGCGCGATACCGTAAACACTCCAGTCACCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAg  <  1:3197695/69‑1 (MQ=255)
                                                  tgtAGCGCGCGATACCGTAAACACTCCAGTCACCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAg  <  1:739508/69‑1 (MQ=255)
                                                  tgtAGCGCGCGATACCGTAAACACTCCAGTCACCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAg  <  1:1885983/69‑1 (MQ=255)
                                                  tgtAGCGCGCAATACCGTAAACACTCCAGTCACCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAg  <  1:715994/69‑1 (MQ=255)
                                                         |                                                             
CCAGGATTTATCCACCATCGGGCTGTCAGTAACTTCATCAGACAGACGGGTGTAGCGTGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAG  >  minE/1182893‑1183011

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: