Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,704,968	C→T	24.1%	Y393Y (TAC→TAT)	damX →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,704,968	0	C	T	24.1%	88.3 / 21.5	70	Y393Y (TAC→TAT)	damX	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (25/28);  new base T (3/14);  total (28/42)
Fisher's exact test for biased strand distribution p-value = 4.55e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.95e-01

TGAAACGACGCGTAATGGTCAGCCGTGGTATGTCCTGGTTTCTGGCGTGTACGCTTCGAAAGAAGAGGCGAAAAAAGCGGTATCTACATTGCCAGCAGATGTCCAGGCCAAAAACCC  >  minE/2704917‑2705033                                                    |                                                                  tGAAACGACGCGTAATGGTCGGCCGTGGTATGTCCTGGTTTCTGGCGTGTACGCTTCGAAAGAAGAgg                                                   >  1:1104352/1‑68 (MQ=255) tGAAACGACGCGTAATGGTCAGCCGTGGTATGTCCTGGTTTCTGGCGTGTACGCTTCGAAAGAAGAgg                                                   >  1:2995122/1‑68 (MQ=255) tGAAACGACGCGTAATGGTCAGCCGTGGTATGTCCTGGTTTCTGGCGTGTACGCTTCGAAAGAAGAgg                                                   >  1:3431056/1‑68 (MQ=255) tGAAACGACGCGTAATGGTCAGCCGTGGTATGTCCTGGTTTCTGGCGTGTACGCTTCGAAAGAAGAgg                                                   >  1:438460/1‑68 (MQ=255) tGAAACGACGCGTAATGGTCAGCCGTGGTATGTCCTGGTTTCTGGCGTGTACGCTTCGAAAGAAGAgg                                                   >  1:1783072/1‑68 (MQ=255) tGAAACGACGCGTAATGGTCAGCCGTGGTATGTCCTGGTTTCTGGCGTGTACGCTTCGAAAGAAGAgg                                                   >  1:53021/1‑68 (MQ=255) tGAAACGACGCGTAATGGTCAGCCGTGGTATGTCCTGGTTTCTGGCGTGTACGCTTCGAAAGAAGAgg                                                   >  1:2420912/1‑68 (MQ=255) tGAAACGACGCGTAATGGTCAGCCGTGGTATGTCCTGGTTTCTGGCGTGTACGCTTCGAAAGAAGAgg                                                   >  1:2494289/1‑68 (MQ=255) tGAAACGACGCGTAATGGTCAGCCGTGGTATGTCCTGGTTTCTGGCGTGTACGCTTCGAAAGAAGAgg                                                   >  1:3083650/1‑68 (MQ=255) tGAAACGACGCGTAATCGTCAGCCGTGGTATGTCCTGGTTTCTGGCGTGTACGCTTCGAAAGAAGAgg                                                   >  1:3226932/1‑68 (MQ=255)       gacgCGTAATGGTCAGCCGTGGTATGTCCTGGTTTCTGGCGTGTATGCTTCGAAAGAAGAGGCGaaaaa                                            <  1:2370109/69‑1 (MQ=255)       gacgCGTAATGGTCAGCCGTGGTATGTCCTGGTTTCTGGCGTGTATGCTTCGAAAGAAGAGGCGaaaaa                                            <  1:3223518/69‑1 (MQ=255)       gacgCGTAATGGTCAGCCGTGGTATGTCCTGGTTTCTGGCGTGTATGCTTCGAAAGAAGAGGCGaaaaa                                            <  1:2773880/69‑1 (MQ=255)       gacgCGTAATGGTCAGCCGTGGTATGTCCTGGTTTCTGGCGTGTATGCTTCGAAAGAAGAGGCGaaaaa                                            <  1:3322036/69‑1 (MQ=255)       gacgCGTAATGGTCAGCCGTGGTATGTCCTGGTTTCTGGCGTGTATGCTTCGAAAGAAGAGGCGaaaaa                                            <  1:2818385/69‑1 (MQ=255)       gacgCGTAATGGTCAGCCGTGGTATGTCCTGGTTTCTGGCGTGTATGCTTCGAAAGAAGAGGCGaaaaa                                            <  1:1043402/69‑1 (MQ=255)       gacgCGTAATGGTCAGCCGTGGTATGTCCTGGTTTCTGGCGTGTATGCTTCGAAAGAAGAGGCGaaaaa                                            <  1:352476/69‑1 (MQ=255)       gacgCGTAATGGTCAGCCGTGGTATGTCCTGGTTTCTGGCGTGTATGCTTCGAAAGAAGAGGCGaaaaa                                            <  1:363664/69‑1 (MQ=255)       gacgCGTAATGGTCAGCCGTGGTATGTCCTGGTTTCTGGCGTGTATGCTTCGAAAGAAGAGGCGaaaaa                                            <  1:1964268/69‑1 (MQ=255)       gacgCGTAATGGTCAGCCGTGGTATGTCCTGGTTTCTGGCGTGTATGCTTCGAAAGAAGAGGCGaaaaa                                            <  1:1784522/69‑1 (MQ=255)       gacgCGTAATGGTCAGCCGTGGTATGTCCTGGTTTCTGGCGTGTATGCTTCGAAAGAAGAGGCGaaaaa                                            <  1:1511291/69‑1 (MQ=255)       gacgCGTAATGGTCAGCCGTGGTATGTCCTGGTTTCTGGCGTGTATGCTTCGAAAGAAGAGGCGaaaaa                                            <  1:100691/69‑1 (MQ=255)       gacgCGTAATGGTCAGCCGTGGTATGTCCTGGTTTCTGGCGTGTATGCTTCGAAAGAAGAGGCGaaaaa                                            <  1:1361603/69‑1 (MQ=255)        acgCGTAATGGTCAGCCGTTGTATGTCCTGGTTTCTGGCGTGTATGCTTCGAAAGAAGAGGCGaaaaa                                            <  1:2563/68‑1 (MQ=255)               aTGGTCAGCCGTGGTATGTCCTGGTTTCTGGCGTGTACGCTTGGAaag                                                         <  1:3276341/48‑1 (MQ=255)               aTGGTCAGCCGTGGTATGTCCTGGTTTCTGGCGTGTACGCTTCGAaag                                                         <  1:664015/48‑1 (MQ=255)                  gTCAGCCGTGGTATGTCCTGGTTTCTGGCGTGTACGCTTCGAAAGAAGAGGCGAAAAAAGCGGTATCTa                                 >  1:1139510/1‑69 (MQ=255)                    cAGCCGTGGTATGTCCTGGTTTCTGGCGTGTACGCTTCGAAAGAAGAGGCGAAAAAAGCGGTATCTAc                                <  1:2876643/68‑1 (MQ=255)                    cAGCCGTGGTATGTCCTGGTTTCTGGCGTGTACGCTTCGAAAGAAGAGGCGAAAAAAGCGGTATCTAc                                <  1:837117/68‑1 (MQ=255)                    cAGCCGTGGTATGTCCTGGTTTCTGGCGTGTACGCTTCGAAAGAAGAGGCGAAAAAAGCGGTATCTAc                                <  1:3075923/68‑1 (MQ=255)                    cAGCCGTGGTATGTCCTGGTTTCTGGCGTGTACGCTTCGAAAGAAGAGGCGAAAAAAGCGGTATCTAc                                <  1:573777/68‑1 (MQ=255)                    cAGCCGTGGTATGTCCTGGTTTCTGGCGTGTACGCTTCGAAAGAAGAGGCGAAAAAAGCGGTATCTAc                                <  1:519626/68‑1 (MQ=255)                    cAGCCGTGGTATGTCCTGGTTTCTGGCGTGTACGCTTCGAAAGAAGAGGCGAAAAAAGCGGTATCTAc                                <  1:3259450/68‑1 (MQ=255)                    cAGCCGTGGTATGTCCTGGTTTCTGGCGTGTACGCTTCGAAAGAAGAGGCGAAAAAAGCGGTATCTAc                                <  1:390124/68‑1 (MQ=255)                    cAGCCGTGGTATGTCCTGGTTTCTGGCGTGTACGCTTCGAAAGAAGAGGCGAAAAAAGCGGTATCTAc                                <  1:2039272/68‑1 (MQ=255)                    cAGCCGTGGTATGTCCTGGTTTCTGGCGTGTACGCTTCGAAAGAAGAGGCGAAAAAAGCGGTATCTAc                                <  1:2872012/68‑1 (MQ=255)                    cAGCCGTGGTATGTCCTGGTTTCTGGCGTGTACGCTTCGAAAGAAGAGGCGAAAAAAGCGGTATCTAc                                <  1:2274857/68‑1 (MQ=255)                    cAGCCGTCGTATGTCCTGGTTTCTGGCGTGTACGCTTCGAAAGAAGAGGCGAAAAAAGCGGTATCTAc                                <  1:2219238/68‑1 (MQ=255)                              aTGTCCTGGTTTCTGGCGTGTACGCTTCGAAAGAAGAGGCGAAAAAAGCGGTATCTACATTGCcagca                      <  1:1048513/68‑1 (MQ=255)                              aTGTCCTGGTTTCTGGCGTGTACGCTTCGAAAGAAGAGGCGAAAAAAGCGGTATCTACATTGCcagca                      <  1:827019/68‑1 (MQ=255)                              aTGTCCTGGTTTCTGGCGTGTACGCTTCGAAAGAAGAGGCGAAAAAAGCGGTATCTACATTGCcagca                      <  1:2923186/68‑1 (MQ=255)                              aTGTCCTGGTTTCTGGCGTGTACGCTTCGAAAGAAGAGGCGAAAAAAGCGGTATCTACATTGCcagca                      <  1:679187/68‑1 (MQ=255)                              aTGTCCTGGTTTCTGGCGTGTACGCTTCGAAAGAAGAGGCGAAAAAAGCGGTATCTACATTGCcagca                      <  1:1052463/68‑1 (MQ=255)                              aTGTCCTGGTTTCTGGCGTGTACGCTTCGAAAGAAGAGGCGAAAAAAGCGGTATCTACATTGCcagca                      <  1:491869/68‑1 (MQ=255)                              aTGTCCTGGTTTCTGGCGTGTACGCTTCGAAAGAAGAGGCGAAAAAAGCGGTATCTACATTGCcagca                      <  1:300789/68‑1 (MQ=255)                              aTGTCCTGGTTTCTGGCGTGTACGCTTCGAAAGAAGAGGCGAAAAAAGCGGTATCTACATTGCcagca                      <  1:3415007/68‑1 (MQ=255)                              aTGTCCTGGTTTCTGGCGTGTACGCTTCGAAAGAAGAGGCGAAAAAAGCGGTATCTACATTGCcagca                      <  1:3308593/68‑1 (MQ=255)                               tGTCCTGGTTTCTGGCGTGTACGCTTCGAAagaaga                                                     >  1:2821030/1‑36 (MQ=255)                               tGTCCTGGTTTCTGGCGTGTACGCTTCGAAAGAAGAGGCGAAAAAAGCGGTATCTACATTGCcagca                      >  1:2234736/1‑67 (MQ=255)                               tGTCCTGGTTTCTGGCGTGTACGCTTCGAAAGAAGAGGCGAAAAAAGCGGTATCTACATTGCcagca                      <  1:2006603/67‑1 (MQ=255)                               tGTCCTGGTTTCTGGCGTGTACGCTTCGAAAGAAGAGGCGAAAAAAGCGGTATCTACATTGCcagca                      >  1:2983604/1‑67 (MQ=255)                                gTCCTGGTTTCTGGCGTGTACGCTTCGAAAGAAGAGGCGAAAAAAGCGGTATCTACATTGCcagca                      <  1:2790990/66‑1 (MQ=255)                                 tCCTTGTTTCTGGCGTGTACGCTTCGAAAGAAGAGGCGAAAAAAGCGGTATCTAc                                <  1:418743/55‑1 (MQ=255)                                 tCCTGGTTTCTGGCGTGTACGCTTCGAAAGAAGAGGCGAAAAAAGCGGTATCt                                  >  1:2243564/1‑53 (MQ=255)                                 tCCTGGTTTCTGGCGTGTACGCTTCGAAAGAAGAGGCGAAAAAAGCGGTATCTa                                 >  1:2249874/1‑54 (MQ=255)                                 tCCTGGTTTCTGGCGTGTACGCTTCGAAAGAAGAGGCGAAAAAAGCGGTATCTa                                 >  1:2305253/1‑54 (MQ=255)                                 tCCTGGTTTCTGGCGTGTACGCTTCGAAAGAAGAGGCGAAAAAAGCGGTATCTa                                 >  1:828079/1‑54 (MQ=255)                                 tCCTGGTTTCTGGCGTGTACGCTTCGAAAGAAGAGGCGAAAAAAGCGGTATCTa                                 >  1:2903182/1‑54 (MQ=255)                                 tCCTGGTTTCTGGCGTGTACGCTTCGAAAGAAGAGGCGAAAAAAGCGGTATCTa                                 >  1:2932879/1‑54 (MQ=255)                                 tCCTGGTTTCTGGCGTGTACGCTTCGAAAGAAGAGGCGAAAAAAGCGGTATCTa                                 >  1:1327569/1‑54 (MQ=255)                                        ttCTGGCGTGTACGCTTCGAAAGAAGAGGCGAAAAAAGCGGTATCTACATTGCCAGCAGATGTCCAg             <  1:3168796/67‑1 (MQ=255)                                        ttCTGGCGTGTACGCTTCGAAAGAAGAGGCGAAAAAAGCGGTATCTACATTGCCAGCAGATGTCCAg             <  1:2232981/67‑1 (MQ=255)                                        ttCTGGCGTGTACGCTTCGAAAGAAGAGGCGAAAAAAGAGGTATCTACATTGCCAGCAGATGTCCAg             <  1:2827190/67‑1 (MQ=255)                                         tCTGGCGTGTACGCTTCGAAAGAAGAGGCGAAAAAAGCGGTATCTACATTGCCAGCAGATGTCCAGGc           >  1:3001947/1‑68 (MQ=255)                                         tCTGGCGTGTACGCTTCGAAAGAAGAGGCGAAAAAAGCGGTATCTACATTGCCAGCAGATGTCCAGGc           >  1:1427481/1‑68 (MQ=255)                                         tCTGGCGTGTACGCTTCGAAAGAAGAGGCGAAAAAAGCGGTATCTACATTGCCAGCAGATGTCCAGGc           >  1:355915/1‑68 (MQ=255)                                         tCTGGCGTGTACGCTTCGAAAGAAGAGGCGAAAAAAGCGGTATCTACATTGCCAGCAGATGTCCAGGc           >  1:2505143/1‑68 (MQ=255)                                                 gtATGCTTCGAAAGAAGAGGCGAAAAAAGCGGTATCTACATTGCCAGCAGATGTCCAGACCAAAAAccc  >  1:1421325/1‑69 (MQ=255)                                                 gtATGCTTCGAAAGAAGAGGCGAAAAAAGCGGTATCTACATTGCCAGCAGATGTCCAGACCAAAAAccc  >  1:541870/1‑69 (MQ=255)                                                 gtATGCTTCGAAAGAAGAGGCGAAAAAAGCGGTATCTACATTGCCAGCAGATGTCCAGACCAAAAAccc  >  1:328365/1‑69 (MQ=255)                                                    cGCTTCGAAAGAAGAGGCGAAAAAAGCGGTATCTACATTGCcagcag                     <  1:2940048/47‑1 (MQ=255)                                                    cGCTTCGAAAGAAGAGGCGAAAAAAGCGGTATCTACATTGCcagcag                     <  1:2027733/47‑1 (MQ=255)                                                    |                                                                  TGAAACGACGCGTAATGGTCAGCCGTGGTATGTCCTGGTTTCTGGCGTGTACGCTTCGAAAGAAGAGGCGAAAAAAGCGGTATCTACATTGCCAGCAGATGTCCAGGCCAAAAACCC  >  minE/2704917‑2705033

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑