Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 385,140 T→C 17.9% M117T (ATG→ACG)  cstA → carbon starvation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE385,1400TC17.9% 144.1 / 25.3 67M117T (ATG→ACG) cstAcarbon starvation protein
Reads supporting (aligned to +/- strand):  ref base T (21/34);  new base C (5/7);  total (26/41)
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.01e-01

CGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTAC  >  minE/385071‑385209
                                                                     |                                                                     
cGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGgatg                                                                      <  1:3329337/71‑1 (MQ=255)
  gCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGgatgat                                                                    <  1:2355728/71‑1 (MQ=255)
    cATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATTGGCTACCTGCCGGGGATGATCt                                                                  <  1:3008884/71‑1 (MQ=255)
    cATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCt                                                                  <  1:3106573/71‑1 (MQ=255)
    cATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCt                                                                  <  1:2075286/71‑1 (MQ=255)
    cATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCt                                                                  >  1:1832261/1‑71 (MQ=255)
            gAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGgctgct                                                           <  1:3736903/70‑1 (MQ=255)
                 ggTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTgg                                                     <  1:2535290/71‑1 (MQ=255)
                  gTCCGCTGGTGGGTCCGGTTCTCGCCGCGCAGATGGGCTACCTGCCTGGCACGCTGTGGCTGCTGGcgg                                                      <  1:2764584/69‑3 (MQ=255)
                  gTCCGCTGGTGGGTCCGGTTCTCGCCGCGCAGATGGGCTACCTGCCTGGCACGCTGTGGCTGCTGGCgggg                                                    >  1:3133188/1‑71 (MQ=255)
                  gTCCGCTGGTGGGTCCGGTTCTCGCCGCGCAGATGGGCTACCTGCCTGGCACGCTGTGGCTGCTGGCgggg                                                    >  1:696918/1‑71 (MQ=255)
                  gTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGc                                                        <  1:450810/67‑1 (MQ=255)
                      gCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGggtgg                                                >  1:2151925/1‑71 (MQ=255)
                       cTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTg                                                      <  1:2541821/64‑1 (MQ=255)
                        tggtggGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGtt                                              >  1:1105524/1‑71 (MQ=255)
                             gggCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGgctgc                                                            <  1:532101/52‑1 (MQ=255)
                             gggCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGc                                         <  1:3899209/71‑1 (MQ=255)
                                ccGGTTCTCGCCGCGCAGATGGGCTACCTGCCTGGCACGCTGTGGCTGCTGGCGGGGGTCGTGCTGGCcgg                                      >  1:1007658/1‑71 (MQ=255)
                                ccGGTTCTCGCCGCGCAGATGGGCTACCTGCCTGGCACGCTGTGGCTGCTGGCGGGGGTCGTGCTGGCcgg                                      >  1:642165/1‑71 (MQ=255)
                                ccGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCcgg                                      <  1:887719/71‑1 (MQ=255)
                                ccGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCcgg                                      <  1:1456049/71‑1 (MQ=255)
                                 cGGTTCTCGCCGCGCAGATGGGCTACCTGCCTGGCACGCTGTGGCTGCTGGCGGGGGTCGTGCTGGCcgg                                      <  1:1106663/70‑1 (MQ=255)
                                   gTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGgct                                                              <  1:3381243/44‑1 (MQ=255)
                                     aCTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTggg                                                    <  1:3049909/52‑1 (MQ=255)
                                     aCTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTggg                                                    <  1:3414683/52‑1 (MQ=255)
                                     aCTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCcggtgcgg                                 >  1:3593323/1‑71 (MQ=255)
                                      cTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGgct                                                              >  1:2079447/1‑41 (MQ=255)
                                      cTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCcggtgcggt                                <  1:2969541/71‑1 (MQ=255)
                                      cTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCcggtgcggt                                >  1:51912/1‑71 (MQ=255)
                                        ggcggcGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCcggtgcggtg                               >  1:2952552/1‑70 (MQ=255)
                                         gcggcgCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCcggtgcggtg                               <  1:920041/69‑1 (MQ=255)
                                         gcggcgCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCcggtgcggtg                               >  1:3180320/1‑69 (MQ=255)
                                          ccgcgcAGATGGGCTACCTGCCTGGCACGCTGTGGCTGCTGGCGGGGGTCGTGCTGGCCGGTGCGGTTCAg                            <  1:2199469/69‑1 (MQ=255)
                                          ccgcgcAGATGGGCTACCTGCCTGGCACGCTGTGGCTGCTGGCGGGGGTCGTGCTGGCCGGTGCGGTTCAg                            <  1:1723176/69‑1 (MQ=255)
                                            gcgcAGATGGGCTACCTGCCTGGCACGCTGTGGCTGCTGGCGGGGGTCGTGCTGGCCGGTGCGGTTCAGGa                          <  1:31942/71‑1 (MQ=255)
                                             cgcAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGa                          <  1:353451/70‑1 (MQ=255)
                                              gcAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTCGTTCTCGCCGGTGCGGTGCAg                            <  1:2034086/67‑1 (MQ=255)
                                               cAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGAttt                       >  1:3917260/1‑71 (MQ=255)
                                                aaaTGGGCTACCTGCCGGGGATGATCTGGGTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGAttt                       >  1:1502356/1‑70 (MQ=255)
                                                aaaTGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCggtgcggtgc                              <  1:1083019/63‑1 (MQ=255)
                                                aaaTGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAg                            <  1:1766098/65‑1 (MQ=255)
                                                aaaTGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGAttt                       >  1:2676202/1‑70 (MQ=255)
                                                 gaTGGGCTACCTGCCTGGCACGCTGTGGCTGCTGGCGGGGGTCGTGCTGGCCGGTGCGGTTCAGGActtta                     <  1:5334/70‑3 (MQ=255)
                                                 aaTGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTtctc                                           >  1:3525463/1‑49 (MQ=255)
                                                  aTGGGCTACCTGCCTGGCACGCTGTGGCTGCTGGCGGGGGTCGTGCTGGCcggtgcggtt                               <  1:3770706/60‑2 (MQ=255)
                                                  aTGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGggtggt                                               <  1:3262511/44‑1 (MQ=255)
                                                  aTGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGggtggt                                               <  1:245569/44‑1 (MQ=255)
                                                  aTGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCAt                    >  1:2791340/1‑71 (MQ=255)
                                                  aTGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCAt                    >  1:221712/1‑71 (MQ=255)
                                                   tGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATg                   >  1:2352343/1‑71 (MQ=255)
                                                      gCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTg                >  1:21204/1‑71 (MQ=255)
                                                       cTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGAtt                        >  1:1403815/1‑62 (MQ=255)
                                                        tACCTGCCTGGCACGCTGTGGCTGCTGGCGGGGGTCGTGCTGGCCGGTGCGGTTCAGGACTTTATGGTga               >  1:1253979/1‑69 (MQ=255)
                                                          ccTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTc                      <  1:2928980/61‑1 (MQ=255)
                                                          ccTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCtgt            >  1:1932713/1‑71 (MQ=255)
                                                          ccTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCt              <  1:2920578/69‑1 (MQ=255)
                                                          ccTGCCGGGGATGATATGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCtgt            >  1:5547/1‑71 (MQ=255)
                                                          ccTGCAGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGAt                         >  1:1525154/1‑58 (MQ=255)
                                                             gCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCtgttt          <  1:1753348/70‑1 (MQ=255)
                                                             gCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCtgttt          <  1:3784738/70‑1 (MQ=255)
                                                              gcGGGGATGATGTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGAttt                       <  1:1342292/55‑1 (MQ=255)
                                                              ccGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCtgtttgt        <  1:2643513/71‑1 (MQ=255)
                                                              ccGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCtgtttgt        >  1:2787022/1‑71 (MQ=255)
                                                              ccGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCtgttt          <  1:2164562/69‑1 (MQ=255)
                                                                ggggATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCcggtgcggtg                               <  1:3786238/46‑1 (MQ=255)
                                                                ggggATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTgtttgttt      <  1:3854765/71‑1 (MQ=255)
                                                                 gggATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTgtttgttt      <  1:3092309/70‑1 (MQ=255)
                                                                  ggATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCcggtgc                                   <  1:2088594/40‑1 (MQ=255)
                                                                  ggATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCcggtgc                                   <  1:1541613/40‑1 (MQ=255)
                                                                    atgatCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTAc  <  1:61701/71‑1 (MQ=255)
                                                                     |                                                                     
CGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTAC  >  minE/385071‑385209

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: