Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,933,771 G→A 9.8% 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,7710GA 9.8% 143.8 / 5.6 51I224I (ATC→ATTlysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base G (24/22);  new base A (3/2);  total (27/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

CAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGC  >  minE/1933702‑1933838
                                                                     |                                                                   
caGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGa                                                                    >  1:446239/1‑71 (MQ=255)
 aGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCAAGCGCGTTATGGTGGGTGa                                                                    >  1:3154480/1‑70 (MQ=255)
  gACGTTTCAGATACAGCTCCGGCGCGATACGCAGATACATATCTAAATCCAGAGCATTATGATGGGTAAt                                                                   >  1:97730/1‑70 (MQ=255)
  gACGTTTCAGATACAGCTCCGGCGCGATACGCAGATACATATCTAAATCCAGAGCATTATGATGGGTAAt                                                                   >  1:649542/1‑70 (MQ=255)
   aCGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATa                                                                  >  1:2535162/1‑70 (MQ=255)
         gAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATaaa                                                                >  1:2801156/1‑66 (MQ=255)
         gAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACgg                                                             <  1:1513571/69‑1 (MQ=255)
         gAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACgg                                                             <  1:1638514/69‑1 (MQ=255)
           ggTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGAcg                                                          >  1:782838/1‑70 (MQ=255)
           gatACAGCTCCGGCGCGATACGCAGATACATATCTAAATCCAGAGCATTATGATGGGTAATAAACGGGcga                                                         <  1:3758614/69‑2 (MQ=255)
              aCAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGAcgcg                                                        >  1:1354442/1‑69 (MQ=255)
               cAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAg                                                     >  1:2361368/1‑71 (MQ=255)
               cAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAg                                                     >  1:1155071/1‑71 (MQ=255)
                aGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGc                                                    <  1:1819896/71‑1 (MQ=255)
                 gTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCg                                                   >  1:1475851/1‑71 (MQ=255)
                 gTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCg                                                   <  1:3654387/71‑1 (MQ=255)
                  ttCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACgg                                                             <  1:1904697/60‑1 (MQ=255)
                   tCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCgg                                                  >  1:3460421/1‑70 (MQ=255)
                      ggCGCGATACGCATGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATaa                                                                 <  1:2995049/52‑1 (MQ=255)
                      ggCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATaa                                                                 <  1:3497882/52‑1 (MQ=255)
                      ggCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACgg                                                             <  1:1727812/56‑1 (MQ=255)
                        cgcgATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATa                                                                  <  1:982465/49‑1 (MQ=255)
                         gcgATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGcc                                           <  1:2126745/71‑1 (MQ=255)
                         gcgATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGcc                                           <  1:2264892/71‑1 (MQ=255)
                         gcgATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGc                                            >  1:3251873/1‑70 (MQ=255)
                           gATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAg                                         <  1:1788250/71‑1 (MQ=255)
                           gATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAg                                         <  1:3706680/71‑1 (MQ=255)
                             tACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAgg                                        >  1:1957852/1‑70 (MQ=255)
                              aCGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGa                                      >  1:453053/1‑71 (MQ=255)
                               cGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGa                                      <  1:880894/70‑1 (MQ=255)
                                gCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGa                                      <  1:3910733/69‑1 (MQ=255)
                                  aGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCAc                                  >  1:516773/1‑71 (MQ=255)
                                    gTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAgg                                        <  1:2381602/63‑1 (MQ=255)
                                    gTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCt                                >  1:2798574/1‑71 (MQ=255)
                                    gTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCt                                >  1:1026421/1‑71 (MQ=255)
                                        aTGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGAc                                                           >  1:1075801/1‑40 (MQ=255)
                                        aTGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATc                                    <  1:3224735/63‑1 (MQ=255)
                                          gTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCa                                                >  1:1257388/1‑49 (MQ=255)
                                            cGAGCTCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCa                                                >  1:431986/1‑47 (MQ=255)
                                            cGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCAcc                                              >  1:1575550/1‑49 (MQ=255)
                                            cGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCatcatc                        >  1:2247221/1‑71 (MQ=255)
                                             gagaTCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCg                       <  1:1550066/71‑1 (MQ=255)
                                                 tCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCg                    <  1:214493/70‑1 (MQ=255)
                                                   cAGAGCATTATGATGGGTAATAAACGGGCGAGCAGATGCCCCACCTGGAATTACCTGCATCATCGGGGttt                 <  1:3709489/71‑1 (MQ=255)
                                                    aGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGttt                 >  1:1721954/1‑70 (MQ=255)
                                                    aGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTc                <  1:3296067/71‑1 (MQ=255)
                                                     gcgcGTTATGGTGGGTGATAAACGGACGCGCAGCg                                                   <  1:3218267/35‑1 (MQ=255)
                                                      cgcgTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCaa              >  1:2246622/1‑71 (MQ=255)
                                                         attATGATGGGTAATAAACGGGCGAGCAGATGCCCCACCTGGAATTACCTGCATCATCGGGGTTTCTACtt           >  1:3654041/2‑71 (MQ=255)
                                                           tATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTc          >  1:3135061/1‑70 (MQ=255)
                                                             tggtggGCGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCAt       <  1:3850904/71‑1 (MQ=255)
                                                             tgatgggtaATAAACGGGCGAGCAGATGCCCCACCTGGAATTACCTGCATCATCGGGGTTTCTACTTcc         >  1:2541611/4‑69 (MQ=255)
                                                               gtgggtgATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATa      >  1:2177956/1‑70 (MQ=255)
                                                               atgggtaATAAACGGGCGAGCAGATGCCCCACCTGGAATTACCTGCATCATCGGGGTTTCTACTTCCATa      <  1:3696467/69‑1 (MQ=255)
                                                               atgggtaATAAACGGGCGAGCAGATGCCCCACCTGGAATTACCTGCATCATCGGGGTTTCTACTTCCATa      <  1:2810000/69‑1 (MQ=255)
                                                                tgggtgATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATa      >  1:1411182/1‑69 (MQ=255)
                                                                 gggtgATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAg   >  1:1414807/1‑71 (MQ=255)
                                                                  ggtgATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAg   <  1:2314223/70‑1 (MQ=255)
                                                                  ggtgATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGc  >  1:3260724/1‑71 (MQ=255)
                                                                   gtgATAAACGGACGCGCAGCGGCTCCGCCAGGGATCACCTGCATCATCGGCGTTTCAACt            >  1:1454886/1‑60 (MQ=255)
                                                                    tgATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGcatcat                         <  1:1989731/46‑1 (MQ=255)
                                                                     |                                                                   
CAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGC  >  minE/1933702‑1933838

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: