Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 385,621 C→A 11.0% A277A (GCC→GCA cstA → carbon starvation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE385,6210CA11.0% 185.9 / 5.1 61A277A (GCC→GCAcstAcarbon starvation protein
Reads supporting (aligned to +/- strand):  ref base C (24/30);  new base A (3/4);  total (27/34)
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

CAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCATTTTGAT  >  minE/385553‑385689
                                                                    |                                                                    
cAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGccc                                                                     <  1:1764602/70‑1 (MQ=255)
 aGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCcg                                                                   <  1:3206180/71‑1 (MQ=255)
   ttGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGATGGCCCcccg                                                                 >  1:852744/1‑71 (MQ=255)
     gACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGt                                                                >  1:1315308/1‑70 (MQ=255)
     gACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGt                                                                >  1:5671959/1‑70 (MQ=255)
       ccTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGAc                                                             >  1:5090822/1‑71 (MQ=255)
          ggATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTAc                                                          >  1:4665770/1‑71 (MQ=255)
          ggATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCACCGCGTg                                                               >  1:5331024/1‑66 (MQ=255)
           gATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGcccc                                                                    >  1:4819517/1‑60 (MQ=255)
           gATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTAcc                                                         >  1:1069388/1‑71 (MQ=255)
           gATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTAcc                                                         >  1:4752566/1‑71 (MQ=255)
           gATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTAc                                                          <  1:4483431/70‑1 (MQ=255)
            aTGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACct                                                        >  1:4723507/1‑71 (MQ=255)
                 ggtggGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTa                                                   >  1:144673/1‑71 (MQ=255)
                 ggtggGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTa                                                   >  1:4706999/1‑71 (MQ=255)
                  gtggGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTAc                                                  >  1:2173621/1‑71 (MQ=255)
                      gTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTc                                              <  1:3594466/71‑1 (MQ=255)
                        tACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCt                                            >  1:556715/1‑71 (MQ=255)
                        tACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCt                                            >  1:4062551/1‑71 (MQ=255)
                         aCGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTAcc                                                         >  1:4165621/1‑57 (MQ=255)
                           ggTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGaa                                         <  1:833257/71‑1 (MQ=255)
                            gTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGaaa                                        <  1:1130967/71‑1 (MQ=255)
                               ttGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTcc                                             <  1:2931034/63‑1 (MQ=255)
                               ttGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATc                                     >  1:3968050/1‑71 (MQ=255)
                                       gcggTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATc                                     <  1:2092252/63‑1 (MQ=255)
                                       gcggTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATc                                     <  1:5671157/63‑1 (MQ=255)
                                       gcggTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATc                                     <  1:3265248/63‑1 (MQ=255)
                                       gcggTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATc                                     <  1:5542898/63‑1 (MQ=255)
                                        cactgctgcCAGTGTGGCTGATCCTCGCACCGCGCGACTATCTGGCAACCTTCCTGaaaa                                       >  1:907676/4‑60 (MQ=255)
                                           tgctgcCAGTGTGGCTGATCCTCGCACCGCGCGACTATCTGGCAACCTTCCTGa                                          <  1:3358470/54‑1 (MQ=255)
                                            gctgcCAGTGTGGCTGATCCTCGCACCGCGCGACTATCTGGCAACCTTCCTGAAAATCGGCGTTAtcgtcg                        >  1:3116366/1‑71 (MQ=255)
                                             ctgcCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTa                              >  1:1879722/1‑64 (MQ=255)
                                             ctgcCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCgg                       >  1:1053986/1‑71 (MQ=255)
                                             ctgcCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCgg                       <  1:4568831/71‑1 (MQ=255)
                                             ctgcCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCgg                       >  1:3163480/1‑71 (MQ=255)
                                             ctgcCAGTGTGGCTGATCCTCGCACCGCGCGACTATCTGGCAACCTTCCTGAAAATCGGCGTTAtcgtcg                        <  1:5292467/70‑1 (MQ=255)
                                             ctgcCAGTGTGGCTGATCCTCGCACCGCGCGACTATCTGGCAACCTTCCTGAAAATCGGCGTTAtcgtc                         <  1:1211863/69‑1 (MQ=255)
                                             ctgcCAGTGTGGCTGATCCTCGCACCGCGCGACTATCTGGCAACCTTCCTGAAAATCGGCGTTAtcgtc                         <  1:5149104/69‑1 (MQ=255)
                                                ccGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGa                                          <  1:1612282/49‑1 (MQ=255)
                                                ccGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCgg                       <  1:1553514/68‑1 (MQ=255)
                                                ccGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCgg                       <  1:5232363/68‑1 (MQ=255)
                                                  ggTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTg                   <  1:1446727/70‑1 (MQ=255)
                                                    tgtgGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGc                 <  1:5483196/70‑1 (MQ=255)
                                                    tgtgGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCg                <  1:457091/71‑1 (MQ=255)
                                                     gtgGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCgg               >  1:2327305/1‑71 (MQ=255)
                                                      tgGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGaa                                         <  1:2177273/44‑1 (MQ=255)
                                                        gTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATGGTCGGt                      <  1:5998/61‑1 (MQ=255)
                                                        gTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCTGTCTGGCGGTAg            <  1:1606371/71‑1 (MQ=255)
                                                        gTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAg            >  1:3797468/1‑71 (MQ=255)
                                                         ttACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAgg           >  1:2483854/1‑71 (MQ=255)
                                                           aCTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATc                                     >  1:5690135/1‑43 (MQ=255)
                                                           aCTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGAc                                >  1:2137469/1‑48 (MQ=255)
                                                           aCTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTg                   <  1:5039293/61‑1 (MQ=255)
                                                             tgctgGCCCCGCGTGACTACCTCTGTACCTTCCTGAAAATCGGGACTATCGTCGGt                      <  1:3512216/56‑1 (MQ=255)
                                                             tgctgGCCCCGCGTGACTACCTCTCTACCTTCCTGaaa                                        <  1:1899124/38‑1 (MQ=255)
                                                             tgctgGCCCCGCGTGACTACCTCTCTACCTTCCTGaaa                                        <  1:1940238/38‑1 (MQ=255)
                                                             tgctgGCCCCGCGTGACTACCTCTCTACCTTCCTGaaa                                        <  1:4448642/38‑1 (MQ=255)
                                                             tgctgGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGTGACTATCGTCGGt                      <  1:1516846/56‑1 (MQ=255)
                                                             tgctgGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGt                      <  1:2109879/56‑1 (MQ=255)
                                                               ctgGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCAtttt     >  1:2013277/1‑71 (MQ=255)
                                                               ctgGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCAtt       <  1:4871700/69‑1 (MQ=255)
                                                                  gCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCATTTTGAt  <  1:613552/71‑1 (MQ=255)
                                                                  gCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCATTTTGAt  <  1:3925491/71‑1 (MQ=255)
                                                                    |                                                                    
CAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGGCATTTTGAT  >  minE/385553‑385689

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: