Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 218,030 A→G 59.2% Y6C (TAC→TGC)  yafF → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE218,0300AG59.2% 43.2 / 83.6 76Y6C (TAC→TGC) yafFconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (19/12);  new base G (30/15);  total (49/27)
Fisher's exact test for biased strand distribution p-value = 6.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.96e-01

GAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTT  >  minE/217964‑218096
                                                                  |                                                                  
gAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCaaa                                                                <  1:3322932/71‑1 (MQ=39)
 aaaTCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCaaaa                                                               >  1:173946/1‑71 (MQ=255)
  aaCCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCaaaa                                                               >  1:726669/1‑70 (MQ=39)
  aaCCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACaaaa                                                               >  1:2559308/1‑70 (MQ=255)
  aaCCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACaaaa                                                               >  1:2347406/1‑70 (MQ=255)
   aCCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACaaaa                                                               >  1:215781/1‑69 (MQ=255)
   aCCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAAt                                                              >  1:3581221/1‑70 (MQ=255)
      aCTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATaag                                                           >  1:1165177/1‑70 (MQ=255)
        tGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACaaaa                                                               >  1:2237008/1‑64 (MQ=255)
         ggCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATaagaa                                                         >  1:3805089/1‑69 (MQ=255)
         ggCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATaagaa                                                         >  1:3800580/1‑69 (MQ=255)
            atgTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAgg                                                     >  1:1617337/3‑70 (MQ=255)
            atgTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGAAGACTGCAAAATAAGAAGAGGa                                                    <  1:1110959/69‑1 (MQ=255)
             tgTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATaag                                                           >  1:2648226/2‑63 (MQ=255)
             tgTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGa                                                    >  1:1323050/2‑70 (MQ=255)
              gTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGaaa                                                  >  1:253522/1‑71 (MQ=255)
                      tAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATgcag                                             >  1:2739297/1‑68 (MQ=255)
                         gCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAat                                       <  1:2540746/71‑1 (MQ=255)
                          cTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAat                                       >  1:1681841/1‑70 (MQ=255)
                           tGCACTGGCGTCTGGACGTGGTAATGAATGAAGCCGACTGCAAAATAAGAAGAGGAAACGCCGCAGAatta                                     >  1:2413801/1‑71 (MQ=255)
                           tGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAatta                                     >  1:3952469/1‑71 (MQ=255)
                           tGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAatt                                      >  1:1169153/1‑70 (MQ=255)
                             cacaGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAattat                                    <  1:616730/70‑1 (MQ=255)
                               ctggCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATaagaag                                                        >  1:395554/3‑48 (MQ=255)
                               ctggCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTAtttt                                 >  1:299572/3‑71 (MQ=255)
                               caGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTAtttt                                 <  1:1367906/71‑1 (MQ=255)
                                aGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGcagcag                                          <  1:3086400/61‑1 (MQ=255)
                                aGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTAtttt                                 >  1:2538903/1‑70 (MQ=255)
                                 ggCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATAagaaga                                                       >  1:2646695/1‑47 (MQ=255)
                                  gCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTAtttt                                 >  1:3619811/1‑68 (MQ=255)
                                  gCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTAtttt                                 >  1:309878/1‑68 (MQ=255)
                                    gTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGaa                                        >  1:3437745/1‑59 (MQ=255)
                                      cTGGACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGa                           >  1:334325/1‑70 (MQ=255)
                                       tGGACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTc                                >  1:2199733/1‑64 (MQ=255)
                                       tGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATa                         >  1:3495409/1‑71 (MQ=255)
                                         gACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATgcagc                                            >  1:3611270/1‑50 (MQ=255)
                                         gACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACg                       >  1:535588/1‑71 (MQ=255)
                                         gACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAat                                       <  1:2514967/55‑1 (MQ=255)
                                         gACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATa                         <  1:1624961/69‑1 (MQ=255)
                                          aCGTGTTAATTAATGAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACgg                      <  1:750625/71‑1 (MQ=255)
                                           cGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCa                               >  1:1290340/1‑61 (MQ=255)
                                             tGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGca                    >  1:3189581/1‑70 (MQ=255)
                                              ggTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGcaca                  >  1:3813831/1‑71 (MQ=255)
                                              ggTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGcaca                  >  1:907077/1‑71 (MQ=255)
                                               gTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACAt                 <  1:41141/71‑1 (MQ=255)
                                                tAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACAtt                >  1:2290569/1‑71 (MQ=255)
                                                aaatgaatgAAGATGCCAGCCGGATAAGAAGAGGAAACgca                                              <  1:1021342/40‑1 (MQ=255)
                                                 aatgaatgAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTg               >  1:3231337/1‑71 (MQ=255)
                                                  atgaatgaAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCa                               <  1:757557/54‑1 (MQ=255)
                                                  atgaatgaAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAg                              >  1:3676027/1‑55 (MQ=255)
                                                  atgaatgaAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGc              <  1:1563068/71‑1 (MQ=255)
                                                   tgaatgaaGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCt             <  1:1235152/71‑1 (MQ=255)
                                                     aatgaaGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAg                              >  1:1358576/1‑52 (MQ=255)
                                                      atgaaGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGAACATTGCTAt           >  1:602004/1‑70 (MQ=255)
                                                      atgaaGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAatta                                     <  1:2802961/44‑1 (MQ=38)
                                                       tgaaGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAg                              <  1:24334/50‑1 (MQ=39)
                                                        gaaGACGACTACAAAATAAGAAGAGGAAACGcagcag                                          <  1:470995/37‑1 (MQ=255)
                                                         aaGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTaa        >  1:133170/1‑70 (MQ=255)
                                                          cgacgacTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGAt                          <  1:567610/50‑1 (MQ=38)
                                                          cgacgacTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTc                                <  1:2962830/44‑1 (MQ=38)
                                                          cgacgacTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATAc                        <  1:1144971/52‑1 (MQ=255)
                                                          aGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAgg                             <  1:814196/48‑1 (MQ=39)
                                                          aGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGATACg                       <  1:3830853/54‑1 (MQ=255)
                                                          aGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAattatt                                   <  1:1573422/42‑1 (MQ=255)
                                                           gacgacTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAatat     >  1:121127/1‑71 (MQ=255)
                                                           gacgacTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAatat     >  1:3281020/1‑71 (MQ=255)
                                                           gacgacTACAAAATAAGAAGAGGGAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAatat     >  1:492523/1‑71 (MQ=255)
                                                            acgacTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTc                                <  1:2079541/43‑1 (MQ=38)
                                                            acgacTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGa                           >  1:204121/1‑48 (MQ=255)
                                                            acgacTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAatat     >  1:1630863/1‑70 (MQ=255)
                                                             cgacTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATAttt   >  1:125265/1‑71 (MQ=255)
                                                             cgacTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCa                               <  1:1548989/43‑1 (MQ=38)
                                                              gacTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATAtttt  >  1:1274865/1‑71 (MQ=255)
                                                              gacTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATAtttt  <  1:3392271/71‑1 (MQ=255)
                                                               acTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATAtttt  <  1:2097155/70‑1 (MQ=255)
                                                                 tACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCa                               >  1:863718/1‑39 (MQ=255)
                                                                  |                                                                  
GAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTT  >  minE/217964‑218096

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: