Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 385,066 T→C 16.4% H92H (CAT→CAC cstA → carbon starvation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE385,0660TC16.4% 173.6 / 25.5 79H92H (CAT→CACcstAcarbon starvation protein
Reads supporting (aligned to +/- strand):  ref base T (19/47);  new base C (5/8);  total (24/55)
Fisher's exact test for biased strand distribution p-value = 5.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

GCCAGCGGTGCGCCATAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGG  >  minE/384997‑385136
                                                                     |                                                                      
gCCAGCGGTGCGCCATAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCAtt                                                                       <  1:506824/71‑1 (MQ=255)
    gCGGTGCGCCATAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGc                                                                   <  1:1271483/71‑1 (MQ=255)
       ctGCGCCATAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGc                                                                <  1:565044/70‑1 (MQ=255)
          cgcCATAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCa                                                              <  1:2759067/70‑1 (MQ=255)
          cgcCATAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCAt                                                             >  1:4927011/1‑71 (MQ=255)
           gcCATAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCAtt                                                            <  1:2971558/71‑1 (MQ=255)
            ccATAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGc                                                                <  1:4333570/66‑1 (MQ=255)
            ccATAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGc                                                                <  1:3791880/66‑1 (MQ=255)
             cATAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGc                                                          <  1:525246/71‑1 (MQ=255)
             cATAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGc                                                          <  1:3267131/71‑1 (MQ=255)
             cATAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGc                                                          <  1:1894336/71‑1 (MQ=255)
             cATAACGACGGGCTGGACTATGTGCCGAAGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGc                                                          <  1:4955664/71‑1 (MQ=255)
              aTAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGcc                                                         <  1:1730019/71‑1 (MQ=255)
               tAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCAtt                                                            <  1:5663818/67‑1 (MQ=255)
               tAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCAtt                                                            <  1:5657472/67‑1 (MQ=255)
               tAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCAtt                                                            <  1:5153350/67‑1 (MQ=255)
               tAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCAtt                                                            <  1:967651/67‑1 (MQ=255)
               tAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCg                                                        <  1:2051506/71‑1 (MQ=255)
                 acgacgGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCAtt                                                                       >  1:1304091/1‑54 (MQ=255)
                 acgacgGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCAtt                                                                       >  1:2071143/1‑54 (MQ=255)
                 acgacgGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCa                                                              >  1:4797923/1‑63 (MQ=255)
                 acgacgGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCggc                                                      >  1:1702531/1‑70 (MQ=255)
                 acgacgGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCgg                                                       <  1:1060206/70‑1 (MQ=255)
                 acgacgGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCgg                                                       <  1:781360/70‑1 (MQ=255)
                 acgacgGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATATTGCGGCCATTGCCg                                                        >  1:370272/1‑69 (MQ=255)
                  cgacgGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAg                                                     >  1:3004488/1‑71 (MQ=255)
                  cgacgGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAg                                                     <  1:2937430/71‑1 (MQ=255)
                  cgacgGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAg                                                     <  1:2894153/71‑1 (MQ=255)
                  cgacgGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAg                                                     <  1:1969869/71‑1 (MQ=255)
                  cgacgGACTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAg                                                     >  1:739857/1‑71 (MQ=255)
                   gacgGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCAtttt                                                                     >  1:5138426/1‑54 (MQ=255)
                   gacgGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGc                                                    <  1:2977286/71‑1 (MQ=255)
                    acgGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCa                                                              >  1:2245817/1‑60 (MQ=255)
                    acgGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAg                                                     <  1:1323330/69‑1 (MQ=255)
                     cgGGCTGGACTATGTGCCGACGGACAGGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAg                                                  <  1:744356/71‑1 (MQ=255)
                     cgGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAg                                                  <  1:5164516/71‑1 (MQ=255)
                     cgGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAg                                                  <  1:5272378/71‑1 (MQ=255)
                      gggCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAgg                                                 <  1:19354/71‑1 (MQ=255)
                      gggCTGGACTATGTGCAGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAgg                                                 <  1:2840582/71‑1 (MQ=255)
                         cTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTcc                                              <  1:3776807/71‑1 (MQ=255)
                          tGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCAt                                                                        >  1:1040396/1‑44 (MQ=255)
                             aCTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGa                                                      <  1:3179177/59‑1 (MQ=255)
                             aCTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCtg                                          <  1:3968233/71‑1 (MQ=255)
                               tATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCtggt                                        >  1:1513531/1‑71 (MQ=255)
                                 tgtgCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTg                                                                    >  1:3319642/1‑41 (MQ=255)
                                 tgtgCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCtggtg                                       <  1:2850864/70‑1 (MQ=255)
                                   tgCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGcc                                                               <  1:3224693/44‑1 (MQ=255)
                                   tgCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGcc                                                               <  1:2117544/44‑1 (MQ=255)
                                   tgCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGcc                                                               <  1:1883475/44‑1 (MQ=255)
                                     ccGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGcc                                  >  1:489017/1‑71 (MQ=255)
                                      cGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCg                                 >  1:4572003/1‑71 (MQ=255)
                                       gACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCgg                                <  1:554019/71‑1 (MQ=255)
                                       gACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCgg                                <  1:1661392/71‑1 (MQ=255)
                                        aCGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCa                                                              >  1:5182327/1‑40 (MQ=255)
                                         ctgACAAGAAAGTGCTGTTCGGGCACCATTTTGCGGCCATTGCCGGAGCAGGGCCGCTGGTGGGGcc                                  <  1:3712220/65‑1 (MQ=255)
                                          ggACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCtg                                          >  1:2584290/1‑58 (MQ=255)
                                          ggACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTAc                             >  1:1046367/1‑71 (MQ=255)
                                             cAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTgg                          <  1:2712877/71‑1 (MQ=255)
                                              aaGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTggc                         >  1:3787673/1‑71 (MQ=255)
                                                    gtgtTGTTTGGTCACCACTTCGCCGCTATCGCCGGTGCTGGTCCGCtggt                                        <  1:1850922/50‑1 (MQ=255)
                                                    gtgtTGTTTGGTCACCACTTCGCCGCTATCGCCGGTGCTGGTCCGCtggt                                        <  1:2651637/50‑1 (MQ=255)
                                                    gTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCa                   <  1:1289279/71‑1 (MQ=255)
                                                    gTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCa                   <  1:1296041/71‑1 (MQ=255)
                                                     tGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCg                                             <  1:4700197/44‑1 (MQ=255)
                                                     tGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCg                                             <  1:33002/44‑1 (MQ=255)
                                                      gttGTTTGGTCACCACTTCGCCGCTATCGCCGGTGCTGGTCCGCTGGTGGGTCCGGTTCTCgccgcg                     <  1:1649463/65‑1 (MQ=255)
                                                       ttGTTTGGTCACCACTTCGCCGCTATCGCCGGTGCTGGTCCGCtggt                                        >  1:4289989/2‑47 (MQ=255)
                                                       ttGTTTGGTCACCACTTCGCCGCTATCGCCGGTGCTGGTCCGCtggt                                        >  1:2677279/2‑47 (MQ=255)
                                                       cTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAAt                >  1:868934/1‑71 (MQ=255)
                                                        tGTTTGGTCACCACTTCGCCGCTATCGCCGGTGCTGGTCCGCTGGTGGGTCCGGTTCTCGCCGCGCa                   >  1:3884430/1‑67 (MQ=255)
                                                        tGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAAt                <  1:5489501/70‑1 (MQ=255)
                                                         gTTTGGTCACCACTTCGCCGCTATCGCCGGTGCTGGTCCGCTGGTGGGTCCGGTTCTCGCCGCGCAGATg               <  1:5612341/70‑1 (MQ=255)
                                                         gTTTGGTCACCACTTCGCCGCTATCGCCGGTGCTGGTCCGCTGGTGGGTCCGGTTCTCGCCGCGCAGATg               <  1:5279289/70‑1 (MQ=255)
                                                           ttggTCACCACTTCGCCGCTATCGCCGGTGCTGGTCCGCTGGTGGGTCCGGTTCTCGCCGCGCAGATGGGc            >  1:230749/3‑71 (MQ=255)
                                                            cGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGc            >  1:2273881/1‑70 (MQ=255)
                                                             ggTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCt           <  1:1631204/70‑1 (MQ=255)
                                                              gTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGc            >  1:204562/1‑68 (MQ=255)
                                                                cACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGt                               <  1:1850941/47‑1 (MQ=255)
                                                                cACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGt                               <  1:625748/47‑1 (MQ=255)
                                                                 aCCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTg      >  1:4949266/1‑71 (MQ=255)
                                                                 aCCACTTCGCCGCTATCGCCGGTGCTGGTCCGCTGGTGGGTCCGGTTCTCGCCGCGCAGATGGGCTACCt       <  1:3629173/70‑1 (MQ=255)
                                                                  ccACTTCGCCGCTATCGCCGGTGCTGGTCCGCTGGTGGGTCCGGTTCTCGCCGCGCAGATGGGCt           <  1:4808523/65‑1 (MQ=255)
                                                                  ccACTTCGCCGCTATCGCCGGTGCTGGTCCGCTGGTGGGTCCGGTTCTCGCCGCGCAGATGGGCt           <  1:273080/65‑1 (MQ=255)
                                                                  ccACTTCGCCGCTATCGCCGGTGCTGGTCCGCTGGTGGGTCCGGTTCTCGCCGCGCAGATGGGCTACCTGc     >  1:3509648/1‑71 (MQ=255)
                                                                     ttttGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCt           <  1:647690/62‑1 (MQ=255)
                                                                     ttttGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCgg  >  1:5385301/1‑71 (MQ=255)
                                                                     ttttGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCgg  >  1:373704/1‑71 (MQ=255)
                                                                     |                                                                      
GCCAGCGGTGCGCCATAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGG  >  minE/384997‑385136

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: