Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,932,982 G→A 12.7% T487T (ACC→ACT lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,932,9820GA12.7% 226.6 / 19.9 87T487T (ACC→ACTlysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base G (44/31);  new base A (5/6);  total (49/38)
Fisher's exact test for biased strand distribution p-value = 5.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.96e-01

AACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTG  >  minE/1932915‑1933052
                                                                   |                                                                      
aaCGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGt                                                                      >  1:2621741/1‑70 (MQ=255)
 aCGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTg                                                                     <  1:2531095/70‑1 (MQ=255)
  cGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGaa                                                                   <  1:2635929/71‑1 (MQ=255)
     aaTGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAg                                                                >  1:645167/1‑71 (MQ=255)
      aTGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAg                                                                >  1:3127157/1‑70 (MQ=255)
      aTGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAg                                                                >  1:3537682/1‑70 (MQ=255)
          ttATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTa                                                           >  1:1149017/1‑71 (MQ=255)
             ttttACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTAcc                                                         <  1:1169249/70‑1 (MQ=255)
             ttttACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCa                                                        >  1:878366/1‑71 (MQ=255)
             ttttACCGGACGCATCGCCGCGAACAGCATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCa                                                        >  1:2314682/1‑71 (MQ=255)
               ttACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCAt                                                             >  1:3650037/1‑64 (MQ=255)
                tACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATa                                                      <  1:1280495/70‑1 (MQ=255)
                tACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATAc                                                     >  1:2324401/1‑71 (MQ=255)
                  ccGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCGGTTGGTGAACAGCATTACCATACg                                                    <  1:3049851/70‑1 (MQ=255)
                    ggACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTc                                                 <  1:1995375/71‑1 (MQ=255)
                    ggACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTc                                                 <  1:142934/71‑1 (MQ=255)
                    ggACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTc                                                 <  1:3554281/71‑1 (MQ=255)
                     gACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTc                                                 >  1:2499762/1‑70 (MQ=255)
                        gCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCAt                                                             >  1:3378994/1‑55 (MQ=255)
                        gCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCAt                                                             >  1:1704730/1‑55 (MQ=255)
                         cATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTg                                                                        <  1:2212840/43‑1 (MQ=255)
                         cATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTg                                                                        <  1:255271/43‑1 (MQ=255)
                         cATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTg                                                                        <  1:2511229/43‑1 (MQ=255)
                          aTCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCa                                                        >  1:342351/1‑58 (MQ=255)
                          aTCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATAc                                            >  1:1589153/1‑70 (MQ=255)
                          aTCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATAc                                            >  1:399179/1‑70 (MQ=255)
                            cGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTgg                                                                       >  1:1545448/1‑41 (MQ=255)
                             gCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTaaa                                                                   >  1:2507471/1‑44 (MQ=21)
                             gCCGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCaa                                         >  1:389074/1‑70 (MQ=255)
                              ccGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATa                                             <  1:2249842/65‑1 (MQ=255)
                              ccGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAAta                                       <  1:1387639/71‑2 (MQ=255)
                              ccGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAAt                                        >  1:2309855/1‑70 (MQ=255)
                               cGGGAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAAta                                       <  1:3545116/70‑2 (MQ=255)
                                gggAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATAcc                                     <  1:2036599/71‑1 (MQ=255)
                                 ggAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCAt                                                       <  1:2573551/52‑1 (MQ=255)
                                 ggAACAGAATAACGTCGCGGATGGTATGGCTGTTTGTGat                                                                   <  1:113343/40‑2 (MQ=39)
                                 ggAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGa                                               <  1:957812/60‑1 (MQ=255)
                                 ggAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTcc                                     >  1:2420740/1‑70 (MQ=255)
                                  gAAGAGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCa                                   <  1:3836490/71‑1 (MQ=255)
                                  gAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTcc                                     >  1:446491/1‑69 (MQ=255)
                                  gAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTcc                                     >  1:51047/1‑69 (MQ=255)
                                   aaCAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTcc                                     >  1:2017540/1‑68 (MQ=255)
                                    aCAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGa                                 >  1:2874504/1‑71 (MQ=255)
                                     gaGAATAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGAc                                >  1:2195732/2‑71 (MQ=255)
                                      aGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCAtt                                                            <  1:3511196/42‑1 (MQ=255)
                                        aaTAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTa                                                           >  1:1159958/1‑41 (MQ=255)
                                        aaTAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATa                                                      >  1:2073444/1‑46 (MQ=255)
                                        aaTAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATa                                                      >  1:2444327/1‑46 (MQ=255)
                                        aaTAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACctg                             >  1:3221528/1‑71 (MQ=255)
                                        aaTAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACctg                             >  1:729160/1‑71 (MQ=255)
                                        aaTAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACctg                             <  1:2702998/71‑1 (MQ=255)
                                        aaTAACGTCGCGGATGGTATGGCTGTTGGTGAACAGAATTACCATa                                                      >  1:2595537/1‑46 (MQ=255)
                                        aaTAACGTCGCGAATAGTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGAcca                              >  1:439714/1‑69 (MQ=255)
                                         aTAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACctgc                            >  1:265374/1‑71 (MQ=255)
                                          tAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACctgct                           <  1:2236446/71‑1 (MQ=255)
                                           aaCGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACctgctg                          >  1:881339/1‑71 (MQ=255)
                                            aCGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGt                         <  1:2107835/71‑1 (MQ=255)
                                            aCGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGt                         <  1:1883498/71‑1 (MQ=255)
                                              gTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTc                                                 >  1:3812308/1‑45 (MQ=255)
                                              gTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACctgct                           <  1:1849638/67‑1 (MQ=255)
                                               tCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGAc                                <  1:1516799/61‑1 (MQ=255)
                                               tCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACctgct                           <  1:2579557/66‑1 (MQ=255)
                                               tCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTcg                       >  1:3709333/1‑70 (MQ=255)
                                                  cgGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTcggcgg                   <  1:949646/71‑1 (MQ=255)
                                                   ggATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACctgct                           <  1:596495/62‑1 (MQ=255)
                                                   ggATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTcggcgg                   >  1:891360/1‑70 (MQ=255)
                                                     atagTATGGCTGTTAGTAAACAGCATAATCATTCGGTCGATACCAATACCCAGACCAGCGGTtggcgg                   >  1:3952323/4‑68 (MQ=255)
                                                     aTGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGt                  >  1:622335/1‑69 (MQ=255)
                                                      tGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAg               <  1:3250687/71‑1 (MQ=255)
                                                       ggTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGc              >  1:143374/1‑71 (MQ=255)
                                                       ggTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGc              >  1:1002513/1‑71 (MQ=255)
                                                       ggTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGc              >  1:2167885/1‑71 (MQ=255)
                                                          aTGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTcgg                      <  1:1447979/60‑1 (MQ=255)
                                                          aTGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTcgg                      <  1:1840271/60‑1 (MQ=255)
                                                           tGGCTGTTGGTGAACAGCATTATCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATg          >  1:1916970/1‑71 (MQ=255)
                                                           tGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTcgg                      >  1:2915610/1‑59 (MQ=255)
                                                           tGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATg          >  1:205013/1‑71 (MQ=255)
                                                           tGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCACAGACCTGCTGTCGGCGGTAAGCCATg          >  1:1699821/1‑71 (MQ=255)
                                                            ggCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACctgct                           <  1:2707742/53‑1 (MQ=255)
                                                            ggCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCa            <  1:1754461/68‑1 (MQ=255)
                                                             gCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGcc             >  1:1147443/1‑66 (MQ=255)
                                                             gCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGt         <  1:1736708/70‑1 (MQ=255)
                                                              cTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTc       <  1:1806520/71‑1 (MQ=255)
                                                              cTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTc       <  1:1157779/71‑1 (MQ=255)
                                                               tGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCAt           >  1:500213/1‑66 (MQ=255)
                                                               tGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTc       >  1:608211/1‑70 (MQ=255)
                                                               tGTTGGTGAACAGCATTACCATACCGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCAt           >  1:389763/1‑66 (MQ=255)
                                                                 ttGGTGAACAGCATTACCATACGGTCGATACCAATTccc                                    <  1:314631/39‑1 (MQ=255)
                                                                 ttGGTGAACAGCATTACCATACGGTCGATACCAATTccc                                    <  1:3230337/39‑1 (MQ=255)
                                                                 ttGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCa     <  1:3740287/70‑1 (MQ=255)
                                                                  tGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACctgctg                          <  1:3630332/48‑1 (MQ=255)
                                                                   ggTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCa     >  1:429589/1‑68 (MQ=255)
                                                                   ggTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGt   <  1:610347/70‑1 (MQ=255)
                                                                   ggTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTg  >  1:1557379/1‑71 (MQ=255)
                                                                   |                                                                      
AACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTG  >  minE/1932915‑1933052

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: