Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 385,180 G→T 23.2% V130V (GTG→GTT cstA → carbon starvation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE385,1800GT23.2% 124.2 / 22.8 56V130V (GTG→GTTcstAcarbon starvation protein
Reads supporting (aligned to +/- strand):  ref base G (20/23);  new base T (3/10);  total (23/33)
Fisher's exact test for biased strand distribution p-value = 2.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

GCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGTCGCTCGCTGGGTGAGCTGGTCAAAGA  >  minE/385112‑385248
                                                                    |                                                                    
gcggcgCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCcggtgcggtg                                                                      >  1:1625793/1‑69 (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:5612495/70‑3 (MQ=255)
        gaTGGGCTACCTGCCTGGCACGCTGTGGCTGCTGGCGGGGGTCGTGCTGGCCGGTGCGGTTCAGGActtta                                                            >  1:3683789/2‑69 (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)
                     ccGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGAttt                                                              <  1:808785/56‑1 (MQ=255)
                     ccGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCtgtttgt                                               >  1:4970713/1‑71 (MQ=255)
                       ggggATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCa                                                            >  1:3114377/1‑56 (MQ=255)
                       ggggATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCa                                                            >  1:890358/1‑56 (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:2458152/70‑1 (MQ=255)
                                   gctgctCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGt                                  <  1:4418362/70‑1 (MQ=255)
                                   gctgctCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGt                                  <  1:3200192/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:1874767/1‑70 (MQ=255)
                                                     tctcGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGTCGCTCGCTGGGt                >  1:4414246/1‑70 (MQ=255)
                                                     tctcGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGTCGCTCGCTGGGt                >  1:4562062/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)
                                                         gCCGGTGCGGTGCAGGATTTCATGGTGCtgtttgtt                                              >  1:372742/1‑36 (MQ=255)
                                                         gCCGGTGCGGTGCAGGATTTCATGGTGCTgtttgttt                                             >  1:2593068/1‑37 (MQ=255)
                                                         gCCGGTGCGGTGCAGGATTTCATGGTGCTgtttgttt                                             >  1:400031/1‑37 (MQ=255)
                                                         gCCGGTGCGGTGCAGGATTTCATGGTGCTgtttgttt                                             >  1:2460449/1‑37 (MQ=255)
                                                         gCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTc                                            <  1:5469044/38‑1 (MQ=255)
                                                         gCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGTCGCTCGCTGGGTGAGc            <  1:2879691/70‑1 (MQ=255)
                                                         gCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGTCGCTCGCTGGGTGAGCt           >  1:2012251/1‑71 (MQ=255)
                                                          ccGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGTCGCTCGCTGGGTGAGCTg          <  1:5132589/71‑1 (MQ=255)
                                                            ggtgcggtgcAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGTCGCTCGCTGGGTGAGCTGGt        >  1:3763399/1‑71 (MQ=255)
                                                             gtgcggtgcAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGTCGCTCGCTGGGTGAGCTGGTc       <  1:1909713/71‑1 (MQ=255)
                                                              tgcggtgcAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGACgg                             <  1:2536175/48‑1 (MQ=255)
                                                                 ggtgcAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGTCGCTCGCTGGGTGAGCTGGTCAaag   <  1:1055205/71‑1 (MQ=255)
                                                                  gtgcAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGACGCTCGCTCGCTGGGTGAGCTGGTCAaaga  <  1:1892611/71‑1 (MQ=255)
                                                                    |                                                                    
GCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGTCGCTCGCTGGGTGAGCTGGTCAAAGA  >  minE/385112‑385248

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: