Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,933,720 T→C 30.3% E241E (GAA→GAG lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,933,7200TC30.3% 155.4 / 91.9 122E241E (GAA→GAGlysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base T (38/47);  new base C (25/12);  total (63/59)
Fisher's exact test for biased strand distribution p-value = 2.96e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

Only 100 of 123 total aligned reads displayed.

GTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCG  >  minE/1933651‑1933789
                                                                     |                                                                     
attACGGAAGTTACGGTTGATTTCGAATACCCGTTCAAAACCGCCTAAACCAGACGTTTCAGATACAGCTc                                                                     >  1:2459265/2‑71 (MQ=255)
   aCGGAAGTTACGGTTGATTTCGAATACCCGTTCAAAACCGCCTACAACCAGACGTTTCAGATACAGCTcc                                                                    >  1:5351536/1‑70 (MQ=255)
   aCGGAAGTTACGGTTGATTTCGAATACCCGTTCAAAACCGCCTACAACCAGACGTTTCAGATACAGCTCCg                                                                   >  1:5008196/1‑71 (MQ=255)
   aCGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTcc                                                                    <  1:2260078/70‑1 (MQ=255)
   aCGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCg                                                                   >  1:3371058/1‑71 (MQ=255)
       aaGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGgcgc                                                               <  1:2875600/71‑1 (MQ=255)
         gTTACGGTTGATTTCGAATACCCGTTCAAAACCGCCTACAACCAGACGTTTCAGATACAGCTCCGGCGCGa                                                             <  1:4795269/71‑1 (MQ=255)
           tACGGTTGATTTCGAATACCCGTTCAAAACCGCCTACAACCAGACGTTTCAGATACAGCTCCGGCGCGATa                                                           >  1:2469760/1‑71 (MQ=255)
            aCGGTTGATTTCGAATACCCGTTCAAAACCGCCTACAACCAGACGTTTCAGATACAGCTCCGGCGCGATa                                                           >  1:4185953/1‑70 (MQ=255)
            aCGGTTGATTTCGAATACCCGTTCAAAACCGCCTACAACCAGACGTTTCAGATACAGCTCCGGCGCGATAc                                                          >  1:4862806/1‑71 (MQ=255)
            aCGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGgcgc                                                               >  1:5387034/1‑66 (MQ=255)
             cGGTTGATTTCGAATACCCGTTCAAAACCGCCTACAACCAGACGTTTCAGATACAGCTCCGGCGCGATACg                                                         >  1:4656395/1‑71 (MQ=255)
              ggTTGATTTCGAATACCCGTTCAAAACCGCCTACAACCAGACGTTCAGATACAGCTCCGGCGCGATACGCa                                                       <  1:201603/71‑1 (MQ=255)
              ggTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGc                                                        <  1:3248809/71‑1 (MQ=255)
              ggTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGc                                                        <  1:1084355/71‑1 (MQ=255)
              ggTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGc                                                        <  1:3931902/71‑1 (MQ=255)
              ggTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGc                                                        <  1:5289292/71‑1 (MQ=255)
                ttGATTTCGAATACCCGTTCAAAACCGCCTACAACCAGACGTTTCAGATACAGCTCCGGCGCGATACGCa                                                       <  1:665385/70‑1 (MQ=255)
                ttGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAg                                                      >  1:3631369/1‑71 (MQ=255)
                   aTTTCGAATACCCGTTCAAAACCGCCTACAACCAGACGTTTCAGATACAGCTCCGGCGCGATACGCAGATa                                                   <  1:4404712/71‑1 (MQ=255)
                     ttCGAATACCCGTTCAAAACCGCCTACAACCAGACGTTTCAGATACAGCTCCGGCGCGATACGCAGATACa                                                 >  1:2757169/1‑71 (MQ=255)
                      tCGAATACCCGTTCAAAACCGCCTACAACCAGACGTTTCAGATACAGCTCCGGCGCGATACGCAGATACa                                                 >  1:3019741/1‑70 (MQ=255)
                       cGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATg                                               >  1:2840604/1‑71 (MQ=255)
                             cacGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTc                                             >  1:5660354/1‑67 (MQ=255)
                              acGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGAt                                        <  1:2723919/71‑1 (MQ=255)
                               cGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTAc                                                  >  1:2664894/1‑60 (MQ=255)
                               cGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTAc                                                  >  1:1057842/1‑60 (MQ=255)
                               cGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTAc                                                  >  1:1303601/1‑60 (MQ=255)
                               cGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCgaga                                         <  1:3263218/69‑1 (MQ=255)
                               cGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATc                                       <  1:2306737/71‑1 (MQ=255)
                               cGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATc                                       <  1:5599401/71‑1 (MQ=255)
                                gCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATcc                                      <  1:3558803/71‑1 (MQ=255)
                                gCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATcc                                      <  1:925502/71‑1 (MQ=255)
                                 cTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGc                                                        >  1:2788999/1‑52 (MQ=255)
                                 cTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGc                                                        >  1:1380682/1‑52 (MQ=255)
                                 cTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGc                                                        >  1:1373156/1‑52 (MQ=255)
                                 cTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACaa                                                >  1:362631/1‑59 (MQ=255)
                                  tCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTc                                                                     <  1:1940562/38‑1 (MQ=255)
                                   cGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACa                                                 >  1:935990/1‑57 (MQ=255)
                                    gAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCg                                            <  1:446745/61‑1 (MQ=255)
                                      aGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGcgcg                                >  1:3754388/1‑71 (MQ=255)
                                       gCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGt                               >  1:1466483/1‑71 (MQ=255)
                                       gCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGt                               <  1:5160679/71‑1 (MQ=255)
                                        ccaccaACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCga                                           <  1:307660/58‑1 (MQ=255)
                                        ccaccaACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCga                                           <  1:2140913/58‑1 (MQ=255)
                                         caccaACCACCAGACGCTTGAGGTACAGTTCCGgcgc                                                               >  1:2939465/1‑37 (MQ=255)
                                         caccaACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTa                                                   >  1:252778/1‑49 (MQ=255)
                                          accaaccaCCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTAt                            <  1:1891689/71‑1 (MQ=255)
                                          accaaccaCCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTAt                            >  1:2121995/1‑71 (MQ=255)
                                             aaccaCCAGACGCTTGAGGTACAGTTCCGGCGCGa                                                             >  1:5087035/1‑35 (MQ=255)
                                             aaccaCCAGACGCTTGAGGTACAGTTCCGGCGCGATACg                                                         <  1:4778200/39‑1 (MQ=255)
                                             aaccaCCAGACGCTTGAGGTACAGTTCCGGCGCGATACg                                                         <  1:149809/39‑1 (MQ=255)
                                             aaccaCCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAgc                                   >  1:1710353/1‑61 (MQ=255)
                                             aaccaCCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTAtggt                         >  1:768292/1‑71 (MQ=255)
                                              accaccAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTAtggtg                        <  1:4447946/71‑1 (MQ=255)
                                              accaccAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTAtgg                          >  1:1444606/1‑69 (MQ=255)
                                              acaaccaGACGTTTCAGATACAGCTCCGGCGCGATACGCAGATACATATCTAAATCCAGAGCATTatgatg                        >  1:4548195/4‑71 (MQ=255)
                                              acaaccaGACGTTTCAGATACAGCTCCGGCGCGATACGCAGATACATATCTAAATCCAGAGCATTatgatg                        >  1:4321830/4‑71 (MQ=255)
                                               caaccaGACGTTTCAGATACAGCTCCGGCGCGATACGCAGATACATATCTAAATCCAGAGCATTAtgatgg                       >  1:1068133/3‑71 (MQ=255)
                                                caccaGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCga                                           >  1:1355163/1‑50 (MQ=255)
                                                aaccaGACGTTTCAGATACAGCTCCGGCGCGATACGCAGATACa                                                 >  1:287815/2‑44 (MQ=19)
                                                 accaGACGTTTCAGATACAGCTCCGGTGCGATACGCAGATACATATCTAAATCCAGAGCATTATgatggg                      >  1:415251/1‑70 (MQ=255)
                                                 accaGACGTTTCAGATACAGCTCCGGCGCGATACGCAGATACGTATCTAAATCCAGAGCATTATgatggg                      >  1:3721280/1‑70 (MQ=255)
                                                 accaGACGTTTCAGATACAGCTCCGGCGCGATACGCAGATACATATCTAAATCCAGAGCATTATgatggg                      >  1:829866/1‑70 (MQ=255)
                                                 accaGACGTTTCAGATACAGCTCCGGCGCGATACGCAGATACATATCTAAATCCAGAGCATTATgatggg                      >  1:3398683/1‑70 (MQ=255)
                                                 accaGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATggtgggt                     >  1:5677978/1‑71 (MQ=255)
                                                  ccaGACGTTTCAGATACAGCTCCGGCGCGATACGCAGATACATATCTAAATCCAGAGCATTATgatgggt                     >  1:739177/1‑70 (MQ=255)
                                                  ccaGACGTTTCAGATACAGCTCCGGCGCGATACGCAGATACATATCTAAATCCAGAGCATTATgatgggt                     >  1:1133008/1‑70 (MQ=255)
                                                  ccaGACGTTTCAGATACAGCTCCGGCGCGATACGCAGATACATATCTAAATCCAGAGCATTATgatggg                      >  1:4618509/1‑69 (MQ=255)
                                                    aGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGa                   >  1:1790977/1‑70 (MQ=255)
                                                     gACGTTTCAGATACAGCTCCGGCGCGATACGCAGATACATATCTAAATCCAGAGCATTATGATGGGTAAt                  >  1:3131719/1‑70 (MQ=255)
                                                     gACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATa                 <  1:3421319/71‑1 (MQ=255)
                                                      aCGTTTCAGATACAGCTCCGGCGCGATACGCAGATa                                                   <  1:3392935/36‑1 (MQ=16)
                                                      aCGTTTCAGATACAGCTCCGGCGCGATACGCAGATACATATCTAAATCCAGAGCATTATGATGGGtaataa                <  1:1714656/71‑1 (MQ=255)
                                                        gCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAAc              <  1:4374755/71‑1 (MQ=255)
                                                        gCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAAc              >  1:5113055/1‑71 (MQ=255)
                                                            gAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACgg            <  1:5619305/69‑1 (MQ=255)
                                                            gAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACgg            <  1:743625/69‑1 (MQ=255)
                                                             agatACAGCTCCGGCGCGATACGCAGATACATATCTAAATCCAGAGCATTATGATGGGTAATAAACGGGcg         <  1:4996120/68‑1 (MQ=255)
                                                               atACAGCTCCGGCGCGATACGCATATACATATCTAAATCCAGAGCATTAt                            <  1:280149/49‑1 (MQ=255)
                                                               atACAGCTCCGGCGCGATACGCAGATACATATCTAAATCCAGAGCATTAt                            <  1:4654202/49‑1 (MQ=255)
                                                               atACAGCTCCGGCGCGATACGCAGATACATATCTAAATCCAGAGCATTATGATGGGTAATAAACGGGcgag       >  1:2999489/2‑69 (MQ=255)
                                                               atACAGCTCCGGCGCGATACGCAGATACATATCTAAATCCAGAGCATTATGATGGGTAATAAACGGGcgag       >  1:5531915/2‑69 (MQ=255)
                                                               atACAGCTCCGGCGCGATACGCAGATACATATCTAAATCCAGAGCATTATGATGGGTAATAAACGGGcgag       >  1:4868777/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)
                                                                 aCAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTAtgg                          >  1:1514538/1‑50 (MQ=255)
                                                                 aCAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGAcgcg       >  1:4957127/1‑69 (MQ=255)
                                                                  cAGTTCCGGTGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCa     <  1:2054571/70‑1 (MQ=255)
                                                                  cAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCa                                     <  1:5547392/38‑1 (MQ=255)
                                                                  cAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCa                                     <  1:1493471/38‑1 (MQ=255)
                                                                  cAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATaaa               >  1:5344057/1‑60 (MQ=255)
                                                                  cAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAg    <  1:4187216/71‑1 (MQ=255)
                                                                  cAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAg    >  1:730371/1‑71 (MQ=255)
                                                                  cAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAg    <  1:2744063/71‑1 (MQ=255)
                                                                    gTTCCGGCGCGATACGCAGGTACATGTCGAGATCCa                                     <  1:5223067/36‑1 (MQ=255)
                                                                    gTTCCGGCGCGATACGCAGGTACATGTCGAGATCCa                                     <  1:2269693/36‑1 (MQ=255)
                                                                    gTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCg  <  1:1366105/71‑1 (MQ=255)
                                                                    gTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCg  >  1:2763614/1‑71 (MQ=255)
                                                                     |                                                                     
GTTACGGAAGTTACGGTTGATTTCGAATACACGCTCGAAGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCG  >  minE/1933651‑1933789

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: