Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 385,072 G→C 11.6% A94A (GCG→GCC cstA → carbon starvation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE385,0720GC11.6% 161.2 / 11.1 60A94A (GCG→GCCcstAcarbon starvation protein
Reads supporting (aligned to +/- strand):  ref base G (20/33);  new base C (2/5);  total (22/38)
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.61e-01

CGCCATAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATG  >  minE/385007‑385141
                                                                 |                                                                     
cgcCATAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTTCTGTTCGGTCACCATTTTGCGGCCa                                                                   <  1:1465084/70‑1 (MQ=255)
cgcCATAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCg                                                                       >  1:2037583/1‑66 (MQ=255)
cgcCATAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGc                                                                     <  1:1564063/68‑1 (MQ=255)
cgcCATAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGc                                                                     <  1:543015/68‑1 (MQ=255)
 gcCATAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCAtt                                                                 <  1:41723/71‑1 (MQ=255)
 gcCATAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCAtt                                                                 <  1:3732035/71‑1 (MQ=255)
 gcCATAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCAtt                                                                 <  1:1943463/71‑1 (MQ=255)
  ccATAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTg                                                                >  1:2884745/1‑71 (MQ=255)
  ccATAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTg                                                                <  1:3269649/71‑1 (MQ=255)
    aTAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGcc                                                              <  1:3678626/71‑1 (MQ=255)
     tAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCg                                                                       >  1:3654825/1‑61 (MQ=255)
      aaCGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGc                                                                     >  1:1783773/1‑62 (MQ=255)
      aaCGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGc                                                                     >  1:3851878/1‑62 (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)
              gCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCgg                                                                      >  1:2801818/1‑53 (MQ=255)
              gCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCg                                                                       >  1:3700973/1‑52 (MQ=255)
              gCTGGACTATGTGACGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCgg                                                                      >  1:2326731/1‑53 (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)
                                                                 |                                                                     
CGCCATAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATG  >  minE/385007‑385141

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: