Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,000,186 G→A 43.6% G179D (GGT→GAT)  ydgD → predicted peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,000,1860GA43.6% 16.1 / 47.0 39G179D (GGT→GAT) ydgDpredicted peptidase
Reads supporting (aligned to +/- strand):  ref base G (11/11);  new base A (8/9);  total (19/20)
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.99e-01

TACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTG  >  minE/1000124‑1000254
                                                              |                                                                    
tACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAg                                                              >  1:1270537/1‑71 (MQ=255)
 aCGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGt                                                             >  1:243246/1‑71 (MQ=255)
 aCGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGt                                                             >  1:1186220/1‑71 (MQ=255)
  cGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGt                                                             >  1:1046773/1‑70 (MQ=255)
  cGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGt                                                             <  1:967119/70‑1 (MQ=255)
  cGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGt                                                             <  1:12396/70‑1 (MQ=255)
     cgttgccgttATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGt                                                                  >  1:447974/1‑62 (MQ=255)
         gccgttATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTcagg                                                     >  1:1919279/1‑71 (MQ=255)
         gccgttATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTcagg                                                     >  1:150109/1‑71 (MQ=255)
           cgttATTTGAGGGAGATAAAGCCGCACTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTcagg                                                     <  1:696163/69‑1 (MQ=255)
               aTTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTaaa                                                               >  1:710732/1‑55 (MQ=255)
                tttGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGTGACTCAGGCAGGCt                                               <  1:1692717/70‑1 (MQ=255)
                     gggAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTg                                         <  1:2197076/71‑1 (MQ=255)
                          aTAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGTGACTCAGGCAGGCTACCCTg                                         <  1:1142170/66‑1 (MQ=255)
                          aTAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGTGACTCAGGCAGGCTACCCTg                                         <  1:873216/66‑1 (MQ=255)
                          aTAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGTGACTCAGGCAGGCTACCCTg                                         <  1:1139666/66‑1 (MQ=255)
                          aTAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGTGACTCAGGCAGGCTACCCTGAAGa                                     <  1:76860/70‑1 (MQ=255)
                           tAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTa                                                                 <  1:625698/41‑1 (MQ=255)
                             aaGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGt                                                                  >  1:443191/1‑38 (MQ=255)
                                       tACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGt                       >  1:1408210/1‑71 (MQ=255)
                                        aCCGCCGCATTAAAAGCGGCAGATCGTAAAGTGACTCAGGCAGGCTACCCTGAAGatcatc                                >  1:937225/1‑61 (MQ=255)
                                             cgcATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACa                 <  1:1308357/71‑1 (MQ=255)
                                                 ttAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCAtctc                              <  1:1258127/54‑1 (MQ=255)
                                                 ttAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCAtctc                              <  1:693027/54‑1 (MQ=255)
                                                 ttAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCAtctc                              <  1:513329/54‑1 (MQ=255)
                                                 ttAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCAtctc                              <  1:350799/54‑1 (MQ=255)
                                                 ttAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGt                       <  1:841917/61‑1 (MQ=255)
                                                  tAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGtcat            >  1:2048686/1‑71 (MQ=255)
                                                  tAAAAGCGGCAGATCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGtcat            >  1:1841889/1‑71 (MQ=255)
                                                   aaaaGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGa                                     >  1:1129324/1‑45 (MQ=255)
                                                   aaaaGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTg                     >  1:1130554/1‑61 (MQ=255)
                                                   aaaaGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGtcat            >  1:1063910/1‑70 (MQ=255)
                                                   aaaaGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGtc              >  1:1670967/1‑68 (MQ=255)
                                                    aaaGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGtcatc           <  1:1305553/70‑1 (MQ=255)
                                                    aaaGCGGCAGATCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGtcat            >  1:1823916/1‑69 (MQ=255)
                                                     aaGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCaa         <  1:614382/71‑1 (MQ=255)
                                                     aaGCGGCAGATCGTAAAGTGACTCAGGCAGGCTAcc                                            <  1:320041/36‑1 (MQ=255)
                                                     aaGCGGCAGATCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCaa         <  1:1981377/71‑1 (MQ=255)
                                                      aGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTAc                  >  1:1055433/1‑61 (MQ=255)
                                                              gTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACg                        <  1:834900/47‑1 (MQ=255)
                                                              gTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACg                        <  1:26686/47‑1 (MQ=255)
                                                              gTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACg                        <  1:174296/47‑1 (MQ=255)
                                                              gTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACtgtg  >  1:1783845/1‑69 (MQ=255)
                                                              |                                                                    
TACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTG  >  minE/1000124‑1000254

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: