Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,801,633 C→T 43.4% intergenic (+53/‑81) yjbN → / → yjbO tRNA‑dihydrouridine synthase A/phage shock protein G

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,801,6330CT43.4% 19.2 / 50.9 69intergenic (+53/‑81)yjbN/yjbOtRNA‑dihydrouridine synthase A/phage shock protein G
Reads supporting (aligned to +/- strand):  ref base C (13/26);  new base T (13/17);  total (26/43)
Fisher's exact test for biased strand distribution p-value = 4.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAAGCGTTAACTTTTCACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTCGCGGGGCGTTTTGCTGTTAAATCAATAGATTATTTTTGGCATGATTCTTGTAATGCCAGCAAGAGATTTCA  >  minE/2801571‑2801697
                                                              |                                                                
tAAGCGTTAACTTTTCACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTCGCGGGGCGTTTTg                                                            >  1:3356151/1‑69 (MQ=255)
tAAGCGTTAACTTTTCACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTCGCGGGGCGTTTTg                                                            >  1:3071926/1‑69 (MQ=255)
tAAGCGTTAACTTTTCACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTCGCGGGGCGTTTTg                                                            >  1:2153423/1‑69 (MQ=255)
tAAGCGTTAACTTTTCACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTCGCGGGGCGTTTTg                                                            >  1:2149573/1‑69 (MQ=255)
             ttCACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTCGCGGGGTGtttt                                                             >  1:1965645/1‑55 (MQ=255)
             ttCACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTCGCGGGGTGTTTTg                                                            >  1:2279331/1‑56 (MQ=255)
             ttCACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTCGCGGGGTGTTTTGCTGTTAAAt                                                   >  1:2449176/1‑65 (MQ=255)
             ttCACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTCGCGGGGTGTTTTGCTGTTAAATc                                                  >  1:2405131/1‑66 (MQ=255)
             ttCACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTCGCGGGGTGTTTTGCTGTTAAATc                                                  >  1:813851/1‑66 (MQ=255)
             ttCACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTCGCGGGGTGTTTTGCTGTTAAATc                                                  >  1:652812/1‑66 (MQ=255)
             ttCACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTCGCGGGGTGTTTTGCTGTTAAATc                                                  >  1:544533/1‑66 (MQ=255)
             ttCACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTCGCGGGGTGTTTTGCTGTTAAATc                                                  >  1:522700/1‑66 (MQ=255)
             ttCACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTCGCGGGGTGTTTTGCTGTTAAATc                                                  >  1:40684/1‑66 (MQ=255)
             ttCACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTCGCGGGGTGTTTTGCTGTTAAATc                                                  >  1:2832610/1‑66 (MQ=255)
             ttCACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTCGCGGGGTGTTTTGCTGTTAAATc                                                  >  1:1894052/1‑66 (MQ=255)
             ttCACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTCGCGGGGTGTTTTGCTGTTAAATc                                                  >  1:1051420/1‑66 (MQ=255)
             ttCACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTCGCGGGGTGTTTTGCTGTTAAATc                                                  >  1:1089970/1‑66 (MQ=255)
               cACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTCGCGGGGTGTTTTGCTGTTAAATCaaca                                              <  1:1792299/68‑3 (MQ=255)
               cACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTCGCGGGGTGTTTTGCTGTTAAATCaaca                                              <  1:516220/68‑3 (MQ=255)
               cACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTCGCGGGGTGTTTTGCTGTTAAATCaaca                                              <  1:1249532/68‑3 (MQ=255)
               cACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTCGCGGGGTGTTTTGCTGTTAAATCaaca                                              <  1:3142905/68‑3 (MQ=255)
               cACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTCGCGGGGTGTTTTGCTGTTAAATCaaca                                              <  1:3130979/68‑3 (MQ=255)
               cACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTCGCGGGGTGTTTTGCTGTTAAATCaaca                                              <  1:1333230/68‑3 (MQ=255)
               cACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTCGCGGGGTGTTTTGCTGTTAAATCaaca                                              <  1:302260/68‑3 (MQ=255)
               cACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTCGCGGGGTGTTTTGCTGTTAAATCaaca                                              <  1:298721/68‑3 (MQ=255)
               cACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTCGCGGGGTGTTTTGCTGTTAAATCaaca                                              <  1:2491660/68‑3 (MQ=255)
               cACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTCGCGGGGTGTTTTGCTGTTAAATCaaca                                              <  1:1417274/68‑3 (MQ=255)
               cACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTCGCGGGGTGTTTTGCTGTTAAATCaaca                                              <  1:2315945/68‑3 (MQ=255)
               cACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTCGCGGGGTGTTTTGCTGTTAAATCaaca                                              <  1:2200932/68‑3 (MQ=255)
               cACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTCGCGGGGTGTTTTGCTGTTAAATCaaca                                              <  1:2013072/68‑3 (MQ=255)
               cACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTCGCGGGGTGTTTTGCTGTTAAATCaaca                                              <  1:2072277/68‑3 (MQ=255)
               cACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTCGCGGGGTGTTTTGCTGTTAAATCaaca                                              <  1:1568680/68‑3 (MQ=255)
               cACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTCGCGGGGTGTTTTGCTGTTAAATCaaca                                              <  1:1959456/68‑3 (MQ=255)
                 ccAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTCGCGGGGTGTTTTGCTGTTAAATCaaca                                              <  1:355088/66‑3 (MQ=255)
                  cAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTCGCGGGGCGTTTTGCTGTTAAATCAATAGatt                                          >  1:199435/1‑69 (MQ=255)
                  cAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTCGCGGGGCGTTTTGCTGTTAAATCAATAGatt                                          >  1:3313363/1‑69 (MQ=255)
                  cAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTCGCGGGGCGTTTTGCTGTTAAATCAATAGat                                           <  1:1723260/68‑1 (MQ=255)
                  cAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTCGCGGGGCGTTTTGCTGTTAAATCAATAGat                                           <  1:3217003/68‑1 (MQ=255)
                  cAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTCGCGGGGCGTTTTGCTGTTAAATCAATAGat                                           >  1:1712508/1‑68 (MQ=255)
                  cAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTCGCGGGGCGTTTTGCTGTTAAATCAATAGat                                           <  1:708604/68‑1 (MQ=255)
                         tagtCAAATTCACCACGCCCTGCGCACCGTCGCGGGGCGTTTTGCTGTTAAATCAATAGATTAttttt                                    >  1:1228282/1‑68 (MQ=255)
                             cAAATTCACCACGCCCTGCGCACCGTCGCGGGGCGTTTTGCTGTTAAATCAATAGATTATTTTTGGCAt                               <  1:2681680/69‑1 (MQ=255)
                             cAAATTCACCACGCCCTGCGCACCGTCGCGGGGCGTTTTGCTGTTAAATCAATAGATTATTTTTGGCAt                               <  1:278396/69‑1 (MQ=255)
                              aaaTTCACCACGCCCTGCGCACCGTCGCGGGGCGTTTTGCTGTTAAATCAATAGATTATTTTTGGCATg                              <  1:702311/69‑1 (MQ=255)
                                 ttCACCACGCCCTGCGCACCGTCGCGGGGCGTTTTGCTGTTAAATCAATAGATTATTTTTGGCATGAtt                           <  1:2115251/69‑1 (MQ=255)
                                 ttCACCACGCCCTGCGCACCGTCGCGGGGCGTTTTGCTGTTAAATCAATAGATTATTTTTGGCATGAtt                           <  1:132026/69‑1 (MQ=255)
                                 ttCACCACGCCCTGCGCACCGTCGCGGGGCGTTTTGCTGTTAAATCAATAGATTATTTTTGGCATGAtt                           <  1:2103581/69‑1 (MQ=255)
                                      cacGCCCTGCGCACCGTCGCGGGGCGTTTTGCTGTTAAATCAATAGATTATTTTTGGCATGATTCTTGt                      <  1:3180549/69‑1 (MQ=255)
                                      cacGCCCTGCGCACCGTCGCGGGGCGTTTTGCTGTTAAATCAATAGATTATTTTTGGCATGATTCTTGt                      <  1:706699/69‑1 (MQ=255)
                                      cacGCCCTGCGCACCGTCGCGGGGCGTTTTGCTGTTAAATCAATAGATTATTTTTGGCATGATTCTTGt                      <  1:2208800/69‑1 (MQ=255)
                                      cacGCCCTGCGCACCGTCGCGGGGCGTTTTGCTGTTAAATCAATAGATTATTTTTGGCATGATTCTTGt                      <  1:140368/69‑1 (MQ=255)
                                      cacGCCCTGCGCACCGTCGCGGGGCGTTTTGCTGTTAAATCAATAGATTATTTTTGGCATGATTCTTGt                      <  1:2828380/69‑1 (MQ=255)
                                      cacGCCCTGCGCACCGTCGCGGGGCGTTTTGCTGTTAAATCAATAGATTATTTTTGGCATGATTCTTGt                      <  1:2692440/69‑1 (MQ=255)
                                      cacGCCCTGCGCACCGTCGCGGGGCGTTTTGCTGTTAAATCAATAGATTATTTTTGGCATGATTCTTGt                      <  1:2568924/69‑1 (MQ=255)
                                      cacGCCCTGCGCACCGTCGCGGGGCGTTTTGCTGTTAAATCAATAGATTATTTTTGGCATGATTCTTGt                      <  1:2489062/69‑1 (MQ=255)
                                       acGCCCTGCGCACCGTCGCGGGGCGTTTTGCTGTTAAATCAATAGa                                            >  1:3250039/1‑46 (MQ=255)
                                       acGCCCTGCGCACCGTCGCGGGGCGTTTTGCTGTTAAATCAATAGATTATTTTTGGCATGATTCTTGTa                     >  1:1004713/1‑69 (MQ=255)
                                       acGCCCTGCGCACCGTCGCGGGGCGTTTTGCTGTTAAATCAATAGATTATTTTTGGCATGATTCTTGTa                     >  1:2298423/1‑69 (MQ=255)
                                       acGCCCGGCGCACCGTCGCGGGGCGTTTTGCTGTTAAATCAATAGATTATTTTTGGCATGATTCTTGTa                     >  1:4235/1‑69 (MQ=255)
                                       acGCCCGGCGCACCGTCGCGGGGCGTTTTGCTGTTAAATCAATAGATTATTTTTGGCATGATTCTTGTa                     >  1:434867/1‑69 (MQ=255)
                                         gCCCTGCGCACCGTCGCGGGGCGTTTTGCTGTTAAATCAATAGATTATTTTTGGCATGATTCTTTTAAt                   <  1:2244046/69‑1 (MQ=255)
                                         gCCCTGCGCACCGTCGCGGGGCGTTTTGCTGTTAAATCAATAGATTATTTTTGGCATGATTCTTGTAAt                   <  1:1828967/69‑1 (MQ=255)
                                         gCCCTGCGCACCGTCGCGGGGCGTTTTGCTGTTAAATCAATAGATTATTTTTGGCATGATTCTTGTAAt                   <  1:3035290/69‑1 (MQ=255)
                                         gCCCTGCGCACCGTCGCGGGGCGTTTTGCTGTTAAATCAATAGATTATTTTTGGCATGATTCTTGTAAt                   <  1:1553482/69‑1 (MQ=255)
                                         gCCCTGCGCACCGTCGCGGGGCGTTTTGCTGTTAAATCAATAGATTATTTTTGGCATGATTCTTGTAAt                   <  1:1063069/69‑1 (MQ=255)
                                                   ccGTCGCGGGGCGTTTTGCTGTTAAATCAATAGATTATTTTTGGCATGATTCTTGTAATGCCAGCAaga         <  1:647706/69‑1 (MQ=255)
                                                           gggCGTTTTGCTGTTAAATCAATAGATTATTTTTGGCATg                              <  1:1906376/40‑1 (MQ=255)
                                                           gggCGTTTTGCTGTTAAATCAATAGATTATTTTTGGCATGATTCTTGTAATGCCAGCAAGAGATTTCa  <  1:546634/68‑1 (MQ=255)
                                                            ggCGTTTTGCTGTTAAATCAATAGATTATTTTTGGCATGATTCTTGt                      <  1:1955883/47‑1 (MQ=255)
                                                              |                                                                
TAAGCGTTAACTTTTCACCAAAAAGTAGTCAAATTCACCACGCCCTGCGCACCGTCGCGGGGCGTTTTGCTGTTAAATCAATAGATTATTTTTGGCATGATTCTTGTAATGCCAGCAAGAGATTTCA  >  minE/2801571‑2801697

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: