Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 385,948 T→A 10.8% D386E (GAT→GAA cstA → carbon starvation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE385,9480TA10.8% 174.6 / 12.3 65D386E (GAT→GAAcstAcarbon starvation protein
Reads supporting (aligned to +/- strand):  ref base T (31/27);  new base A (4/3);  total (35/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

TCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGT  >  minE/385878‑386016
                                                                      |                                                                    
tCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGAt                                                                      >  1:2968835/1‑71 (MQ=255)
tCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTACGCCTGTATCATCGAt                                                                      <  1:174871/71‑1 (MQ=255)
 cGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGAt                                                                      <  1:1027440/70‑1 (MQ=255)
  ggCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATc                                                                     <  1:2534966/70‑1 (MQ=255)
  ggCTACGGCGCAATGCTGATGGAGTCCTTCGTGGCGATTATGGCGCTGGTTGCTGCGTCCATCATCGAAcc                                                                    <  1:1272316/71‑1 (MQ=255)
    cTACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCgg                                                                  >  1:1846836/1‑71 (MQ=255)
       cGGTGGGATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCg                                                               <  1:2697598/71‑1 (MQ=255)
            ggATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCgg                                                                  <  1:2117396/63‑1 (MQ=255)
              aTGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTAttt                                                        >  1:2733142/1‑71 (MQ=255)
               tGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTAtttt                                                       <  1:3409863/71‑1 (MQ=255)
               tGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTCTCCGCCTGTATCATCGATcc                                                                    >  1:743209/1‑58 (MQ=255)
                gTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCggg                                                                 <  1:1559909/60‑1 (MQ=255)
                gTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTAttt                                                        >  1:3677136/1‑69 (MQ=255)
                   aaTGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTAttt                                                        >  1:617591/1‑66 (MQ=255)
                   aaTGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTAttt                                                        <  1:1664047/66‑1 (MQ=255)
                   aaTGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGcc                                                    >  1:406668/1‑70 (MQ=255)
                    aTGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCAt                                                  >  1:2331529/1‑71 (MQ=255)
                     tGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGc                                                                >  1:2330920/1‑56 (MQ=255)
                     tGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATg                                                 <  1:3649642/71‑1 (MQ=255)
                     tGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATg                                                 <  1:1408146/71‑1 (MQ=255)
                         gtCCTTCGTGGCGATTATGGCGCTGGTTGCTGCGTCCATCATCGAACCGGGTCTTTACTTCGCGATGAACa                                             <  1:2450192/70‑1 (MQ=255)
                         aTCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATg                                                 >  1:1329927/1‑67 (MQ=255)
                         aTCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATGa                                                >  1:1494920/1‑68 (MQ=255)
                         aTCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATGAACa                                             <  1:2018197/71‑1 (MQ=255)
                          tCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGc                                                                >  1:3828925/1‑51 (MQ=255)
                          tCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTg                                                      <  1:2107254/61‑1 (MQ=255)
                            cTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATGAACAGcc                                          <  1:1425091/71‑1 (MQ=255)
                            cTTCGTGGCGATTATGGCGCTGGTTGCTGCGTCCATCATCGAACCGGGTCTTTACTTCGCGATGAACac                                            >  1:2351686/1‑68 (MQ=255)
                              tCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCgt                                                              >  1:3333672/1‑49 (MQ=255)
                              tCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATGAACAGccc                                         >  1:1677542/1‑70 (MQ=255)
                               cGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATGAACAGccc                                         <  1:3012109/69‑1 (MQ=255)
                               cGTGGCGATTATGGCGCTGGTTGCTGCGTCCATCATCGAACCGGGTCTTTACTTCGCGATGAACACCCCgc                                       >  1:3528104/1‑70 (MQ=255)
                                gTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTAtt                                                         <  1:3694265/52‑1 (MQ=255)
                                gTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATg                                                 <  1:3506746/60‑1 (MQ=255)
                                       ttATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGt                                >  1:2256305/1‑70 (MQ=255)
                                       ttATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGt                                >  1:3021694/1‑70 (MQ=255)
                                       ttATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTg                               <  1:14396/71‑1 (MQ=255)
                                       ttATGGCGCTGGTTGCTGCGTCCATCATCGAACCGGGTCTTTACTTCGCGATGaa                                               >  1:1355696/1‑55 (MQ=255)
                                       ttATGGCGCTGGTTGCTGCGTCCATCATCGAACCGGGTCTTTACTTCGCGATGaa                                               >  1:3722328/1‑55 (MQ=255)
                                        tATGGCGCTGGTTGCTGCGTCCATCATCGAACCGGGTCTTTACTTCGCGATGAAc                                              <  1:2112680/55‑1 (MQ=255)
                                         aTGGCGCTGGTTTCCGCCTGTATCATCGATCCggg                                                                 >  1:3884738/1‑35 (MQ=255)
                                         aTGGCGCTGGTTTCCGCCTGTATCATCGATCCggg                                                                 >  1:808085/1‑35 (MQ=255)
                                         aTGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATGa                                                <  1:1431758/52‑1 (MQ=255)
                                          tGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGt                                <  1:2764697/67‑1 (MQ=255)
                                             cgcTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATGAGCAGCCCGATGGCGGTGCTGGCt                         <  1:724819/71‑1 (MQ=255)
                                             cgcTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCt                         <  1:449226/71‑1 (MQ=255)
                                              gcTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTc                        >  1:1405036/1‑71 (MQ=255)
                                               cTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATg                                                 <  1:110534/45‑1 (MQ=255)
                                               cTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCAAGAACAGcc                                          >  1:2128327/1‑52 (MQ=255)
                                                  gTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCgg                     <  1:3141410/70‑1 (MQ=255)
                                                     tCCGCCTGTATCATCGATCCGGGCGTGTATTTTGcc                                                    >  1:1259441/1‑36 (MQ=255)
                                                      ccGCCTGTATCATCGATCCGGGCGTGTATTTTGcc                                                    >  1:375767/1‑35 (MQ=255)
                                                      ccGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAgg                 >  1:2988646/1‑70 (MQ=255)
                                                       cGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAggg                <  1:1597973/70‑1 (MQ=255)
                                                       cGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGa               >  1:254284/1‑71 (MQ=255)
                                                       cGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGa               >  1:3278000/1‑71 (MQ=255)
                                                         ccTGTATCATCGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGAcg             <  1:285236/71‑1 (MQ=255)
                                                            gTATCATCGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGAcggcg          >  1:3698815/1‑71 (MQ=255)
                                                             tATCATCGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGc                          <  1:2039867/54‑1 (MQ=255)
                                                             tATCATCGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGc                          <  1:1562283/54‑1 (MQ=255)
                                                             tATCATCGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGAcggcg          <  1:2980094/70‑1 (MQ=255)
                                                                catcGATCCGTGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGAtg      >  1:2472677/1‑71 (MQ=255)
                                                                 atcGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCgg                     >  1:2122624/1‑55 (MQ=255)
                                                                   cGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTgg   >  1:3672266/1‑71 (MQ=255)
                                                                   cGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTgg   >  1:1207560/1‑71 (MQ=255)
                                                                    gATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGc                    >  1:2911890/1‑53 (MQ=255)
                                                                    gATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCa                   >  1:547902/1‑54 (MQ=255)
                                                                    gATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTgg   <  1:3684144/70‑1 (MQ=255)
                                                                      tCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGt  <  1:3016322/69‑1 (MQ=255)
                                                                      |                                                                    
TCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGT  >  minE/385878‑386016

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: