Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,911,643 A→G 57.4% L184L (TTG→CTG)  thyA ← thymidylate synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,911,6430AG57.4% 25.3 / 54.3 75L184L (TTG→CTG) thyAthymidylate synthetase
Reads supporting (aligned to +/- strand):  ref base A (7/25);  new base G (26/17);  total (33/42)
Fisher's exact test for biased strand distribution p-value = 1.04e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.30e-01

TCGCCACCGGTCCAGACAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAATAACGCGTAGCTGGCAATGTTGAACGGCAGGCCGAGGAAGACGTCACAGGAGCGCTG  >  minE/1911579‑1911700
                                                                |                                                         
tCGCCACCGGTCCAGACAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAATAAc                                                       >  1:2263060/1‑69 (MQ=255)
           ccAGACAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAATAACGCGTAGCTGGc                                            <  1:1055930/69‑1 (MQ=255)
           ccAGACAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAATAACGCGTAGCTGGc                                            <  1:991476/69‑1 (MQ=255)
           ccAGACAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAATAACGCGTAGCTGGc                                            <  1:887183/69‑1 (MQ=255)
           ccAGACAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAATAACGCGTAGCTGGc                                            <  1:855503/69‑1 (MQ=255)
           ccAGACAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAATAACGCGTAGCTGGc                                            <  1:581154/69‑1 (MQ=255)
           ccAGACAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAATAACGCGTAGCTGGc                                            <  1:554783/69‑1 (MQ=255)
           ccAGACAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAATAACGCGTAGCTGGc                                            <  1:389823/69‑1 (MQ=255)
           ccAGACAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAATAACGCGTAGCTGGc                                            <  1:3403750/69‑1 (MQ=255)
           ccAGACAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAATAACGCGTAGCTGGc                                            <  1:329587/69‑1 (MQ=255)
           ccAGACAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAATAACGCGTAGCTGGc                                            <  1:3186362/69‑1 (MQ=255)
           ccAGACAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAATAACGCGTAGCTGGc                                            <  1:3059240/69‑1 (MQ=255)
           ccAGACAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAATAACGCGTAGCTGGc                                            <  1:2996489/69‑1 (MQ=255)
           ccAGACAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAATAACGCGTAGCTGGc                                            <  1:2439183/69‑1 (MQ=255)
           ccAGACAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAATAACGCGTAGCTGGc                                            <  1:1933146/69‑1 (MQ=255)
           ccAGACAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAATAACGCGTAGCTGGc                                            <  1:1186381/69‑1 (MQ=255)
           ccAGACAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAATAACGCGTAGCTGGc                                            <  1:1856346/69‑1 (MQ=255)
           ccAGACAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAATAACGCGTAGCTGGc                                            <  1:1089503/69‑1 (MQ=255)
           ccAGACAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAATAACGCGTAGCTGGc                                            <  1:124583/69‑1 (MQ=255)
           ccAGACAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAATAACGCGTAGCTGGc                                            <  1:1744828/69‑1 (MQ=255)
                cAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAGTAAcgc                                                     >  1:53653/1‑55 (MQ=255)
                cAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAGTAACGCGTAGCt                                               >  1:145855/1‑61 (MQ=255)
                cAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAGTAACGCGTAGCTGGCAAt                                         >  1:3038999/1‑67 (MQ=255)
                cAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAGTAACGCGTAGCTGGCAAt                                         >  1:2654296/1‑67 (MQ=255)
                cAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAGTAACGCGTAGCTGGCAAt                                         >  1:515943/1‑67 (MQ=255)
                cAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAGTAACGCGTAGCTGGCAAt                                         >  1:2885685/1‑67 (MQ=255)
                cAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAGTAACGCGTAGCTGGCAAt                                         >  1:438246/1‑67 (MQ=255)
                cAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAGTAACGCGTAGCTGGCAAt                                         >  1:1742349/1‑67 (MQ=255)
                cAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAGTAACGCGTAGCTGGCAAt                                         >  1:3003066/1‑67 (MQ=255)
                cAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAGTAACGCGTAGCTGGCAAt                                         >  1:3344377/1‑67 (MQ=255)
                cAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAGTAACGCGTAGCTGGCAAt                                         >  1:3399407/1‑67 (MQ=255)
                cAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAGTAACGCGTAGCTGGCAAt                                         >  1:3156390/1‑67 (MQ=255)
                cAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAGTAACGCGTAGCTGGCAAt                                         >  1:3170550/1‑67 (MQ=255)
                cAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAGTAACGCGTAGCTGGCAAt                                         >  1:3288111/1‑67 (MQ=255)
                cAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAGTAACGCGTAGCTGGCAAt                                         >  1:57879/1‑67 (MQ=255)
                cAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAGTAACGCGTAGCTGGCAAt                                         >  1:1022199/1‑67 (MQ=255)
                cAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAGTAACGCGTAGCTGGCAAt                                         >  1:1175127/1‑67 (MQ=255)
                cAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAGTAACGCGTAGCTGGCAAt                                         >  1:2490957/1‑67 (MQ=255)
                cAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAGTAACGCGTAGCTGGCAAt                                         >  1:7202/1‑67 (MQ=255)
                cAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAGTAACGCGTAGCTGGCAAt                                         >  1:2313026/1‑67 (MQ=255)
                cAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAGTAACGCGTAGCTGGCAAt                                         >  1:2278318/1‑67 (MQ=255)
                cAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAGTAACGCGTAGCTGGCAAt                                         >  1:1075516/1‑67 (MQ=255)
                cAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAGTAACGCGTAGCTGGCAAt                                         >  1:22576/1‑67 (MQ=255)
                cAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAGTAACGCGTAGCTGGCAAt                                         >  1:2123292/1‑67 (MQ=255)
                cAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAGTAACGCGTAGCTGGCAAt                                         >  1:2081541/1‑67 (MQ=255)
                cAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAGTAACGCGTAGCTGGCAAt                                         >  1:2038104/1‑67 (MQ=255)
                                  gATCGCACTGCTGCGCCATCATATGCACCAATAACGCGTAGCTGGCAATGTTGAACGGCAGGCCGAgg                      <  1:260018/68‑1 (MQ=255)
                                  gATCGCACTGCTGCGCCATCATATGCACCAATAACGCGTAGCTGGCAATGTTGAACGGCAGGCCGAgg                      <  1:1858705/68‑1 (MQ=255)
                                  gATCGCACTGCTGCGCCATCATATGCACCAATAACGCGTAGCTGGCAATGTTGAACGGCAGGCCGAgg                      <  1:3175695/68‑1 (MQ=255)
                                  gATCGCACTGCTGCGCCATCATATGCACCAATAACGCGTAGCTGGCAATGTTGAACGGCAGGCCGAgg                      <  1:1798452/68‑1 (MQ=255)
                                  gATCGCACTGCTGCGCCATCATATGCACCAATAACGCGTAGCTGGCAATGTTGAACGGCAGGCCGAgg                      <  1:2748990/68‑1 (MQ=255)
                                   aTCGCACTGCTGCGCCATCATATGCACCAATAACGCGTAGCTGGCAATGTTGAACGGCAGGCCGAGGaa                    <  1:2328934/69‑1 (MQ=255)
                                         ctgctgCGCCATCATATGCACCAGTAACGCGTAGCTGGCAATGTTGAACGGCAGGCCGAGGAAGACGTc              <  1:522314/69‑1 (MQ=255)
                                         ctgctgCGCCATCATATGCACCAGTAACGCGTAGCTGGCAATGTTGAACGGCAGGCCGAGGAAGACGTc              <  1:3066113/69‑1 (MQ=255)
                                         ctgctgCGCCATCATATGCACCAGTAACGCGTAGCTGGCAATGTTGAACGGCAGGCCGAGGAAGACGTc              <  1:1426607/69‑1 (MQ=255)
                                         ctgctgCGCCATCATATGCACCAGTAACGCGTAGCTGGCAATGTTGAACGGCAGGCCGAGGAAGACGTc              <  1:185740/69‑1 (MQ=255)
                                         ctgctgCGCCATCATATGCACCAGTAACGCGTAGCTGGCAATGTTGAACGGCAGGCCGAGGAAGACGTc              <  1:181693/69‑1 (MQ=255)
                                         ctgctgCGCCATCATATGCACCAGTAACGCGTAGCTGGCAATGTTGAACGGCAGGCCGAGGAAGACGTc              <  1:2136197/69‑1 (MQ=255)
                                         ctgctgCGCCATCATATGCACCAGTAACGCGTAGCTGGCAATGTTGAACGGCAGGCCGAGGAAGACGTc              <  1:249539/69‑1 (MQ=255)
                                         ctgctgCGCCATCATATGCACCAGTAACGCGTAGCTGGCAATGTTGAACGGCAGGCCGAGGAAGACGTc              <  1:2610362/69‑1 (MQ=255)
                                         ctgctgCGCCATCATATGCACCAGTAACGCGTAGCTGGCAATGTTGAACGGCAGGCCGAGGAAGACGTc              <  1:2775021/69‑1 (MQ=255)
                                         ctgctgCGCCATCATATGCACCAGTAACGCGTAGCTGGCAATGTTGAACGGCAGGCCGAGGAAGACGTc              <  1:2824048/69‑1 (MQ=255)
                                         ctgctgCGCCATCATATGCACCAGTAACGCGTAGCTGGCAATGTTGAACGGCAGGCCGAGGAAGACGTc              <  1:2868652/69‑1 (MQ=255)
                                         ctgctgCGCCATCATATGCACCAGTAACGCGTAGCTGGCAATGTTGAACGGCAGGCCGAGGAAGACGTc              <  1:2936344/69‑1 (MQ=255)
                                         ctgctgCGCCATCATATGCACCAGTAACGCGTAGCTGGCAATGTTGAACGGCAGGCCGAGGAAGACGTc              <  1:3101642/69‑1 (MQ=255)
                                         ctgctgCGCCATCATATGCACCAGTAACGCGTAGCTGGCAATGTTGAACGGCAGGCCGAGGAAGACGTc              <  1:3315550/69‑1 (MQ=255)
                                         ctgctgCGCCATCATATGCACCAATAACGCGTATCTGGCAATGTTGAACGGCAGGCCGAGGAAGACGt               >  1:543863/1‑68 (MQ=255)
                                          tgctgcGCCATCATATGCACCAGTAACGCGTAGCTGGCAATGTTGAACGGCAGGCCGAGGAAGACGTc              <  1:526036/68‑1 (MQ=255)
                                          tgctgcGCCATCATATGCACCAGTAACGCGTAGCTGGCAATGTTGAACGGCAGGCCGAGGAAGACGTc              <  1:2454815/68‑1 (MQ=255)
                                                      atatGCACCAGTAACGCGTAGCTGGCAATGTTGAACGGCAGGCCGAGGAAGACGTc              <  1:817953/56‑1 (MQ=255)
                                                      atatGCACCAATAACGCGTAGCTGGCAATGTTGAACGGCAGGCCGAGGAAGACGTCACAGGAGCGCTg  >  1:530141/1‑68 (MQ=255)
                                                      atatGCACCAATAACGCGTAGCTGGCAATGTTGAACGGCAGGCCGAGGAAGACGTCACAGGAGCGCTg  >  1:496938/1‑68 (MQ=255)
                                                      atatGCACCAATAACGCGTAGCTGGCAATGTTGAACGGCAGGCCGAGGAAGACGTCACAGGAGCGCTg  >  1:566015/1‑68 (MQ=255)
                                                      atatGCACCAATAACGCGTAGCTGGCAATGTTGAACGGCAGGCCGAGGAAGACGTCACAGGAGCGCTg  >  1:619990/1‑68 (MQ=255)
                                                      atatGCACCAATAACGCGTAGCTGGCAATGTTGAACGGCAGGCCGAGGAAGACGTCACAGGAGCGCTg  >  1:1409552/1‑68 (MQ=255)
                                                               aaTAACGCGTAGCTGGCAATGTTGAACGGCAGGCCGAGGAAg                   <  1:551801/42‑1 (MQ=255)
                                                               aaTAACGCGTAGCTGGCAATGTTGAACGGCAGGCCGAGGAAg                   <  1:2654113/42‑1 (MQ=255)
                                                               aaTAACGCGTAGCTGGCAATGTTGAACGGCAGGCCGAGGAAg                   <  1:1612755/42‑1 (MQ=255)
                                                                |                                                         
TCGCCACCGGTCCAGACAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAATAACGCGTAGCTGGCAATGTTGAACGGCAGGCCGAGGAAGACGTCACAGGAGCGCTG  >  minE/1911579‑1911700

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: