Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,933,249 A→G 46.4% S398S (TCT→TCC lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,933,2490AG46.4% 20.8 / 105.0 84S398S (TCT→TCClysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base A (17/28);  new base G (22/17);  total (39/45)
Fisher's exact test for biased strand distribution p-value = 1.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGACCACCAATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCG  >  minE/1933181‑1933315
                                                                    |                                                                  
cGACCACCAATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGGCAGCGgag                                                                   >  1:4861846/1‑70 (MQ=255)
cGACCACCAATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGgag                                                                   >  1:2485440/1‑70 (MQ=255)
 gACCACCGATGAAGAATTCAAAACGGTCGGTGATTTCCGGGTTAACATCATTACGGCTTGCCAGCGgggac                                                                 <  1:4978718/71‑2 (MQ=255)
   ccaccaATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAAc                                                               >  1:1245625/1‑71 (MQ=255)
   ccaccaATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAAc                                                               >  1:1928819/1‑71 (MQ=255)
     accaATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGaa                                                                 <  1:5295457/67‑1 (MQ=255)
        gaTGAAGAATTCAAAACGGTCGGTGATTTCCGGGTTAACATCATTACGGCGTGCCAGCGGGGACACTTCTg                                                          <  1:5022301/70‑1 (MQ=255)
        actGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCt                                                           <  1:5183383/68‑1 (MQ=255)
         aTGAAGAATTCAAAACGGGCGGTGATTTCCGGGTTAACATCATTACGGCGTGCCAGCGGGGACACTTCTGc                                                         <  1:162148/71‑1 (MQ=255)
          tGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCt                                                           <  1:3702050/68‑1 (MQ=255)
           gaagaaCTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAAc                                                               <  1:1030113/63‑1 (MQ=255)
           gaagaaCTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAAc                                                               <  1:3077125/63‑1 (MQ=255)
            aagaaTTCAAAACGGTCGGTGATTTCCGGGTTAACATCATTACGGCGTGCCAGCGGGGACACTTCTGCCg                                                       <  1:4965131/70‑1 (MQ=255)
            aagaaCTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCgg                                                      <  1:3482038/71‑1 (MQ=255)
            aagaaCTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCg                                                       >  1:2914009/1‑70 (MQ=255)
             agaaCTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAAc                                                               <  1:4604216/61‑1 (MQ=255)
              gaaTTCAAAACGGTCGGTGATTTCCGGGTTAACATCATTACGGCGTGCCAGCGGGGACACTTCTGCCGGa                                                     >  1:2674907/1‑70 (MQ=255)
                  tCAAAACGGTCGGTGATTTCCGGGTTAACATCATTACGGCGTGCCAGCGGGGACACTTCTGCCg                                                       >  1:2266241/1‑64 (MQ=255)
                  tCAAAACGGTCGGTGATTTCCGGGTTAACATCATTACGGCGTGCCAGCGGGGACACTTCTGCCg                                                       >  1:4740757/1‑64 (MQ=255)
                  tCAAAACGGTCGGTGATTTCCGGGTTAACATCATTACGGCGTGCCAGCGGGGACACTTCTGCCGGATATTc                                                >  1:5697636/1‑71 (MQ=255)
                   cAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATAtt                                                 >  1:4524972/1‑69 (MQ=255)
                    aaaaCGGTCGGTGATTTCCGGGTTAACATCATTACGGCGTGCCAGCGGGGACACTTCTGCCGGATATTccg                                              >  1:4684276/1‑69 (MQ=255)
                    aaaGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGa                                                     >  1:1058383/1‑64 (MQ=255)
                      aacggtcggtGATTTCCGGGTTAACATCATTACGGCGTGCCAGCGGGGACACTTCTGCCGGATATTCCGTa                                            <  1:5613074/69‑1 (MQ=255)
                      aacggtcggtGATTTCCGGGTTAACATCATTACGGCGTGCCAGCGGGGACACTTCTGCCGGATATTCCGTa                                            <  1:2695751/69‑1 (MQ=255)
                      aGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTa                                            <  1:4494449/71‑1 (MQ=255)
                      aGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTa                                            <  1:5057228/71‑1 (MQ=255)
                       gCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTaa                                           >  1:2623724/1‑71 (MQ=255)
                        cggtcggtGATTTCCGGGTTAACATCATTACGGCGTGCCAGCGggga                                                                  >  1:3525716/1‑47 (MQ=19)
                        cGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAAt                                          <  1:2902075/71‑1 (MQ=255)
                         ggTCGGTGATTTCCGGGTTAACATCATTACGGCGTGCCAGCGGGGACACTTCTGCCGGata                                                   >  1:3379239/1‑61 (MQ=255)
                             tgtgATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACtt                                                             <  1:3115275/47‑1 (MQ=255)
                             ggtgATTTCCGGGTTAACATCATTACGGCGTGCCAGCGGGGACACTTCTGCCGGa                                                     >  1:546466/2‑55 (MQ=255)
                              gtgATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGt                                    >  1:4739304/1‑71 (MQ=255)
                              gtgATTTCCGGGTTAACATCATTACGGCGTGCCAGCGGGGACACTTCTGcc                                                        <  1:3685336/51‑1 (MQ=255)
                              gtgATTTCCGGGTTAACATCATTACGGCGTGCCAGCGGGGACACTTCTGcc                                                        <  1:4504936/51‑1 (MQ=255)
                               tgATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGatat                                                  <  1:2685242/56‑1 (MQ=255)
                               tgATTTCCGGGTTAACGTCGTTACGACGCGCAAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTc                                   >  1:4597413/1‑71 (MQ=255)
                               tgATTTCCGGGTTAACATCATTACGGCGTGCCAGCGGGGACACTTCTGCCGGATATTCCGTAATAAAGGtt                                   <  1:310383/71‑2 (MQ=255)
                                 aTTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACtt                                                             >  1:1742510/1‑43 (MQ=255)
                                     ccGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAAc                                                               <  1:645268/37‑1 (MQ=255)
                                      cGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGaa                            >  1:239370/1‑71 (MQ=255)
                                       gggTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGa                             <  1:3645299/69‑1 (MQ=255)
                                       gggTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGa                             <  1:2620717/69‑1 (MQ=255)
                                         gTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAAt                           >  1:5534655/1‑69 (MQ=255)
                                         gTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCa                         <  1:4157015/71‑1 (MQ=255)
                                         gTTAACATCATTACGGCGTGCCAGCGGGGACACTTCTGCCGGATATTCCGTAATAAAGGTTGGCTGAATCa                         <  1:5486266/71‑1 (MQ=255)
                                          ttAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAg                        <  1:5676097/71‑1 (MQ=255)
                                          ttAACATCATTACGGCGTGCCAGCGGGGACACTTCTGCCGGATATTCCGTAATAAAGGTTGGCTGAATCAg                        <  1:2243293/71‑1 (MQ=255)
                                            aaCGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGa                       <  1:1599594/70‑1 (MQ=255)
                                            aaCATCATTACGGCGTGCCAGCGGGGACACTTCTGCCGGATATTCCGTAATAAAGGTTGGCTGaa                            >  1:1530850/1‑65 (MQ=255)
                                            aaCATCATTACGGCGTGCCAGCGGGGACACTTCTGCCGGATATTCCGTAATAAAGGTTGGCTGaa                            >  1:5617751/1‑65 (MQ=255)
                                            aaCATCATTACGGCGTGCCAGCGGGGACACTTCTGCCGGATATTCCGTAATAAAGGTTGGCTGAATCAGAt                      >  1:4055516/1‑71 (MQ=255)
                                              cgtcgtTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATg                                         <  1:3767432/50‑1 (MQ=255)
                                              cgtcgtTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATg                                         <  1:2597515/50‑1 (MQ=255)
                                                tcgtTACGACGCGCCAGCGGAGAAACTTCTGCCGGa                                                     <  1:3596984/36‑1 (MQ=255)
                                                tcattACGGCGTGCCAGCGGGGACACTTCTGCCGGATATTCCGTAATAAAGGTTGGCTGAATCAGa                       >  1:5360930/4‑66 (MQ=255)
                                                tcattACGGCGTGCCAGCGGGGACACTTCTGCCGGATATTCCGTAATAAAGGTTGGCTGAATCAGAtg                     <  1:2922824/65‑1 (MQ=255)
                                                  gtTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCtt                >  1:2449454/1‑71 (MQ=255)
                                                  gtTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCtt                <  1:3358170/71‑1 (MQ=255)
                                                  attACGGCGTGCCAGCGGGGACACTTCTGCCGGATATTCCGTAATAAAg                                      >  1:580606/2‑49 (MQ=19)
                                                  attACGGCGTGCCAGCGGGGACACTTCTGCCGGATATTCCGTAATAAAg                                      >  1:4785117/2‑49 (MQ=19)
                                                  attACGGCGTGCCAGCGGGGACACTTCTGCCGGATATTCCGTAATAAAGGtt                                   >  1:4030066/2‑51 (MQ=19)
                                                  attACGGCGTGCCAGCGGGGACACTTCTGCCGGATATTCCGTAATAAAGGtt                                   >  1:2834892/2‑51 (MQ=19)
                                                  attACGGCGTGCCAGCGGGGACACTTCTGCCGGATATTCCGTAATAAAGGTTGGCTGAATCAGATGTGCtt                <  1:187248/70‑1 (MQ=255)
                                                   ttACGGCGTGCCAGCGGGGACACTTCTGCCGGATATTCCGTAATaaa                                       <  1:2350342/47‑1 (MQ=19)
                                                   ttACGGCGTGCCAGCGGGGACACTTCTGCCGGATATTCCGTAATAAAGGttg                                  <  1:3570228/52‑3 (MQ=19)
                                                   ttACGGCGTGCCAGCGGGGACACTTCTGCCGGATATTCCGTAATAAAGGTTGGCTGAATCAGATGTGCTTc               >  1:2088010/1‑71 (MQ=255)
                                                    tACGGCGTGCCAGCGGGGACACTTCTGCCGGATATTCCGTAATAAAGGTTGGCTGAATCa                         <  1:4351435/60‑1 (MQ=255)
                                                    tACGGCGTGCCAGCGGGGACACTTCTGCCGGATATTCCGTAATAAAGGTTGGCTGAATCa                         <  1:40475/60‑1 (MQ=255)
                                                     acggcgTGCCAGCGGGGACACTTCTGCCGGATATTCCGTAATAAAGGttg                                  >  1:2387713/1‑48 (MQ=19)
                                                     acggcgTGCCAGCGGGGACACTTCTGCCGGATATTCCGTAATAAAGGttg                                  >  1:638248/1‑48 (MQ=19)
                                                     acggcgTGCCAGCGGGGACACTTCTGCCGGATATTCCGTAATAAAGGTTGGCTGAATCAGATGTGCTTCTg             >  1:3141325/1‑71 (MQ=255)
                                                       gacgCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGcc           <  1:3442001/71‑1 (MQ=255)
                                                         cgcgCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCAc         <  1:2829183/71‑1 (MQ=255)
                                                          gcgcCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAAt                           <  1:5121648/52‑1 (MQ=255)
                                                            gcCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCActtc      >  1:4819275/1‑71 (MQ=255)
                                                             ccAGCGGGGACACTTCTGCCGGATATTCCGTAATAAAGGTTGGCTGAATCAGATGTGCTTCTGCCActtca     >  1:5224330/1‑70 (MQ=255)
                                                             ccAGCGGGGACACTTCTGCCGGAGATTCCGTAATAAAGGTTGGCTg                              >  1:2226820/1‑46 (MQ=255)
                                                             ccAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCActtct     >  1:243279/1‑71 (MQ=255)
                                                              cAGCGGGGACACTTCTGCCGGATATTCCGTAATAAAGGTTGGCTGAATCAGATGTGCtt                <  1:366448/59‑1 (MQ=255)
                                                                gCGGGGACACTTCTGCCGGATATTCCGTAATAAAGGTTGGCTGAATCAGATGTGCTTCTGCCActtcata   >  1:4625207/1‑67 (MQ=255)
                                                                gCGGAGAAACTTCTGCCGGATATTCAGTAATGAAgg                                     >  1:2318334/1‑36 (MQ=255)
                                                                gCGGAGAAACTTCTGCCGGATATTCAGTAATGAAgg                                     >  1:3758930/1‑36 (MQ=255)
                                                                gCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACttcttc   >  1:137987/1‑70 (MQ=255)
                                                                gCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCg  <  1:4607680/71‑1 (MQ=255)
                                                                gCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCg  <  1:934750/71‑1 (MQ=255)
                                                                 cGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGcc           <  1:336913/61‑1 (MQ=255)
                                                                   gagaAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCAc         >  1:2806325/1‑61 (MQ=255)
                                                                   gagaAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCAc         >  1:3008522/1‑61 (MQ=255)
                                                                    |                                                                  
CGACCACCAATGAAGAACTCAAAGCGGTCTGTGATTTCCGGGTTAACGTCGTTACGACGCGCCAGCGGAGAAACTTCTGCCGGATATTCAGTAATGAAGGTCGGCTGAATCAGATGTGCTTCTGCCACTTCTTCG  >  minE/1933181‑1933315

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: