Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,327,112 C→T 40.6% intergenic (‑102/+174) wcaK ← / ← wzxC predicted pyruvyl transferase/colanic acid exporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,327,1120CT40.6% 4.0 / 18.9 32intergenic (‑102/+174)wcaK/wzxCpredicted pyruvyl transferase/colanic acid exporter
Reads supporting (aligned to +/- strand):  ref base C (13/6);  new base T (5/8);  total (18/14)
Fisher's exact test for biased strand distribution p-value = 1.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGA  >  minE/1327049‑1327173
                                                               |                                                              
tgttttgttCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAg                                                           >  1:2646659/1‑69 (MQ=255)
tgttttgttCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAg                                                           >  1:2461383/1‑69 (MQ=255)
tgttttgttCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAg                                                           >  1:1884126/1‑69 (MQ=255)
                                     tGACGCGTCTTATCAGGCCTACAAGATCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCa                     >  1:3024215/1‑69 (MQ=255)
                                     tGACGCGTCTTATCAGGCCTACAAGATCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCa                     >  1:3194158/1‑69 (MQ=255)
                                     tGACGCGTCTTATCAGGCCTACAAGATCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCa                     >  1:3320817/1‑69 (MQ=255)
                                     tGACGCGTCTTATCAGGCCTACAAGATCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCa                     >  1:1701513/1‑69 (MQ=255)
                                     tGACGCGTCTTATCAGGCCTACAAGATCCGAGCACAGAACCGTAGGACGCATAAGGCGTTCACGCCGCa                     >  1:732058/1‑69 (MQ=17)
                                            tCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCa              <  1:1396881/69‑1 (MQ=255)
                                                  cAGGCCTACAAGACCCGAGCACAGAACCGTAGTACGGATAAGGCGTTCACGCCGCATCCgg                <  1:3380825/61‑1 (MQ=255)
                                                  cAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCgg                <  1:1986961/61‑1 (MQ=255)
                                                  cAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCgg                <  1:2043437/61‑1 (MQ=255)
                                                  cAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTc           <  1:739906/66‑1 (MQ=255)
                                                   aGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCgg                <  1:2483893/60‑1 (MQ=255)
                                                       cTACAAGACCCGAGCACAGAGCCGTAGGACGGATAAGGCGTTTTACGCCGCATCCGGCATTCAGTGCCt    >  1:646099/1‑69 (MQ=14)
                                                       cTACAAGACCCGAGCACAGAGCCGTAGGACGGATAAGGCGTTTTACGCCGCATCCGGCATTCAGTGCCt    >  1:62155/1‑69 (MQ=14)
                                                       cTACAAGACCCGAGCACAGAGCCGTAGGACGGATAAGGCGTTTTACGCCGCATCCGGCATTCAGTGCCt    >  1:87635/1‑69 (MQ=14)
                                                       cTACAAGACCCGAGCACAGAGCCGTAGGACGGATAAGGCGTTTTACGCCGCATCCGGCATTCAGTGCCt    >  1:920557/1‑69 (MQ=14)
                                                       cTACAAGACCCGAGCACAGAGCCGTAGGACGGATAAGGCGTTTTACGCCGCATCCGGCATTCAGTGCCt    >  1:2562143/1‑69 (MQ=14)
                                                       cTACAAGACCCGAGCACAGAGCCGTAGGACGGATAAGGCGTTTTACGCCGCATCCGGCATTCAGTGCCt    >  1:2889840/1‑69 (MQ=14)
                                                       cTACAAGACCCGAGCACAGAGCCGTAGGACGGATAAGGCGTTTTACGCCGCATCCGGCATTCAGTGCCt    >  1:2768978/1‑69 (MQ=14)
                                                       cTACAAGACCCGAGCACAGAGCCGTAGGACGGATAAGGCGTTTTACGCCGCATCCGGCATTCAGTGCCt    >  1:1074142/1‑69 (MQ=14)
                                                       cTACAAGACCCGAGCACAGAGCCGTAGGACGGATAAGGCGTTTTACGCCGCATCCGGCATTCAGTGCCt    >  1:2527897/1‑69 (MQ=14)
                                                       cTACAAGACCCGAGCACAGAGCCGTAGGACGGATAAGGCGTTTTACGCCGCATCCGGCATTCAGTGCCt    >  1:149458/1‑69 (MQ=14)
                                                        tACAAGATCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGa  <  1:3135804/69‑1 (MQ=255)
                                                        tACAAGATCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGa  <  1:2561272/69‑1 (MQ=255)
                                                        tACAAGATCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGa  <  1:60414/69‑1 (MQ=255)
                                                        tACAAGATCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGa  <  1:1877043/69‑1 (MQ=255)
                                                        tACAAGATCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGa  <  1:1821625/69‑1 (MQ=255)
                                                        tACAAGATCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGa  <  1:1767146/69‑1 (MQ=255)
                                                        tACAAGATCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGa  <  1:164430/69‑1 (MQ=255)
                                                        tACAAGATCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGa  <  1:1151189/69‑1 (MQ=255)
                                                               |                                                              
TGTTTTGTTCGTTGCCGCTTTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACAAGACCCGAGCACAGAACCGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTCAGTGCCTGA  >  minE/1327049‑1327173

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: