Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 909,663 A→G 46.6% I614V (ATC→GTC)  ycjT → predicted hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE909,6630AG46.6% 5.6 / 80.2 58I614V (ATC→GTC) ycjTpredicted hydrolase
Reads supporting (aligned to +/- strand):  ref base A (18/13);  new base G (14/13);  total (32/26)
Fisher's exact test for biased strand distribution p-value = 7.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCGTGTCTTGCCAATCTGCAATTTTATGAACCGCGCACTATTCACGACTCGTCATTAAGTAAAGCAATCCACGGCATTGTTGCCGCACGCTGTGGCCTGCTGACCCAAAGTTATCAGTTCTGGCGCGAGGGGAC  >  minE/909596‑909730
                                                                   |                                                                   
aTCGTGTCTTGCCAATCTGCATTTTTATGAACCGCGCACTATTCACGACTCGTCATTAAGTAAAGCAGTc                                                                   <  1:593456/70‑1 (MQ=255)
 tCGTGTCTTGCCAATCTGCAATTTTATGAACCGCGCACTATTCACGACTCGTCATTAAGTAAAGCAATCCa                                                                 <  1:1494266/71‑1 (MQ=255)
   gtgtCTTGCCAATCTGCAATTTTATGAACCGCGCACTATTCACGACTCGTCATTAAGTAAAGCAGTCCACg                                                               >  1:3081615/1‑71 (MQ=255)
    tgtCTTGCCAATCTGCAATTTTATGAACCGCGCACTATTCACGACTCGTCATTAAGTAAAGCAATCCACgg                                                              <  1:3918250/71‑1 (MQ=255)
        ttGCCAATCTGCAATTTTATGAACCGCGCACTATTCACGACTCGTCATTAAGTAAAGCAGTCCACGGCAtt                                                          >  1:379199/1‑71 (MQ=255)
        ttGCCAATCTGCAATTTTATGAACCGCGCACTATTCACGACTCGTCATTAAGTAAAGCAGTCCACGGCAt                                                           >  1:1087167/1‑70 (MQ=255)
          gCCAATCTGCAATTTTATGAACCGCGCACTATTCACGACTCGTCATTAAGTAAAGCAGTCCACGGCAttgt                                                        >  1:1859680/1‑71 (MQ=255)
          gCCAATCTGCAATTTTATGAACCGCGCACTATTCACGACTCGTCATTAAGTAAAGCAGTCCACGGCAt                                                           <  1:1016130/68‑1 (MQ=255)
                 tGCAATTTTATGAACCGCTCACTATTCACGACTCGTCATTAAGTAAAGCAATCCACg                                                               >  1:2591768/1‑57 (MQ=255)
                 tGCAATTTTATGAACCGCGCACTATTCACGACTCGTCATTAAGTAAAGCAATCCACg                                                               >  1:3118027/1‑57 (MQ=255)
                 tGCAATTTTATGAACCGCGCACTATTCACGACTCGTCATTAAGTAAAGCAATCCACg                                                               >  1:3556163/1‑57 (MQ=255)
                    aaTTTTATGAACCGCGCACTATTCACGACTCGTCATTAAGTAAAGCAATCCACGGCATTGTTGCCg                                                   <  1:348424/66‑1 (MQ=255)
                      ttatATGAACCGCGCACTATTCACGACTCGTCATTAAGTAAAGCAGTCCACGGCATTGTTGCCGCACGCtg                                            >  1:1941286/4‑71 (MQ=255)
                       tttATGAACCGCGCACTATTCACGACTCGTCATTAAGTAAAGCAATCCACGGCATTGTTGCCGCACGCtgt                                           >  1:3713370/1‑71 (MQ=255)
                        ttATGAACCGCGCACTATTCACGACTCGTCATTAAGTAAAGCAGTCCACGGCATTGTTGCCGCACGCtgtg                                          >  1:476332/1‑71 (MQ=255)
                        ttATGAACCGCGCACTATTCACGACTCGTCATTAAGTAAAGCAATCCACGGCATTGTTGCCGCACGCtgtg                                          >  1:755221/1‑71 (MQ=255)
                         tATGAACCGCGCACTATTCACGACTCGTCATTAAGTAAAGCAGTCCACGGCATTGTTGCCGCACGCtgtg                                          <  1:3564915/70‑1 (MQ=255)
                              aCCGCGCACTATTCACGACTCGTCATTAAGTAAAGCAGTCCACGGCAttgt                                                        >  1:2340308/1‑51 (MQ=255)
                              aCCGCGCACTATTCACGACTCGTCATTAAGTAAAGCAGTCCACGGCAttgt                                                        >  1:1645990/1‑51 (MQ=255)
                              aCCGCGCACTATTCACGACTCGTCATTAAGTAAAGCAGTCCACGGCATTGTTGCCGCACGCTGTGGCctgc                                    >  1:1714678/1‑71 (MQ=255)
                              aCCGCGCACTATTCACGACTCGTCATTAAGTAAAGCAATCCACGGCATTGTTGCCGCACGCTGTGGCctg                                     <  1:3376218/70‑1 (MQ=255)
                                cgcgCACTATTCACGACTCGTCATTAAGTAAAGCAGTCCACGGCATTGTTGCCGCACGCTGTGGCctgct                                   <  1:3891500/70‑1 (MQ=255)
                                    cACTATTCACGACTCGTCATTAAGTAAAGCAATCCACGGCATTGTTGCCGCACGCTGTGGCct                                      <  1:626519/63‑1 (MQ=255)
                                     aCTATTCACGACTCGTCATTAAGTAAAGCAATCCACGGCAtt                                                          <  1:1400316/42‑1 (MQ=255)
                                     aCTATTCACGACTCGTCATTAAGTAAAGCAATCCACGGCATTGTTGCCGCACGCTGTGGCctg                                     <  1:1239426/63‑1 (MQ=255)
                                       tATTCACGACTCGTCATTAAGTAAAGCAGTCCACGGCATTGTTGCCGCACGCtg                                            <  1:1496020/54‑1 (MQ=255)
                                       tATTCACGACTCGTCATTAAGTAAAGCAGTCCACGGCATTGTTGCCGCACGCt                                             <  1:4086859/53‑1 (MQ=255)
                                       tATTCACGACTCGTCATTAAGTAAAGCAGTCCACGGCATTGTTGCCGCACGCt                                             <  1:3699642/53‑1 (MQ=255)
                                       tATTCACGACTCGTCATTAAGTAAAGCAATCCACGGCAttgtt                                                       <  1:23556/43‑1 (MQ=255)
                                        aTTCACGACTCGTCATTAAGTAAAGCAGTCCACGGCATTGTTGCCGCACGc                                              <  1:1375207/51‑1 (MQ=255)
                                        aTTCACGACTCGTCATTAAGTAAAGCAATCCACGGCAttgtt                                                       <  1:2592789/42‑1 (MQ=255)
                                        aTTCACGACTCGTCATTAAGTAAAGCAATCCACGGCAttgtt                                                       <  1:1814806/42‑1 (MQ=255)
                                        aTTCACGACTCGTCATTAAGTAAAGCAATCCACGGCATTGTTGCCGCACGCTGTGGCCTGCTGACCCAAAg                          >  1:771028/1‑71 (MQ=255)
                                        aTTCACGACTCGTCATTAAGTAAAGCAATCCACGGCATTGTTGCCGCACGCTGTGGCCTGCTGACCCAAAg                          >  1:849353/1‑71 (MQ=255)
                                         ttCACGACTCGTCATTAAGTAAAGCAGTCCACGGCATTGTTGCCGCACGCtgt                                           >  1:3415257/1‑53 (MQ=255)
                                         ttCACGACTCGTCATTAAGTAAAGCAGTCCACGGCATTGTTGCCGCACGCTGTGGCCTGCTGCCCCAAAGt                         >  1:3879099/1‑71 (MQ=255)
                                         ttCACGACTCGTCATTAAGTAAAGCAATCCACGGCAt                                                           >  1:3185560/1‑37 (MQ=255)
                                          tCACGACTCGTCATTAAGTAAAGCAATCCACGGCATTGTTGCCGCACGCTGTGGCCTGCTGa                                 <  1:1367784/62‑1 (MQ=255)
                                           cACGACTCGTCATTAAGTAAAGCAATCCACGGCATTGTTGCCGCACGCTGTGGCctgct                                   <  1:624701/59‑1 (MQ=255)
                                           cACGACTCGTCATTAAGTAAAGCAATCCACGGCATTGTTGCCGCACGCTGTGGCctgc                                    >  1:3280730/1‑58 (MQ=255)
                                           cACGACTCGTCATTAAGTAAAGCAATCCACGGCATTGTTGCCGCACGCTGTGGCctgc                                    >  1:3050673/1‑58 (MQ=255)
                                           cACGACTCGTCATTAAGTAAAGCAATCCACGGCATTGTTGCCGCACGCTGTGGCctgc                                    >  1:462422/1‑58 (MQ=255)
                                           cACGACTCGTCATTAAGTAAAGCAATCCACGGCATTGTTGCCGCACGCTGTGGCCTGCTGACCCAAAGTTa                       >  1:713545/1‑71 (MQ=255)
                                           cACGACTCGTCATTAAGTAAAGCAATCCACGGCATTGTTGCCGCACGCTGTGGCCTGCTGACCCAAAGTTa                       >  1:2802119/1‑71 (MQ=255)
                                            aCGACTCGTCATTAAGTAAAGCAGTCCACGGCATTGTTGCCGCACGCTGTGGCCTGCTGACCCAAAGTTAt                      >  1:2525375/1‑71 (MQ=255)
                                            aCGACTCGTCATTAAGTAAAGCAATCCACGGCATTGTTGCCGCACGCTGTGGCCTGCTGACCCAAAGTTAt                      >  1:3082669/1‑71 (MQ=255)
                                             cGACTCGTCATTAAGTAAAGCAATCCACGGCATTGTTGCCg                                                   <  1:3582292/41‑1 (MQ=255)
                                              gACTCGTCATTAAGTAAAGCAATCCACGGCATTGTTGCCGCACGCTGTGGCCTGCTGACCCAAAGTTAt                      >  1:3927270/1‑69 (MQ=255)
                                               aCTCGTCATTAAGTAAAGCAGTCCACGGCATTGTTGCCGCACGCTGTGGCCTGCTGACCCAAAGTTATCAg                   >  1:1402486/1‑71 (MQ=255)
                                                   gTCATTAAGTAAAGCAGTCCACGGCATTGTTGCCGCACGCTGTGGCCTGCTGACCCAAAGTTATCAGTTCt               >  1:2800658/1‑71 (MQ=255)
                                                   gTCATTAAGTAAAGCAATCCACGGCATTGTTGCCGCACGCTGTGGCCTGCTGACCCAAAGTTATCAGTTCt               >  1:3662017/1‑71 (MQ=255)
                                                      aTTAAGTAAAGCAGTCCACGGCATTGTTGCCGCACGCTGTGGCCTGCTGACCCAAAGTTATCAGTTCTgg             <  1:2836073/70‑1 (MQ=255)
                                                        tAAGTAAAGCAGTCCACGGCATTGTTGCCGCACGCt                                             <  1:599457/36‑1 (MQ=255)
                                                        tAAGTAAAGCAGTCCACGGCATTGTTGCCGCACGCt                                             <  1:1663665/36‑1 (MQ=255)
                                                           gTAAAGCAGTCCACGGCATTGTTGCCGCACGCTGTGGCCTGCTGACCCAAAGTTATCAGTTCTGGCGCGAg       <  1:1793540/71‑1 (MQ=255)
                                                           gTAAAGCAATCCACGGCATTGTTGCCGCACGCTGTGGCCTGCTGACCCAAAGTTATCAGTTCTGGCGCCAg       >  1:1863911/1‑71 (MQ=255)
                                                             aaaGCAATCCACGGCATTGTTGCCGCACGCTGTGGCCTGCTGACCCAAAGTTATCAGTTCTg              >  1:3870361/1‑62 (MQ=255)
                                                                gCAGTCCACGGCATTGTTGCCGCACGCTGTGGCCTGCTGACCCAAAGTTATCAGTTCTGGCGCGAGGGGAc  <  1:2281977/71‑1 (MQ=255)
                                                                   |                                                                   
ATCGTGTCTTGCCAATCTGCAATTTTATGAACCGCGCACTATTCACGACTCGTCATTAAGTAAAGCAATCCACGGCATTGTTGCCGCACGCTGTGGCCTGCTGACCCAAAGTTATCAGTTCTGGCGCGAGGGGAC  >  minE/909596‑909730

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: