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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE385,8160TC20.3% 116.8 / 25.9 59H342H (CAT→CACcstAcarbon starvation protein
Reads supporting (aligned to +/- strand):  ref base T (27/20);  new base C (3/9);  total (30/29)
Fisher's exact test for biased strand distribution p-value = 5.75e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.21e-01

GGTATGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGC  >  minE/385747‑385882
                                                                     |                                                                  
ggTATGGACCGGTAACCTGTTCCCGTTCCTGTTTATCAGCATCGCCTGTGGCGCGGTGTCTGGCTTCCATg                                                                   >  1:3508245/1‑71 (MQ=255)
    tGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCt                                                               <  1:1213685/71‑1 (MQ=255)
     ggACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCt                                                               >  1:2096827/1‑70 (MQ=255)
       aCCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGAt                                                            <  1:551934/71‑1 (MQ=255)
        ccGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCAt                                                                    >  1:3712824/1‑62 (MQ=255)
            tAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCtctt                                                       <  1:2811760/71‑1 (MQ=255)
             aaCCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCtcttc                                                      >  1:1538696/1‑71 (MQ=255)
             aaCCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCtcttc                                                      >  1:589601/1‑71 (MQ=255)
             aaCCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCtcttc                                                      >  1:1879872/1‑71 (MQ=255)
                   ttCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGAt                                                            <  1:3291672/59‑1 (MQ=255)
                   ttCCCGTTCCTGTTCATCACCATCGCCTGTGGTGCGGTATCTGGCTTCCACGCGCTGATCTCTTCCggta                                                 <  1:516376/70‑3 (MQ=255)
                    tCCCGTTCATGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGAc                                                >  1:2241936/1‑70 (MQ=255)
                      ccGTTCCTGTTCATCACCATCGCCTGTGGTGCGGTATCTGGCTTCCACGCGCTGATCTCTTCCGGTAcga                                              <  1:198973/70‑3 (MQ=255)
                      ccGTTCCTGTTCATCACCATCGCCTGTGGTGCGGTATCTGGCTTCCACGCGCTGATCTCTTCCGGTAcga                                              <  1:2624561/70‑3 (MQ=255)
                       cGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCAc                                             >  1:3610705/1‑70 (MQ=255)
                       cGTTCCTGTTGATCACCATCGCCTGTGGTGCGGTATCTGGCTTCCACGCGCTGATCTCTTCCGGTacgacg                                            >  1:3111843/1‑71 (MQ=255)
                            cTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCAGGCCGaa                                       <  1:513094/71‑1 (MQ=255)
                              gTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGAtctc                                                         >  1:2085774/1‑51 (MQ=255)
                              gTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCg                                         >  1:871837/1‑67 (MQ=255)
                              gTTCATCACCATCGCCTGTGGTGCGGTATCTGGCTTCCACGCGCTGATCTCTTCCGGTAcga                                              <  1:345865/62‑3 (MQ=255)
                               tttATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATgca                                                                 >  1:2946110/1‑41 (MQ=255)
                                  aTCACCATCGCCTGTGGTGCGGTATCTGGCTTCCACGCGCTGATCTCTTCCGGTACGACGCCAAAACTGCt                                 <  1:2748866/71‑1 (MQ=255)
                                  aTCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGc                                  <  1:3476404/70‑1 (MQ=255)
                                  aTCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCt                                 >  1:793095/1‑71 (MQ=255)
                                   tCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTg                                <  1:1647770/71‑1 (MQ=255)
                                   tCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTg                                <  1:652499/71‑1 (MQ=255)
                                     aCCATCGCCTGTGGTGCGGTATCTGGCTTCCACGCGCTGAt                                                            >  1:301806/1‑41 (MQ=255)
                                      ccATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGc                                  <  1:3153357/66‑1 (MQ=255)
                                       cATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGGTGATCTCTTCTGGGACCACGc                                           >  1:973459/1‑56 (MQ=255)
                                        aTCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGcgaa                           <  1:983120/71‑1 (MQ=255)
                                        aTCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGcgaa                           <  1:609374/71‑1 (MQ=255)
                                         tCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGcgaac                          >  1:3689463/1‑71 (MQ=255)
                                           gCCTGTGGTGCGGTATCTGGCTTCCACGCGCTGATCTCTTCCGGTACGACGCCAAAACTGCTGGctaatga                        <  1:2135801/71‑1 (MQ=255)
                                           gCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGcgaacga                        <  1:1849348/71‑1 (MQ=255)
                                            ccTGTGGTGCGGTATCTGGCTTCCACGCGCTGATCTCTTCCGGTACGACGCCAAAACTGCTGGctaatgaa                       <  1:2030088/71‑1 (MQ=255)
                                            ccTGTGGTGCGGTATCTGGCTTCCACGCGCTGATCTCTTCCGGTACGACGCCAAAACTGCTGGctaatgaa                       <  1:2278416/71‑1 (MQ=255)
                                             cTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAg                      <  1:1613947/71‑1 (MQ=255)
                                             cTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAg                      >  1:2297925/1‑71 (MQ=255)
                                                tgGTGCGGTATCTGGCTTCCACGCGCTGATCTCTTCCGGTACGACGCCAAAACTGCTGGctaatgaaacc                    >  1:3047714/1‑67 (MQ=255)
                                                tgGTGCGGTATCTGGCTTCCACGCGCTGATCTCTTCCGGTACGACGCCAAAACTGCTGGctaatgaaac                     <  1:3812158/69‑3 (MQ=255)
                                                tgGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGc                                           <  1:277114/47‑1 (MQ=255)
                                                tgGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGc                   <  1:2970180/71‑1 (MQ=255)
                                                  gcgcGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGcga                            <  1:910831/60‑1 (MQ=255)
                                                  gcgcGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAg                 <  1:3158300/71‑1 (MQ=255)
                                                       gtgtCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGt             >  1:2231017/1‑70 (MQ=255)
                                                          tCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCtt         >  1:2903705/1‑71 (MQ=255)
                                                          tCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCtt         >  1:2899356/1‑71 (MQ=255)
                                                          tCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCtt         >  1:969233/1‑71 (MQ=255)
                                                           cTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCtt         >  1:3010349/1‑70 (MQ=255)
                                                            tGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTg            <  1:1400676/66‑1 (MQ=255)
                                                            tGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTa       >  1:611994/1‑71 (MQ=255)
                                                              gCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTa       <  1:3108896/69‑1 (MQ=255)
                                                              gCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATc     >  1:2382780/1‑71 (MQ=255)
                                                               cTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTg            >  1:3858316/1‑63 (MQ=255)
                                                               cTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCg    >  1:805865/1‑71 (MQ=255)
                                                                ttCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGcg                             >  1:1392380/1‑45 (MQ=255)
                                                                ttCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCgg   >  1:1162794/1‑71 (MQ=255)
                                                                  ccATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCTAACGAAGGGCAGGCGTGCTTTAt      >  1:1428709/1‑66 (MQ=255)
                                                                  ccATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAgg                     <  1:230622/51‑1 (MQ=255)
                                                                  ccATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGc  >  1:819680/1‑70 (MQ=255)
                                                                     |                                                                  
GGTATGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGC  >  minE/385747‑385882

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: