Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 217,996 A→T 9.7% intergenic (‑88/‑18) yafU ← / → yafF predicted inner membrane protein/conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE217,9960AT 9.7% 172.3 / 9.7 62intergenic (‑88/‑18)yafU/yafFpredicted inner membrane protein/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (36/20);  new base T (4/2);  total (40/22)
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.46e-01

GCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTT  >  minE/217927‑218065
                                                                     |                                                                     
gctgATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAg                                                                      >  1:1222233/1‑71 (MQ=255)
 ctgATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAgg                                                                     >  1:215928/1‑71 (MQ=255)
 ctgATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAg                                                                      <  1:1920259/70‑1 (MQ=255)
 ctgATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTAcaca                                                                       <  1:1886206/69‑1 (MQ=255)
    aTTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCg                                                                   <  1:155539/70‑1 (MQ=255)
    aTTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGt                                                                  <  1:1578321/71‑1 (MQ=255)
     tttAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTc                                                                 >  1:2123645/1‑71 (MQ=255)
     tttAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTc                                                                 <  1:1775655/71‑1 (MQ=255)
       tAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCt                                                                <  1:3709668/70‑1 (MQ=255)
       tAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCt                                                                <  1:2880065/70‑1 (MQ=255)
          ccGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGAc                                                            <  1:1749249/71‑1 (MQ=255)
             cAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACg                                                           >  1:1165205/1‑69 (MQ=255)
                agaAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTa                                                      <  1:915251/71‑1 (MQ=255)
                agaAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTa                                                      <  1:930339/71‑1 (MQ=255)
                      tCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTaatgaat                                                >  1:3099758/1‑71 (MQ=255)
                       cGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAgg                                                                     >  1:2683932/1‑49 (MQ=255)
                       cGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAgg                                                                     >  1:992667/1‑49 (MQ=255)
                        gCCTCAGCAATCCGAAACCGCTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTaatgaatg                                               <  1:1040046/70‑1 (MQ=255)
                         ccACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTaatgaatg                                               >  1:2317458/1‑69 (MQ=255)
                          cacaGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTaatgaatg                                               >  1:835169/1‑68 (MQ=255)
                           acaGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGc                                                                    >  1:676341/1‑46 (MQ=255)
                           acaGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAga                                           >  1:568172/1‑71 (MQ=255)
                           acaGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAga                                           >  1:2418661/1‑71 (MQ=255)
                           acaGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAga                                           >  1:1916021/1‑71 (MQ=255)
                           acaGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAga                                           >  1:964734/1‑71 (MQ=255)
                           acaGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAga                                           >  1:2312940/1‑71 (MQ=255)
                            caGCGATCCGATATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAgac                                          >  1:1024586/1‑71 (MQ=255)
                            caGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGTTAATGAATGAAgac                                          <  1:2919111/71‑1 (MQ=255)
                            caGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAgac                                          >  1:2118770/1‑71 (MQ=255)
                            caGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAgac                                          >  1:2757676/1‑71 (MQ=255)
                            caGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAgac                                          >  1:3698327/1‑71 (MQ=255)
                            caGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAgac                                          >  1:2099471/1‑71 (MQ=255)
                               cAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAgacgac                                       >  1:2616744/1‑71 (MQ=255)
                                 aTCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTa                                     >  1:383285/1‑71 (MQ=255)
                                 aTCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTa                                     >  1:3480905/1‑71 (MQ=255)
                                   ccGAAACCACTGGCACGTGGAGAATAAGCTGcaca                                                                       >  1:767015/1‑35 (MQ=255)
                                   ccGAAACCACTGGCACGTGGAGAATAAGCTGcaca                                                                       >  1:389507/1‑35 (MQ=255)
                                   ccGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAgac                                          <  1:2786662/64‑1 (MQ=255)
                                      aaaCCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTa                                                      >  1:485452/1‑49 (MQ=255)
                                      aaaCCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAAtgaatgaa                                             >  1:516419/1‑58 (MQ=255)
                                       aaCCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGTACGTGGTAATGAATGAAGACGACTACaaaa                                >  1:2242711/1‑70 (MQ=255)
                                       aaCCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACaaaa                                <  1:2723146/70‑1 (MQ=255)
                                       aaCCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACaaaa                                <  1:3385195/70‑1 (MQ=255)
                                       aaCCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACaaaa                                >  1:1292249/1‑70 (MQ=255)
                                         tcACTGGCACGTGGAGAATAAGCTGCACTGGCGTCt                                                                <  1:3583781/35‑1 (MQ=39)
                                         tcACTGGCACGTGGAGAATAAGCTGCACTGGCGTCt                                                                <  1:1239757/35‑1 (MQ=39)
                                          cACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAAt                               >  1:2429495/1‑68 (MQ=255)
                                           aCTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATaaga                           <  1:1578328/71‑1 (MQ=255)
                                            cTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGCCGACTACAAAATaag                            >  1:349178/1‑69 (MQ=255)
                                             tGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAgacgac                                       >  1:1752202/1‑57 (MQ=255)
                                             tGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTa                                     >  1:544938/1‑59 (MQ=255)
                                             tGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTa                                     >  1:2053138/1‑59 (MQ=255)
                                               gCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAagaaga                        <  1:3267037/70‑1 (MQ=255)
                                               gCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAagag                       >  1:3335552/1‑71 (MQ=255)
                                                cACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAagag                       >  1:2942886/1‑70 (MQ=255)
                                                cACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAagag                       >  1:1323149/1‑70 (MQ=255)
                                                cACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAgg                      <  1:25306/71‑1 (MQ=255)
                                                    tGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTa                                     >  1:2129699/1‑52 (MQ=255)
                                                      gagaATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAgac                                          >  1:3220188/1‑45 (MQ=255)
                                                      gagaATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAgac                                          >  1:2115603/1‑45 (MQ=255)
                                                        gaATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGaa                    >  1:2931290/1‑65 (MQ=255)
                                                        gaATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACgca               <  1:302336/70‑1 (MQ=255)
                                                           tAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGcagcag           <  1:1231078/71‑1 (MQ=255)
                                                              gCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAat        <  1:880816/71‑1 (MQ=255)
                                                               cTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACg                 >  1:2291139/1‑61 (MQ=255)
                                                                     aGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTAtttt  >  1:533800/1‑70 (MQ=255)
                                                                     aGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTAtttt  >  1:530657/1‑70 (MQ=255)
                                                                     |                                                                     
GCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTT  >  minE/217927‑218065

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: