Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,933,825 A→T 10.5% V206V (GTT→GTA lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,933,8250AT10.5% 149.8 / 8.9 57V206V (GTT→GTAlysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base A (28/23);  new base T (3/3);  total (31/26)
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.64e-01

CGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGTTCACCATGAACTGGCGAATACCAGAGAGGATCTGCGAGCGCACTTTAAAG  >  minE/1933756‑1933894
                                                                     |                                                                     
cgcgTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCaa                                                                      >  1:2246622/1‑71 (MQ=255)
   attATGATGGGTAATAAACGGGCGAGCAGATGCCCCACCTGGAATTACCTGCATCATCGGGGTTTCTACtt                                                                   >  1:3654041/2‑71 (MQ=255)
     tATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTc                                                                  >  1:3135061/1‑70 (MQ=255)
       tggtggGCGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCAt                                                               <  1:3850904/71‑1 (MQ=255)
       tgatgggtaATAAACGGGCGAGCAGATGCCCCACCTGGAATTACCTGCATCATCGGGGTTTCTACTTcc                                                                 >  1:2541611/4‑69 (MQ=255)
         gtgggtgATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATa                                                              >  1:2177956/1‑70 (MQ=255)
         atgggtaATAAACGGGCGAGCAGATGCCCCACCTGGAATTACCTGCATCATCGGGGTTTCTACTTCCATa                                                              <  1:2810000/69‑1 (MQ=255)
         atgggtaATAAACGGGCGAGCAGATGCCCCACCTGGAATTACCTGCATCATCGGGGTTTCTACTTCCATa                                                              <  1:3696467/69‑1 (MQ=255)
          tgggtgATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATa                                                              >  1:1411182/1‑69 (MQ=255)
           gggtgATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAg                                                           >  1:1414807/1‑71 (MQ=255)
            ggtgATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAg                                                           <  1:2314223/70‑1 (MQ=255)
            ggtgATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGc                                                          >  1:3260724/1‑71 (MQ=255)
             gtgATAAACGGACGCGCAGCGGCTCCGCCAGGGATCACCTGCATCATCGGCGTTTCAACt                                                                    >  1:1454886/1‑60 (MQ=255)
                aTAAACGGGCGAGCAGATGCCCCACCTGGAATTACCTGCATCATCGGGGTTTCTACTTCCATAAAGCcgc                                                       <  1:3420029/70‑1 (MQ=255)
                 tAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCgg                                                     >  1:3834172/1‑71 (MQ=255)
                 tAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCgg                                                     >  1:2073661/1‑71 (MQ=255)
                  aaaCGGGCGAGCAGATGCCCCACCTGGAATTACCTGCATCATCGGGGTTTCTACTTCCATAAAGCcgcgcg                                                    >  1:269631/1‑69 (MQ=255)
                  aaaCGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGt                                                    >  1:333916/1‑71 (MQ=255)
                   aaCGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGtt                                                   <  1:3052754/71‑1 (MQ=255)
                       gACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGTTCAcc                                               <  1:1859456/71‑1 (MQ=255)
                       gACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGTTCAcc                                               <  1:1877772/71‑1 (MQ=255)
                       gACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGTTCAcc                                               <  1:3683957/71‑1 (MQ=255)
                       gACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGTTCAcc                                               <  1:1816317/71‑1 (MQ=255)
                       gACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGTTCAcc                                               <  1:1664655/71‑1 (MQ=255)
                        aCGCGCAGCGGCACCGCCAGGGATCACGTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGTTCACCa                                              >  1:1597963/1‑71 (MQ=255)
                         cgcgCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTcc                                                                 >  1:2770204/1‑51 (MQ=255)
                          gcgcAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGTTCACCAt                                             <  1:1095457/70‑1 (MQ=255)
                               gCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGTTCACCATGAACTg                                       <  1:550273/71‑1 (MQ=255)
                                  gCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGTTCACCATGAACTGGc                                     <  1:1325918/70‑1 (MQ=255)
                                  gCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGTTCACCATGAACTGGCg                                    <  1:2587624/71‑1 (MQ=255)
                                  gCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGTTCACCATGAACTGGCg                                    <  1:565835/71‑1 (MQ=255)
                                     ccGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGTTCACCATGa                                           >  1:585602/1‑61 (MQ=255)
                                      cGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGTTCACCATGAACTGGCGAATa                                >  1:625992/1‑71 (MQ=255)
                                       gCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGTTCACCATGaa                                          >  1:1156056/1‑60 (MQ=255)
                                       gCCAGGGATCACCTGCATCATCGGCGTTTCAACATCCATAAAGCCGCGGTTCACCATGaa                                          >  1:1745595/1‑60 (MQ=255)
                                          aGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGTTCACCATGAACTGGCGAATAc                               >  1:2608047/1‑68 (MQ=255)
                                             gATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGTTCACCATGAACTGGCg                                    <  1:154930/60‑1 (MQ=255)
                                             gATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGTTCACCATGAACTGGCGAATACCagaga                         <  1:1872459/71‑1 (MQ=255)
                                              aTCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGTTCAcc                                               >  1:3484387/1‑48 (MQ=255)
                                              aTCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGTTCACCATGAACTGGCGa                                   >  1:854294/1‑60 (MQ=255)
                                              aTCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGTTCACCATGAACTGGCGAATACCagagag                        <  1:3835040/71‑1 (MQ=255)
                                              aTCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGTTCACCATGAACTGGCGAATACCagagag                        <  1:3152919/71‑1 (MQ=255)
                                               tCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGTTCACCATGAACTGGCGAATACCagagag                        >  1:1154825/1‑70 (MQ=255)
                                               tCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGTTCACCATGAACTGGCGAATACCAGAGAgg                       <  1:2071003/71‑1 (MQ=255)
                                                 aCCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGTTCAc                                                >  1:989205/1‑44 (MQ=255)
                                                 aCCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGTTCACCATGAACTGGCGAATACCAGAGAGGAt                     >  1:2978667/1‑71 (MQ=255)
                                                    tGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGTTCACCATGAACTGGCGAATACCAGAGAGGATCTg                  >  1:3429732/1‑71 (MQ=255)
                                                      catcatCGGCGTTTCAACTTCCATAAAGCCGCGGTTCACCATGAACTgg                                      >  1:2739514/1‑49 (MQ=255)
                                                      catcatCGGCGTTTCAACTTCCATAAAGCCGCGGTTCACCATGAACTGGCGAATACCAGAGAGGATCTGc                 >  1:14197/1‑70 (MQ=255)
                                                           tcGGCGTTTCAACTTCCATAAAGCCGCGGTTCACCATGAACTGGCGAATACCAGAGAGGATCTGCGAgcgc           <  1:3439968/71‑1 (MQ=255)
                                                            cGGCGTTTCAACTTCCATAAAGCCGCGGTTCACCATGAACTGGCGAATACCAGAGAGGATCTGCGAGCGCa          >  1:2287857/1‑71 (MQ=255)
                                                             ggCGTTTCAACTTCCATAAAGCCGCGGTTCACCATGAACTGGCGAATACCAGAGAGGATCTGCGAGCGCa          <  1:2028368/70‑1 (MQ=255)
                                                                 tttCAACTTCCATAAAGCCGCGGTTCACCATGAACTGGCGAATACCAGAGAGGATCt                   <  1:702976/57‑1 (MQ=255)
                                                                 tttCAACTTCAATAAAGCCGCGGTTCACCATGAACTGGCGAATACCAGAGAGGATCTGCGAGCGCACTTTa     <  1:749636/71‑1 (MQ=255)
                                                                  ttCAACTTCCATAAAGCCGCGGTTCACCATGAACTGGCGAATACCa                             <  1:1858755/46‑1 (MQ=255)
                                                                   tCAACTTCCATAAAGCCGCGGTTCACCATGAACTGGCGAATACCAGAGAGGAt                     >  1:619693/1‑53 (MQ=255)
                                                                   tCAACTTCCATAAAGCCGCGGTTCACCATGAACTGGCGAATACCAGAGAGGATCTGCGAg              >  1:2394206/1‑60 (MQ=255)
                                                                    cAACTTCCATAAAGCCGCGGTTCACCATGAACTGGCGAATACCAGAGAGGATCTGCGAGCGCACTTTAAAg  >  1:2238388/1‑71 (MQ=255)
                                                                     aaCTTCCATAAAGCCGCGGTTCACCATGAACTGGCGAATa                                >  1:2267633/1‑40 (MQ=255)
                                                                     aaCTTCCATAAAGCCGCGGTTCACCATGAACTGGCGAATACCAGAGAGGATCTGCGAgcgc           <  1:3488199/61‑1 (MQ=255)
                                                                     |                                                                     
CGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGTTCACCATGAACTGGCGAATACCAGAGAGGATCTGCGAGCGCACTTTAAAG  >  minE/1933756‑1933894

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: