Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 385,072 G→C 15.8% A94A (GCG→GCC cstA → carbon starvation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE385,0720GC15.8% 195.6 / 25.3 82A94A (GCG→GCCcstAcarbon starvation protein
Reads supporting (aligned to +/- strand):  ref base G (22/47);  new base C (5/8);  total (27/55)
Fisher's exact test for biased strand distribution p-value = 7.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GTGCGCCATAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGG  >  minE/385004‑385137
                                                                    |                                                                 
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:3791880/66‑1 (MQ=255)
     ccATAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGc                                                                 <  1:4333570/66‑1 (MQ=255)
      cATAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGc                                                           <  1:3267131/71‑1 (MQ=255)
      cATAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGc                                                           <  1:525246/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:5153350/67‑1 (MQ=255)
        tAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCAtt                                                             <  1:5657472/67‑1 (MQ=255)
        tAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCAtt                                                             <  1:5663818/67‑1 (MQ=255)
        tAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCAtt                                                             <  1:967651/67‑1 (MQ=255)
        tAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCg                                                         <  1:2051506/71‑1 (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)
            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)
                      aCTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGa                                                       <  1:3179177/59‑1 (MQ=255)
                      aCTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCtg                                           <  1:3968233/71‑1 (MQ=255)
                        tATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCtggt                                         >  1:1513531/1‑71 (MQ=255)
                          tgtgCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCtggtg                                        <  1:2850864/70‑1 (MQ=255)
                            tgCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGcc                                                                <  1:3224693/44‑1 (MQ=255)
                            tgCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGcc                                                                <  1:1883475/44‑1 (MQ=255)
                            tgCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGcc                                                                <  1:2117544/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:2651637/50‑1 (MQ=255)
                                             gtgtTGTTTGGTCACCACTTCGCCGCTATCGCCGGTGCTGGTCCGCtggt                                         <  1:1850922/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:5279289/70‑1 (MQ=255)
                                                  gTTTGGTCACCACTTCGCCGCTATCGCCGGTGCTGGTCCGCTGGTGGGTCCGGTTCTCGCCGCGCAGATg                <  1:5612341/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:625748/47‑1 (MQ=255)
                                                         cACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGt                                <  1:1850941/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)
                                                                ttGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCggg  >  1:283420/1‑70 (MQ=255)
                                                                ttGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCggg  >  1:2082745/1‑70 (MQ=255)
                                                                    |                                                                 
GTGCGCCATAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGG  >  minE/385004‑385137

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: