Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,932,985 G→A 11.1% F486F (TTC→TTT lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,932,9850GA11.1% 247.1 / 17.0 90F486F (TTC→TTTlysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base G (45/35);  new base A (4/6);  total (49/41)
Fisher's exact test for biased strand distribution p-value = 5.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.35e-01

ACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCG  >  minE/1932916‑1933054
                                                                     |                                                                     
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:3554281/71‑1 (MQ=255)
                   ggACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTc                                                   <  1:1995375/71‑1 (MQ=255)
                   ggACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTc                                                   <  1:142934/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)
                         aTCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCa                                                          >  1:342351/1‑58 (MQ=255)
                         aTCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATAc                                              >  1:1589153/1‑70 (MQ=255)
                         aTCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATAc                                              >  1:399179/1‑70 (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:2444327/1‑46 (MQ=255)
                                       aaTAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATa                                                        >  1:2073444/1‑46 (MQ=255)
                                       aaTAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACctg                               >  1:3221528/1‑71 (MQ=255)
                                       aaTAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACctg                               <  1:2702998/71‑1 (MQ=255)
                                       aaTAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACctg                               >  1:729160/1‑71 (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:1002513/1‑71 (MQ=255)
                                                      ggTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGc                >  1:143374/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:3230337/39‑1 (MQ=255)
                                                                ttGGTGAACAGCATTACCATACGGTCGATACCAATTccc                                      <  1:314631/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)
                                                                   gTGAACAGCATTACCATACTGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTg    <  1:81702/70‑1 (MQ=255)
                                                                   gTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTg    <  1:1431549/70‑1 (MQ=255)
                                                                    tGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCg  <  1:126706/71‑1 (MQ=255)
                                                                     |                                                                     
ACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAAGCCATGTTCCAGTGCG  >  minE/1932916‑1933054

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: