Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 385,613 C→A 11.0% L275M (CTG→ATG) ‡ cstA → carbon starvation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE385,6130CA11.0% 159.5 / 4.6 53L275M (CTG→ATG) ‡cstAcarbon starvation protein
Reads supporting (aligned to +/- strand):  ref base C (25/22);  new base A (2/4);  total (27/26)
Fisher's exact test for biased strand distribution p-value = 4.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGG  >  minE/385544‑385680
                                                                     |                                                                   
aCCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTAct                                                                    >  1:196598/1‑71 (MQ=255)
    gCGTGCAGTTGACCTGGATGCTGGTGGGTTACTGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTActgctg                                                                >  1:3156345/1‑71 (MQ=255)
    gCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTActgctg                                                                >  1:3949890/1‑71 (MQ=255)
         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)
                 cTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTAc                                                                     <  1:851523/53‑1 (MQ=255)
                 cTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTAc                                                                     <  1:311343/53‑1 (MQ=255)
                 cTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTAc                                                                     <  1:1617529/53‑1 (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:4752566/1‑71 (MQ=255)
                    gATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTAcc                                                >  1:1069388/1‑71 (MQ=255)
                    gATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTAc                                                 <  1:4483431/70‑1 (MQ=255)
                     aTGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACct                                               >  1:4723507/1‑71 (MQ=255)
                          ggtggGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTa                                          >  1:4706999/1‑71 (MQ=255)
                          ggtggGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTa                                          >  1:144673/1‑71 (MQ=255)
                           gtggGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTAc                                         >  1:2173621/1‑71 (MQ=255)
                            tggGTTACGGTTTTGTGGCAGCGGTGCTGCCGGTGTGGTTActgctg                                                                <  1:2753046/47‑1 (MQ=255)
                               gTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTAct                                                                    <  1:1890492/40‑1 (MQ=255)
                               gTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTc                                     <  1:3594466/71‑1 (MQ=255)
                                 tACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCt                                   >  1:4062551/1‑71 (MQ=255)
                                 tACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCt                                   >  1:556715/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:3265248/63‑1 (MQ=255)
                                                gcggTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATc                            <  1:5671157/63‑1 (MQ=255)
                                                gcggTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATc                            <  1:5542898/63‑1 (MQ=255)
                                                gcggTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATc                            <  1:2092252/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:3163480/1‑71 (MQ=255)
                                                      ctgcCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCgg              >  1:1053986/1‑71 (MQ=255)
                                                      ctgcCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCgg              <  1:4568831/71‑1 (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:5232363/68‑1 (MQ=255)
                                                         ccGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCgg              <  1:1553514/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)
                                                                     |                                                                   
ACCGGCGTGCAGTTGACCTGGATGCTGGTGGGTTACGGTTTTGTGGCGGCGGTGCTGCCGGTGTGGTTACTGCTGGCCCCGCGTGACTACCTCTCTACCTTCCTGAAAATCGGGACTATCGTCGGTCTGGCGGTAGG  >  minE/385544‑385680

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: