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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE385,8160TC 7.0% 252.4 / 7.4 86H342H (CAT→CACcstAcarbon starvation protein
Reads supporting (aligned to +/- strand):  ref base T (32/48);  new base C (2/4);  total (34/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.78e-01

CGGTATGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTAC  >  minE/385746‑385885
                                                                      |                                                                     
cGGTATGGACCGGTAACCTGTTCCCGTTCATGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCAt                                                                       >  1:4729444/1‑71 (MQ=255)
 ggTATGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATg                                                                      <  1:4664225/71‑1 (MQ=255)
    aTGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATgcgc                                                                   >  1:3233171/1‑71 (MQ=255)
    aTGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATgcgc                                                                   >  1:5124392/1‑71 (MQ=255)
    aTGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATgcgc                                                                   >  1:338770/1‑71 (MQ=255)
    aTGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATgcgc                                                                   >  1:4593394/1‑71 (MQ=255)
    aTGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATgcgc                                                                   >  1:2293883/1‑71 (MQ=255)
     tGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATg                                                                      <  1:1916734/67‑1 (MQ=255)
     tGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATg                                                                      <  1:3573833/67‑1 (MQ=255)
     tGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATg                                                                      <  1:2172157/67‑1 (MQ=255)
     tGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATg                                                                      <  1:1708200/67‑1 (MQ=255)
     tGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATg                                                                      <  1:3341666/67‑1 (MQ=255)
     tGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATg                                                                      <  1:3253962/67‑1 (MQ=255)
     tGGACCGGTAACATGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATg                                                                      <  1:2504617/67‑1 (MQ=255)
      ggACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATg                                                                      <  1:5199712/66‑1 (MQ=255)
       gACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATgcg                                                                    >  1:3877731/1‑67 (MQ=255)
       gACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTg                                                                 <  1:3259093/70‑1 (MQ=255)
        aCCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGa                                                                >  1:2782091/1‑70 (MQ=255)
        aCCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGAt                                                               >  1:279950/1‑71 (MQ=255)
        aCCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGAt                                                               >  1:549291/1‑71 (MQ=255)
        aCCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGAt                                                               >  1:2420182/1‑71 (MQ=255)
           ggTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGAtctc                                                            <  1:2254526/71‑1 (MQ=255)
            gTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGAtctc                                                            <  1:647515/70‑1 (MQ=255)
            gTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATctct                                                           >  1:1815246/1‑71 (MQ=255)
            gTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATctct                                                           >  1:4263547/1‑71 (MQ=255)
            gTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATctct                                                           >  1:4124511/1‑71 (MQ=255)
                         cttCCTGTTCATCACCATCGCCTGCGGCGCCGTCTCGGGCTTCCACGCGCTGATCTCGTCGGGCACCACGc                                              >  1:2451785/2‑71 (MQ=255)
                         cttCCTGTTCATCACCATCGCCTGCGGCGCCGTCTCGGGCTTCCACGCGCTGATCTCGTCGGGCACCACGc                                              >  1:4431645/2‑71 (MQ=255)
                          ttCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATgcg                                                                    <  1:4726845/48‑1 (MQ=255)
                          ttCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATgcg                                                                    <  1:382769/48‑1 (MQ=255)
                          ttCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATgcg                                                                    <  1:2896013/48‑1 (MQ=255)
                          ttCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATgcg                                                                    <  1:679302/48‑1 (MQ=255)
                          ttCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATgcg                                                                    <  1:4589741/48‑1 (MQ=255)
                            ccTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCg                                            >  1:3400471/1‑70 (MQ=255)
                            ccTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCg                                            >  1:3810118/1‑70 (MQ=255)
                            ccTGTTCATCACCATCGCCTGCGGCGCCGTCTCGGGCTTCCACGCGCTGAt                                                               <  1:3124818/51‑1 (MQ=19)
                            ccTGTTCATCACCATCGCCTGCGGCGCCGTCTCGGGCTTCCACGCGCTGAt                                                               <  1:5014864/51‑1 (MQ=19)
                            ccTGTTCATCACCATCGCCTGCGGCGCCGTCTCGGGCTTCCACGCGCTGAt                                                               <  1:1897151/51‑1 (MQ=19)
                              tGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCg                                            >  1:4142249/1‑68 (MQ=255)
                                tttATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTgg                                                      <  1:3288923/56‑1 (MQ=255)
                                  tATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTgg                                                      >  1:1055302/1‑54 (MQ=255)
                                  tATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGa                                                    >  1:1016631/1‑56 (MQ=255)
                                  tATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGa                                                    >  1:1204589/1‑56 (MQ=255)
                                  tATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGa                                                    >  1:2066375/1‑56 (MQ=255)
                                     cACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTgg                                  >  1:1807834/1‑71 (MQ=255)
                                      aCCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGc                                 >  1:3757847/1‑71 (MQ=255)
                                      aCCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGc                                 >  1:2577401/1‑71 (MQ=255)
                                       ccATCGCCTGCGGCGCCGTCTCGGGCTTCCACGCGCTGATCTCGTCGGGCACCACGCCCAAGATGCtcg                                  <  1:754585/69‑3 (MQ=255)
                                        cATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTg                                                                 <  1:519833/37‑1 (MQ=255)
                                          tCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGcgaac                             >  1:4112289/1‑71 (MQ=255)
                                            gCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGcgaacg                            <  1:3601303/70‑1 (MQ=255)
                                            gCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGcgaacg                            <  1:3289528/70‑1 (MQ=255)
                                            gCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGcgaacg                            <  1:4240537/70‑1 (MQ=255)
                                            gCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGcgaacg                            <  1:4238302/70‑1 (MQ=255)
                                            gCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGcgaacg                            <  1:5012806/70‑1 (MQ=255)
                                            gCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGcgaacg                            <  1:390373/70‑1 (MQ=255)
                                            gCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGcgaacg                            >  1:3509276/1‑70 (MQ=255)
                                            gCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGcgaacg                            <  1:5030359/70‑1 (MQ=255)
                                            gCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGcgaacg                            <  1:2187567/70‑1 (MQ=255)
                                                       ggTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGt                >  1:1695415/1‑71 (MQ=255)
                                                          gtCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCt             >  1:1330233/1‑71 (MQ=255)
                                                           tCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGaa                                          <  1:982607/41‑1 (MQ=255)
                                                           tCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGaa                                          <  1:3567546/41‑1 (MQ=255)
                                                           tCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGaa                                          <  1:3547780/41‑1 (MQ=255)
                                                           tCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGaa                                          <  1:3415445/41‑1 (MQ=255)
                                                           tCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGaa                                          <  1:23700/41‑1 (MQ=255)
                                                           tCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGaa                                          <  1:2073813/41‑1 (MQ=255)
                                                           tCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGc                      <  1:4076970/61‑1 (MQ=255)
                                                           tCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGc                      <  1:914255/61‑1 (MQ=255)
                                                           tCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGc                      <  1:4973737/61‑1 (MQ=255)
                                                           tCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGc                      <  1:1735831/61‑1 (MQ=255)
                                                           tCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCGGGCGAACGAAGGGc                      <  1:717418/61‑1 (MQ=255)
                                                               gCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTAt         >  1:5191099/1‑70 (MQ=255)
                                                               gCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTAt         >  1:1123056/1‑70 (MQ=255)
                                                                  tCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGcgaacga                           <  1:156570/49‑1 (MQ=255)
                                                                  tCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGcgaacga                           <  1:537286/49‑1 (MQ=255)
                                                                  tCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGcgaacga                           <  1:1180339/49‑1 (MQ=255)
                                                                  tCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGcgaacga                           <  1:116797/49‑1 (MQ=255)
                                                                  tCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGcgaacga                           <  1:494122/49‑1 (MQ=255)
                                                                  tCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGcgaacga                           <  1:722011/49‑1 (MQ=255)
                                                                  tCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGcgaacga                           <  1:3974681/49‑1 (MQ=255)
                                                                   ccATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGcgaacga                           <  1:3146232/48‑1 (MQ=255)
                                                                   ccATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGc     >  1:4794289/1‑70 (MQ=255)
                                                                   ccATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGc     >  1:1099885/1‑70 (MQ=255)
                                                                   ccATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAAAGAAGGGCAGGCGTGCTTTATCGGc     >  1:647684/1‑70 (MQ=255)
                                                                     aTGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTAc  <  1:4485932/71‑1 (MQ=255)
                                                                     aTGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTAc  <  1:1527584/71‑1 (MQ=255)
                                                                      |                                                                     
CGGTATGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTAC  >  minE/385746‑385885

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: