Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,933,771 G→A 15.4% I224I (ATC→ATT lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,933,7710GA15.4% 216.4 / 26.3 91I224I (ATC→ATTlysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base G (34/43);  new base A (9/5);  total (43/48)
Fisher's exact test for biased strand distribution p-value = 2.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

AGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGC  >  minE/1933703‑1933841
                                                                    |                                                                      
aGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGa                                                                       >  1:1790977/1‑70 (MQ=255)
 gACGTTTCAGATACAGCTCCGGCGCGATACGCAGATACATATCTAAATCCAGAGCATTATGATGGGTAAt                                                                      >  1:3131719/1‑70 (MQ=255)
 gACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATa                                                                     <  1:585195/71‑1 (MQ=255)
 gACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATa                                                                     <  1:3421319/71‑1 (MQ=255)
  aCGTTTCAGATACAGCTCCGGCGCGATACGCAGATACATATCTAAATCCAGAGCATTATGATGGGtaataa                                                                    <  1:1714656/71‑1 (MQ=255)
    gCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAAc                                                                  <  1:4374755/71‑1 (MQ=255)
    gCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAAc                                                                  >  1:5113055/1‑71 (MQ=255)
        gAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACgg                                                                <  1:743625/69‑1 (MQ=255)
        gAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACgg                                                                <  1:4481914/69‑1 (MQ=255)
        gAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACgg                                                                <  1:5619305/69‑1 (MQ=255)
         agatACAGCTCCGGCGCGATACGCAGATACATATCTAAATCCAGAGCATTATGATGGGTAATAAACGGGcg                                                             <  1:4996120/68‑1 (MQ=255)
           atACAGCTCCGGCGCGATACGCAGATACATATCTAAATCCAGAGCATTATGATGGGTAATAAACGGGcgag                                                           >  1:2999489/2‑69 (MQ=255)
           atACAGCTCCGGCGCGATACGCAGATACATATCTAAATCCAGAGCATTATGATGGGTAATAAACGGGcgag                                                           >  1:4868777/2‑69 (MQ=255)
           atACAGCTCCGGCGCGATACGCAGATACATATCTAAATCCAGAGCATTATGATGGGTAATAAACGGGcgag                                                           >  1:5531915/2‑69 (MQ=255)
            tACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGAcgcg                                                           >  1:5297841/1‑70 (MQ=255)
            tACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGAcgcg                                                           >  1:1524639/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)
              cAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATaaa                                                                   >  1:5344057/1‑60 (MQ=255)
              cAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAg                                                        <  1:4187216/71‑1 (MQ=255)
              cAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAg                                                        <  1:2744063/71‑1 (MQ=255)
              cAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAg                                                        >  1:730371/1‑71 (MQ=255)
              cAGCTCCGGCGCGATACGCAGATACATATCTAAATCCAGAGCATTATGATGGGTAATAAACGGGCGAGCa                                                         <  1:3098146/70‑1 (MQ=255)
              cAGCTCCGGCGCGATACGCAGATACATATCTAAATCCAGAGCATTATGATGGGTAATAAACGGGCGAGCa                                                         <  1:1520136/70‑1 (MQ=255)
                gTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCg                                                      >  1:2763614/1‑71 (MQ=255)
                gTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCg                                                      <  1:1366105/71‑1 (MQ=255)
                 ttCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATa                                                                     <  1:3277419/55‑1 (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)
                           aTACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATa                                                                     <  1:1623944/45‑1 (MQ=255)
                           aTACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACg                                                                 <  1:886589/49‑1 (MQ=255)
                           aTACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACg                                                                 <  1:2549212/49‑1 (MQ=255)
                           aTACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACgcgc                                                          >  1:2779726/1‑56 (MQ=255)
                           aTACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAgg                                           >  1:4463565/1‑71 (MQ=255)
                           aTACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAgg                                           >  1:4510568/1‑71 (MQ=255)
                           aTACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAgg                                           >  1:3536593/1‑71 (MQ=255)
                           aTACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAgg                                           >  1:2634568/1‑71 (MQ=255)
                            tACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGa                                                                       >  1:5393793/1‑42 (MQ=255)
                            tACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGa                                                                       >  1:4060995/1‑42 (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)
                                  ggTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACg                                                                 <  1:167878/42‑1 (MQ=255)
                                   gTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCt                                   <  1:2411358/71‑1 (MQ=255)
                                    tACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACgcgc                                                          <  1:4856365/47‑1 (MQ=255)
                                    tACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTg                                  >  1:3909127/1‑71 (MQ=255)
                                      cATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAAc                                                                  <  1:196368/37‑1 (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:4157284/71‑1 (MQ=255)
                                                       cgTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACt               >  1:3994413/1‑71 (MQ=255)
                                                        gTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACtt              >  1:4633813/1‑71 (MQ=255)
                                                         ttATGATGGGTAATAAACGGGCGAGCAGATGCCCCACCTGGAATTACCTGCATCATCGGGGTTTCTACTTc             >  1:4341644/1‑71 (MQ=255)
                                                         ttATGATGGGTAATAAACGGGCGAGCAGATGCCCCACCTGGAATTACCTGCATCATCGGGGTTTCTACTTc             >  1:1266663/1‑71 (MQ=255)
                                                          tATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTc             >  1:2907183/1‑70 (MQ=255)
                                                          tATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTc             >  1:2900545/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:1764973/71‑1 (MQ=255)
                                                                 ggtgATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGc     <  1:2354082/71‑1 (MQ=255)
                                                                 ggtaATAAACGGGCGAGCAGATGCCCCACCTGGAATTACCTGCATCATCGGGGTTTCTACTTCCATAAAg      <  1:1906686/70‑1 (MQ=255)
                                                                    gATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCcgc  >  1:3084545/1‑71 (MQ=255)
                                                                    |                                                                      
AGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGC  >  minE/1933703‑1933841

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: