Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 385,847 A→C 10.6% M353L (ATG→CTG)  cstA → carbon starvation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE385,8470AC10.6% 147.5 / 14.5 66M353L (ATG→CTG) cstAcarbon starvation protein
Reads supporting (aligned to +/- strand):  ref base A (30/29);  new base C (2/5);  total (32/34)
Fisher's exact test for biased strand distribution p-value = 4.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.79e-01

ATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGA  >  minE/385781‑385916
                                                                  |                                                                     
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)
    ccATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGc                                                                    <  1:3153357/66‑1 (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:2278416/71‑1 (MQ=255)
          ccTGTGGTGCGGTATCTGGCTTCCACGCGCTGATCTCTTCCGGTACGACGCCAAAACTGCTGGctaatgaa                                                         <  1:2030088/71‑1 (MQ=255)
           cTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAg                                                        >  1:2297925/1‑71 (MQ=255)
           cTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAg                                                        <  1:1613947/71‑1 (MQ=255)
              tgGTGCGGTATCTGGCTTCCACGCGCTGATCTCTTCCGGTACGACGCCAAAACTGCTGGctaatgaaacc                                                      >  1:3047714/1‑67 (MQ=255)
              tgGTGCGGTATCTGGCTTCCACGCGCTGATCTCTTCCGGTACGACGCCAAAACTGCTGGctaatgaaac                                                       <  1:3812158/69‑3 (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:2899356/1‑71 (MQ=255)
                        tCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCtt                                           >  1:2903705/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)
                                    gcgcTGATCTCTTCCGGTACGACGCCAAAACTGCTGGCTAATGAAACCGACGCGCGTTTCATCGGCTACgg                               >  1:3022534/1‑71 (MQ=255)
                                      gcTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTg                             <  1:3509131/71‑1 (MQ=255)
                                       cTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTgg                            <  1:1662811/71‑1 (MQ=255)
                                       cTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTgg                            >  1:3016323/1‑71 (MQ=255)
                                       cTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTg                             <  1:1583949/70‑1 (MQ=255)
                                        tGATCTCTTCTGGGACCACTCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCttt                                          <  1:2659821/56‑1 (MQ=255)
                                        tGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCttt                                          <  1:1582497/56‑1 (MQ=255)
                                        tGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTgg                            >  1:374929/1‑70 (MQ=255)
                                         gATCTCTTCTGGGACCACGCCGAAGATGCTGGcgaacga                                                          <  1:88656/39‑1 (MQ=255)
                                         gATCTCTTCTGGGACCACGCCGAAGATGCTGGcgaacga                                                          <  1:1591502/39‑1 (MQ=255)
                                          aTCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTggg                           <  1:3076262/69‑1 (MQ=255)
                                           tctcTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTg                             >  1:696242/1‑66 (MQ=255)
                                               ttctGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTaa                    >  1:3139172/1‑71 (MQ=255)
                                                tctGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGt                       >  1:515839/1‑67 (MQ=255)
                                                 ctGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGc                                             >  1:3580757/1‑44 (MQ=255)
                                                 ctGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATg                  >  1:2695806/1‑71 (MQ=255)
                                                   gggACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATgg                 >  1:1923530/1‑70 (MQ=255)
                                                    ggACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATgg                 <  1:2062311/69‑1 (MQ=255)
                                                    ggACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGa                <  1:3520558/70‑1 (MQ=255)
                                                        cACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCt           <  1:3416980/71‑1 (MQ=255)
                                                         aCGCCGAAGATGGTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGa                          >  1:1521606/1‑55 (MQ=255)
                                                          cGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTc         >  1:424096/1‑71 (MQ=255)
                                                          cGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTc         >  1:2964384/1‑71 (MQ=255)
                                                            ccGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGa                          <  1:46536/52‑1 (MQ=255)
                                                            ccGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCg        <  1:3240639/70‑1 (MQ=255)
                                                            ccGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCGt       >  1:3810750/1‑71 (MQ=255)
                                                              gAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATgg                 <  1:1053513/59‑1 (MQ=255)
                                                                aGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGc    >  1:961064/1‑70 (MQ=255)
                                                                aGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCg   <  1:1370659/71‑1 (MQ=255)
                                                                 gATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGa  >  1:2975520/1‑71 (MQ=255)
                                                                  aTGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCtt          <  1:2360372/62‑1 (MQ=255)
                                                                  |                                                                     
ATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGA  >  minE/385781‑385916

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: