Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,933,756 C→A 7.5% A229A (GCG→GCT lysS ← lysine tRNA synthetase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,933,7560CA 7.5% 215.3 / 2.6 66A229A (GCG→GCTlysSlysine tRNA synthetase, constitutive
Reads supporting (aligned to +/- strand):  ref base C (28/33);  new base A (3/2);  total (31/35)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.53e-01

AGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAA  >  minE/1933689‑1933826
                                                                   |                                                                      
aGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGcgcg                                                                     >  1:2184303/1‑71 (MQ=255)
 gCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGt                                                                    <  1:2949407/71‑1 (MQ=255)
  ccaccaACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGt                                                                    >  1:3802675/1‑70 (MQ=255)
      caaccaCCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTAtgg                                                               <  1:2314149/71‑1 (MQ=255)
      caaccaCCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTAtgg                                                               >  1:1888811/1‑71 (MQ=255)
       aaccaCCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAAGTACATGTCGAGATCCAGCGCGTTAtggt                                                              >  1:2978053/1‑71 (MQ=255)
        accaccAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTAtg                                                                >  1:135577/1‑68 (MQ=255)
         ccaccaGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTAtggtgg                                                            <  1:1134267/71‑1 (MQ=255)
         ccaccaGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTAtggtgg                                                            <  1:2897533/71‑1 (MQ=255)
         ccaccaGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTAtggtgg                                                            <  1:2061418/71‑1 (MQ=255)
         caaccaGACGTTTCAGATACAGCTCCGGCGCGATACGCAGATACATATCTAAATCCAGAGCATTAtgatgg                                                            >  1:1558106/3‑71 (MQ=255)
            ccaGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATggtgggt                                                          >  1:3239038/1‑70 (MQ=255)
            ccaGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATggtgggt                                                          >  1:3460855/1‑70 (MQ=255)
             caGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGa                                                        >  1:446239/1‑71 (MQ=255)
              aGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCAAGCGCGTTATGGTGGGTGa                                                        >  1:3154480/1‑70 (MQ=255)
               gACGTTTCAGATACAGCTCCGGCGCGATACGCAGATACATATCTAAATCCAGAGCATTATGATGGGTAAt                                                       >  1:649542/1‑70 (MQ=255)
               gACGTTTCAGATACAGCTCCGGCGCGATACGCAGATACATATCTAAATCCAGAGCATTATGATGGGTAAt                                                       >  1:97730/1‑70 (MQ=255)
                aCGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATa                                                      >  1:2535162/1‑70 (MQ=255)
                  gCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTAt                                                                 <  1:3047046/57‑1 (MQ=255)
                  gCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTAt                                                                 <  1:3002532/57‑1 (MQ=255)
                  gCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTAt                                                                 <  1:654660/57‑1 (MQ=255)
                  gCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTAt                                                                 <  1:3645636/57‑1 (MQ=255)
                      gAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATaaa                                                    >  1:2801156/1‑66 (MQ=255)
                      gAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACgg                                                 <  1:1638514/69‑1 (MQ=255)
                      gAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACgg                                                 <  1:1513571/69‑1 (MQ=255)
                        ggTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGAcg                                              >  1:782838/1‑70 (MQ=255)
                        gatACAGCTCCGGCGCGATACGCAGATACATATCTAAATCCAGAGCATTATGATGGGTAATAAACGGGcga                                             <  1:3758614/69‑2 (MQ=255)
                         gTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGt                                                                    <  1:3457676/47‑1 (MQ=255)
                         gTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGt                                                                    <  1:1438502/47‑1 (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:3654387/71‑1 (MQ=255)
                              gTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCg                                       >  1:1475851/1‑71 (MQ=255)
                               ttCCGGCGCGATACGCAGGTACATGTCGAGATCCAgcg                                                                       <  1:596014/38‑1 (MQ=255)
                               ttCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACgg                                                 <  1:1904697/60‑1 (MQ=255)
                                tCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCgg                                      >  1:3460421/1‑70 (MQ=255)
                                 ccGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATggtggg                                                           <  1:2377535/48‑1 (MQ=255)
                                   ggCGCGATACGCATGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATaa                                                     <  1:2995049/52‑1 (MQ=255)
                                   ggCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATaa                                                     <  1:3497882/52‑1 (MQ=255)
                                   ggCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACgg                                                 <  1:1727812/56‑1 (MQ=255)
                                    gcgcGATACGCAGGTACATGTCGAGATCCAGCGCGTTAtggt                                                              <  1:610453/42‑1 (MQ=255)
                                     cgcgATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATa                                                      <  1:982465/49‑1 (MQ=255)
                                      gcgATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGcc                               <  1:2264892/71‑1 (MQ=255)
                                      gcgATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGcc                               <  1:2126745/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:1026421/1‑71 (MQ=255)
                                                 gTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCt                    >  1:2798574/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)
                                                                   |                                                                      
AGCCACCAACCACCAGACGCTTGAGGTACAGTTCCGGCGCGATACGCAGGTACATGTCGAGATCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAA  >  minE/1933689‑1933826

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: