Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 385,156 T→G 31.1% A122A (GCT→GCG cstA → carbon starvation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE385,1560TG31.1% 50.0 / 53.1 61A122A (GCT→GCGcstAcarbon starvation protein
Reads supporting (aligned to +/- strand):  ref base T (16/26);  new base G (6/13);  total (22/39)
Fisher's exact test for biased strand distribution p-value = 7.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGTCG  >  minE/385086‑385224
                                                                      |                                                                    
cAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCt                                                                      <  1:1548247/71‑1 (MQ=255)
      cGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTAGCTGGggt                                                                 <  1:4945212/70‑1 (MQ=255)
       gCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGggtgg                                                               <  1:2738382/71‑1 (MQ=255)
        cTGGTGGGTCCGGTTCTCGCCGCGCAGATGGGCTACCTGCCTGGCACGCTGTGGCTGCTGGCGGGggtcgt                                                              <  1:3437747/71‑1 (MQ=255)
           gtggGTCCGGTTCTCGCCGCGCAGATGGGCTACCTGCCTGGCACGCTGTGGCTGCTGGCGGGGGTCGTgct                                                           >  1:3439234/1‑71 (MQ=255)
           gtggGTCCGGTTCTCGCCGCGCAGATGGGCTACCTGCCTGGCACGCTGTGGCTGCTGGCGGGGGTCGTgct                                                           >  1:449019/1‑71 (MQ=255)
              ggtccGGTTCTCGCCGCGCAGATGGGCTACCTGCCTGGCACGCTGTGGCTGCTGGCGGGGGTCGTGCTGGc                                                        >  1:3274463/4‑71 (MQ=255)
               ggCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGcc                                                       <  1:5507184/71‑1 (MQ=255)
                     tACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGggtgg                                                               <  1:5656408/57‑1 (MQ=255)
                     tACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGCCTGCTCGCTGGggtgg                                                               <  1:3690331/57‑1 (MQ=255)
                      aCTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTgg                                                                    >  1:4122070/1‑51 (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:1569900/65‑1 (MQ=255)
                          gccgcgcAGATGGGCTACCTGCCTGGCACGCTGTGGCTGCTGGCGGGGGTCGTGCTGGCcggtgcggt                                               <  1:12417/65‑1 (MQ=255)
                             gcgcAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAg                                           <  1:2642605/69‑1 (MQ=255)
                              cgcAGATGGGCTACCTGCCTGGCACGCTGTGGCTGCTGGCGGGGGTCGTGCTGGCCGGTGCGGTTCAGGac                                        <  1:3892571/71‑2 (MQ=255)
                              cgcAGATGGGCTACCTGCCTGGCACGCTGTGGCTGCTGGCGGGGGTCGTGCTGGCCGGTGCGGTTCAGGac                                        <  1:4336456/71‑2 (MQ=255)
                               gcAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGAtt                                       >  1:5187398/1‑71 (MQ=255)
                                cagaTGGGCTACCTGCCTGGCACGCTGTGGCTGCTGGCGGGGGTCGTGCTGGCCGGTGCGGTTCAGGActt                                      <  1:4055455/68‑1 (MQ=255)
                                cAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTggg                                                                   <  1:1156919/42‑1 (MQ=255)
                                cAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTggg                                                                   <  1:4314084/42‑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)
                                  aaTGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGtggtgg                                                               <  1:4952236/44‑1 (MQ=255)
                                  aaTGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGggtgg                                                               <  1:890215/44‑1 (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:4303072/40‑1 (MQ=255)
                                                  gggATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCcggtg                                                   <  1:247791/40‑1 (MQ=255)
                                                  gggATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCcggtg                                                   <  1:4647609/40‑1 (MQ=255)
                                                   ggATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCt                   <  1:416204/71‑1 (MQ=255)
                                                   ggATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCt                   <  1:3221073/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)
                                                                      |                                                                    
CAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGTCG  >  minE/385086‑385224

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: