Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 385,075 C→T 12.3% A95A (GCC→GCT cstA → carbon starvation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE385,0750CT12.3% 149.7 / 10.7 57A95A (GCC→GCTcstAcarbon starvation protein
Reads supporting (aligned to +/- strand):  ref base C (18/32);  new base T (2/5);  total (20/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.72e-01

CGCCATAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCT  >  minE/385007‑385145
                                                                    |                                                                      
cgcCATAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTTCTGTTCGGTCACCATTTTGCGGCCa                                                                       <  1:1465084/70‑1 (MQ=255)
 gcCATAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCAtt                                                                     <  1:41723/71‑1 (MQ=255)
 gcCATAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCAtt                                                                     <  1:1943463/71‑1 (MQ=255)
 gcCATAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCAtt                                                                     <  1:3732035/71‑1 (MQ=255)
  ccATAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTg                                                                    <  1:3269649/71‑1 (MQ=255)
  ccATAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTg                                                                    >  1:2884745/1‑71 (MQ=255)
    aTAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGcc                                                                  <  1:3678626/71‑1 (MQ=255)
        cgacgGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAg                                                              >  1:2086705/1‑71 (MQ=255)
         gacgGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGc                                                             <  1:3376753/71‑1 (MQ=255)
          aggggCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGGCACCATTTTGCGGCCATTGCCgg                                                                <  1:3757942/65‑1 (MQ=255)
           cgGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGc                                                                   <  1:1576451/63‑1 (MQ=255)
           cgGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCa                                                            <  1:3265336/70‑1 (MQ=255)
               cttgACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGTCCATTGCCGGAGCAGGTcc                                                       <  1:1972803/68‑1 (MQ=255)
                  gACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGc                                                     <  1:3340742/70‑1 (MQ=255)
                   aCTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCtg                                                   >  1:1799774/1‑71 (MQ=255)
                   aCTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCtg                                                   >  1:3081190/1‑71 (MQ=255)
                    cTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCtgg                                                  <  1:329394/71‑1 (MQ=255)
                    cTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCtgg                                                  <  1:784927/71‑1 (MQ=255)
                     tATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCtggt                                                 >  1:3048428/1‑71 (MQ=255)
                       tgtgCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCtggtgg                                               <  1:3058147/71‑1 (MQ=255)
                       tgtgCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCtggt                                                 <  1:3521552/69‑1 (MQ=255)
                        gtgCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGcc                                                                  <  1:2067321/51‑1 (MQ=255)
                        gtgCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCg                                                      <  1:435790/63‑1 (MQ=255)
                         tgCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTgggg                                             <  1:3075079/71‑1 (MQ=255)
                         tgCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGc                                            <  1:525590/71‑1 (MQ=255)
                          gCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCAt                                                                      >  1:128646/1‑45 (MQ=255)
                           ccGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGc                                                                   <  1:700735/47‑1 (MQ=255)
                            cGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTg                                                                    <  1:3448280/45‑1 (MQ=255)
                              aCGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAg                                                              >  1:666472/1‑49 (MQ=255)
                              aCGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCgg                                         >  1:2809642/1‑70 (MQ=255)
                                ggACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTAc                                      >  1:1463868/1‑71 (MQ=255)
                                   cAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCg                                                                 <  1:2910785/41‑1 (MQ=255)
                                   cAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTgg                                   <  1:29184/71‑1 (MQ=255)
                                    aaGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCt                                                    >  1:2987707/1‑53 (MQ=255)
                                    aaGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCt                                                    >  1:2453953/1‑53 (MQ=255)
                                    aaGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTggc                                  >  1:1271111/1‑71 (MQ=255)
                                    aaGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTgg                                   >  1:3370545/1‑70 (MQ=255)
                                        acgtgtTGTTTGGTCACCACTTCGCCGCTATCGCCGGTGCTGGTCCGCTGGTGGGTCCGGTTCTCgccgcg                              >  1:2778902/3‑71 (MQ=255)
                                        aaGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTggcggc                               <  1:952885/70‑1 (MQ=255)
                                         aGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTAc                                      <  1:1889097/62‑1 (MQ=255)
                                             ttGTTTGGTCACCACTTCGCCGCTATCGCCGGTGCTGGTCCGCTGGTGGGTCCGGTTCTCGCCGCGCAGAt                         <  1:2429710/70‑1 (MQ=255)
                                             cTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCaaa                          >  1:1152441/1‑70 (MQ=255)
                                                ttCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATggg                      >  1:491594/1‑71 (MQ=255)
                                                ttCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATg                        <  1:2390870/69‑1 (MQ=255)
                                                    gTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTAc                  <  1:788121/71‑1 (MQ=255)
                                                      caccacTTCGCCGCTATCGCCGGTGCTGGTCCGCTGGTGGGTCCGGTTCtcg                                   >  1:553342/1‑50 (MQ=255)
                                                       aCCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTg               >  1:633270/1‑71 (MQ=255)
                                                       aCCACTTCGCCGCTATCGCCGGTGCTGGTCCGCTGGTGGGTCCGGTTCTCGCCGCGCAGATGGGCTACCTg               <  1:3351358/71‑1 (MQ=255)
                                                        ccATTTTGCGGCCATTGCCCGAGCAGGACCGCTGGTGGGGCCgct                                        >  1:3282226/1‑43 (MQ=255)
                                                        ccACTTCGCCGCTATCGCCGGTGCTGGTTCGCTGGTGGGTCCGGTTCTCGCCGCGCAGATGGGCTACCTGc              <  1:3578533/71‑1 (MQ=255)
                                                        ccACTTCGCCGCTATCGCCGGTGCTGGTCCGCTGGTGGGTCCGGTTCTCGCCGCGCAGATGGGCTACCTGc              <  1:2743310/71‑1 (MQ=255)
                                                        ccACTTCGCCGCTATCGCCGGTGCTGGTCCGCTGGTGGGTCCGGTTCTCGCCGCGCAGATAGGCTACCTGc              <  1:1676088/71‑1 (MQ=255)
                                                            tttGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCggg          >  1:2037630/1‑71 (MQ=255)
                                                             ttGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCgg           <  1:2008806/69‑1 (MQ=255)
                                                             ttGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGAGGCGCAAATGGGCTACCTGCCggg          <  1:938978/70‑1 (MQ=255)
                                                                cGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGgatg      <  1:3329337/71‑1 (MQ=255)
                                                                  gCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGgatgat    <  1:2355728/71‑1 (MQ=255)
                                                                   ccATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGcc             <  1:1011411/61‑1 (MQ=255)
                                                                    cATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATTGGCTACCTGCCGGGGATGATCt  <  1:3008884/71‑1 (MQ=255)
                                                                    cATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCt  <  1:2075286/71‑1 (MQ=255)
                                                                    cATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCt  >  1:1832261/1‑71 (MQ=255)
                                                                    cATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCt  <  1:3106573/71‑1 (MQ=255)
                                                                    |                                                                      
CGCCATAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCT  >  minE/385007‑385145

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: