Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,103,636	T→C	C208C (TGT→TGC)	mutY →	adenine DNA glycosylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,103,636	0	T	C	94.4%	53.5 / ‑4.1	18	C208C (TGT→TGC)	mutY	adenine DNA glycosylase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (9/8);  total (10/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.89e-01

AGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGATGGATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGAAATGTTCGCTCTGTCCGCTACAAAACGGATGTATTGCCGCCGCCAACAATAGCTGGGCGCTTTATCCGGGCAAAAAACCGAAACAGACGCTGCCGGAGCGCACCGGCTACTT  >  NC_000913/3103508‑3103716                                                                                                                                 |                                                                                 aGTTTGGGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGATGGATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGAAATGTTCGCTCTGTCCGCTACAAAACGGATGCATTgccgccg                                                                        <  2:542671/139‑1 (MQ=255)     tGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGATGGATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGAAATGTTCGCTCTGTCCGCTACAAAACGGATGCATTGCCGCCGCcaa                                                                    >  1:39251/1‑139 (MQ=255)     tGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGATGGATTTGGGTGCGATGATTAGTACGCGCTTGAAACCGAAATGTTCGCTCTGTCCGCTACAAAACGGATGTATTGCCGCCGCcaa                                                                    >  2:412034/1‑139 (MQ=255)                      cccGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGATGGATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGAAATGTTCGCTCTGTCCGCTACAAAACGGATGCATTGCCGCCGCCAACAATAGCTGGGCGCttt                                                   >  2:749071/1‑139 (MQ=255)                        cgcgGTTGGCGTGGAACGGTTTAATCAGGCGATGATGGATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGAAATGTTCGCTCTGTCCGCTACAAAACGGATGCATTGCCGCCGCCAACAATAGCTGGGCGCTTTAt                                                 <  1:749071/139‑1 (MQ=255)                                              aaTCAGGCGATGATGGATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGAAATGTTCGCTCTGTCCGCTACAAAACGGATGCATTGCCGCCGCCAACAATAGCTGGGCGCTTTATCCGGGCAAAAAACCGAAACAGa                           <  2:213896/139‑1 (MQ=255)                                              aaTCAGGCGATGATGGATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGAAATGTTCGCTCTGTCCGCTACAAAACGGATGCATTGCCGCCGCCAACAATAGCTGGGCGCTTTATCCGGGCAAAAAACCGAAACAGa                           <  2:39251/139‑1 (MQ=255)                                                   ggCGATGATGGATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGAAATGTTCGCTCTGTCCGCTACAAAACGGATGCATTgccg                                                                           <  1:879872/86‑1 (MQ=255)                                                   ggCGATGATGGATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGAAATGTTCGCTCTGTCCGCTACAAAACGGATGCATTgccg                                                                           >  2:879872/1‑86 (MQ=255)                                                    gCGATGATGGATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGAAATGTTCGCTCTGTCCGCTACAAAACGGATGCATTGCCGCCGCCAACAATAGCTGGGCGCTTTATCCGGGCAAAAAACCGAAACAGACGCTGc                     >  2:102217/1‑139 (MQ=255)                                                    gCGATGATGGATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGAAATGTTCGCTCTGTCCGCTACAAAACGGATGCATTGCCGCCGCCAACAATAGCTGGGCGCTTTATCCGGGCAAAAAACCGAAACAGACGCTGc                     >  2:972073/1‑139 (MQ=255)                                                            ggATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGAAATGTTCGCTCTGTCCGCTACAAAACGGATGCATTGCCGCCGCCAACAATAGCTGGGCGCTTTATCCGGGCAAAAAACCGAAACAGACGCTGCCGGAgcgc             <  1:102217/139‑1 (MQ=255)                                                                  gggTGCGATGATTTGTACGCGCTCGAAACCGAAATGTTCGCTCTGTCCGCTACAAAACGGATGCATTGCCGCCGCCAACAATAGCTGGGCGCTTTATCCGGGCAAAAAACCGAAACAGACGCTGCCGGAGCGCACCGGc       >  2:594204/1‑139 (MQ=255)                                                                    gTGCGATGATTTGTACGCGCTCGAAACCGAAATGTTCGCTCTGTCCGCTACAAAACGGATGCATTGCCGCCGCCAACAATAGCTGGGCGCTTTATCCGGGCAAAAAACCGAAACAGACGCTGCCGGAGCGCACCGGCTa     >  2:689885/1‑139 (MQ=255)                                                                       cGATGATTTGTACGCGCTCGAAACCGAAATGTTCGCTCTGTCCGCTACAAAACGGATGCATTGCCGCCGCCAACAATAGCTGGGCGCTTTATCCGGGCAAAAAACCGAAACAGACGCTGCCGGAGCGCACCGGCTACtt  >  1:920417/1‑139 (MQ=255)                                                                       cGATGATTTGTACGCGCTCGAAACCGAAATGTTCGCTCTGTCCGCTACAAAACGGATGCATTGCCGCCGCCAACAATAGCTGGGCGCTTTATCCGGGCAAAAAACCGAAACAGACGCTGCCGGAGCGCACCGGCTACtt  <  2:920417/139‑1 (MQ=255)                                                                        gatgatTTGTACGCGCTCGAAACCGAAATGTTCGCTCTGTCCGCTACAAAACGGATGCATTGCCGCCGCCAACAATAGCTGGGCGCTTTATcc                                               >  1:96711/1‑93 (MQ=255)                                                                        gatgatTTGTACGCGCTCGAAACCGAAATGTTCGCTCTGTCCGCTACAAAACGGATGCATTGCCGCCGCCAACAATAGCTGGGCGCTTTATcc                                               <  2:96711/93‑1 (MQ=255)                                                                                                                                 |                                                                                 AGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGTGGAACGGTTTAATCAGGCGATGATGGATTTGGGTGCGATGATTTGTACGCGCTCGAAACCGAAATGTTCGCTCTGTCCGCTACAAAACGGATGTATTGCCGCCGCCAACAATAGCTGGGCGCTTTATCCGGGCAAAAAACCGAAACAGACGCTGCCGGAGCGCACCGGCTACTT  >  NC_000913/3103508‑3103716

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