Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,708,236 A→G 8.1% E260G (GAA→GGA)  trpS → tryptophanyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,708,2360AG 8.1% 210.4 / 8.6 74E260G (GAA→GGA) trpStryptophanyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base A (34/34);  new base G (4/2);  total (38/36)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

TGCAGAACAAAGCGGGCGTTTCCAACCTGTTGGATATCCTTTCAGCGGTAACGGGCCAGAGCATCCCAGAACTGGAAAAACAGTTCGAAGGCAAGATGTATGGTCATCTGAAAGGTGAAGTGGCTGATGCCGTTTCCGG  >  minE/2708167‑2708305
                                                                     |                                                                     
tgCAGAACAAAGCGGGCGTTTCCAACCTGTTGGATATCCTTTCAGCGGTAACGGGCCAGAGCATCCCAGaa                                                                      <  1:1213054/71‑1 (MQ=255)
tgCAGAACAAAGCGGGCGTTTCCAACCTGTTGGATATCCTTTCAGCGGTAACGGGCCAGAGCATCCCAGa                                                                       <  1:130615/70‑1 (MQ=255)
  cAGAACAAAGCGGGCGTTTCCAACCTGTTGGATATCCTTTCAGCGGTAACGGGCCAGAGCATCCCAGAAc                                                                     <  1:1034607/70‑1 (MQ=255)
   aGAACAAAGCGGGCGTTTCCAACCTGTTGGATATCCTTTCAGCGGTAACGGGCCAGAGCATCCCAGAACTg                                                                   >  1:983569/1‑71 (MQ=255)
   aGAACAAAGCGGGCGTTTCCAACCTGTTGGATATCCTTTCAGCGGTAACGGGCCAGAGCATCCCAGAACTg                                                                   >  1:3145158/1‑71 (MQ=255)
   aGAACAAAGCGGGCGTTTCCAACCTGTTGGATATCCTTTCAGCGGTAACGGGCCAGAGCATCCCAGAACTg                                                                   <  1:1901613/71‑1 (MQ=255)
        aaaGCGGGCGTTTCCAACCTGTTGGATATCCTTTCAGCGGTAACGGGCCAGAGCATCCCAGAACTGGaaa                                                               >  1:3290736/1‑70 (MQ=255)
        aaaGCGGGCGTTTCCAACCTGTTGGATATCCTTTCAGCGGTAACGGGCCAGAGCATCCCAGAACTGGaaa                                                               >  1:2149361/1‑70 (MQ=255)
             gggCGTTTCCAACCTGTTGGATATCCTTTCAGCGGTAACGGGCCAGAGCATCCCAGAACTGGAAAAACAGt                                                         >  1:1665954/1‑71 (MQ=255)
                 gcttCCAACCTGTTGGATATCCTTTCAGCGGTAACGGGCCAGAGCATCCCAGAACTGGAAAAACAGTTCGa                                                     <  1:2795609/69‑1 (MQ=255)
                 gTTTCCAACCTGTTGGATATCCTTTCAGCGGTAACGGGCCAGAGCATCCCAGAACTGGAAAAACAGTTCGa                                                     <  1:2764709/71‑1 (MQ=255)
                 gTTTCCAACCTGTTGGATATCCTTTCAGCGGTAACGGGCCAGAGCATCCCAGAACTGGAAAAACAGTTCGa                                                     <  1:1494866/71‑1 (MQ=255)
                  tttCCAACCTGTTGGATATCCTTTCAGCGGTAACGGGCCAGAGCATCCCAGAACTGGAAAAACAGTTCg                                                      >  1:479958/1‑69 (MQ=255)
                  tttCCAACCTGTTGGATATCCTTTCAGCGGTAACGGGCCAGAGCATCCCAGAACTGGAAAAACAGTTCGa                                                     <  1:509267/70‑1 (MQ=255)
                    tCCAACCTGTTGGATATCCTTTCAGCGGTAACGGGCCAGAGCATCCCAGAACTGGAAAAACAGTTCGAAgg                                                  >  1:3459032/1‑71 (MQ=255)
                    tCCAACCTGTTGGATATCCTTTCAGCGGTAACGGGCCAGAGCATCCCAGAACTGGAAAAACAGTTCGAAgg                                                  >  1:1357172/1‑71 (MQ=255)
                     ccAACCTGTTGGATATCCTTTCAGCGGTAACGGGCCAGAGCATCCCAGAACTGGAAAAACAGTTCGAAgg                                                  >  1:2518154/1‑70 (MQ=255)
                     ccAACCTGTTGGATATCCTTTCAGCGGTAACGGGCCAGAGCATCCCAGAACTGGAAAAACAGTTCGAAgg                                                  >  1:693640/1‑70 (MQ=255)
                     ccAACCTGTTGGATATCCTTTCAGCGGTAACGGGCCAGAGCATCCCAGAACTGGAAAAACAGTTCGAAgg                                                  >  1:1825316/1‑70 (MQ=255)
                     ccAACCTGTTGGATATCCTTTCAGCGGTAACGGGCCAGAGCATCCCAGAACTGGAAAAACAGTTCGAAgg                                                  >  1:1708233/1‑70 (MQ=255)
                      cAACCTGTTGGATATCCTTTCAGCGGTAACGGGCCAGAGCATCCCAGAACTGGAAAAACAGTTCg                                                      <  1:2044377/65‑1 (MQ=255)
                      cAACCTGTTGGATATCCTTTCAGCGGTAACGGGCCAGAGCATCCCAGAACTGGAAAAACAGTTCg                                                      <  1:584052/65‑1 (MQ=255)
                      cAACCTGTTGGATATCCTTTCAGCGGTAACGGGCCAGAGCATCCCAGAACTGGAAAAACAGTTCGa                                                     <  1:1709842/66‑1 (MQ=255)
                      cAACCTGTTGGATATCCTTTCAGCGGTAACGGGCCAGAGCATCCCAGAACTGGAAAAAAAGTTCg                                                      <  1:475840/65‑1 (MQ=255)
                       aaCCTGTTGGATATCCTTTCAGCGGTAACGGGCCAGAGCATCCCAGAACTGGAAAAACAGTTCGAAGGc                                                 <  1:2953155/69‑1 (MQ=255)
                         ccTGTTGGATATCCTTTCAGCGGTAACGGGCCAGAGCATCCCAGAACTGGAAAAACAGTTCg                                                      >  1:1788080/1‑62 (MQ=255)
                         ccTGTTGGATATCCTTTCAGCGGTAACGGGCCAGAGCATCCCAGAACTGGAAAAACAGTTCg                                                      >  1:838232/1‑62 (MQ=255)
                         ccTGTTGGATATCCTTTCAGCGGTAACGGGCCAGAGCATCCCAGAACTGGAAAAACAGTTCg                                                      >  1:56788/1‑62 (MQ=255)
                          cTGTTGGATATCCTTTCAGCGGTAACGGGCCAGAGCATCCCAGAACTGGAAAAACAGTTCGAAGGCAAGAt                                            <  1:1823834/71‑1 (MQ=255)
                             ttGGATATCCTTTCAGCGGTAACGGGCCAGAGCATCCCAGAACTGGAAAAACAGTTCGAAGGCAAGATGTa                                         <  1:495728/71‑1 (MQ=255)
                               ggATATCCTTTCAGCGGTAACTGGCCAGAGCATCCCAGAACTGGAAAAACAGTTCGAAGGCAAGATg                                           <  1:942717/67‑1 (MQ=255)
                               ggATATCCTTTCAGCGGTAACGGGCCAGAGCATCCCAGGACTGGAAAAACAGTTCGAAGGCAAGATGTATg                                       >  1:500508/1‑71 (MQ=255)
                                 atatCCTTTCAGCGGTAACGGGCCAGAGCATCCCAGAACTGGAAAAACAGTTCGAAGGCAAGATGTATGGt                                     >  1:283013/1‑71 (MQ=255)
                                  tatCCTTTCAGCGGTAACGGGCCAGAGCATCCCAGAACTGGAAAAACAGTTCGa                                                     <  1:2290305/54‑1 (MQ=255)
                                   atCCTTTCAGCGGTAACGGGCCAGAGCATCCCAGaa                                                                      <  1:1996312/36‑1 (MQ=255)
                                   atCCTTTCAGCGGTAACGGGCCAGAGCATCCCAGaa                                                                      <  1:1252592/36‑1 (MQ=255)
                                   atCCTTTCAGCGGTAACGGGCCAGAGCATCCCAGaa                                                                      <  1:719729/36‑1 (MQ=255)
                                   atCCTTTCAGCGGTAACGGGCCAGAGCATCCCAGaa                                                                      <  1:500451/36‑1 (MQ=255)
                                   atCCTTTCAGCGGTAACGGGCCAGAGCATCCCAGAACTGGAAAAACAGt                                                         <  1:1339735/49‑1 (MQ=255)
                                     ccTTTCAGCGGTAACGGGCCAGAGCATCCCAGAACt                                                                    >  1:3436270/1‑36 (MQ=255)
                                     ccTTTCAGCGGTAACGGGCCAGAGCATCCCAGAACt                                                                    >  1:3152181/1‑36 (MQ=255)
                                     ccTTTCAGCGGTAACGGGCCAGAGCATCCCAGAACt                                                                    >  1:3032863/1‑36 (MQ=255)
                                     ccTTTCAGCGGTAACGGGCCAGAGCATCCCAGAACTGGAAAAACAGTTCGAAGGCAAGATGTa                                         <  1:2518012/63‑1 (MQ=255)
                                       tttCAGCGGTAACGGGCCAGAGCATCCCAGAACTGGaaa                                                               <  1:3232510/39‑1 (MQ=255)
                                       tttCAGCGGTAACGGGCCAGAGCATCCCAGAACTGGaaa                                                               <  1:3460926/39‑1 (MQ=255)
                                         tCAGCGGTAACGGGCCAGAGCATCCCAGAACTGGAAAAACAGTTCGAAGGCAAGa                                             >  1:3434862/1‑55 (MQ=255)
                                         tCAGCGGTAACGGGCCAGAGCATCCCAGAACTGGAAAAACAGTTCGAAGGCAAGa                                             >  1:1904373/1‑55 (MQ=255)
                                         tCAGCGGTAACGGGCCAGAGCATCCCAGAACTGGAAAAACAGTTCGAAGGCAAGATGTATGGTCATCTGa                              >  1:185323/1‑70 (MQ=255)
                                          cAGCGGTAACGGGCCAGAGCATCCCAGGACTGGAAAAACAGTTCGAAGGCAAGATGTATGGTc                                    <  1:91488/63‑1 (MQ=255)
                                            gCGGTAACGGGCCAGAGCATCCCAGGACTGGAAAAACAGTTCGAAGGCAAGATGTATGGTCAt                                  >  1:3251881/1‑63 (MQ=255)
                                              ggTAACGGGCCAGAGCATCCCAGAACTGGAAAAACAGTTCGAAGGCAAGATGTATGGTCATCTGAAAGGTg                        <  1:1994403/71‑1 (MQ=255)
                                                tAACGGGCCAGAGCATCCCAGAACTGGAAAAACAGTTCGa                                                     <  1:1617215/40‑1 (MQ=255)
                                                tAACGGGCCAGAGCATCCCAGAACTGGAAAAACAGTTCGAAGGCAAGATGTATGGTCATCTGa                              <  1:3594495/63‑1 (MQ=255)
                                                tAACGGGCCAGAGCATCCCAGAACTGGAAAAACAGTTCGAAGGCAAGATGTATGGTCATCTGa                              <  1:2446589/63‑1 (MQ=255)
                                                tAACGGGCCAGAGCATCCCAGAACTGGAAAAACAGTTCGAAGGCAAGATGTATGGTCATCTGAAAGGTg                        <  1:3617810/69‑1 (MQ=255)
                                                tAACGGGCCAGAGCATCCCAGAACTGGAAAAACAGTTCGAAGGCAAGATGTATGGTCATCTGAAAGGTGa                       >  1:3017153/1‑70 (MQ=255)
                                                tAACGGGCCAGAGCATCCCAGAACTGGAAAAACAGTTCGAAGGCAAGATGTATGGTCATCTGAAAGGTGa                       >  1:133534/1‑70 (MQ=255)
                                                 aaCGGGCCAGAGCATCCCAGAACTGGAAAAACAGTTCGAAGGCAAGATGTATGGTCATCTGa                              <  1:2161763/62‑1 (MQ=255)
                                                  aCGGGCCAGAGCATCCCAGAACTGGAAAAACAGTTCGAAGGCAAGATGTATGGTCATCTGAAAGGTGAAGt                    >  1:753183/1‑71 (MQ=255)
                                                    gggCCAGAGCATCCCAGAACTGGAAAAACAGTTCGAAGGCAAGATGTATGGTCATCTGAAAGGTGAAGTgg                  >  1:1230172/1‑71 (MQ=255)
                                                     ggCCAGAGCATCCCAGAACTGGAAAAACAGTTCGAAGGCAAGATGTATGGTCATc                                 >  1:832847/1‑55 (MQ=255)
                                                     ggCCAGAGCATCCCAGAACTGGAAAAACAGTTCGAAGGCAAGATGTATGGTCATCTGAAAGGTGAAGTGGc                 <  1:2947250/71‑1 (MQ=255)
                                                      gCCAGAGCATCCCAGAACTGGAAAAACAGTTCGAAGGCAAGATGTATGGTCATCTGAAAGGTGAAGTGGc                 >  1:1264733/1‑70 (MQ=255)
                                                      gCCAGAGCATCCCAGAACTGGAAAAACAGTTCGAAGGCAAGATGTATGGTCATCTGAAAGGTGAAGTGGc                 >  1:593155/1‑70 (MQ=255)
                                                       ccAGAGCATCCCAGAACTGGAAAAACAGTTCGAAGGCAAGATGTa                                         >  1:1893730/1‑45 (MQ=255)
                                                        cAGAGCATCCCAGAACTGGAAAAACAGTTCGAAGGCAAGATGTATGGTCATCt                                <  1:377297/53‑1 (MQ=255)
                                                        cAGAGCATCCCAGAACTGGAAAAACAGTTCGAAGGCAAGATGTATGGTCATCt                                <  1:2353415/53‑1 (MQ=255)
                                                        cAGAGCATCCCAGAACTGGAAAAACAGTTCGAAGGCAAGATGTATGGTCATCTGAAAGGTGAAGTGGCTGa              >  1:3276232/1‑71 (MQ=255)
                                                        cAGAGCATCCCAGAACTGGAAAAACAGTTCGAAGGCAAGATGTATGGTCATCTGAAAGGTGAAGAGGCTGa              <  1:481372/71‑1 (MQ=255)
                                                         agagCATCCCAGGACTGGAAAAACAGTTCGAAGGCAAGATGTATGGTCATCTGAAAGGTGAAGTGGCTGAt             >  1:1920253/1‑71 (MQ=255)
                                                          gagCATCCCAGGACTGGAAAAACAGTTCGAAGGCAAGATGTATGGTCATCTGAAAGGTGAAGTGGCTGATg            >  1:292972/1‑71 (MQ=255)
                                                          gagCATCCCAGAACTGGAAAAACAGTTCGAAGGCAAGATGTATGGTCATCTGAAAGGTGAAGTGGCTGATg            >  1:1184822/1‑71 (MQ=255)
                                                           agCATCCCAGAACTGGAAAAACAGTTCGAAGGCAAGATGTATGGTCATCTGAAAGGTGAAGTGGCTGATGc           >  1:1100736/1‑71 (MQ=255)
                                                                cccAGAACTGGAAAAACAGTTCGAAGGCAAGATGTATGGTCATCTGAAAGGTg                        <  1:1116429/53‑1 (MQ=255)
                                                                 ccAGAACTGGAAAAACAGTTCGAAGGCAAGATGTATGGTCATCTGAAAGGTGAAGTGGCTGATGCCGTTTc     <  1:1907137/71‑1 (MQ=255)
                                                                 ccAGAACTGGAAAAACAGTTCGAAGGCAAGATGTATGGTCATCTGAAAGGTGAAGTGGCTGATGCCGTTTc     <  1:1191610/71‑1 (MQ=255)
                                                                 ccAGAACTGGAAAAACAGTTCGAAGGCAAGATGTATGGTCATCTGAAAGGTGAAGTGGCTGATGCCGTTTc     <  1:84073/71‑1 (MQ=255)
                                                                 ccAGAACTGGAAAAACAGTTCGAAGGCAAGATGTATGGTCATCTGAAAGGTGAAGTGGCTGATGCCGTTTc     <  1:1496687/71‑1 (MQ=255)
                                                                  cAGGACTGGAAAAACAGTTCGAAGGCAAGATGTATGGTCATCTGAAAGGTGAAGTGGCTGATGCCGttt      <  1:3435081/69‑1 (MQ=255)
                                                                  cAGAACTGGAAAAACAGTTCGAAGGCAAGATGTATGGTCATCTGAAAGGTGAAGTGGCTGATGc           >  1:3055167/1‑64 (MQ=255)
                                                                     aaCTGGAAAAACAGTTCGAAGGCAAGATGTATGGTCATCTGAAAGGTGAAGTGGCTGATGCCGTTTc     >  1:2112057/1‑67 (MQ=255)
                                                                     aaCTGGAAAAACAGTTCGAAGGCAAGATGTATGGTCATCTGAAAGGTGAAGTGGCTGATGCCGTTTc     >  1:3737484/1‑67 (MQ=255)
                                                                     aaCTGGAAAAACAGTTCGAAGGCAAGATGTATGGTCATCTGAAAGGTGAAGTGGCTGATGCCGTTTCCgg  >  1:142688/1‑70 (MQ=255)
                                                                     |                                                                     
TGCAGAACAAAGCGGGCGTTTCCAACCTGTTGGATATCCTTTCAGCGGTAACGGGCCAGAGCATCCCAGAACTGGAAAAACAGTTCGAAGGCAAGATGTATGGTCATCTGAAAGGTGAAGTGGCTGATGCCGTTTCCGG  >  minE/2708167‑2708305

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: