Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 106,166 C→T 31.2% Q42Q (CAG→CAA aroP ← aromatic amino acid transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE106,1660CT31.2% 45.2 / 21.2 51Q42Q (CAG→CAAaroParomatic amino acid transporter
Reads supporting (aligned to +/- strand):  ref base C (18/17);  new base T (9/7);  total (27/24)
Fisher's exact test for biased strand distribution p-value = 7.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAA  >  minE/106098‑106233
                                                                    |                                                                   
tGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGAc                                                                     <  1:1274692/69‑1 (MQ=255)
tGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGAc                                                                     <  1:147252/69‑1 (MQ=255)
    gCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTg                                                                   <  1:3049280/67‑1 (MQ=255)
    gCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTg                                                                   <  1:2448169/67‑1 (MQ=255)
    gCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTg                                                                   <  1:355016/67‑1 (MQ=255)
    gCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTg                                                                   <  1:455473/67‑1 (MQ=255)
    gCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTg                                                                   <  1:18215/67‑1 (MQ=255)
                  gctaTAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGATTGTATTACGGAGGCGCTa                                                   >  1:2136181/4‑69 (MQ=255)
                  gCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGATTGTATTACGGAGGCGCTa                                                   >  1:899767/1‑69 (MQ=255)
                  gCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGATTGTATTACGGAGGCGCTa                                                   >  1:3184020/1‑69 (MQ=255)
                  gCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGATTGTATTACGGAGGCGCTa                                                   >  1:2265750/1‑69 (MQ=255)
                  gCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGATTGTATTACGGAGGCGCTa                                                   >  1:532440/1‑69 (MQ=255)
                  gCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGATTGTATTACGGAGGCGCTa                                                   >  1:567347/1‑69 (MQ=255)
                  gCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGATTGTATTACGGAGGCGCTa                                                   >  1:1947104/1‑69 (MQ=255)
                  gCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGATTGTATTACGGAGGCGCTa                                                   >  1:1687740/1‑69 (MQ=255)
                  gCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGATTGTATTACGGAGGCGCTa                                                   >  1:1452527/1‑69 (MQ=255)
                             aGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTACGGAGGCGCTACCCAGGAATaa                                        <  1:1761498/69‑1 (MQ=255)
                             aGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTACGGAGGCGCTACCCAGGAATaa                                        <  1:2122008/69‑1 (MQ=255)
                             aGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTACGGAGGCGCTACCCAGGAATaa                                        <  1:1236755/69‑1 (MQ=255)
                                aaTGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTACGGAGGCGCTACCCAGGAATAAccc                                     >  1:2425662/1‑69 (MQ=255)
                                   ggCGTAACCCAGGATAATCCCTGGCCCTGCGGTCTg                                                                   <  1:853366/36‑1 (MQ=255)
                                   ggCGTAACCCAGGATAATCCCTGGCCCTGCGGACTg                                                                   <  1:3104866/36‑1 (MQ=255)
                                   ggCGTAACCCAGGATAATCCCTGGCCCTGCGGACTg                                                                   <  1:1529843/36‑1 (MQ=255)
                                   ggCGTAACCCAGGATAATCCCTGGCCCTGCGGACTg                                                                   <  1:886368/36‑1 (MQ=255)
                                   ggCGTAACCCAGGATAATCCCTGGCCCTGCGGACTg                                                                   <  1:1205991/36‑1 (MQ=255)
                                   ggCGTAACCCAGGATAATCCCTGGCCCTCCGGACTg                                                                   <  1:2420230/36‑1 (MQ=255)
                                    gCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTACGGAGGCGCTACCCAGGAATAACCCGGt                                  >  1:2687180/1‑68 (MQ=255)
                                    gCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTc                                 >  1:841509/1‑69 (MQ=255)
                                    gCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTc                                 >  1:2506651/1‑69 (MQ=255)
                                    gCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTc                                 >  1:1767134/1‑69 (MQ=255)
                                    gCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTACGGAGGCGCTACCCAGGAATAACCAgg                                   >  1:850214/1‑67 (MQ=255)
                                         aCCCAGGATAATCCCTGGCCCTGCGGACTGTATTACGGa                                                          >  1:944706/1‑39 (MQ=255)
                                         aCCCAGGATAATCCCTGGCCCTGCGGACTGTATTACGGa                                                          >  1:2925383/1‑39 (MQ=255)
                                         aCCCAGGATAATCCCTGGCCCTGCGGACTGTATTACGGa                                                          >  1:2696364/1‑39 (MQ=255)
                                         aCCCAGGATAATCCCTGGCCCTGCGGACTGTATTACGGa                                                          >  1:122908/1‑39 (MQ=255)
                                         aCCCAGGATAATCCCTGGCCCTGCGGACTGTATTACGGa                                                          >  1:2025663/1‑39 (MQ=255)
                                                     cccTGGCCCTGCGGACTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCa                 >  1:1093255/1‑68 (MQ=255)
                                                            ccTGCGGATTGTATTACGGAGGCGCTACCCAGGAATAGCCCTGTCCCTATCGCGCCACCCAGCGCGATa         <  1:2289666/69‑1 (MQ=255)
                                                            ccTGCGGATTGTATTACGGAGGCGCTACCCAGGAATAGCCCTGTCCCTATCGCGCCACCCAGCGCGATa         <  1:1002337/69‑1 (MQ=255)
                                                            ccTGCGGATTGTATTACGGAGGCGCTACCCAGGAATAGCCCTGTCCCTATCGCGCCACCCAGCGCGATa         <  1:1430094/69‑1 (MQ=255)
                                                            ccTGCGGATTGTATTACGGAGGCGCTACCCAGGAATAGCCCTGTCCCTATCGCGCCACCCAGCGCGATa         <  1:1867741/69‑1 (MQ=255)
                                                            ccTGCGGATTGTATTACGGAGGCGCTACCCAGGAATAGCCCTGTCCCTATCGCGCCACCCAGCGCGATa         <  1:3220486/69‑1 (MQ=255)
                                                            ccTGCGGATTGTATTACGGAGGCGCTACCCAGGAATAGCCCTGTCCCTATCGCGCCACCCAGCGCGATa         <  1:2649859/69‑1 (MQ=255)
                                                            ccTGCGGATTGTATTACGGAGGCGCTACCCAGGAATAGCCCTGTCCCTATCGCGCCACCCAGAGCGATa         <  1:2657/69‑1 (MQ=255)
                                                                cGGACTGTATTACGGAGGCGCTACCCAGGAATAAccc                                     >  1:3033097/1‑37 (MQ=255)
                                                                cGGACTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGAt          >  1:2198099/1‑64 (MQ=255)
                                                                cGGACTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGc      >  1:3358502/1‑68 (MQ=255)
                                                                cGGACTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGc      >  1:2225960/1‑68 (MQ=255)
                                                                cGGACTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGc      >  1:1743358/1‑68 (MQ=255)
                                                                cGGACTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGc      >  1:2784717/1‑68 (MQ=255)
                                                                   aCTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGaa  <  1:329739/69‑1 (MQ=255)
                                                                   aCTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGaa  <  1:279906/69‑1 (MQ=255)
                                                                   aCTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGaa  <  1:3215153/69‑1 (MQ=255)
                                                                    |                                                                   
TGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAA  >  minE/106098‑106233

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: