Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,741,082 T→C 12.9% K120R (AAA→AGA)  rseB ← anti‑sigma factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,741,0820TC12.9% 136.5 / 10.8 54K120R (AAA→AGA) rseBanti‑sigma factor
Reads supporting (aligned to +/- strand):  ref base T (22/25);  new base C (3/4);  total (25/29)
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

ACTTCGCAAAGACGATCAGCAATACGCGTGCGCCCGACGGAGATAAAGTCGTAGTAAGGAGAAAGGCGTTTGAAATCGGTATAGATAAGCGATGGCAGAGAATCAACAATGTAATCGCCATTAAGCGTGAACGGTTC  >  minE/1741013‑1741149
                                                                     |                                                                   
aCTTCGCAAAGACGATCAGCAATACGCGTGCGCCCGACGGAGATAAAGTCGTAGTAAGGAGAAAGGCGtt                                                                     <  1:3090482/70‑1 (MQ=255)
   tCGCAAAGACGATCAGCAATACGCGTGCGCCCGACGGAGATAAAGTCGTAGTAAGGAGAAAGGCGTTTg                                                                   >  1:1477595/1‑69 (MQ=255)
     gCAAAGACGATCAGCAATACGCGTGCGCCCGACGGAGATAAAGTCGTAGTAAGGAGAAAGGCGTTTGAAAt                                                               <  1:2433553/71‑1 (MQ=255)
     gCAAAGACGATCAGCAATACGCGTGCGCCCGACGGAGATAAAGTCGTAGTAAGGAGAAAGGCGTTTGAAAt                                                               <  1:3431616/71‑1 (MQ=255)
      cAAAGACGATCAGCAATACGCGTGCGCCCGACGGAGATAAAGTCGTAGTAAGGAGAAAGGCGTTTGAAATc                                                              <  1:1957553/71‑1 (MQ=255)
       aaaGACGATCAGCAATACGCGTGCGCCCGACGGAGATAAAGTCGTAGTAAGGAGAAAGGCGTTTGAAATCg                                                             <  1:1400593/71‑1 (MQ=255)
        aaGACGATCAGCAATACGCGTGCGCCCGACGGAGATAAAGTCGTAGTAAGGAGAAAGGCGTTTGAAATCgg                                                            <  1:1327730/71‑1 (MQ=255)
         aGACGATCAGCAATACGCGTGCGCCCGACGGAGATAAAGTCGTAGTAAGGAGAAAGGCGTTTGAAATCgg                                                            >  1:3254961/1‑70 (MQ=255)
          gACGATCAGCAATACGCGTGCGCCCGACGGAGATAAAGTCGTAGTAAGGAGAAAGGCGTTTGAAATCGGta                                                          <  1:653058/71‑1 (MQ=255)
          gACGATCAGCAATACGCGTGCGCCCGACGGAGATAAAGTCGTAGTAAGGAGAAAGGCGTTTGAAATCGGt                                                           <  1:2586232/70‑1 (MQ=255)
             gATCAGCAATACGCGTGCGCCCGACGGAGATAAAGTCGTAGTAAGGAGAAAGGCGTTTGAAATCgg                                                            <  1:2262664/66‑1 (MQ=255)
                   cAATACGCGTGCGCCCGACGGAGATAAAGTCGTAGTAAGGAGAAAGGCGTTTGAAATCGGTATAGATAAGc                                                 >  1:2606985/1‑71 (MQ=255)
                   cAATACGCGTGCGCCCGACGGAGATAAAGTCGTAGTAAGGAGAAAGGCGTCTGAAATCgg                                                            <  1:1554617/60‑1 (MQ=255)
                     aTACGCGTGCGCCCGACGGAGATAAAGTCGTAGTAAGGAGAAAGGCGTTTGAAATCGGTATAGa                                                      >  1:3507136/1‑64 (MQ=255)
                     aTACGCGTGCGCCCGACGGAGATAAAGTCGTAGTAAGGAGAAAGGCGTTTGAAATCGGTATAGATAAGCGa                                               >  1:1855693/1‑71 (MQ=255)
                     aTACGCGTGCGCCCGACGGAGATAAAGTCGTAGTAAGGAGAAAGGCGTCTGAAATCGGTATAGATAAGCGa                                               >  1:2537463/1‑71 (MQ=255)
                       aCGCGTGCGCCCGACGGAGATAAAGTCGTAGTAAGGAGAAAGGCGTTTGAAATCGGTATAGATAAGCGATg                                             >  1:782441/1‑71 (MQ=255)
                         gcgTGCGCCCGACGGAGATAAAGTCGTAGTAAGGAGAAAGGCGTTTGaa                                                                 <  1:1202523/49‑1 (MQ=255)
                          cgTGCGCCCGACGGAGATAAAGTCGTAGTAAGGAGAAAGGCGTTTGAAATCGGTATAGATAAGCGATGGCa                                          <  1:1813191/71‑1 (MQ=255)
                           gTGCGCCCGACGGAGATAAAGTCGTAGTAAGGAGAAAGGCGTTTGAAAt                                                               >  1:2248375/1‑49 (MQ=255)
                            tGCGCCCGACGGAGATAAAGTCGTAGTAAGGAGAAAGGCGTTTGAAATCGGTATAGATAAGCGATGGCa                                          >  1:3724271/1‑69 (MQ=255)
                                 ccGACGGAGATAAAGTCGTAGTAAGGAGAAAGGCGTTTGAAATCGGTATAGATAAGCGATGGCAGAGAATc                                   <  1:956968/71‑1 (MQ=255)
                                   gACGGAGATAAAGTCGTAGTAAGGAGAAAGGCGTTTGAAATCGGTATAGATAAGCGATGGCaga                                        <  1:968188/64‑1 (MQ=255)
                                   gACGGAGATAAAGTCGTAGTAAGGAGAAAGGCGTTTGAAATCGGTATAGATAAGCGATGGCAGAGAATca                                  >  1:839486/1‑70 (MQ=255)
                                     cGGAGATAAAGTCGTAGTAAGGAGAAAGGCGTTTGAAATCGGTATAGATAAGCg                                                >  1:3713255/1‑54 (MQ=255)
                                     cGGAGATAAAGTCGTAGTAAGGAGAAAGGCGTTTGAAATCGGTATAGATAAGCGATGGCAGAGAATcaaca                               <  1:1875596/71‑1 (MQ=255)
                                     cGGAGATAAAGTCGTAGTAAGGAGAAAGGCGTCTGAAATCGGTATAGATAAGCGATGGCAGAGAAt                                    >  1:1487141/1‑66 (MQ=255)
                                      ggAGATAAAGTCGTAGTAAGGAGAAAGGCGTTTGAAATCGGTATAGATAAGCGATGGCAGAGAATcaacaa                              >  1:2228769/1‑71 (MQ=255)
                                       gagaTAAAGTCGTAGTAAGGAGAAAGGCGTCTGaaa                                                                <  1:2237457/36‑1 (MQ=255)
                                           tAAAGTCGTAGTAAGGAGAAAGGCGTTTGAAATCGGTATAGATAAGCGATGGCAGAGAATCAACAATGTaa                         <  1:3471726/71‑1 (MQ=255)
                                            aaaGTCGTAGTAAGGAGAAAGGCGTTTGAAATCGGTATAGATAAGCGATGGCAGAGAATcaaca                               <  1:212005/64‑1 (MQ=255)
                                             aaGTCGTAGTAAGGAGAAAGGCGTCTGAAATCGGTATAGATAAGCGATGGCAGAGAATCAACAATGTAATc                       <  1:2928538/71‑1 (MQ=255)
                                              aGTCGTAGTAAGGAGAAAGGCGTTTGAAATCGGTATAGATAAGCGATg                                             >  1:3850231/1‑48 (MQ=255)
                                              aGTCGTAGTAAGGAGAAAGGCGTTTGAAATCGGTATAGATAAGCGATg                                             >  1:3342363/1‑48 (MQ=255)
                                              aGTCGTAGTAAGGAGAAAGGCGTTTGAAATCGGTATAGATAAGCGATGGCAGAGAATCAACAATGTAATCg                      >  1:3260402/1‑71 (MQ=255)
                                              aGTCGTAGTAAGGAGAAAGGCGTCTGAAATCGGTATAGATAAGCGATGGCAGAGAATCAACAATGTAATCg                      >  1:1606518/1‑71 (MQ=255)
                                                tCGTAGTAAGGAGAAAGGCGTTTGAAATCGGTATAGATAAGCGATGGCAGAGAATCAACAATGTAATCGcc                    <  1:2780309/71‑1 (MQ=255)
                                                   tagtaAGGAGAAAGGCGTTTGAAATCGGTATAGATAAGCGATGGCAGAGAATCAACAATGTAATCGCCAtt                 <  1:3858838/71‑1 (MQ=255)
                                                   tagtaAGGAGAAAGGCGTTTGAAATCGGTATAGATAAGCGATGGCAGAGAATCAACAATGTAATCGCCAtt                 <  1:3786410/71‑1 (MQ=255)
                                                     gtaAGGAGAAAGGCGTTTGAAATCGGTATAGATAAGCGATGGCAGAGAATCAACAATGTAATCGCCAtt                 >  1:2133534/1‑69 (MQ=255)
                                                     gtaAGGAGAAAGGCGTTTGAAATCGGTATAGATAAGCGATGGCAGAGAATCAACAATGTAATCGCCATTaa               <  1:3191121/71‑1 (MQ=255)
                                                           agaAAGGCGTTTGAAATCGGTATAGATAAGCGATGGCAGAGAATCAACAATGTaa                         <  1:3777908/55‑1 (MQ=255)
                                                            gaAAGGCGTTTGAAATCGGTATAGATAAGCGATGGCAGAGAATCAACAATGTAATCGCCa                   >  1:227924/1‑60 (MQ=255)
                                                            gaAAGGCGTTTGAAATCGGTATAGATAAGCGATGGCAGAGAATCAACAATGTAATCGCCATTAAGCGTGaa        >  1:2594020/1‑71 (MQ=255)
                                                            gaAAGGCGTTTGAAATCGGTATAGATAAGCGATGGCAGAGAATAAACAATGTAATCGCCa                   >  1:381665/1‑60 (MQ=255)
                                                            caaaGGCGTTTGAAATCGGTATAGATAAGCGATGGCAGAGAATCAACAATGTAATCGCCATTAAGc             <  1:3436107/65‑1 (MQ=255)
                                                              aaGGCGTTTGAAATCGGTATAGATAAGCGATGGCAGAGaa                                     >  1:187415/1‑40 (MQ=255)
                                                              aaGGCGTTTGAAATCGGTATAGATAAGCGATGGCAGAGAATCAACAATGTAATCGCCATTAAGCGTGAACg      <  1:1821275/71‑1 (MQ=255)
                                                               aGGCGTTTGAAATCGGTATAGATAAGCGATGGCAGAGAATCAACAATGTAATCGCCATTAAGTGTGAACgg     >  1:1941989/1‑71 (MQ=255)
                                                               aGGCGTTTGAAATCGGTATAGATAAGCGATGGCAGAGAATCAACAATGTAATCGCCATTAAGCGTGAACgg     >  1:1302682/1‑71 (MQ=255)
                                                               aGGCGTTTGAAATCGGTATAGATAAGCGATGGCAGAGAATCAACAATGTAATCGCCATTAAGCGTGAACgg     >  1:2755785/1‑71 (MQ=255)
                                                                ggCGTTTGAAATCGGTATAGATAAGCGATGGCAGAGAATCAACAATGTAATCGCCATTAAGCGTGAACGGt    <  1:3672361/71‑1 (MQ=255)
                                                                 gCGTTTGAAATCGGTATAGATAAGCGATGGCAGAGAATCAACAATGTAATCGCCATTAAGCGTGAACGGtt   <  1:2503925/71‑1 (MQ=255)
                                                                  cGTTTGAAATCGGTATAGATAAGCGATGGCAGAGAATCAACAATGTAATCGCCATTAAGCGTGAACGGTTc  <  1:1124097/71‑1 (MQ=255)
                                                                  cGTCTGAAATCGGTATAGATAAGCGATGGCAGAGAATCAACAATGTAATCGCCATTAAGCGTGAACGGTTc  <  1:3816489/71‑1 (MQ=255)
                                                                    tttGAAATCGGTATAGATAAGCGATGGCAGAGAATcaaca                               <  1:3814715/40‑1 (MQ=255)
                                                                     |                                                                   
ACTTCGCAAAGACGATCAGCAATACGCGTGCGCCCGACGGAGATAAAGTCGTAGTAAGGAGAAAGGCGTTTGAAATCGGTATAGATAAGCGATGGCAGAGAATCAACAATGTAATCGCCATTAAGCGTGAACGGTTC  >  minE/1741013‑1741149

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑