Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 385,974 G→C 11.3% S395T (AGC→ACC)  cstA → carbon starvation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE385,9740GC11.3% 114.4 / 6.2 44S395T (AGC→ACC) cstAcarbon starvation protein
Reads supporting (aligned to +/- strand):  ref base G (12/27);  new base C (3/2);  total (15/29)
Fisher's exact test for biased strand distribution p-value = 3.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.58e-01

CTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGCTTCTGCCGCGCAGGTGGTGAGTAG  >  minE/385906‑386043
                                                                    |                                                                     
cTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATGAACAGcc                                                                     >  1:144419/1‑71 (MQ=255)
cTTCGTGGCGATTATGGCGCTGGTTGCTGCGTCCATCATCGAACCGGGTCTTTACTTCGCGATGAACAccc                                                                     >  1:5577617/1‑71 (MQ=255)
  tCGTGGCGATTATGGCGCTGGTTGCTGCGTCCATCATCGAACCGGGTCTTTACTTCGCGATGAACACCCCg                                                                   >  1:1968933/1‑71 (MQ=255)
   cGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATGAACAGccc                                                                    >  1:4800594/1‑69 (MQ=255)
   cGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATGAACAGccc                                                                    <  1:5523451/69‑1 (MQ=255)
   cGTGGCGATTATGGCGCTGGTTGCTGCGTCCATCATCGAACCGGGTCTTTACTTCGCGATGAACACCCCgc                                                                  >  1:2396515/1‑70 (MQ=255)
   cGTGGCGATTATGGCGCTGGTTGCTGCGTCCATCATCGAACCGGGTCTTTACTTCGCGATGAACACCCCgc                                                                  >  1:1954198/1‑70 (MQ=255)
   cGTGGCGATTACGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATGAACAGccc                                                                    <  1:1744653/69‑1 (MQ=255)
       gCGATTATGGCGCTGGTTGCTGCGTCCATCATCGAACCGGGTCTTTACTTCGCGATGAACACCCCgc                                                                  <  1:1177082/67‑2 (MQ=255)
       gCGATTATGGCGCTGGTTGCTGCGTCCATCATCGAACCGGGTCTTTACTTCGCGATGAACACCCCgc                                                                  <  1:1170488/67‑2 (MQ=255)
          aTTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGt                                                           <  1:5356536/71‑1 (MQ=255)
            tATGGCGCTGGTTGCTGCGTCCATCATCGAACCGGGTCTTTACTTCGCGATGAACAccc                                                                     >  1:2232338/1‑59 (MQ=255)
            tATGGCGCTGGTTGCTGCGTCCATCATCGAACCGGGTCTTTACTTCGCGATGAACAccc                                                                     >  1:122627/1‑59 (MQ=255)
              tGGCGCTGGTTTCCGCCTGTATCATCGAGCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCt                                                        <  1:2370925/70‑1 (MQ=255)
                   cTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATg                                                                >  1:4033427/1‑57 (MQ=255)
                     ggTTTCCGCCTGTATCATCGATCCTGGCGTGTATTTTGCCATGAACAGccc                                                                    <  1:46340/51‑1 (MQ=255)
                     ggTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATGAACAGccc                                                                    <  1:1721911/51‑1 (MQ=255)
                             ccTGTATCATCGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGAcg                                        >  1:3031944/1‑71 (MQ=255)
                                 tATCATCGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATg                                                                <  1:764686/43‑1 (MQ=255)
                                  atcatcGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATg                                                                <  1:3170218/42‑1 (MQ=255)
                                   tcatcGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGAt                                  >  1:3087739/1‑71 (MQ=255)
                                          tCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCg                           <  1:3520844/71‑1 (MQ=255)
                                          tCCGGGCGTGTATTTTACCATTAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCg                           <  1:1715959/71‑1 (MQ=255)
                                           ccGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGc                          >  1:5260507/1‑71 (MQ=255)
                                            cGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGCt                         <  1:4650254/71‑1 (MQ=255)
                                              ggCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGCTTc                       <  1:336755/71‑1 (MQ=255)
                                               gCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGCTTCt                      <  1:3744493/71‑1 (MQ=255)
                                                cGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGCTTCTg                     <  1:3751628/71‑1 (MQ=255)
                                                cGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGCTTCTg                     <  1:48219/71‑1 (MQ=255)
                                                    tATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGCTTCTGCcgc                 >  1:3079273/1‑71 (MQ=255)
                                                    tATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGCTTCTGCcg                  >  1:963691/1‑70 (MQ=255)
                                                      ttttGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGTCGGATGTGGTCGCTTCTGCcgcg                <  1:1938472/70‑1 (MQ=255)
                                                      ttttGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGc                          <  1:2033383/60‑1 (MQ=255)
                                                      ttttGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGCTTCTGCcgcg                <  1:365381/70‑1 (MQ=255)
                                                      ttttGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGCTTCTGCcgcg                >  1:19775/1‑70 (MQ=255)
                                                      ttttGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGCTTCTGCcgcg                <  1:4090815/70‑1 (MQ=255)
                                                      ttttGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGCTTCTGCCgcgc               <  1:5174180/71‑1 (MQ=255)
                                                        ttGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGCTTCTGCCGCGCAg             >  1:3913300/1‑71 (MQ=255)
                                                           ccATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGCTTCTGCCGCGCAggtg          <  1:3772223/71‑1 (MQ=255)
                                                           ccATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGCTTCTGCCGCGCAggt           <  1:213990/70‑1 (MQ=255)
                                                           ccATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGCTTCTGCCGCGCAggt           <  1:3610177/70‑1 (MQ=255)
                                                             aTGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGCTTCTGCCGCGCAggtggt        <  1:478261/71‑1 (MQ=255)
                                                             aTGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGCTTCTGCCGCGCAggtggt        <  1:3162243/71‑1 (MQ=255)
                                                              tGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGCTTCTGCCGCGCAggtggt        >  1:2724767/1‑70 (MQ=255)
                                                              tGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGCTTCTGCCGCGCAGgtggtg       <  1:216487/71‑1 (MQ=255)
                                                                   aGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGCTTCTGCCGCGCAGGTGGTGAGTAg  >  1:4632892/1‑71 (MQ=255)
                                                                   aGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGCTTCTGCCGCGCAGGTGGTGAGTAg  <  1:3763759/71‑1 (MQ=255)
                                                                    gCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGCTTCTGCCGCGCa              <  1:2846716/58‑1 (MQ=255)
                                                                    |                                                                     
CTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGCTTCTGCCGCGCAGGTGGTGAGTAG  >  minE/385906‑386043

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: