Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 199,135 T→A 35.4% S207C (AGC→TGC)  metQ ← DL‑methionine transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE199,1350TA35.4% 32.1 / 25.8 31S207C (AGC→TGC) metQDL‑methionine transporter subunit
Reads supporting (aligned to +/- strand):  ref base T (16/4);  new base A (7/4);  total (23/8)
Fisher's exact test for biased strand distribution p-value = 4.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTTACGTACGGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCATAGGTGGTATTGATAACTGCCAGAGCGATTTGCGCGTCGTCCAGAGAACGCGGCAGTTG  >  minE/199067‑199198
                                                                    |                                                               
gTTTACGTACGGTGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATctggctgg                                                               <  1:380518/71‑1 (MQ=255)
gTTTACGTACGGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATctggctgg                                                               <  1:848924/71‑1 (MQ=255)
        acgGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCATAggtg                                                       >  1:2727212/1‑71 (MQ=255)
               ctctTTATCTTCAACAAAGATACCGTCTTTCGCCGTAGTCAGGCCAATCTGGCTGGCATAGGTGGTATTGa                                                >  1:3365771/1‑71 (MQ=255)
               ctctTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCATAGGTGGTATTGa                                                >  1:3751760/1‑71 (MQ=255)
               ctctTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCATAGGTGGTATTGa                                                >  1:2823891/1‑71 (MQ=255)
               ctctTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCAGGCATAGGTGGTATTGa                                                <  1:4029020/71‑1 (MQ=255)
                  tttATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATctggcagg                                                               <  1:2561386/53‑1 (MQ=255)
                  tttATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCATAGGTGGTATTGATaa                                             >  1:3800005/1‑71 (MQ=255)
                   ttATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCa                                                             >  1:4026601/1‑54 (MQ=255)
                    tATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCATAGGTGGTATTGATAACt                                           >  1:3792517/1‑71 (MQ=255)
                    tATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCAGGCATAGGTGGTATTGATAACt                                           >  1:3146616/1‑71 (MQ=255)
                    tATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCAGGCATAGGTGGTATTGATAACt                                           >  1:3263347/1‑71 (MQ=255)
                      tCTTCAACAAAGATACCGTCTTTAGCCGGAGTCAGGCCAATctggctg                                                                >  1:1030356/1‑48 (MQ=255)
                        ttCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCAGGCATAGGTGGTATTGAt                                               <  1:444026/63‑1 (MQ=255)
                        ttCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCAGGCATAGGTGGTATTGAt                                               <  1:505073/63‑1 (MQ=255)
                          caacaaAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCATAGGTGGTATTGATAACTg                                          >  1:2469795/1‑66 (MQ=255)
                          caacaaAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCATAGGTGGTATTGATAACTGc                                         <  1:3833012/67‑1 (MQ=255)
                              aaaGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCATAGGTGGTATTGATAACTGCCAGAGCGa                                 >  1:2832995/1‑71 (MQ=255)
                               aaGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCATAGGTGGTATTGATAACTGCCAGAGCGAt                                >  1:1273133/1‑71 (MQ=255)
                                    aCCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCATAg                                                          >  1:2770637/1‑40 (MQ=255)
                                            tCGCCGGAGTCAGGCCAATCTGGCAGGCATAGGTGGTATTGATAACTGCCAGAGCGAtt                               >  1:2271349/1‑59 (MQ=255)
                                            tCGCCGGAGTCAGGCCAATCTGGCAGGCATAGGTGGTATTGATAACTGCCAGAGCGATTTGCGCgtcgtc                    >  1:1425883/1‑70 (MQ=255)
                                            tCGCCGGAGTCAGGCCAATCTGGCAGGCATAGGTGGTATTGATAACTGCCAGAGCGATTTGCGCgtcgtc                    >  1:3843395/1‑70 (MQ=255)
                                             cGCCGGAGTCAGGCCAATCTGGCAGGCATAGGTGGTATTGATAACTGCCAGAGCGATTTGCGCgtcgtc                    >  1:2807230/1‑69 (MQ=255)
                                                 ggAGTCAGGCCAATCTGGCTGGCATAGGTGGTATTGATAACTGCCAGAGCGATTTGCGCGTCGTCCagag               >  1:3603288/1‑70 (MQ=255)
                                                 ggAGTCAGGCCAATCTGGCTGGCATAGGTGGTATTGATAACTGCCAGAGCGATTTGCGCGTCGTCCagag               >  1:718361/1‑70 (MQ=255)
                                                 ggAGTCAGGCCAATCTGGCAGGCATAGGTGGTATTGATAACTGCCAGAGCGATTTGCGCGTCGTCCagag               >  1:2288593/1‑70 (MQ=255)
                                                    gTCAGGCCAATCTGGCTGGCATAGGTGGTATTGATAACTGCCAGAGCGATTTGCGCGTCGTCCAGAGAAcg           >  1:2388588/1‑71 (MQ=255)
                                                    gTCAGGCCAATCTGGCTGGCATAGGTGGTATTGATAACTGCCAGAGCGATTTGCGCGTCGTCCAGAGAAcg           >  1:2904717/1‑71 (MQ=255)
                                                    gTCAGGCCAATCTGGCTGGCATAGGTGGTATTGATAACTGCCAGAGCGATTTGCGCGTAGTCCAGAGAAcg           >  1:1930932/1‑71 (MQ=255)
                                                    gTCAGGCCAATCTGGCAGGCATAGGTGGTATTGATAACTGCCAGAGCGATTTGCGCGTCGTCCAGAGAAcg           <  1:3567915/71‑1 (MQ=255)
                                                         gCCAATCTGGCTGGCATAGGTGGTATTGATAACTGCCAGAGCGATTTGCGCGTCGTCCAGAGAACGCGGCa      <  1:3271733/71‑1 (MQ=255)
                                                              tCTGGCTGGCATAGGTGGTATTGATAACTGCCAGAGCGATTTGCGCGTCGTCCAGAGAACGCGGCAGTTg  <  1:765048/70‑1 (MQ=255)
                                                              tCTGGCTGGCATAGGTGGTATTGATAACTGCCAGAGCGATTTGCGCGTCGTCCAGAGAACGCGGCAGTTg  <  1:3768004/70‑1 (MQ=255)
                                                               ctggctggCATAGGTGGTATTGATAACTGCCAGAGCGATTTGCGCGTCGTCCaga                >  1:1759224/1‑55 (MQ=255)
                                                               ctggctggCATAGGTGGTATTGATAACTGCCAGAGCGATTTGCGCGTCGTCCaga                >  1:3002659/1‑55 (MQ=255)
                                                               ctggctggCATAGGTGGTATTGATAACTGCCAGAGCGATTTGCGCGTCGTCCaga                >  1:25348/1‑55 (MQ=255)
                                                               ctggctggCAT‑‑GTGGTATTGATAACTGCCAGAGCGATTTGCGCGTCGTCCaga                >  1:1192729/1‑53 (MQ=255)
                                                                    |                                                               
GTTTACGTACGGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCATAGGTGGTATTGATAACTGCCAGAGCGATTTGCGCGTCGTCCAGAGAACGCGGCAGTTG  >  minE/199067‑199198

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: