Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,933,780 A→G 21.6% R221R (CGT→CGC lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,933,7800AG21.6% 167.3 / 41.2 88R221R (CGT→CGClysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base A (32/37);  new base G (12/7);  total (44/44)
Fisher's exact test for biased strand distribution p-value = 3.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.45e-01

AGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGTTCAC  >  minE/1933712‑1933848
                                                                    |                                                                    
agatACAGCTCCGGCGCGATACGCAGATACATATCTAAATCCAGAGCATTATGATGGGTAATAAACGGGcg                                                                    <  1:4996120/68‑1 (MQ=255)
  atACAGCTCCGGCGCGATACGCAGATACATATCTAAATCCAGAGCATTATGATGGGTAATAAACGGGcgag                                                                  >  1:5531915/2‑69 (MQ=255)
  atACAGCTCCGGCGCGATACGCAGATACATATCTAAATCCAGAGCATTATGATGGGTAATAAACGGGcgag                                                                  >  1:2999489/2‑69 (MQ=255)
  atACAGCTCCGGCGCGATACGCAGATACATATCTAAATCCAGAGCATTATGATGGGTAATAAACGGGcgag                                                                  >  1:4868777/2‑69 (MQ=255)
   tACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGAcgcg                                                                  >  1:1524639/1‑70 (MQ=255)
   tACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGAcgcg                                                                  >  1:5297841/1‑70 (MQ=255)
   tACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACgcgc                                                                 <  1:3213693/71‑1 (MQ=255)
    aCAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGAcgcg                                                                  >  1:4957127/1‑69 (MQ=255)
     cAGTTCCGGTGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCa                                                                <  1:2054571/70‑1 (MQ=255)
     cAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAg                                                               >  1:730371/1‑71 (MQ=255)
     cAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAg                                                               <  1:2744063/71‑1 (MQ=255)
     cAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAg                                                               <  1:4187216/71‑1 (MQ=255)
     cAGCTCCGGCGCGATACGCAGATACATATCTAAATCCAGAGCATTATGATGGGTAATAAACGGGCGAGCa                                                                <  1:1520136/70‑1 (MQ=255)
     cAGCTCCGGCGCGATACGCAGATACATATCTAAATCCAGAGCATTATGATGGGTAATAAACGGGCGAGCa                                                                <  1:3098146/70‑1 (MQ=255)
       gTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCg                                                             <  1:1366105/71‑1 (MQ=255)
       gTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCg                                                             >  1:2763614/1‑71 (MQ=255)
        ctCCGGCGCGATACGCAGATACATATCTAAATCCAGAGCATTATGATGGGTAATAAACGGGCGAGCa                                                                <  1:2434266/66‑1 (MQ=255)
         tCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCgg                                                            >  1:1682187/1‑70 (MQ=255)
         tCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCgg                                                            >  1:3631966/1‑70 (MQ=255)
          ccGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGc                                                           <  1:3232740/70‑1 (MQ=255)
          ccGGCGCGATACGCAGATACATATCTAAATCCAGAGCATTATGATGGGTAATAAACGGGCGAGCAga                                                              >  1:2668014/1‑66 (MQ=255)
           cGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCAc                                                         >  1:37593/1‑71 (MQ=255)
           cGGCGCGATACGCAGATACATATCTAAATCCAGAGCATTATGATGGGTAATAAACGGGcga                                                                   <  1:4585500/61‑2 (MQ=255)
             gcgcGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGAc                                                                     <  1:2283399/57‑1 (MQ=255)
              cgcgATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGCTAAACGGACGCGCAGCGGCACCGc                                                      <  1:5365137/71‑1 (MQ=255)
               gcgATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGcc                                                     <  1:993792/71‑1 (MQ=255)
                 gATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAg                                                   <  1:1173550/71‑1 (MQ=255)
                 gATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAg                                                   <  1:1631914/71‑1 (MQ=255)
                 gATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAg                                                   <  1:1766010/71‑1 (MQ=255)
                  aTACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACgcgc                                                                 >  1:2779726/1‑56 (MQ=255)
                  aTACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAgg                                                  >  1:2634568/1‑71 (MQ=255)
                  aTACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAgg                                                  >  1:4510568/1‑71 (MQ=255)
                  aTACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAgg                                                  >  1:4463565/1‑71 (MQ=255)
                  aTACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAgg                                                  >  1:3536593/1‑71 (MQ=255)
                   tACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGcc                                                     <  1:3841223/67‑1 (MQ=255)
                   tACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAggg                                                 <  1:182418/71‑1 (MQ=255)
                     cGCAGGTACATGTCGAGATCCAGCGCGTTATTGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGa                                                <  1:827530/70‑1 (MQ=255)
                     cGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCa                                                    >  1:2929089/1‑66 (MQ=255)
                     cGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGa                                                <  1:3943651/70‑1 (MQ=255)
                      gCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCa                                                          <  1:211022/59‑1 (MQ=255)
                       cAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCa                                             <  1:2229378/71‑1 (MQ=255)
                        aGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCAc                                            <  1:1902217/71‑1 (MQ=255)
                          gTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCt                                          <  1:2411358/71‑1 (MQ=255)
                           tACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACgcgc                                                                 <  1:4856365/47‑1 (MQ=255)
                           tACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTg                                         >  1:3909127/1‑71 (MQ=255)
                              aTGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGcat                                      >  1:3289528/1‑71 (MQ=255)
                              aTGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGca                                       >  1:5083439/1‑70 (MQ=255)
                               tGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGAcgcg                                                                  <  1:1338097/42‑1 (MQ=255)
                               tGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGcatc                                     <  1:2953586/71‑1 (MQ=255)
                                 tCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGcc                                                     <  1:2705047/53‑1 (MQ=255)
                                 tCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGcatca                                    <  1:3846328/70‑1 (MQ=255)
                                   gagaTCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCt                                          >  1:2118789/1‑62 (MQ=255)
                                   gagaTCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCg                                 >  1:1910079/1‑71 (MQ=255)
                                    agaTCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGc                                                           >  1:755035/1‑44 (MQ=255)
                                    agaTCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGc                                                           >  1:5679920/1‑44 (MQ=255)
                                      aTCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGc                               >  1:218407/1‑70 (MQ=255)
                                        ccAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCAcc                                           >  1:4806698/1‑56 (MQ=255)
                                          aGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTc                          <  1:299525/71‑1 (MQ=255)
                                           gcgcGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTc                          >  1:1235064/1‑70 (MQ=255)
                                           gcgcGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCa                         <  1:2873253/71‑1 (MQ=255)
                                             gcgTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCaa                        <  1:4238244/70‑1 (MQ=255)
                                              cgTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACt                      >  1:3994413/1‑71 (MQ=255)
                                              cgTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACt                      <  1:4157284/71‑1 (MQ=255)
                                               gTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACtt                     >  1:4633813/1‑71 (MQ=255)
                                                ttATGATGGGTAATAAACGGGCGAGCAGATGCCCCACCTGGAATTACCTGCATCATCGGGGTTTCTACTTc                    >  1:1266663/1‑71 (MQ=255)
                                                ttATGATGGGTAATAAACGGGCGAGCAGATGCCCCACCTGGAATTACCTGCATCATCGGGGTTTCTACTTc                    >  1:4341644/1‑71 (MQ=255)
                                                 tATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTc                    >  1:2900545/1‑70 (MQ=255)
                                                 tATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTc                    >  1:2907183/1‑70 (MQ=255)
                                                  atgatgGGTAATAAACGGGCGAGCAGATGCCCCACCTGGAATTACCTGCATCATCGGGGTTTCTACTTCCa                  >  1:5466831/1‑71 (MQ=255)
                                                  atgatgGGTAATAAACGGGCGAGCAGATGCCCCACCTGGAATTACCTGCATCATCGGGGTTTCTACTTCCa                  >  1:1629431/1‑71 (MQ=255)
                                                  aTGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCa                  >  1:2968749/1‑71 (MQ=255)
                                                   tggtggGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCa                  <  1:768147/70‑1 (MQ=255)
                                                   tggtggGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCAt                 <  1:5063368/71‑1 (MQ=255)
                                                    gatgggtaATAAACGGGCGAGCAGATGCCCCACCTGGAATTACCTGCATCATCGGGGTTTCTACTACCAt                 >  1:5275157/3‑70 (MQ=255)
                                                     gtgggtgATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATaa               <  1:2920968/71‑1 (MQ=255)
                                                        ggtgATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGc            <  1:2354082/71‑1 (MQ=255)
                                                        ggtgATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGc            <  1:1764973/71‑1 (MQ=255)
                                                        ggtaATAAACGGGCGAGCAGATGCCCCACCTGGAATTACCTGCATCATCGGGGTTTCTACTTCCATAAAg             <  1:1906686/70‑1 (MQ=255)
                                                         gtaaTAAACGGGCGAGCAGATGCCCCACCTGGAATTACCTGCATCATCGGGGTTTCTACTTCCATAAAGcc           <  1:4942852/68‑1 (MQ=255)
                                                           gATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCcgc         >  1:3084545/1‑71 (MQ=255)
                                                            aTAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCcgcg        >  1:1694715/1‑71 (MQ=255)
                                                             tAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCgg       <  1:4098648/71‑1 (MQ=255)
                                                              aaaCGGGCGAGCAGATGCCCCACCTGGAATTACCTGCATCATCGGGGTTTCTACTTCCATaaa              >  1:810487/1‑63 (MQ=12)
                                                              aaaCGGGCGAGCAGATGCCCCACCTGGAATTACCTGCATCATCGGGGTTTCTACTTCCATAAAGCcgcgc       >  1:980354/1‑69 (MQ=255)
                                                               aaCGGGCGAGCAGATGCCCCACCTGGAATTACCTGCATCATCGGGGTTTCTACTTCCATAAAGCcgcgcgc     >  1:2829337/1‑68 (MQ=255)
                                                                  ggACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCcg          <  1:1974073/63‑1 (MQ=255)
                                                                  ggACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGTTCAc  <  1:3778846/71‑1 (MQ=255)
                                                                   gACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGTTCAc  >  1:5549767/1‑70 (MQ=255)
                                                                    |                                                                    
AGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGTTCAC  >  minE/1933712‑1933848

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: