Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,089,817 C→T 16.6% V99V (GTC→GTT yhaK → predicted pirin‑related protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,089,8170CT16.6% 154.1 / 21.7 66V99V (GTC→GTTyhaKpredicted pirin‑related protein
Reads supporting (aligned to +/- strand):  ref base C (24/31);  new base T (5/6);  total (29/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

CGATAGCGAAGGCAATCATGTTCAGGCCAGCGCCGGAGAGGCGTTGCTGCTCTCTACCCAGCCGGGTGTCAGCTATAGCGAACACAATCTCAGCAAAGACAAACCGTTAACGCGAATGCAGCTTTGGCTGGACGCCTGCC  >  minE/2089748‑2089887
                                                                     |                                                                      
cgATAGCGAAGGCAATCATGTTCAGGCCAGCGCCGGAGAGGCGTTGCTGCTCTCTACCCAGCCGGGTGTCa                                                                       <  1:2993848/71‑1 (MQ=255)
 gATAGCGAAGGCAATCATGTTCAGGCCAGCGCCGGAGAGGCGTTGCTGCTCTCTACCCAGCCGGGTGTCAg                                                                      >  1:103253/1‑71 (MQ=255)
  aTAGCGAAGGCAATCATGTTCAGGCCAGCGCCGGAGAGGCGTTGCTGCTCTCTACCCAGCCGGGTGTCAGc                                                                     <  1:3580230/71‑1 (MQ=255)
  aTAGCGAAGGCAATCATGTTCAGGCCAGCGCCGGAGAGGCGTTGATGCTCTCTACCCAGCCGGGTGTc                                                                        <  1:3327351/68‑1 (MQ=255)
      cGAAGGCAATCATGTTCAGGCCAGCGCCGGAGAGGCGTTGCTGCTCTCTACCCAGCCGGGTGTCAGCta                                                                   >  1:3050564/1‑69 (MQ=255)
        aaGGCAATCATGTTCAGGCCAGCGCCGGAGAGGCGTTGCTGCTCTCTACCCAGCCGGGTGTCAGCTATAGc                                                               >  1:2895826/1‑71 (MQ=255)
         aGGCAATCATGTTCAGGCCAGCGCCGGAGAGGCGTTGCTGCTCTCTACCCAGCCGGGTGTCAGCTATAGCg                                                              <  1:947625/71‑1 (MQ=255)
         aGGCAATCATGTTCAGGCCAGCGCCGGAGAGGCGTTGCTGCTCTCTACCCAGCCGGGTGTCAGCTATAGCg                                                              <  1:891320/71‑1 (MQ=255)
          ggCAATCATGTTCAGGCCAGCGCCGGAGAGGCGTTGCTGCTCTCTACCCAGCCGGGTGTCAGCTATAGCg                                                              >  1:2888190/1‑70 (MQ=255)
          ggCAATCATGTTCAGGCCAGCGCCGGAGAGGCGTTGCTGCTCTCTACCCAGCCGGGTGTCAGCTATAGCg                                                              >  1:2073582/1‑70 (MQ=255)
             aaTCATGTTCAGGCCAGCGCCGGAGAGGCGTTGCTGCTCTCTACCCAGCCGGGTGTCAGCTATAGCGAaca                                                          >  1:2397963/1‑71 (MQ=255)
             aaTCATGTTCAGGCCAGCGCCGGAGAGGCGTTGCTGCTCTCTACCCAGCCGGGTGTCAGCTATAGCGAaca                                                          <  1:2614307/71‑1 (MQ=255)
                 aTGTTCAGGCCAGCGCCGGAGAGGCGTTGCTGCTCTCTACCCAGCCGGGTGTCAGCTATAGCGAACACAAt                                                      >  1:3905007/1‑71 (MQ=255)
                 aTGTTCAGGCCAGCGCCGGAGAGGCGTTGCTGCTCTCTACCCAGCCGGGTGTCAGCTATAGCGAACACAAt                                                      >  1:3401915/1‑71 (MQ=255)
                     tCAGGCCAGCGCCGGAGAGGCGTTGCTGCTCTCTACCCAGCCGGGTGTCAGCTATAGCGAACACAAtct                                                    >  1:666983/1‑69 (MQ=255)
                     tCAGGCCAGCGCCGGAGAGGCGTTGCTGCTCTCTACCCAGCCGGGTGTCAGCTATAGCGAACACAAtct                                                    >  1:1618357/1‑69 (MQ=255)
                          ccAGCGCCGGAGAGGCGTTGCTGCTCTCTACCCAGCCGGGTGTCAGCTATAGCGAACACAATCTCAGCaaa                                             <  1:104719/71‑1 (MQ=255)
                          ccAGCGCCGGAGAGGCGTTGCTGCTCTCTACCCAGCCGGGTGTCAGCTATAGCGAACACAATCTCAGCaaa                                             <  1:2464578/71‑1 (MQ=255)
                              cgcCGGAGAGGCGTTGCTGCTCTCTACCCAGCCGGGTGTCAGCTATAGCGAACACAATCTCAGCAAAGACa                                         <  1:3382813/71‑1 (MQ=255)
                              cgcCGGAGAAGCGTTGCTGCTCTCTACCCAGCCGGGTGTCAGCTATAGCGAACACAAACTCAGCAAAGACa                                         <  1:1897389/71‑1 (MQ=255)
                                 cGGAGAGGCGTTGCTGCTCTCTACCCAGCCGGGTGTCAGCTATAGCGAACACAATCTCAGCAAAGACAAAc                                      >  1:2399519/1‑71 (MQ=255)
                                  ggAGAGGCGTTGCTGCTCTCTACCCAGCCGGGTGTTAGCTATAGCGAACACAATCTCAGCAAAGACAAAcc                                     >  1:32083/1‑71 (MQ=255)
                                   gagaGGCGTTGCTGCTCTCTACCCAGCCGGGTGTTAGCTATAGCGAACACAATCTCAGCAAAGACAAACCg                                    >  1:2049807/1‑71 (MQ=255)
                                   gagaGGCGTTGCTGCTCTCTACCCAGCCGGGTGTCAGCTATAGCGAACACAATCTCAGCAAAGACAAACCg                                    >  1:397789/1‑71 (MQ=255)
                                    agagGCGTTGCTGCTCTCTACCCAGCCGGGTGTCAGCTATAGCGAACACAATCTCAGCAAAGACAAACCGt                                   <  1:2475581/71‑1 (MQ=255)
                                    agagGCGTTGCTGCTCTCTACCCAGCCGGGTGTCAGCTATAGCGAACACAATCTCAGCAAAGACAAACCGt                                   <  1:3416831/71‑1 (MQ=255)
                                    agagGCGTTGCTGCTCTCTACCCAGCCGGGTGTCAGCTATAGCGAACACAATCTCAGCAAAGACAAACCGt                                   <  1:2213957/71‑1 (MQ=255)
                                     gagGCGTTGCTGCTCTCTACCCAGCCGGGTGTCAGCTATAGCGAACACAATCTCAGCAAAGACAAACCGtt                                  <  1:1466268/71‑1 (MQ=255)
                                     gagGCGTTGCTGCTCTCTACCCAGCCGGGTGTCAGCTATAGCGAACACAATCTCAGCAAAGACAAACCGt                                   >  1:3417615/1‑70 (MQ=255)
                                      agGCGTTGCTGCTCTCTACCCAGACGGGTGTCAGCTATAGCGAACACAATCTCAGCAAAGACAAACCGtt                                  <  1:2901492/70‑1 (MQ=255)
                                      agGCGTTGCTGCTCTCAACCCAGCCGGGTGTTAGCTATAGCGAAcaca                                                        >  1:2544001/1‑48 (MQ=255)
                                       ggCGTTGCTGCTCTCTACCCAGCCGGGTGTTAGCTATAGCGAACACAATCTCAGCAAAGACAAACCGTTaa                                <  1:3391440/71‑1 (MQ=255)
                                        gCGTTGCTGCTCTCTACCCAGCCGGGTGTTAGCTATAGCGAACACAATCTCAGCAAAGACAAACCGTTAAc                               >  1:3507646/1‑71 (MQ=255)
                                        gCGTTGCTGCTCTCTACCCAGCCGGGTGTCAGCTATAGCGAACACAATCTCAGCAAAGACAAAcc                                     <  1:36612/65‑1 (MQ=255)
                                            tgctgcTCTCTACCCAGCCGGGTGTCAGCTATAGCGAACACAATCTCAGCAAAGACAAACCATTAAcgcg                            >  1:3556392/1‑70 (MQ=255)
                                              ctgctCTCTACCCAGCCGGGTGTCAGCTATAGCGAACACAATCTCAGCAAAGACAAACCGTTAACGCGaa                          <  1:1658135/70‑1 (MQ=255)
                                               tgctCTCTACCCAGCCGGGTGTCAGCTATAGCGAacac                                                         <  1:2826845/38‑1 (MQ=255)
                                                 ctctctACCCAGCCGGGTGTCAGCTATAGCGAACACAATCTCAGCAAAGACAAACCGTTAAc                               >  1:343327/1‑62 (MQ=255)
                                                 ctctctACCCAGCCGGGTGTCAGCTATAGCGAACACAATCTCAGCAAAGACAAACCGTTAACGCGAATGCa                      <  1:2932554/71‑1 (MQ=255)
                                                  tctctACCCAGCCGGGTGTCAGCTATAGCGAACACAATCTCAGCAAAGACAAACCGTTAACGCGAATg                        <  1:2926509/68‑1 (MQ=255)
                                                    tctACCCAGCCTGGTGTCAGCTATAGCGAACACAATCTCAGCAAAGACAAAc                                      <  1:2187721/52‑1 (MQ=255)
                                                    tctACCCAGCCGGGTGTCAGCTATAGCGAACACAATCTCAGCAAAGACAAACCGTTAACGCGAATGCAGCt                   >  1:3782505/1‑71 (MQ=255)
                                                     ctACGCAGCCGGGTGTCAGCTATAGCGAACACAATCTCAGCAAAGACAAACCGTTAACGCGAATGCAGCtt                  >  1:2946742/1‑71 (MQ=255)
                                                      tACCGAGCCGGGTGTCAGCTATAGCGAACACAATCTCAGCAAAGACAAACCGTTAACGCGAATGCAGCttt                 >  1:3833780/1‑71 (MQ=255)
                                                       agccAGCCGGGTGTCAGCTATAGCGAACACAATCTCAGCAAAGACAAACCGTTAACGCGAATGCAGCt                   <  1:2420439/66‑1 (MQ=255)
                                                          cAGCCGGGTGTCAGCTATAGCGAACACAATCTCAGCaaa                                             >  1:453900/1‑39 (MQ=255)
                                                          cAGCCGGGTGTCAGCTATAGCGAACACAATCTCAGCaaa                                             >  1:1819527/1‑39 (MQ=255)
                                                           aGCCGGGTGTTAGCTATAGCGAACACAATCTCAGCAAAGACAAACCGTTAACGCGAATGCAGCTTTGGCTg            <  1:3785620/71‑1 (MQ=255)
                                                           aGCCGGGTGTTAGCTATAGCGAACACAATCTCAGCAAAGACAAACCGTTAACGCGAATGCAGCTTTGGCTg            <  1:3255109/71‑1 (MQ=255)
                                                           aGCCGGGTGTCAGCTATAGCGAACACAATCTCAGCAAAGACaa                                        <  1:99827/43‑1 (MQ=255)
                                                           aGCCGGGTGTCAGCTATAGCGAACACAATCTCAGCAAAGACaa                                        <  1:16923/43‑1 (MQ=255)
                                                           aGCCGGGTGTCAGCTATAGCGAACACAATCTCAGCAAAGACAAACCGTTAACGCGAATGCAGCTTTGGCTg            >  1:990427/1‑71 (MQ=255)
                                                            gCCGGGTGTCAGCTATAGCGAACACAATCTCAGCAAAGACaaa                                       <  1:2544334/43‑1 (MQ=255)
                                                            gCCGGGTGTCAGCTATAGCGAACACAATCTCAGCAAAGACAAACCGTTAACGCGAATGCAGCTTTGGCTg            >  1:2685777/1‑70 (MQ=255)
                                                              cGGGTGTTAGCTATAGCGAACACAATCTCAGCAAAGACAAACCGTTAACGCGAATGCAGCt                   <  1:394880/61‑1 (MQ=255)
                                                              cGGGTGTTAGCTATAGCGAACACAATCTCAGCAAAGACAAACCGTTAACGCGAATGCAGCTTTGGCTGGAc         <  1:2647643/71‑1 (MQ=255)
                                                              cGGGTGTTAGCTATAGCGAACACAATCTCAGCAAAGACAAACCGTTAACGCGAATGCAGCTTTGGCTGGAc         <  1:529779/71‑1 (MQ=255)
                                                              cGGGTGTCAGCTATAGCGAACACAATCTCAGCAAAGACAAACCGTTaa                                >  1:1169207/1‑48 (MQ=255)
                                                              cGGGTGTCAGCTATAGCGAACACAATCTCAGCAAAGACAAACCGTTaa                                >  1:2560408/1‑48 (MQ=255)
                                                              cGGGTGTCAGCTATAGCGAACACAATCTCAGCAAAGACAAACCGTTAACGCGAATGCAGCTTTGGCTGGAc         <  1:2416198/71‑1 (MQ=255)
                                                                ggTGTTAGCTATAGCGAACACAATCTCAGCAAAGACAAAc                                      >  1:3891675/1‑40 (MQ=255)
                                                                  tgtCAGCTATAGCGAACACAATCTCAGCAAAGACAAACCGTTAACGCGAATGCAGCTTTGGCTGGACGc       <  1:181834/69‑1 (MQ=255)
                                                                    tCAGCTATAGCGAACACAATCTCAGCAAAGACAAAc                                      <  1:3147385/36‑1 (MQ=255)
                                                                    tCAGCTATAGCGAACACAATCTCAGCAAAGACAAAc                                      <  1:2798741/36‑1 (MQ=255)
                                                                    tCAGCTATAGCGAACACAATCTCAGCAAAGACAAACCGTTa                                 <  1:1762658/41‑1 (MQ=255)
                                                                    tCAGCTATAGCGAACACAATCTCAGCAAAGACAAACCGTTAACGCGAATGCAGCTTTGGCTGGACGCCTGc   <  1:666722/71‑1 (MQ=255)
                                                                    tCAGCTATAGCGAACACAATCTCAGCAAAGACAAACCGTTAACGCGAATGCAGCTTTGGCTGGACGCCTGc   <  1:1395614/71‑1 (MQ=255)
                                                                     cAGCTATAGCGAACACAATCTCAGCAAAGACAAACCGTTAACGCGAAt                         >  1:544477/1‑48 (MQ=255)
                                                                     cAGCTATAGCGAACACAATCTCAGCAAAGACAAACCGTTAACGCGAATGCAGCTTTGGCTGGACGCCTGcc  >  1:636301/1‑71 (MQ=255)
                                                                     cAGCTATAGCGAACACAATCTCAGCAAAGACAAACCGTTAACGAGAATGCAGCTTTGGCTGGACGCCTGcc  >  1:2503837/1‑71 (MQ=255)
                                                                     |                                                                      
CGATAGCGAAGGCAATCATGTTCAGGCCAGCGCCGGAGAGGCGTTGCTGCTCTCTACCCAGCCGGGTGTCAGCTATAGCGAACACAATCTCAGCAAAGACAAACCGTTAACGCGAATGCAGCTTTGGCTGGACGCCTGCC  >  minE/2089748‑2089887

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: