Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 385,084 A→T 15.3% G98G (GGA→GGT cstA → carbon starvation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE385,0840AT15.3% 129.3 / 14.1 52G98G (GGA→GGTcstAcarbon starvation protein
Reads supporting (aligned to +/- strand):  ref base A (18/26);  new base T (2/6);  total (20/32)
Fisher's exact test for biased strand distribution p-value = 4.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCT  >  minE/385015‑385152
                                                                     |                                                                    
cgacgGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAg                                                                     >  1:2086705/1‑71 (MQ=255)
 gacgGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGc                                                                    <  1:3376753/71‑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: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)
                gtgCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCg                                                             <  1:435790/63‑1 (MQ=255)
                 tgCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTgggg                                                    <  1:3075079/71‑1 (MQ=255)
                 tgCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGc                                                   <  1:525590/71‑1 (MQ=255)
                      aCGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAg                                                                     >  1:666472/1‑49 (MQ=255)
                      aCGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCgg                                                >  1:2809642/1‑70 (MQ=255)
                        ggACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTAc                                             >  1:1463868/1‑71 (MQ=255)
                           cAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTgg                                          <  1:29184/71‑1 (MQ=255)
                            aaGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCt                                                           >  1:2453953/1‑53 (MQ=255)
                            aaGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCt                                                           >  1:2987707/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:3106573/71‑1 (MQ=255)
                                                            cATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCt         >  1:1832261/1‑71 (MQ=255)
                                                                 ccGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTAc                         >  1:3760422/1‑50 (MQ=255)
                                                                  cGGTGCTGGTCCGCTGGTGGGTCCGGTTCTCGCCGCGCag                                  <  1:1982117/40‑2 (MQ=12)
                                                                    gAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGcc                    <  1:961999/52‑1 (MQ=255)
                                                                    gAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGgctgct  <  1:3736903/70‑1 (MQ=255)
                                                                     |                                                                    
CGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCT  >  minE/385015‑385152

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: