Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 385,807 T→G 6.8% S339S (TCT→TCG cstA → carbon starvation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE385,8070TG 6.8% 235.4 / 8.6 87S339S (TCT→TCGcstAcarbon starvation protein
Reads supporting (aligned to +/- strand):  ref base T (34/47);  new base G (2/4);  total (36/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCCGGTATGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCT  >  minE/385743‑385874
                                                                |                                                                   
gCCCGGTATGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCt                                                                 <  1:4652632/69‑1 (MQ=255)
gCCCGGTATGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCt                                                                 <  1:586759/69‑1 (MQ=255)
  ccGGTATGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCTGTGTCTGGCTTCCa                                                             <  1:1446556/71‑1 (MQ=255)
  ccGGTATGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCa                                                             <  1:1427593/71‑1 (MQ=255)
  ccGGTATGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCa                                                             <  1:4212524/71‑1 (MQ=255)
   cGGTATGGACCGGTAACCTGTTCCCGTTCATGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCAt                                                            >  1:4729444/1‑71 (MQ=255)
    ggTATGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATg                                                           <  1:4664225/71‑1 (MQ=255)
       aTGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATgcgc                                                        >  1:5124392/1‑71 (MQ=255)
       aTGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATgcgc                                                        >  1:4593394/1‑71 (MQ=255)
       aTGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATgcgc                                                        >  1:2293883/1‑71 (MQ=255)
       aTGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATgcgc                                                        >  1:338770/1‑71 (MQ=255)
       aTGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATgcgc                                                        >  1:3233171/1‑71 (MQ=255)
        tGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATg                                                           <  1:2172157/67‑1 (MQ=255)
        tGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATg                                                           <  1:1916734/67‑1 (MQ=255)
        tGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATg                                                           <  1:3253962/67‑1 (MQ=255)
        tGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATg                                                           <  1:3341666/67‑1 (MQ=255)
        tGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATg                                                           <  1:3573833/67‑1 (MQ=255)
        tGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATg                                                           <  1:1708200/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:2420182/1‑71 (MQ=255)
           aCCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGAt                                                    >  1:549291/1‑71 (MQ=255)
           aCCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGAt                                                    >  1:279950/1‑71 (MQ=255)
             cGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCTCGGTGTCTGGCTTCCa                                                             <  1:3458365/60‑1 (MQ=255)
              ggTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGAtctc                                                 <  1:2254526/71‑1 (MQ=255)
               gTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGAtctc                                                 <  1:647515/70‑1 (MQ=255)
               gTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATctct                                                >  1:4124511/1‑71 (MQ=255)
               gTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATctct                                                >  1:1815246/1‑71 (MQ=255)
               gTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATctct                                                >  1:4263547/1‑71 (MQ=255)
                           cGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGc                                                                  <  1:4435936/41‑1 (MQ=255)
                           cGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGc                                                                  <  1:1670049/41‑1 (MQ=255)
                           cGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGc                                                                  <  1:527102/41‑1 (MQ=255)
                            cttCCTGTTCATCACCATCGCCTGCGGCGCCGTCTCGGGCTTCCACGCGCTGATCTCGTCGGGCACCACGc                                   >  1:4431645/2‑71 (MQ=255)
                            cttCCTGTTCATCACCATCGCCTGCGGCGCCGTCTCGGGCTTCCACGCGCTGATCTCGTCGGGCACCACGc                                   >  1:2451785/2‑71 (MQ=255)
                             ttCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCt                                                                     >  1:2542423/1‑36 (MQ=255)
                             ttCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCtt                                                                >  1:1908714/1‑41 (MQ=255)
                             ttCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCtt                                                                >  1:4088896/1‑41 (MQ=255)
                             ttCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCtt                                                                >  1:891711/1‑41 (MQ=255)
                             ttCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCtt                                                                >  1:3847013/1‑41 (MQ=255)
                             ttCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATgcg                                                         <  1:4589741/48‑1 (MQ=255)
                             ttCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATgcg                                                         <  1:2896013/48‑1 (MQ=255)
                             ttCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATgcg                                                         <  1:679302/48‑1 (MQ=255)
                             ttCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATgcg                                                         <  1:4726845/48‑1 (MQ=255)
                             ttCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATgcg                                                         <  1:382769/48‑1 (MQ=255)
                               ccTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCg                                 >  1:3810118/1‑70 (MQ=255)
                               ccTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCg                                 >  1:3400471/1‑70 (MQ=255)
                               ccTGTTCATCACCATCGCCTGCGGCGCCGTCTCGGGCTTCCACGCGCTGAt                                                    <  1:3124818/51‑1 (MQ=19)
                               ccTGTTCATCACCATCGCCTGCGGCGCCGTCTCGGGCTTCCACGCGCTGAt                                                    <  1:1897151/51‑1 (MQ=19)
                               ccTGTTCATCACCATCGCCTGCGGCGCCGTCTCGGGCTTCCACGCGCTGAt                                                    <  1:5014864/51‑1 (MQ=19)
                                 tGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCg                                 >  1:4142249/1‑68 (MQ=255)
                                   tttATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTgg                                           <  1:3288923/56‑1 (MQ=255)
                                     tATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCa                                                             >  1:1367795/1‑36 (MQ=255)
                                     tATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCa                                                             >  1:1730582/1‑36 (MQ=255)
                                     tATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTgg                                           >  1:1055302/1‑54 (MQ=255)
                                     tATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGa                                         >  1:2066375/1‑56 (MQ=255)
                                     tATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGa                                         >  1:1204589/1‑56 (MQ=255)
                                     tATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGa                                         >  1:1016631/1‑56 (MQ=255)
                                        cACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTgg                       >  1:1807834/1‑71 (MQ=255)
                                         aCCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGc                      >  1:2577401/1‑71 (MQ=255)
                                         aCCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGc                      >  1:3757847/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:5012806/70‑1 (MQ=255)
                                               gCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGcgaacg                 <  1:5030359/70‑1 (MQ=255)
                                               gCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGcgaacg                 <  1:2187567/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:3509276/1‑70 (MQ=255)
                                               gCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGcgaacg                 <  1:390373/70‑1 (MQ=255)
                                               gCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGcgaacg                 <  1:3601303/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:4973737/61‑1 (MQ=255)
                                                              tCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGc           <  1:4076970/61‑1 (MQ=255)
                                                              tCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGc           <  1:914255/61‑1 (MQ=255)
                                                              tCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGc           <  1:1735831/61‑1 (MQ=255)
                                                              tCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCGGGCGAACGAAGGGc           <  1:717418/61‑1 (MQ=255)
                                                                |                                                                   
GCCCGGTATGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCT  >  minE/385743‑385874

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: