Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 385,087 A→T 16.7% A99A (GCA→GCT cstA → carbon starvation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE385,0870AT16.7% 117.5 / 16.5 48A99A (GCA→GCTcstAcarbon starvation protein
Reads supporting (aligned to +/- strand):  ref base A (16/24);  new base T (2/6);  total (18/30)
Fisher's exact test for biased strand distribution p-value = 6.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.12e-01

CGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCT  >  minE/385018‑385152
                                                                     |                                                                 
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:784927/71‑1 (MQ=255)
         cTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCtgg                                                         <  1:329394/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)
                   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:1832261/1‑71 (MQ=255)
                                                         cATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCt         <  1:3106573/71‑1 (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)
                                                                     |                                                                 
CGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCT  >  minE/385018‑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: