Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 385,168 C→G 31.1% L126L (CTC→CTG cstA → carbon starvation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE385,1680CG31.1% 60.6 / 38.2 52L126L (CTC→CTGcstAcarbon starvation protein
Reads supporting (aligned to +/- strand):  ref base C (13/22);  new base G (4/12);  total (18/34)
Fisher's exact test for biased strand distribution p-value = 5.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGTCGCTCGCTGGGTG  >  minE/385100‑385235
                                                                    |                                                                   
ggtccGGTTCTCGCCGCGCAGATGGGCTACCTGCCTGGCACGCTGTGGCTGCTGGCGGGGGTCGTGCTGGc                                                                   >  1:3274463/4‑71 (MQ=255)
 ggCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGcc                                                                  <  1:5507184/71‑1 (MQ=255)
        aCTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCcggtgcg                                                            >  1:4901219/1‑70 (MQ=255)
         cTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTtctc                                                                     >  1:2271265/1‑60 (MQ=255)
           ggcggcGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTtctc                                                                     >  1:3012489/1‑58 (MQ=255)
            gcggcgCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCcggtgcggtg                                                         >  1:1625793/1‑69 (MQ=255)
            gccgcgcAGATGGGCTACCTGCCTGGCACGCTGTGGCTGCTGGCGGGGGTCGTGCTGGCcggtgcggt                                                          <  1:12417/65‑1 (MQ=255)
            gccgcgcAGATGGGCTACCTGCCTGGCACGCTGTGGCTGCTGGCGGGGGTCGTGCTGGCcggtgcggt                                                          <  1:1569900/65‑1 (MQ=255)
               gcgcAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAg                                                      <  1:2642605/69‑1 (MQ=255)
                cgcAGATGGGCTACCTGCCTGGCACGCTGTGGCTGCTGGCGGGGGTCGTGCTGGCCGGTGCGGTTCAGGac                                                   <  1:4336456/71‑2 (MQ=255)
                cgcAGATGGGCTACCTGCCTGGCACGCTGTGGCTGCTGGCGGGGGTCGTGCTGGCCGGTGCGGTTCAGGac                                                   <  1:3892571/71‑2 (MQ=255)
                 gcAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGAtt                                                  >  1:5187398/1‑71 (MQ=255)
                  cagaTGGGCTACCTGCCTGGCACGCTGTGGCTGCTGGCGGGGGTCGTGCTGGCCGGTGCGGTTCAGGActt                                                 <  1:4055455/68‑1 (MQ=255)
                  cAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGAttt                                                 <  1:1112480/71‑1 (MQ=255)
                    gaTGGGCTACCTGCCTGGCACGCTGTGGCTGCTGGCGGGGGTCGTGCTGGCCGGTGCGGTTCAGGActtta                                               >  1:3683789/2‑69 (MQ=255)
                    gaTGGGCTACCTGCCTGGCACGCTGTGGCTGCTGGCGGGGGTCGTGCTGGCCGGTGCGGTTCAGGActtta                                               <  1:5612495/70‑3 (MQ=255)
                    gaTGGGCTACCTGCCTGGCACGCTGTGGCTGCTGGCGGGGGTCGTGCTGGCCGGTGCGGTTCAGGActtta                                               <  1:1917284/70‑3 (MQ=255)
                      tGGGCTACCTGCCTGGCACGCTGTGGCTGCTGGCGGGGGTCGTGCTGGCCGGTGCGGTTCAGGACTTTATg                                             <  1:4973445/71‑1 (MQ=255)
                      tGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTAGCCGGTGCGGTGCAGGATTTCATg                                             >  1:5038970/1‑71 (MQ=255)
                           tACCTGCCTGGCACGCTGTGGCTGCTGGCGGGGGTCGTGCTGGCCGGTGCGGTTCAGGACTTTATGGTGCt                                        >  1:268300/1‑71 (MQ=255)
                           tACCTGCCTGGCACGCTGTGGCTGCTGGCGGGGGTCGTGCTGGCCGGTGCGGTTCAGGACTTTATGGTGCt                                        >  1:2429487/1‑71 (MQ=255)
                             ccTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCtgt                                      <  1:3102020/71‑1 (MQ=255)
                              cTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCtgtt                                     >  1:4242997/1‑71 (MQ=255)
                               tGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCcg                                                                 >  1:5497278/1‑42 (MQ=255)
                               tGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCcg                                                                 >  1:3477614/1‑42 (MQ=255)
                                 ccGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGAttt                                                 <  1:808785/56‑1 (MQ=255)
                                 ccGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCtgtttgt                                  >  1:4970713/1‑71 (MQ=255)
                                   ggggATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCa                                               >  1:890358/1‑56 (MQ=255)
                                   ggggATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCa                                               >  1:3114377/1‑56 (MQ=255)
                                    gggATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCcggtg                                                              <  1:247791/40‑1 (MQ=255)
                                    gggATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCcggtg                                                              <  1:4647609/40‑1 (MQ=255)
                                    gggATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCcggtg                                                              <  1:4303072/40‑1 (MQ=255)
                                     ggATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCt                              <  1:3221073/71‑1 (MQ=255)
                                     ggATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCt                              <  1:416204/71‑1 (MQ=255)
                                      gatgatCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCt                              <  1:356398/70‑1 (MQ=255)
                                         gatCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTAcgc                          >  1:4146158/1‑71 (MQ=255)
                                          atCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGAttt                                                 <  1:3860205/47‑1 (MQ=255)
                                          atCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGAttt                                                 <  1:3721332/47‑1 (MQ=255)
                                             tGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGTTGCAGGAtt                                                  <  1:3479728/43‑1 (MQ=255)
                                               gctgctGGCGGGGGTCGTGCTGGCCGGTGCGGTTCAGGACTTTATGGTGCTGTTTATCTCCTCTCGCCGt                     <  1:1542498/70‑1 (MQ=255)
                                               gctgctGGCGGGGGTCGTGCTGGCCGGTGCGGTTCAGGACTTTATGGTGCTGTTTATCTCCTCTCGCCGt                     <  1:3090677/70‑1 (MQ=255)
                                               gctgctGGCGGGGGGCGTGCTGGCCGGTGCGGTTCAGGACTTTATGGTGCTGTTTATCTCCTCTCGCCGt                     <  1:500199/70‑1 (MQ=255)
                                               gctgctCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGcc                       <  1:4305117/68‑1 (MQ=255)
                                               gctgctCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGt                     <  1:3200192/70‑1 (MQ=255)
                                               gctgctCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGt                     <  1:4418362/70‑1 (MQ=255)
                                               gctgctCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGt                     <  1:2458152/70‑1 (MQ=255)
                                                ctgctgGCGGGGGTCGTGCTGGCCGGTGCGGTTCAGGACTTTATGGTGCtgttt                                    <  1:1922410/54‑1 (MQ=255)
                                                     cGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGtc              >  1:5490728/1‑71 (MQ=255)
                                                      gCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGtcg             >  1:3046235/1‑71 (MQ=255)
                                                         ggggTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGt               <  1:5310234/66‑1 (MQ=255)
                                                          gggTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGtcgctc          <  1:1595512/70‑1 (MQ=255)
                                                           ggtggtTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGtcgctcgc        >  1:3698662/1‑71 (MQ=255)
                                                            gtggtTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGTcgctcgct       <  1:1056214/71‑1 (MQ=255)
                                                               gtTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTa                             <  1:1259021/46‑1 (MQ=255)
                                                                 tctcGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGTCGCTCGCTGGGt   >  1:4562062/1‑70 (MQ=255)
                                                                 tctcGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGTCGCTCGCTGGGt   >  1:1874767/1‑70 (MQ=255)
                                                                 tctcGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGTCGCTCGCTGGGt   >  1:4414246/1‑70 (MQ=255)
                                                                  ctcGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTa                             >  1:5014995/1‑43 (MQ=255)
                                                                  ctcGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGTCGCTCGCTGGGt   >  1:2785190/1‑69 (MQ=255)
                                                                  ctcGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGTCGCTCGCTGGGTg  <  1:614695/70‑1 (MQ=255)
                                                                   tcGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTc                               <  1:365446/40‑1 (MQ=255)
                                                                    |                                                                   
GGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGTCGCTCGCTGGGTG  >  minE/385100‑385235

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: