Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,697,632 A→G 17.6% Y195C (TAT→TGT)  yfhR → predicted peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,697,6320AG17.6% 111.7 / 20.0 51Y195C (TAT→TGT) yfhRpredicted peptidase
Reads supporting (aligned to +/- strand):  ref base A (20/22);  new base G (4/5);  total (24/27)
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

CTGGATGTTATTGGTCGGGGTGATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATCGCCAACCAAATGATCCCCGGCAGTGGCTACTTACTTGATGAGAGTTACAGCGGCGAAAAT  >  minE/1697562‑1697702
                                                                      |                                                                      
cTGGATGTTATTGGTCGGGGTGATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTa                                                                        >  1:1155465/1‑71 (MQ=255)
cTGGATGTTATTGGTCGGGGTGATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTa                                                                        >  1:1958783/1‑71 (MQ=255)
cTGGATGTTATTGGTCGGGCTGATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTa                                                                        >  1:1952377/1‑71 (MQ=255)
  ggATGTTATTGGTCGGGGTGATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTAt                                                                       >  1:1899678/1‑70 (MQ=255)
          ttGGTCGGGGTGATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTGTGCAAcc                                                                 <  1:1421448/68‑1 (MQ=255)
                ggggTGATCGTGAATGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATCGccaac                                                        >  1:2016473/1‑71 (MQ=255)
                ggggTGATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATCGccaac                                                        >  1:520267/1‑71 (MQ=255)
                  ggTGATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATCGccaacca                                                      <  1:1947677/71‑1 (MQ=255)
                  ggTGATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATCGccaacca                                                      <  1:1723325/71‑1 (MQ=255)
                    tGATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATCGCCAACCaaa                                                    <  1:2111981/71‑1 (MQ=255)
                      aTCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATCGccaaccaa                                                     <  1:1046107/68‑1 (MQ=255)
                      aTCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATCGccaaccaa                                                     <  1:476003/68‑1 (MQ=255)
                      aTCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATCGCCAACCAAATg                                                  <  1:1052318/71‑1 (MQ=255)
                       tCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTGTGCAACCATCGCCAACCaaa                                                    <  1:2130633/68‑1 (MQ=255)
                       tCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGc                                                                     <  1:2180468/51‑1 (MQ=255)
                       tCGTGAAAGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTGTGCAACCATCGCCAACCaaa                                                    <  1:230904/68‑1 (MQ=255)
                        cGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATCGCCAACCAAATGAt                                                >  1:1310970/1‑71 (MQ=255)
                            aaGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCAt                                                               >  1:1565124/1‑52 (MQ=255)
                            aaGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCAt                                                               >  1:714537/1‑52 (MQ=255)
                             aGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCAt                                                               >  1:1246717/1‑51 (MQ=255)
                             aGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCAt                                                               >  1:847013/1‑51 (MQ=255)
                              ggCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATCGCCAACCAAATGATCCCCgg                                          <  1:1005538/71‑1 (MQ=255)
                               gCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATCGCCAACCAAATGATCCCCGGc                                         <  1:879349/71‑1 (MQ=255)
                                 aTACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATCGCCAACCAAATGATCCCCGGCAg                                       <  1:745853/71‑1 (MQ=255)
                                 aTACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATCGCCAACCAAATGATCCCCGGCAg                                       <  1:1258869/71‑1 (MQ=255)
                                  tACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCAt                                                               <  1:1542108/46‑1 (MQ=255)
                                    cGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATCGCCAACCAAATGATCCCCGGCAg                                       >  1:250138/1‑68 (MQ=255)
                                    cGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATCGCCAACCAAATGATCCCCGGCAGTgg                                    >  1:1869539/1‑71 (MQ=255)
                                    cGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATCGCCAACCAAATGATCCCCGGCAGTgg                                    >  1:1866996/1‑71 (MQ=255)
                                    cGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATAGccaacca                                                      >  1:900866/1‑53 (MQ=255)
                                         ggTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATCGccaacca                                                      <  1:1254246/48‑1 (MQ=255)
                                         ggTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATCGccaacca                                                      <  1:746517/48‑1 (MQ=255)
                                         ggTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATCGccaacca                                                      <  1:892667/48‑1 (MQ=255)
                                          gTGATCCTCGACTCCACATTTGCCTCTTGTGCAACCATCGCCAACCAAATGATCCCCGGCa                                        >  1:460587/1‑61 (MQ=255)
                                           tGATCCTCGACTCCACATTTGCCTCTTGTGCAACCATCGCCAACCAAATGATCCCCGGCAGTGGCtactta                             >  1:1821073/1‑71 (MQ=255)
                                              tCCTCGACTCCACATTTGCCTCTTATGCAACCATCGccaacca                                                      <  1:1450717/43‑1 (MQ=255)
                                              tCCTCGACTCCACATTTGCCTCTTATGCAACCATCGCCAACCAAATGATCCCCGGCAGTGGCTacttactt                          <  1:1593014/71‑1 (MQ=255)
                                                cTCGACTCCACATTTGCCTCTTGTGCAACCATCGccaaccaa                                                     <  1:85092/42‑1 (MQ=255)
                                                  cGACTCCACATTTGCCTCTTATGCAACCATCGCCAACCAAATGATCCCCGGCAGTGGCTACTTACTtgat                       >  1:1155572/1‑70 (MQ=255)
                                                   gACTCCACATTTGCCTCTTGTGCAACCATCGCCAACCAAATGATCCCCGGCAGTGGCTACTTACTtgat                       >  1:265719/1‑69 (MQ=255)
                                                   gACTCCACATTTGCCTCTTGTGCAACCATCGCCAACCAAATGATCCCCGGCAGTGGCTACTTACTtgat                       >  1:462372/1‑69 (MQ=255)
                                                   gACTCCACATTTGCCTCTTATGCAACCATCGCCAACCAAATGATCCCCGGCAGTGGCTACTTACTtgatga                     <  1:231959/71‑1 (MQ=255)
                                                   gACTCCACATTTGCCTCTTATGCAACCATCGCCAACCAAATGATCCCCGGCAGTGGCTACTTACTtgatga                     <  1:1065747/71‑1 (MQ=255)
                                                   gACTCCACATTTGCCTCTTATGAACCATCGCCAACCAAATGATCCCCGGCAGTGGCTACTTACTtgatga                     <  1:1340543/70‑1 (MQ=255)
                                                    aCTCCACATTTGCCTCTTATGCAACCATCGCCAACCAAATGATCCCCGGCAGTg                                     <  1:962116/54‑1 (MQ=255)
                                                      tCCACATTTGCCTCTTGTGCAACCATCGCCAACCAAATg                                                  <  1:307240/39‑1 (MQ=255)
                                                      tCCACATTTGCCTCTTATGCAACCATCGCCAACCAAATGATCCCCGGCAGTGGCTACTTACTtgatg                      >  1:1744803/1‑67 (MQ=255)
                                                        cacaTTTGCCTCTTATGCAACCATCGCCAACCAAATGATCCCCGGCAGTGGCtacttac                            >  1:780131/1‑59 (MQ=255)
                                                        cacaTTTGCCTCTTATGCAACCATCGCCAACCAAATGATCCCCGGCAGTGGCTacttactt                          >  1:545347/1‑61 (MQ=255)
                                                         acaTTTGCCTCTTATGCAACCATCGCCAACCAAATGATCCCCGGCAGTGGCTACTTACTTGATGAGAGTTa               <  1:284642/71‑1 (MQ=255)
                                                            tttGCCTCTTATGCAACCATCGCCAACCAAATGATCCCCgg                                          >  1:1855531/1‑41 (MQ=255)
                                                            tttGCCTCTTATGCAACCATCGCCAACCAAATGATCCCCgg                                          >  1:1556266/1‑41 (MQ=255)
                                                                  tctTATGCAACCATCGCCAACCAAATGATCCCCGGCAGTGGCTACTTACTTGATgaga                   >  1:431660/1‑58 (MQ=255)
                                                                   ctTATGCAACCATCGCCAACCAAATGATCCCCGGCAGTGGCTACTTACTTGATGAGAGTTACAgcggcg       >  1:75044/1‑69 (MQ=255)
                                                                      aTGCAACCATCGCCAACCAAATGATCCCCGGCAGTGGCTACTTACTtga                        >  1:689597/1‑49 (MQ=255)
                                                                      aTGCAACCATCGCCAACCAAATGATCCCCGGCAGTGGCTACTTACTTGATGAGAGTTACAGCGGCGaaaa   >  1:1908281/1‑70 (MQ=255)
                                                                      aTGCAACCATCGCCAACCAAATGATCCCCGGCAGTGGCTACTTACTTGATGAGAGTTACAGCGGCGAAAAt  >  1:1580448/1‑71 (MQ=255)
                                                                      |                                                                      
CTGGATGTTATTGGTCGGGGTGATCGTGAAGGCATACGTGCGGTGATCCTCGACTCCACATTTGCCTCTTATGCAACCATCGCCAACCAAATGATCCCCGGCAGTGGCTACTTACTTGATGAGAGTTACAGCGGCGAAAAT  >  minE/1697562‑1697702

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: