Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 385,069 T→C 16.4% F93F (TTT→TTC cstA → carbon starvation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE385,0690TC16.4% 177.7 / 24.5 79F93F (TTT→TTCcstAcarbon 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 = 1.00e+00

GCGGTGCGCCATAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGG  >  minE/385001‑385137
                                                                    |                                                                    
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: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:967651/67‑1 (MQ=255)
           tAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCAtt                                                             <  1:5153350/67‑1 (MQ=255)
           tAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCAtt                                                             <  1:5657472/67‑1 (MQ=255)
           tAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCAtt                                                             <  1:5663818/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:781360/70‑1 (MQ=255)
             acgacgGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCgg                                                        <  1:1060206/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)
                         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: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: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:273080/65‑1 (MQ=255)
                                                              ccACTTCGCCGCTATCGCCGGTGCTGGTCCGCTGGTGGGTCCGGTTCTCGCCGCGCAGATGGGCt            <  1:4808523/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)
                                                                    |                                                                    
GCGGTGCGCCATAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGG  >  minE/385001‑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: