Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 894,873 C→T 5.5% E398K (GAA→AAA)  sapA ← predicted antimicrobial peptide transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE894,8730CT 5.5% 243.3 / 2.7 72E398K (GAA→AAA) sapApredicted antimicrobial peptide transporter subunit
Reads supporting (aligned to +/- strand):  ref base C (32/36);  new base T (2/2);  total (34/38)
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.93e-01

GAAAGCGACCTTCTACCGGCACAATCACCACTTTTACGCCAACCTGCGCCATATCCGCCTGAATCAGTTCGGCAGTTTTCAGTGGACTGGGGTTCCACGCCTGCGAACGTGTGGGCACCCACAGTTTCAGCGTTAAATTT  >  minE/894804‑894943
                                                                     |                                                                      
gAAAGCGACCTTCTACCGGCACAATCACCACTTTTACGCCAACCTGCGCCATATCCGCCTGAATCAGTTCg                                                                       >  1:1407785/1‑71 (MQ=255)
 aaaGCGACCTTCTACCGGCACAATCACCACTTTTACGCCAACCTGCGCCATATCCGCCTGAATCAGTTCgg                                                                      >  1:1078706/1‑71 (MQ=255)
 aaaGCGACCTTCTACCGGCACAATCACCACTTTTACGCCAACCTGCGCCATATCCGCCTGAATCAGTTCgg                                                                      <  1:891540/71‑1 (MQ=255)
 aaaGCGACCTTCTACCGGCACAATCACCACTTTTACGCCAACCTGCGCCATATCCGCCTGAATCAGTTCgg                                                                      >  1:864189/1‑71 (MQ=255)
 aaaGCGACCTTCTACCGGCACAATCACCACTTTTACGCCAACCTGCGCCATATCCGCCTGAATCAGTTCgg                                                                      >  1:972216/1‑71 (MQ=255)
      gACCTTCTACCGGCACAATCACCACTTTTACGCCAACCTGCGCCATATCCGCCTGAATCAGTTCGGCAGtt                                                                 <  1:38721/71‑1 (MQ=255)
       aCCTTCTACCGGCACAATCACCACTTTTACGCCAACCTGCGCCATATCCGCCTGAATCAGTTCGGCAGtt                                                                 <  1:2138786/70‑1 (MQ=255)
        ccTTCTACCGGCACAATCACCACTTTTACGCCAACCTGCGCCATATCCGCCTGAATCAGTTTGGCAGtttt                                                               >  1:1312544/1‑71 (MQ=255)
        ccTTCTACCGGCACAATCACCACTTTTACGCCAACCTGCGCCATATCCGCCTGAATCAGTTCGGCAGtttt                                                               <  1:2133403/71‑1 (MQ=255)
         cTTCTACCGGCACAATCACCACTTTTACGCCAACCTGCGCCATATCCGCCTGAATCAGTTCGGCAGTTTTc                                                              <  1:2200253/71‑1 (MQ=255)
         cTTCTACCGGCACAATCACCACTTTTACGCCAACCTGCGCCATATCCGCCTGAATCAGTTCGGCAGTTTTc                                                              <  1:2586021/71‑1 (MQ=255)
         cTTCTACCGGCACAATCACCACTTTTACGCCAACCTGCGCCATATCCGCCTGAATCAGTTCGGCAGTTTTc                                                              <  1:2682694/71‑1 (MQ=255)
         cTTCTACCGGCACAATCACCACTTTTACGCCAACCTGCGCCATATCCGCCTGAATCAGTTCGGCAGTTTTc                                                              <  1:1752008/71‑1 (MQ=255)
           tCTACCGGCACAATCACCACTTTTACGCCAACCTGCGCCATATCCGCCTGAATCAGTTCGGCAGtt                                                                 <  1:3601928/66‑1 (MQ=255)
           tCTACCGGCACAATCACCACTTTTACGCCAACCTGCGCCATATCCGCCTGAATCAGTTCGGCAGtt                                                                 <  1:1831274/66‑1 (MQ=255)
            cTACCGGCACAATCACCACTTTTACGCCAACCTGCGCCATATCCGCCTGAATCAGTTCGGCAGTTTTCa                                                             <  1:1596790/69‑1 (MQ=255)
            cTACCGGCACAATCACCACTTTTACGCCAACCTGCGCCATATCCGCCTGAATCAGTTCGGCAGTTTTCa                                                             <  1:3253237/69‑1 (MQ=255)
             tACCGGCACAATCACCACTTTTACGCCAACCTGCGCCATATCCGCCTGAATCAGTTCGGCAGtt                                                                 <  1:2562322/64‑1 (MQ=255)
             tACCGGCACAATCACCACTTTTACGCCAACCTGCGCCATATCCGCCTGAATCAGTTCGGCAGTTTTCa                                                             <  1:1798989/68‑1 (MQ=255)
             tACCGGCACAATCACCACTTTTACGCCAACCTGCGCCATATCCGCCTGAATCAGTTCGGCAGTTTTCAGTg                                                          >  1:3671287/1‑71 (MQ=255)
               ccGGCACAATCACCACTTTTACGCCAACCTGCGCCATATCCGCCTGAATCAGTTCGGCAGTTTTCAGTg                                                          <  1:2185743/69‑1 (MQ=255)
                    acaATCACCACTTTTACGCCAACCTGCGCCATATCCGCCTGAATCAGTTCGGCAGTTTTCAGTGGACTggg                                                   >  1:1712019/1‑71 (MQ=255)
                       ctCACCACTTTTACGCCAACCTGCGCCATATCCGCCTGAATCAGTTCGGCAGTTTTCAg                                                            <  1:1377424/58‑1 (MQ=255)
                       aTCACCACTTTTACGCCAACCTGCGCCATATCCGCCTGAATCAGTTCGGCAGTTTTCAg                                                            <  1:2115434/59‑1 (MQ=255)
                       aTCACCACTTTTACGCCAACCTGCGCCATATCCGCCTGAATCAGTTCGGCAGTTTTCAg                                                            <  1:3313818/59‑1 (MQ=255)
                       aTCACCACTTTTACGCCAACCTGCGCCATATCCGCCTGAATCAGTTCGGCAGTTTTCAg                                                            <  1:2799399/59‑1 (MQ=255)
                       aTCACCACTTTTACGCCAACCTGCGCCATATCCGCCTGAATCAGTTCGGCAGTTTTCAGt                                                           >  1:3130319/1‑60 (MQ=255)
                       aTCACCACTTTTACGCCAACCTGCGCCATATCCGCCTGAATCAGTTCGGCAGTTTTCAGt                                                           >  1:13864/1‑60 (MQ=255)
                       aTCACCACTTTTACGCCAACCTGCGCCATATCCGCCTGAATCAGTTCGGCAGTTTTCAGt                                                           >  1:1302124/1‑60 (MQ=255)
                       aTCACCACTTTTACGCCAACCTGCGCCATATCCGCCTGAATCAGTTCGGCAGTTTTCAGt                                                           >  1:1766184/1‑60 (MQ=255)
                       aTCACCACTTTTACGCCAACCTGCGCCATATCCGCCTGAATCAGTTCGGCAGTTTTCAGt                                                           >  1:1635194/1‑60 (MQ=255)
                        tCACCACTTTTACGCCAACCTGCGCCATATCCGCCTGAATCAGTTCGGCAGtttt                                                               <  1:2437188/55‑1 (MQ=255)
                        tCACCACTTTTACGCCAACCTGCGCCATATCCGCCTGAATCAGTTCGGCAGtttt                                                               <  1:920641/55‑1 (MQ=255)
                        tCACCACTTTTACGCCAACCTGCGCCATATCCGCCTGAATCAGTTCGGCAGtttt                                                               <  1:832555/55‑1 (MQ=255)
                        tCACCACTTTTACGCCAACCTGCGCCATATCCGCCTGAATCAGTTCGGCAGtttt                                                               <  1:2002974/55‑1 (MQ=255)
                         caccacTTTTACGCCAACCTGCGCCATATCCGCCTGAATCAGTTCGGCAGTTTTCAg                                                            <  1:2113932/57‑1 (MQ=255)
                         caccacTTTTACGCCAACCTGCGCCATATCCGCCTGAATCAGTTCGGCAGTTTTCAGTGGACTGGGGtt                                                >  1:1779428/1‑69 (MQ=255)
                          accacTTTTACGCCAACCTGCGCCATATCCGCCTGAATCAGTTCGGCAGTTTTCAGTGGACTGGGGTTcc                                              <  1:746007/70‑1 (MQ=255)
                             acTTTTACGCCAACCTGCGCCATATCCGCCTGAATCAGTTTGGCAGTTTTCAGTGGACTGGGGTTCCACGc                                          >  1:1299016/1‑71 (MQ=255)
                              cTTTTACGCCAACCTGCGCCATATCCGCCTGAATCAGTTCGGCAGtttt                                                               <  1:3471637/49‑1 (MQ=255)
                              cTTTTACGCCAACCTGCGCCATATCCGCCTGAATCAGTTCGGCAGtttt                                                               <  1:3564933/49‑1 (MQ=255)
                              cTTTTACGCCAACCTGCGCCATATCCGCCTGAATCAGTTCGGCAGtttt                                                               <  1:3603868/49‑1 (MQ=255)
                              cTTTTACGCCAACCTGCGCCATATCCGCCTGAATCAGTTCGGCAGtttt                                                               <  1:878177/49‑1 (MQ=255)
                              cTTTTACGCCAACCTGCGCCATATCCGCCTGAATCAGTTCGGCAGtttt                                                               <  1:2149249/49‑1 (MQ=255)
                              cTTTTACGCCAACCTGCGCCATATCCGCCTGAATCAGTTCGGCAGTTTTCAGTGGACTGGGGTTCCACGcc                                         >  1:2402544/1‑71 (MQ=255)
                               ttttACGCCAACCTGCGCCATATCCGCCTGAATCAGTTCGGCAGTTTTCa                                                             <  1:1949774/50‑1 (MQ=255)
                               ttttACGCCAACCTGCGCCATATCCGCCTGAATCAGTTCGGCAGTTTTCAGTGGACTGGGGTTCCACGCCt                                        >  1:1153066/1‑71 (MQ=255)
                                 ttACGCCAACCTGCGCCATATCCGCCTGAATCAGTTCGGCAGt                                                                  >  1:2976899/1‑43 (MQ=255)
                                 ttACGCCAACCTGCGCCATATCCGCCTGAATCAGTTCGGCAGTTTTCAGTGGACTGGGGTTCCACGCCTGc                                      >  1:731686/1‑71 (MQ=255)
                                 ttACGCCAACCTGCGCCATATCCGCCTGAATCAGTTCGGCAGTTTTCAGTGGACTGGGGTTCCACGCCTGc                                      >  1:423516/1‑71 (MQ=255)
                                 ttACGCCAACCTGCGCCATATCCGCCTGAATCAGTTCGGCAGTTTTCAGTGGACTGGGGTTCCACGCCTGc                                      >  1:3463486/1‑71 (MQ=255)
                                 ttACGCCAACCTGCGCCATATCCGCCTGAATCAGTTCGGCAGTTTTCAGTGGACTGGGGTTCCACGCCTGc                                      >  1:1516031/1‑71 (MQ=255)
                                 ttACGCCAACCTGCGCCATATCCGCCTGAATCAGTTCGGCAGTTTTCAGTGGACTGGGGTTCCACGCCTGc                                      >  1:1611831/1‑71 (MQ=255)
                                  tACGCCAACCTGCGCCATATCCGCCTGAATCAGTTCGGCAGtttt                                                               <  1:956050/45‑1 (MQ=255)
                                  tACGCCAACCTGCGCCATATCCGCCTGAATCAGTTCGGCAGTTTTCa                                                             >  1:1084059/1‑47 (MQ=255)
                                   aCGCCAACCTGCGCCATATCCGCCTGAATCAGTTTGGCAGTTTTCAGTGGACTGGGGTTCCACGCCTGCg                                     <  1:1042125/70‑1 (MQ=255)
                                   aCGCCAACCTGCGCCATATCCGCCTGAATCAGTTCGGCAGTTTTCAGTGGACTGGGGTTCCACGCCTGCg                                     >  1:230448/1‑70 (MQ=255)
                                   aCGCCAACCTGCGCCATATCCGCCTGAATCAGTTCGGCAGTTTTCAGTGGACTGGGGTTCCACGCCTGCg                                     >  1:1394646/1‑70 (MQ=255)
                                     gCCAACCTGCGCCATATCCGCCTGAATCAGTTCGGCAGTTTTCAGTGGACTGGGGTTCCACGCCTGCGAAc                                  <  1:2888965/71‑1 (MQ=255)
                                        aaCCTGCGCCATATCCGCCTGAATCAGTTCGGCAGTTTTCAGTGGACTGGGGTTCCACGCCTGCGAACgtg                               <  1:436741/71‑1 (MQ=255)
                                        aaCCTGCGCCATATCCGCCTGAATCAGTTCGGCAGTCTTCAGTGGACTGGGGTTCCACGCCTGCGAACgtg                               >  1:358924/1‑71 (MQ=255)
                                         aCCTGCGCCATATCCGCCTGAATCAGTTCGGCAGTTTTCAGTGGACTGGGGTTCCACGCCTGCGAACgtgt                              >  1:1220175/1‑71 (MQ=255)
                                           cTGCGCCATATCCGCCTGAATCAGTTCGGCAGTTTTCAGTGGACTGGGGTTCCACGCCTGCGAACGTGTgg                            >  1:2980539/1‑71 (MQ=255)
                                              cgcCATATCCGCCTGAATCAGTTCGGCAGTTTTCAGTGGACTGGGGTTCCACGCCTGCGAACGTGTGGGCa                         >  1:627020/1‑71 (MQ=255)
                                                ccATATCCGCCTGAATCAGTTTGGCAGTTTTCAGTGGACTGGGGTTCCACGCCTGCGAACGTGTGGGCAc                        <  1:565219/70‑1 (MQ=255)
                                                    atCCGCCTGAATCAGTTCGGCAGTTTTCAGTGGACTGGGGTTCCACGCCTGCGAACgtg                               <  1:1063284/59‑1 (MQ=255)
                                                     tCCGCCTGAATCAGTTCGGCAGTTTTCAGTGGACTGGGGTTCCACGCCTGCGAACGTGTGGGCACCCACAg                  >  1:2414209/1‑71 (MQ=255)
                                                     tCCGCCTGAATCAGTTCGGCAGTTTTCAGTGGACTGGGGTTCCACGCCTGCGAACGTGTGGGCACCCACAg                  >  1:214128/1‑71 (MQ=255)
                                                       cGCCTGAATCAGTTCGGCAGTTTTCAGTGGACTGGGGTTCCACGCCTGCGAACGTGTgg                            >  1:3526737/1‑59 (MQ=255)
                                                       cGCCTGAATCAGTTCGGCAGTTTTCAGTGGACTGGGGTTCCACGCCTGCGAACGTGTgg                            >  1:984247/1‑59 (MQ=255)
                                                       cGCCTGAATCAGTCCGGCAGTTTTCAGTGGACTGGGGTTCCACGCCTGCGAACGTGTGGGCACCCACAGtt                >  1:211865/1‑71 (MQ=255)
                                                                 aGTTCGGCAGTTTTCAGTGGACTGGGGTTCCACGCCTGCGAACGTGTGGGCACCCACAGTTTCAGCGTTa       <  1:2379289/70‑1 (MQ=255)
                                                                     cGGCAGTTTTCAGTGGACTGGGGTTCCACGCCTGCGAACGTGTGGGCACCCACAGTTTCAGCGTTAAAttt  >  1:261688/1‑71 (MQ=255)
                                                                     |                                                                      
GAAAGCGACCTTCTACCGGCACAATCACCACTTTTACGCCAACCTGCGCCATATCCGCCTGAATCAGTTCGGCAGTTTTCAGTGGACTGGGGTTCCACGCCTGCGAACGTGTGGGCACCCACAGTTTCAGCGTTAAATTT  >  minE/894804‑894943

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: