Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 385,069 T→C 10.9% F93F (TTT→TTC cstA → carbon starvation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE385,0690TC10.9% 169.4 / 10.8 64F93F (TTT→TTCcstAcarbon starvation protein
Reads supporting (aligned to +/- strand):  ref base T (25/32);  new base C (2/5);  total (27/37)
Fisher's exact test for biased strand distribution p-value = 6.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.88e-01

CAGCGGTGCGCCATAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGG  >  minE/384999‑385137
                                                                      |                                                                    
cAGCGGTGCGCCATAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCAtttt                                                                      >  1:2700143/1‑71 (MQ=255)
 aGCGGTGCGCCATAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTg                                                                     >  1:1512428/1‑71 (MQ=255)
  gCGGTGCGCCATAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGc                                                                    <  1:3435647/71‑1 (MQ=255)
  gCGGTGCGCCATAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGc                                                                    <  1:2525464/71‑1 (MQ=255)
        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:3732035/71‑1 (MQ=255)
         gcCATAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCAtt                                                             <  1:41723/71‑1 (MQ=255)
         gcCATAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCAtt                                                             <  1:1943463/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)
             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:3081190/1‑71 (MQ=255)
                           aCTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCtg                                           >  1:1799774/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)
                                                                      |                                                                    
CAGCGGTGCGCCATAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGG  >  minE/384999‑385137

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: