Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,531,297 C→A 17.8% L118L (CTG→CTT coaD ← pantetheine‑phosphate adenylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,531,2970CA17.8% 127.0 / 14.2 61L118L (CTG→CTTcoaDpantetheine‑phosphate adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base C (40/10);  new base A (7/4);  total (47/14)
Fisher's exact test for biased strand distribution p-value = 2.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCCCTGATGGCGCGCCACCTCTTTCACCAACGATGAAGAGATAAACGACCACTCTTTCGACGGCATCAGAAACACACTTTCCAGTTCCGGCATTAAGTGGCGATTCATATGCGCCAGCTGCATTTCATA  >  minE/2531230‑2531359
                                                                   |                                                              
cGCCCTGATGGCGCGCCACCTCTTTCACCAACGATGAAGAGATAAACGACCACTCTTTCGACGGCATc                                                                >  1:1052998/1‑68 (MQ=255)
cGCCCTGATGGCGCGCCACCTCTTTCACCAACGATGAAGAGATAAACGACCACTCTTTCGACGGCATc                                                                >  1:199241/1‑68 (MQ=255)
agCCCTGATGGCGCGACACCTCTTTCACCAACGATGAAGAGATAAACGACCACTCTTTCGACGGCATc                                                                >  1:1201753/2‑68 (MQ=255)
  cccTGATGGCGCGCCACCTCTTTCACCAACGATGAAGAGATAAACGACCACTCTTTCGACGGCATAAGa                                                             >  1:48283/1‑69 (MQ=255)
  cccTGATGGCGCGCCACCTCTTTCACCAACGATGAAGAGATAAACGACCACTCTTTCGACGGCATAAGa                                                             >  1:2801780/1‑69 (MQ=255)
  cccTGATGGCGCGCCACCTCTTTCACCAACGATGAAGAGATAAACGACCACTCTTTCGACGGCATAAGa                                                             >  1:694513/1‑69 (MQ=255)
  cccTGATGGCGCGCCACCTCTTTCACCAACGATGAAGAGATAAACGACCACTCTTTCGACGGCATAAGa                                                             >  1:1729513/1‑69 (MQ=255)
  cccTGATGGCGCGCCACCTCTTTCACCAACGATGAAGAGATAAACGACCACTCTTTCGACGGCATAAGa                                                             >  1:1363619/1‑69 (MQ=255)
  cccTGATGGCGCGCCACCTCTTTCACCAACGATGAAGAGATAAACGACCACTCTTTCGACGGCATAAGa                                                             >  1:953363/1‑69 (MQ=255)
   ccTGATGGCGCGCCACCTCTTTCACCAACGATGAAGAGATAAACGACCACTCTTTCGACGGCATAAGa                                                             <  1:496882/68‑1 (MQ=255)
   ccTGATGGCGCGCCACCTCTTTCACCAACGATGAAGAGATAAACGACCACTCTTTCGACGGCATAAGa                                                             <  1:1213186/68‑1 (MQ=255)
   ccTGATGGCGCGCCACCTCTTTCACCAACGATGAAGAGATAAACGACCACTCTTTCGACGGCATAAGa                                                             <  1:3054388/68‑1 (MQ=255)
    cTGATGGCGCGCCACCTCTTTCACCAACGATGAAGAGATAAACGACCACTCTTTCGACGGCATCAGaaa                                                           >  1:2347130/1‑69 (MQ=255)
    cTGATGGCGCGCCACCTCTTTCACCAACGATGAAGAGATAAACGACCACTCTTTCGACGGCATCAGaaa                                                           >  1:2492796/1‑69 (MQ=255)
    cTGATGGCGCGCCACCTCTTTCACCAACGATGAAGAGATAAACGACCACTCTTTCGACGGCATCAGaaa                                                           >  1:409674/1‑69 (MQ=255)
    cTGATGGCGCGCCACCTCTTTCACCAACGATGAAGAGATAAACGACCACTCTTTCGACGGCATCAGaaa                                                           >  1:3133636/1‑69 (MQ=255)
    cTGATGGCGCGCCACCTCTTTCACCAACGATGAAGAGATAAACGACCACTCTTTCGACGGCATCAGaaa                                                           >  1:1542782/1‑69 (MQ=255)
    cTGATGGCGCGCCACCTCTTTCACCAACGATGAAGAGATAAACGACCACTCTTTCGACGGCATCAGaaa                                                           >  1:933501/1‑69 (MQ=255)
    cTGATGGCGCGCCACCTCTTTCACCAACGATGAAGAGATAAACGACCACTCTTTCGACGGCATCAGaaa                                                           >  1:3194405/1‑69 (MQ=255)
    cTGATGGCGCGCCACCTCTTTCACCAACGATGAAGAGATAAACGACCACTCTTTCGACGGCATCAGaaa                                                           >  1:1377312/1‑69 (MQ=255)
    cTGATGGCGCACCACCTCTTTCACCAACGATGAAGAGATAAACGACCACTCTTTCGACGGCATCAGaaa                                                           >  1:3289531/1‑69 (MQ=255)
         ggCGCGCCACCTCTTTCACCAACGATGAAGAGATAAACGCCCACTCTTTCGACGGCATCAGAAacacac                                                      >  1:3281771/1‑69 (MQ=255)
         ggCGCGCCACCTCTTTCACCAACGATGAAGAGATAAACGACCACTCTTTCGACGGCATCAGAAacacac                                                      >  1:2981326/1‑69 (MQ=255)
         ggCGCGCCACCTCTTTCACCAACGATGAAGAGATAAACGACCACTCTTTCGACGGCATCAGAAacacac                                                      >  1:2474068/1‑69 (MQ=255)
         ggCGCGCCACCTCTTTCACCAACGATGAAGAGATAAACGACCACTCTTTCGACGGCATCAGAAacacac                                                      >  1:3271532/1‑69 (MQ=255)
         ggCGCGCCACCTCTTTCACCAACGATGAAGAGATAAACGACCACTCTTTCGACGGCATCAGAAacacac                                                      >  1:1867556/1‑69 (MQ=255)
            gcgcCACCTCTTTCACCAACGATGAAGAGATAAACGACCACTCTTTCGACGGCATCAGAAACACACtt                                                    >  1:3010345/1‑68 (MQ=255)
              gcCACCTCTTTCACCAACGATGAAGAGATAAACGACCACTCTTTCGACGGCATCAGAAACACACTTTcc                                                 <  1:1380789/69‑1 (MQ=255)
               ccACCTCTTTCACCAACGATGAAGAGATAAACGACCACTCTTTCGACGGCATCAGAAACACACTTTcc                                                 <  1:1333207/68‑1 (MQ=255)
               ccACCTCTTTCACCAACGATGAAGAGATAAACGACCACTCTTTCGACGGCATCAGAAACACACTTTcc                                                 <  1:1088758/68‑1 (MQ=255)
                     ctTTCACCAACGATGAAGAGATAAACGACCACTCTTTCGACGGCATCAGAAACACACTTTcc                                                 >  1:2161818/1‑62 (MQ=255)
                     ctTTCACCAACGATGAAGAGATAAACGACCACTCTTTCGACGGCATCAGAAACACACTTTcc                                                 >  1:409788/1‑62 (MQ=255)
                     ctTTCACCAACGATGAAGAGATAAACGACCACTCTTTCGACGGCATCAGAAACACACTTTcc                                                 >  1:1397996/1‑62 (MQ=255)
                     ctTTCACCAACGATGAAGAGATAAACGACCACTCTTTCGACGGCATCAGAAACACACTTTcc                                                 >  1:3145191/1‑62 (MQ=255)
                     ctTTCACCAACGATGAAGAGATAAACGACCACTCTTTCGACGGCATCAGAAACACACTTTcc                                                 >  1:516815/1‑62 (MQ=255)
                     ctTTCACCAACGATGAAGAGATAAACGACCACTCTTTCGACGGCATCAGAAACACACTTTcc                                                 >  1:853477/1‑62 (MQ=255)
                     ctTTCACCAACGATGAAGAGATAAACGACCACTCTTTCGACGGCATCAGAAACACACTTTCCAGTTcc                                           >  1:3280778/1‑68 (MQ=255)
                        tCACCAACGATGAAGAGATAAACGACCACTCTTTCGACGGCATCAGAAACACACTTTCCAGt                                              >  1:3180892/1‑62 (MQ=255)
                        tCACCAACGATGAAGAGATAAACGACCACTCTTTCGACGGCATCAGAAACACACTTTCCAGTTCCGGc                                        >  1:630961/1‑68 (MQ=255)
                             aaCGATGAAGAGATAAACGACCACTCTTTCGACGGCATCAGAAACACACTTTCCAGTTCCGGCATTAAg                                  >  1:452059/1‑69 (MQ=255)
                              aCGATGCAGAGATAAACGACCACTCTTTCGACGGCATAAGAAACACACTTTCCAGTTCCGGCATTAAGt                                 >  1:436708/1‑69 (MQ=255)
                              aCGATGAAGAGATAAACGACCACTCTTTCGACGGCATCAGAAACACACTTTCCAGTTCCGGCATTAAg                                  >  1:2776664/1‑68 (MQ=255)
                              aCGATGAAGAGATAAACGACCACTCTTTCGACGGCATCAGAAACACACTTTCCAGTTCCGGCATTAAg                                  >  1:2138776/1‑68 (MQ=255)
                              aCGATGAAGAGATAAACGACCACTCTTTCGACGGCATCAGAAACACACTTTCCAGTTCCGGCATTAAGt                                 >  1:928910/1‑69 (MQ=255)
                              aCGATGAAGAGATAAACGACCACTCTTTCGACGGCATCAGAAACACACTTTCCAGTTCCGGCATTAAGt                                 >  1:108806/1‑69 (MQ=255)
                              aCGATGAAGAGATAAACGACCACTCTTTCGACGGCATCAGAAACACACTTTCCAGTTCCGGCATTAAGt                                 >  1:196934/1‑69 (MQ=255)
                              aCGATGAAGAGATAAACGACCACTCTTTCGACGGCATCAGAAACACACTTTCCAGTTCCGGCATTAAGt                                 >  1:1127299/1‑69 (MQ=255)
                              aCGATGAAGAGATAAACGACCACTCTTTCGACGGCATCAGAAACACACTTTCCAGTTCCGGCATTAAGt                                 >  1:1425807/1‑69 (MQ=255)
                              aCGATGAAGAGATAAACGACCACTCTTTCGACGGCATCAGAAACACACTTTCCAGTTCCGGCATTAAGt                                 >  1:1774427/1‑69 (MQ=255)
                              aCGATGAAGAGATAAACGACCACTCTTTCGACGGCATCAGAAACACACTTTCCAGTTCCGGCATTAAGt                                 >  1:1778407/1‑69 (MQ=255)
                              aCGATGAAGAGATAAACGACCACTCTTTCGACGGCATCAGAAACACACTTTCCAGTTCCGGCATTAAGt                                 >  1:2011377/1‑69 (MQ=255)
                              aCGATGAAGAGATAAACGACCACTCTTTCGACGGCATCAGAAACACACTTTCCAGTTCCGGCATTAAGt                                 >  1:2402019/1‑69 (MQ=255)
                              aCGATGAAGAGATAAACGACCACTCTTTCGACGGCATCAGAAACACACTTTCCAGTTCCGGCATTAAGt                                 >  1:2781015/1‑69 (MQ=255)
                                   gAAGAGATAAACGACCACTCTTTCGACGGCATCAGAAACACACTTTCCAGTTCCGGCATTAAGTGGCGa                            >  1:797363/1‑69 (MQ=255)
                                   gAAGAGATAAACGACCACTCTTTCGACGGCATCAGAAACACACTTTCCAGTTCCGGCATTAAGTGGCGa                            >  1:1747060/1‑69 (MQ=255)
                                    aaGAGATAAACGACCACTCTTTCGACGGCATAAGa                                                             <  1:1901570/35‑1 (MQ=255)
                                          tAAACGACCACTCTTTCGACGGCATCAGAAACACACTATCCAGTTCCGGCATTTAGTGGCGATTCata                      <  1:182108/68‑1 (MQ=255)
                                                    ctctTTCGACGGCATCAGAAACACACTTTCCAGTTCCGGCATTAAGTGGCGATTCATATGCGCCAGCTg           <  1:2291948/69‑1 (MQ=255)
                                                    ctctTTCGACGGCATCAGAAACACACTTTCCAGTTCCGGCATTAAGTGGCGATTCATATGCGCCAGCTg           <  1:2290802/69‑1 (MQ=255)
                                                    ctctTTCGACGGCATCAGAAACACACTTTCCAGTTCCGGCATTAAGTGGCGATTCATATGCGCCAGCTg           <  1:2790332/69‑1 (MQ=255)
                                                    ctctTTCGACGGCATCAGAAACACACTTTCCAGTTCCGGCATTAAGTGGCGATTCATATGCGCCAGCTg           <  1:641499/69‑1 (MQ=255)
                                                    ctctTTCGACGGCATCAGAAACACACTTTCCAGTTCCGGCATTAAGTGGCGATTCATATGCGCCAGCTg           <  1:2528861/69‑1 (MQ=255)
                                                    ctctTTCGACGGCATCAGAAACACACTTTCCAGTTCCGGCATTAAGTGGCGATTCATATGCGCCAGCTg           <  1:3398419/69‑1 (MQ=255)
                                                              ggCATCAGAAACACACTTTCCAGTTCCGGCATTAAGTGGCGATTCATATGCGCCAGCTGCATTTCata  >  1:2714911/1‑68 (MQ=255)
                                                                   |                                                              
CGCCCTGATGGCGCGCCACCTCTTTCACCAACGATGAAGAGATAAACGACCACTCTTTCGACGGCATCAGAAACACACTTTCCAGTTCCGGCATTAAGTGGCGATTCATATGCGCCAGCTGCATTTCATA  >  minE/2531230‑2531359

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: