Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 915,319 T→A 26.4% P421P (CCT→CCA ycjX → conserved hypothetical protein with nucleoside triphosphate hydrolase domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE915,3190TA26.4% 104.8 / 38.6 72P421P (CCT→CCAycjXconserved hypothetical protein with nucleoside triphosphate hydrolase domain
Reads supporting (aligned to +/- strand):  ref base T (7/46);  new base A (6/13);  total (13/59)
Fisher's exact test for biased strand distribution p-value = 9.06e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAATCGACTTAGCGATGGCGCACCGCTCACTGTTTATCCTGGCGAAGTTCCCGCACGTTTGCCTGGTCAGGCGTTCTGGGATAAGCAAGGCTTCCAGTTTGAGGCATTTCGTCCGCAGGTGATGG  >  minE/915256‑915380
                                                               |                                                             
tAATCGACTCAGCGATGGAGCACCGCTCACTGTTTATCCTGGCGAAGTTCCCGCACGTTTGCCAGGTc                                                           >  1:1023472/1‑68 (MQ=255)
tAATCGACTCAGCGATGGAGCACCGCTCACTGTTTATCCTGGCGAAGTTCCCGCACGTTTGCCAGGTc                                                           >  1:39874/1‑68 (MQ=255)
tAATCGACTCAGCGATGGAGCACCGCTCACTGTTTATCCTGGCGAAGTTCCCGCACGTTTGCCAGGTCa                                                          >  1:1233286/1‑69 (MQ=255)
tAATCGACTCAGCGATGGAGCACCGCTCACTGTTTATCCTGGCGAAGTTCCCGCACGTTTGCCAGGTCa                                                          >  1:2399660/1‑69 (MQ=255)
     gACTTAGCGATGGCGCACCGCTCACTGTTTATCCTGGCGAAGTTCCCGCACGTTTGCCTGGTCAGGCGt                                                     <  1:2505303/69‑1 (MQ=255)
     gACTTAGCGATGGCGCACCGCTCACTGTTTATCCTGGCGAAGTTCCCGCACGTTTGCCTGGTCAGGCGt                                                     <  1:2009501/69‑1 (MQ=255)
     gACTTAGCGATGGCGCACCGCTCACTGTTTATCCTGGCGAAGTTCCCGCACGTTTGCCTGGTCAGGCGt                                                     <  1:3262538/69‑1 (MQ=255)
                               gTTTATCCTGGCGAAGTTCCCGCACGTTTGCCTGGTCAGGCGTTCTGGGATAAGCAAGGCTTCCAGt                             >  1:2733528/1‑67 (MQ=255)
                               gTTTATCCTGGCGAAGTTCCCGCACGTTTGCCTGGTCAGGCGTTCTGGGATAAGCAAGGCTTCCAGt                             >  1:2499804/1‑67 (MQ=255)
                               gTTTATCCTGGCGAAGTTCCCGCACGTTTGCCTGGTCAGGCGTTCTGGGATAAGCAAGGCTTCCAGt                             >  1:2918708/1‑67 (MQ=255)
                               gTTTATCCTGGCGAAGTTCCCGCACGTTTGCCTGGTCAGGCGTTCTGGGATAAGCAAGGCTTCCAGt                             >  1:3285987/1‑67 (MQ=255)
                               gTTTATCCTGGCGAAGTTCCCGCACGTTTGCCTGGTCAGGCGTTCTGGGATAAGCAAGGCTTCCAGt                             >  1:36533/1‑67 (MQ=255)
                                tttATCCTGGCGAAGTTCCCGCACGTTTGCCTGGTCAGGCGTTCTGGGATAAGCAAGGCTTCCAGTTTg                          >  1:527330/1‑69 (MQ=255)
                                  tATCCTGGCGAAGTTCCCGCACGTTTGCCTGGTCAGGCGTTCTGGGATAAGCAAGGCTTCCAGTTTGa                         <  1:349613/68‑1 (MQ=255)
                                  tATCCTGGCGAAGTTCCCGCACGTTTGCCTGGTCAGGCGTTCTGGGATAAGCAAGGCTTCCAGTTTGa                         <  1:556013/68‑1 (MQ=255)
                                  tATCCTGGCGAAGTTCCCGCACGTTTGCCTGGTCAGGCGTTCTGGGATAAGCAAGGCTTCCAGTTTGa                         <  1:3154880/68‑1 (MQ=255)
                                  tATCCTGGCGAAGTTCCCGCACGTTTGCCTGGTCAGGCGTTCTGGGATAAGCAAGGCTTCCAGTTTGa                         <  1:960811/68‑1 (MQ=255)
                                  tATCCTGGCGAAGTTCCCGCACGTTTGCCTGGTCAGGCGTTCTGGGATAAGCAAGGCTTCCAGTTTGa                         >  1:2143252/1‑68 (MQ=255)
                                      cTGGCGAAGTTCCCGCACGTTTGCCTGGTCAGGCGTTCTGGGATAAGCAAGGCTTCCAGTTTGAGGCAt                    <  1:358366/69‑1 (MQ=255)
                                      cTGGCGAAGTTCCCGCACGTTTGCCTGGTCAGGCGTTCTGGGATAAGCAAGGCTTCCAGTTTGAGGCAt                    <  1:872505/69‑1 (MQ=255)
                                                tCCCGCACGTTTGCCTGGTCAGGCGTTCTGGGATAAGCAAGGCTTCCAGTTTGAGGCATTTCGTCCGCa          <  1:2977994/69‑1 (MQ=255)
                                                tCCCGCACGTTTGCCTGGTCAGGCGTTCTGGGATAAGCAAGGCTTCCAGTTTGAGGCATTTCGTCCGCa          <  1:852137/69‑1 (MQ=255)
                                                tCCCGCACGTTTGCCTGGTCAGGCGTTCTGGGATAAGCAAGGCTTCCAGTTTGAGGCATTTCGTCCGCa          <  1:275917/69‑1 (MQ=255)
                                                tCCCGCACGTTTGCCTGGTCAGGCGTTCTGGGATAAGCAAGGCTTCCAGTTTGAGGCATTTCGTCCGCa          <  1:3039245/69‑1 (MQ=255)
                                                tCCCGCACGTTTGCCAGGTCAGGCGTTCTGGGATAAGCAAGGCTTCCAGTTTGAGGCGTTTCGTcc             <  1:827031/66‑1 (MQ=255)
                                                tCCCGCACGTTTGCCAGGTCAGGCGTTCTGGGATAAGCAAGGCTTCCAGTTTGAGGCGTTTCGTcc             <  1:2859873/66‑1 (MQ=255)
                                                tCCCGCACGTTTGCCAGGTCAGGCGTTCTGGGATAAGCAAGGCTTCCAGTTTGAGGCGTTTCGTcc             <  1:458830/66‑1 (MQ=255)
                                                tCCCGCACGTTTGCCAGGTCAGGCGTTCTGGGATAAGCAAGGCTTCCAGTTTGAGGCGTTTCGTcc             <  1:751816/66‑1 (MQ=255)
                                                tCCCGCACGTTTGCCAGGTCAGGCGTTCTGGGATAAGCAAGGCTTCCAGTTTGAGGCGTTTCGTcc             <  1:3112420/66‑1 (MQ=255)
                                                tCCCGCACGTTTGCCAGGTCAGGCGTTCTGGGATAAGCAAGGCTTCCAGTTTGAGGCGTTTCGTcc             <  1:976993/66‑1 (MQ=255)
                                                tCCCGCACGTTTGCCAGGTCAGGCGTTCTGGGATAAGCAAGGCTTCCAGTTTGAGGCGTTTCGTcc             <  1:1364723/66‑1 (MQ=255)
                                                tCCCGCACGTTTGCCAGGTCAGGCGTTCTGGGATAAGCAAGGCTTCCAGTTTGAGGCGTTTCGTcc             <  1:1379353/66‑1 (MQ=255)
                                                tCCCGCACGTTTGCCAGGTCAGGCGTTCTGGGATAAGCAAGGCTTCCAGTTTGAGGCGTTTCGTcc             <  1:1624435/66‑1 (MQ=255)
                                                tCCCGCACGTTTGCCAGGTCAGGCGTTCTGGGATAAGCAAGGCTTCCAGTTTGAGGCGTTTCGTcc             <  1:2687971/66‑1 (MQ=255)
                                                tCCCGCACGTTTGCCAGGTCAGGCGTTCTGGGATAAGCAAGGCTTCCAGTTTGAGGCGTTTCGTcc             <  1:186824/66‑1 (MQ=255)
                                                tCCCGCACGTTTGCCAGGTCAGGCGTTCTGGGATAAGCAAGGCTTCCAGTTTGAGGCGTTTCGTcc             <  1:1924692/66‑1 (MQ=255)
                                                tCCCGCACGTTTGCCAGGTCAGGCGTTCTGGGATAAGCAAGGCTTCCAGTTTGAGGCGTTTCGTcc             <  1:1994818/66‑1 (MQ=255)
                                                tCCCGCACGTTTGCCAGGTCAGGCGTTCTGGGATAAGCAAGGCTTCCAGTTTGAGGCGTTTCGTCCGCa          >  1:800615/1‑69 (MQ=255)
                                                tCCCGCACGTTTGCCAGGTCAGGCGTTCTGGGATAAGCAAGGCTTCCAGTTTGAGGCGTTTCGTCCGCa          >  1:1555790/1‑69 (MQ=255)
                                                    gCACGTTTGCCTGGTCAGGCGTTCTGGGATAAGCAAGGCTTCCAGTTTGAGGCATTTCGTCCGGAGGt       <  1:962342/68‑1 (MQ=255)
                                                    gCACGTTTGCCTGGTCAGGCGTTCTGGGATAAGCAAGGCTTCCAGTTTGAGGCATTTCGTCCGCAGGt       <  1:3348317/68‑1 (MQ=255)
                                                    gCACGTTTGCCTGGTCAGGCGTTCTGGGATAAGCAAGGCTTCCAGTTTGAGGCATTTCGTCCGCAGGt       <  1:596902/68‑1 (MQ=255)
                                                    gCACGTTTGCCTGGTCAGGCGTTCTGGGATAAGCAAGGCTTCCAGTTTGAGGCATTTCGTCCGCAGGt       <  1:595173/68‑1 (MQ=255)
                                                    gCACGTTTGCCTGGTCAGGCGTTCTGGGATAAGCAAGGCTTCCAGTTTGAGGCATTTCGTCCGCAGGt       <  1:342062/68‑1 (MQ=255)
                                                    gCACGTTTGCCTGGTCAGGCGTTCTGGGATAAGCAAGGCTTCCAGTTTGAGGCATTTCGTCCGCAGGt       <  1:2929116/68‑1 (MQ=255)
                                                    gCACGTTTGCCTGGTCAGGCGTTCTGGGATAAGCAAGGCTTCCAGTTTGAGGCATTTCGTCCGCAGGt       <  1:2162425/68‑1 (MQ=255)
                                                    gCACGTTTGCCTGGTCAGGCGTTCTGGGATAAGCAAGGCTTCCAGTTTGAGGCATTTCGTCCGCAGGt       <  1:2880410/68‑1 (MQ=255)
                                                    gCACGTTTGCCTGGTCAGGCGTTCTGGGATAAGCAAGGCTTCCAGTTTGAGGCATTTCGTCCGCAGGt       <  1:2988937/68‑1 (MQ=255)
                                                    gCACGTTTGCCTGGTCAGGCGTTCTGGGATAAGCAAGGCTTCCAGTCTGAGGCATTTCGTCCGCAGGt       <  1:2687519/68‑1 (MQ=255)
                                                       cGTTTGCCTGGTCAGGCGTTCTGGGATATGCAAGGCTTCCAGTTTGAGGCATTTCGTCCGCAGGTgatg   <  1:2983450/69‑1 (MQ=255)
                                                       cGTTTGCCTGGTCAGGCGTTCTGGGATAAGCAAGGCTTCCAGTTTGAGGCATTTCGTCCGCAGGTgatg   <  1:1436313/69‑1 (MQ=255)
                                                       cGTTTGCCTGGTCAGGCGTTCTGGGATAAGCAAGGCTTCCAGTTTGAGGCATTTCGTCCGCAGGTgatg   <  1:2950854/69‑1 (MQ=255)
                                                       cGTTTGCCTGGTCAGGCGTTCTGGGATAAGCAAGGCTTCCAGTTTGAGGCATTTCGTCCGCAGGTgatg   <  1:1432385/69‑1 (MQ=255)
                                                        gTTTGCCTGGTCAGGCGTTCTGGGATAAGCAAGGCTTCCAGTTTGAGGCATTTCGTCCGCAGGTgatgg  <  1:3172921/69‑1 (MQ=255)
                                                        gTTTGCCTGGTCAGGCGTTCTGGGATAAGCAAGGCTTCCAGTTTGAGGCATTTCGTCCGCAGGTgatgg  <  1:285602/69‑1 (MQ=255)
                                                        gTTTGCCTGGTCAGGCGTTCTGGGATAAGCAAGGCTTCCAGTTTGAGGCATTTCGTCCGCAGGTgatgg  <  1:1451132/69‑1 (MQ=255)
                                                        gTTTGCCTGGTCAGGCGTTCTGGGATAAGCAAGGCTTCCAGTTTGAGGCATTTCGTCCGCAGGTgatg   <  1:1647191/68‑1 (MQ=255)
                                                         tttGCCTGGTCAGGCGTTCTGGGATAAGCAAGGCTTCCAGTTTGAGGCATTTCGTCCGc           <  1:1384154/59‑1 (MQ=255)
                                                         tttGCCTGGTCAGGCGTTCTGGGATAAGCAAGGCTTCCAGTTTGAGGCATTTCGTCCGc           <  1:844214/59‑1 (MQ=255)
                                                         tttGCCTGGTCAGGCGTTCTGGGATAAGCAAGGCTTCCAGTTTGAGGCATTTCGTCCGc           <  1:1488252/59‑1 (MQ=255)
                                                         tttGCCTGGTCAGGCGTTCTGGGATAAGCAAGGCTTCCAGTTTGAGGCATTTCGTCCGc           <  1:2730110/59‑1 (MQ=255)
                                                         tttGCCTGGTCAGGCGTTCTGGGATAAGCAAGGCTTCCAGTTTGAGGCATTTCGTCCGc           <  1:968703/59‑1 (MQ=255)
                                                         tttGCCTGGTCAGGCGTTCTGGGATAAGCAAGGCTTCCAGTTTGAGGCATTTCGTCCGc           <  1:2725154/59‑1 (MQ=255)
                                                         tttGCCTGGTCAGGCGTTCTGGGATAAGCAAGGCTTCCAGTTTGAGGCATTTCGTCCGc           <  1:2326523/59‑1 (MQ=255)
                                                         tttGCCTGGTCAGGCGTTCTGGGATAAGCAAGGCTTCCAGTTTGAGGCATTTCGTCCGc           <  1:1440726/59‑1 (MQ=255)
                                                         tttGCCTGGTCAGGCGTTCTGGGATAAGCAAGGCTTCCAGTTTGAGGCATTTCGTCCGc           <  1:2278188/59‑1 (MQ=255)
                                                         tttGCCTGGTCAGGCGTTCTGGGATAAGCAAGGCTTCCAGTTTGAGGCATTTCGTCCGc           <  1:1660983/59‑1 (MQ=255)
                                                         tttGCCTGGTCAGGCGTTCTGGGATAAGCAAGGCTTCCAGTTTGAGGCATTTCGTCCGc           <  1:1772923/59‑1 (MQ=255)
                                                         tttGCCTGGTCAGGCGTTCTGGGATAAGCAAGGCTTCCAGTTTGAGGCATTTCGTCCGc           <  1:3324949/59‑1 (MQ=255)
                                                         tttGCCTGGTCAGGCGTTCTGGGATAAGCAAGGCTTCCAGTTTGAGGCATTTCGTCCGc           <  1:1803338/59‑1 (MQ=255)
                                                         tttGCCTGGTCAGGCGTTCTGGGATAAGCAAGGCTTCCAGTTTGAGGCATTTCGTCCGc           <  1:2070250/59‑1 (MQ=255)
                                                         tttGCCTGGTCAGGCGTTCTGGGATAAGCAAGGCTTCCAGTTTGAGGCATTTCGTCCGc           <  1:1954836/59‑1 (MQ=255)
                                                               |                                                             
TAATCGACTTAGCGATGGCGCACCGCTCACTGTTTATCCTGGCGAAGTTCCCGCACGTTTGCCTGGTCAGGCGTTCTGGGATAAGCAAGGCTTCCAGTTTGAGGCATTTCGTCCGCAGGTGATGG  >  minE/915256‑915380

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: