Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 849,820 G→A 100% intergenic (‑1/‑223) yciI ← / → tonB hypothetical protein/membrane spanning protein in TonB‑ExbB‑ExbD complex

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE849,8200GA80.8% 41.3 / 3.8 31intergenic (‑1/‑223)yciI/tonBhypothetical protein/membrane spanning protein in TonB‑ExbB‑ExbD complex
Reads supporting (aligned to +/- strand):  ref base G (4/2);  new base A (0/25);  total (4/27)
Fisher's exact test for biased strand distribution p-value = 4.77e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAA  >  minE/849770‑849882
                                                  |                                                              
cGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTc                                               >  1:615662/1‑68 (MQ=255)
                      tGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTa                         <  1:1309046/68‑1 (MQ=255)
                                 tAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCg             <  1:875364/69‑1 (MQ=255)
                                   aaTAACATACAACACAGGCACAACTCCTTGTTCGTGAAAGTTGTa                                   <  1:1918149/45‑1 (MQ=255)
                                   aaTAACATACAACACAGGCACAACTCCTTGTTCGGGAAAGTTGTa                                   <  1:2427077/45‑1 (MQ=255)
                                   aaTAACATACAACACAGGCACAACTCCTTGTTCGGGAAAGTTGTa                                   <  1:807962/45‑1 (MQ=255)
                                   aaTAACATACAACACAGGCACAACTCCTTGTTCGGGAAAGTTGTa                                   <  1:487225/45‑1 (MQ=255)
                                   aaTAACATACAACACAGGCACAACTCCTTGTTCGGGAAAGTTGTa                                   <  1:456123/45‑1 (MQ=255)
                                   aaTAACATACAACACAGGCACAACTCCTTGTTCGGGAAAGTTGTa                                   <  1:3085071/45‑1 (MQ=255)
                                   aaTAACATACAACACAGGCACAACTCCTTGTTCGGGAAAGTTGTa                                   <  1:3070775/45‑1 (MQ=255)
                                   aaTAACATACAACACAGGCACAACTCCTTGTTCGGGAAAGTTGTa                                   <  1:301997/45‑1 (MQ=255)
                                   aaTAACATACAACACAGGCACAACTCCTTGTTCGGGAAAGTTGTa                                   <  1:2965529/45‑1 (MQ=255)
                                   aaTAACATACAACACAGGCACAACTCCTTGTTCGGGAAAGTTGTa                                   <  1:2869801/45‑1 (MQ=255)
                                   aaTAACATACAACACAGGCACAACTCCTTGTTCGGGAAAGTTGTa                                   <  1:2868371/45‑1 (MQ=255)
                                   aaTAACATACAACACAGGCACAACTCCTTGTTCGGGAAAGTTGTa                                   <  1:2771694/45‑1 (MQ=255)
                                   aaTAACATACAACACAGGCACAACTCCTTGTTCGGGAAAGTTGTa                                   <  1:2403627/45‑1 (MQ=255)
                                   aaTAACATACAACACAGGCACAACTCCTTGTTCGGGAAAGTTGTa                                   <  1:2369221/45‑1 (MQ=255)
                                   aaTAACATACAACACAGGCACAACTCCTTGTTCGGGAAAGTTGTa                                   <  1:2290374/45‑1 (MQ=255)
                                   aaTAACATACAACACAGGCACAACTCCTTGTTCGGGAAAGTTGTa                                   <  1:2246480/45‑1 (MQ=255)
                                   aaTAACATACAACACAGGCACAACTCCTTGTTCGGGAAAGTTGTa                                   <  1:211867/45‑1 (MQ=255)
                                   aaTAACATACAACACAGGCACAACTCCTTGTTCGGGAAAGTTGTa                                   <  1:2076010/45‑1 (MQ=255)
                                   aaTAACATACAACACAGGCACAACTCCTTGTTCGGGAAAGTTGTa                                   <  1:2019925/45‑1 (MQ=255)
                                   aaTAACATACAACACAGGCACAACTCCTTGTTCGGGAAAGTTGTa                                   <  1:1863019/45‑1 (MQ=255)
                                   aaTAACATACAACACAGGCACAACTCCTTGTTCGGGAAAGTTGTa                                   <  1:1466385/45‑1 (MQ=255)
                                   aaTAACATACAACACAGGCACAACTCCTTGTTCGGGAAAGTTGTa                                   <  1:1381144/45‑1 (MQ=255)
                                   aaTAACATACAACACAGGCACAACTCCTTGTTCGGGAAAGTTGTa                                   <  1:1345620/45‑1 (MQ=255)
                                    aTAACATACAACACAGGCACAACTCCTTGTTCGGGAAAGTTGTa                                   <  1:2216484/44‑1 (MQ=255)
                                     tAACATTCAACACAGGCACAACTCCTTGTTCGGGATAGTTGTa                                   <  1:3241899/43‑1 (MQ=255)
                                      aaCATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGAt        >  1:1333860/1‑69 (MQ=255)
                                            caacaCGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCaa  >  1:3239806/1‑69 (MQ=255)
                                            caacaCGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCaa  >  1:776334/1‑69 (MQ=255)
                                                  |                                                              
CGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAA  >  minE/849770‑849882

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: