Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,216,248 A→G 9.2% V82A (GTA→GCA)  htpX ← predicted endopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,216,2480AG 9.2% 235.6 / 15.0 87V82A (GTA→GCA) htpXpredicted endopeptidase
Reads supporting (aligned to +/- strand):  ref base A (35/44);  new base G (5/3);  total (40/47)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

GTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGT  >  minE/1216178‑1216316
                                                                      |                                                                    
gTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTa                                                                      <  1:1043001/71‑1 (MQ=255)
 ttGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTAc                                                                     <  1:1715604/71‑1 (MQ=255)
  tGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCACCGCCTGACGAGCCTGGGTTGCTa                                                                      >  1:1050058/1‑69 (MQ=255)
   gATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTa                                                                      >  1:3597036/1‑68 (MQ=255)
   gATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAg                                                                   >  1:3577732/1‑71 (MQ=255)
    aTGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTGCAg                                                                   >  1:569783/1‑70 (MQ=255)
    aTGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAg                                                                   <  1:2904962/70‑1 (MQ=255)
    aTGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAg                                                                   >  1:2320108/1‑70 (MQ=255)
     tGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTa                                                                 <  1:933740/71‑1 (MQ=255)
     tGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTa                                                                 <  1:881022/71‑1 (MQ=255)
       tCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTAtt                                                               >  1:1924203/1‑71 (MQ=255)
       tCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTAt                                                                <  1:1655986/70‑1 (MQ=255)
         cGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGa                                                             >  1:316028/1‑71 (MQ=255)
                tGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGcc                                                       >  1:3771716/1‑70 (MQ=255)
                tGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGcc                                                       >  1:1763577/1‑70 (MQ=255)
                  gTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAAc                                                    <  1:2634184/71‑1 (MQ=255)
                  gTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAAc                                                    <  1:1648936/71‑1 (MQ=255)
                     gATAGCCACTTGCGGCATAGCGATACCCGCCTGACGAGACTGGGTTGCTGCAGTATTGACCAGCCAACGtt                                                 <  1:1375230/71‑1 (MQ=255)
                       tAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTcc                                               <  1:3858558/71‑1 (MQ=255)
                        aGCCACTTGCGGCATAGCGATCCCCGCCTGAGGAGCCTGGGTTGCTGCAGTATTGACCAGACAACGTTccc                                              >  1:3861846/1‑71 (MQ=255)
                        aGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGttgc                                               >  1:1427232/1‑68 (MQ=255)
                          ccACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCtt                                            >  1:1142474/1‑71 (MQ=255)
                          ccACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCtt                                            <  1:2956723/71‑1 (MQ=255)
                          ccACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCtt                                            >  1:255778/1‑71 (MQ=255)
                          ccACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCtt                                            <  1:1661280/71‑1 (MQ=255)
                            aCTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGcc                                                       <  1:128059/58‑1 (MQ=255)
                             cTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCg                                         <  1:78847/71‑1 (MQ=255)
                             cTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCg                                         <  1:3080393/71‑1 (MQ=255)
                               tGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTGCAGTATTGACCAGCCAACGTTCCCTTTCGt                                        >  1:2408089/1‑70 (MQ=255)
                               tGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGa                                                             >  1:1837314/1‑49 (MQ=255)
                               tGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGa                                                             >  1:3683551/1‑49 (MQ=255)
                               tGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTc                                          <  1:137898/68‑1 (MQ=255)
                               tGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGt                                        >  1:3034303/1‑70 (MQ=255)
                                gCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGt                                        <  1:1050710/69‑1 (MQ=255)
                                gCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTa                                      >  1:3643362/1‑71 (MQ=255)
                                 cGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTa                                      <  1:3066241/70‑1 (MQ=255)
                                    cATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGAc                                                            <  1:653426/45‑1 (MQ=255)
                                     aTAGCGATCCCCGCCTGACGAGCCTGGGTTGCTGCAGTATTGACCAGCCAACGTTCCCtt                                            <  1:1820895/60‑1 (MQ=255)
                                     aTAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCgg                                 <  1:2839157/71‑1 (MQ=255)
                                     aTAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCgg                                 <  1:3628680/71‑1 (MQ=255)
                                      tAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCg                                         <  1:679923/62‑1 (MQ=255)
                                      tAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCgg                                 <  1:61710/70‑1 (MQ=255)
                                       aGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAg                                                         >  1:1286579/1‑45 (MQ=255)
                                       aGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCtt                                            <  1:2046493/58‑1 (MQ=255)
                                       aGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGt                                        <  1:2310627/62‑1 (MQ=255)
                                       aGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGtt                               >  1:1527507/1‑71 (MQ=255)
                                        gCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTg                              >  1:4026907/1‑71 (MQ=255)
                                        gCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTg                              <  1:762208/71‑1 (MQ=255)
                                         cGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTc                                          >  1:2361595/1‑58 (MQ=255)
                                         cGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGc                             <  1:3266132/71‑1 (MQ=255)
                                          gATCCCCGCCTGACGAGCCTGGGTTGCTGCAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCt                            >  1:805567/1‑71 (MQ=255)
                                          gATCCCCGCCTGACGAGCCTGGGTTGCTGCAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCt                            >  1:1206135/1‑71 (MQ=255)
                                          gATCCCCGCCTGACGAGCCTGGGTTGCTGCAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCt                            <  1:1304298/71‑1 (MQ=255)
                                            tCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTc                           >  1:2430558/1‑70 (MQ=255)
                                            tCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTc                           >  1:1321126/1‑70 (MQ=255)
                                             ccccGCCTGACGAGCCTGGGTTGCTACAGTATTGAcc                                                           >  1:79393/1‑37 (MQ=255)
                                              cccGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGa                         <  1:501076/70‑1 (MQ=255)
                                                   cTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACtt                   <  1:3784175/71‑1 (MQ=255)
                                                   cTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACtt                   <  1:3736737/71‑1 (MQ=255)
                                                   cTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACtt                   <  1:1029672/71‑1 (MQ=255)
                                                    tGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGa                         >  1:3207193/1‑64 (MQ=255)
                                                    tGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGa                         >  1:2715208/1‑64 (MQ=255)
                                                    tGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGa                         >  1:3979534/1‑64 (MQ=255)
                                                      aCGAGCCTGGGTTGCTACAGTATTGACCAGCCAACg                                                   >  1:309351/1‑36 (MQ=255)
                                                      aCGAGCCTGGGTTGCTACAGTATTGACCAGCCAACg                                                   >  1:1192699/1‑36 (MQ=255)
                                                      aCGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTATGCGGTTGCTCGATCACTTccc                >  1:3333020/1‑71 (MQ=255)
                                                      aCGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTcc                 <  1:2659770/70‑1 (MQ=255)
                                                       cGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTcccc               <  1:1552236/71‑1 (MQ=255)
                                                        gAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTcc                                               >  1:1820752/1‑38 (MQ=255)
                                                        gAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTcc                                               >  1:1444367/1‑38 (MQ=255)
                                                        gAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCg              >  1:1721802/1‑71 (MQ=255)
                                                        gAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCg              >  1:91150/1‑71 (MQ=255)
                                                        gAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCg              >  1:147271/1‑71 (MQ=255)
                                                            cTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTcc                 >  1:563270/1‑64 (MQ=255)
                                                            cTGGGTTGCTACAGAATTGACCAGCCAACGTTCCCTTTCGTTACGCGCTTGCTCGATCACtt                   >  1:855197/1‑62 (MQ=255)
                                                             tGGGTTGCTACAGTATTGCCCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCaa          >  1:412383/1‑70 (MQ=255)
                                                             tGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGc                             <  1:3596022/51‑1 (MQ=255)
                                                             tGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACt                    >  1:2570874/1‑60 (MQ=255)
                                                               ggTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAg       <  1:3679982/71‑1 (MQ=255)
                                                                    cTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTc                           <  1:1182270/46‑1 (MQ=255)
                                                                    cTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTc                           <  1:646467/46‑1 (MQ=255)
                                                                    cTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTc                           <  1:1698785/46‑1 (MQ=255)
                                                                    cTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGt  <  1:3305404/71‑1 (MQ=255)
                                                                    cTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGt  <  1:1004034/71‑1 (MQ=255)
                                                                    cTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGt  <  1:200160/71‑1 (MQ=255)
                                                                    cTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGt  <  1:4008315/71‑1 (MQ=255)
                                                                    cTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGt  <  1:2243069/71‑1 (MQ=255)
                                                                    cTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGt  <  1:2623415/71‑1 (MQ=255)
                                                                    cTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGt  <  1:3256922/71‑1 (MQ=255)
                                                                     tACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGAt                        <  1:2020761/48‑1 (MQ=255)
                                                                      aCAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGAt                        >  1:1324859/1‑47 (MQ=255)
                                                                      aCAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGAt                        <  1:2203720/47‑1 (MQ=255)
                                                                      aCAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGAt                        >  1:3755423/1‑47 (MQ=255)
                                                                      |                                                                    
GTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCATAGCGATCCCCGCCTGACGAGCCTGGGTTGCTACAGTATTGACCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGT  >  minE/1216178‑1216316

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: ‑