Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 385,078 T→C 12.3% I96I (ATT→ATC cstA → carbon starvation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE385,0780TC12.3% 146.0 / 10.7 57I96I (ATT→ATCcstAcarbon starvation protein
Reads supporting (aligned to +/- strand):  ref base T (18/32);  new base C (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 = 1.00e+00

GCCATAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCT  >  minE/385008‑385145
                                                                      |                                                                   
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)
       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: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)
                       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)
                          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)
                                                                      |                                                                   
GCCATAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCT  >  minE/385008‑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: