Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	92,975	T→C	K299K (AAA→AAG)	unknown ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	92,975	0	T	C	94.7%	43.8 / ‑6.5	19	K299K (AAA→AAG)	unknown	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (18/0);  total (18/1)
Fisher's exact test for biased strand distribution p-value = 5.26e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCCGGCTGCAATCGGATAAACGGTGGTTAAAGTGCGGTCATAATTGCCCACGACCGCGCGTTTGAGCTTGATTTTGTTGGATTCAAGCGCACTAATCTCGACCATTTCTTTGTTTTTTCCTGTAATCAAAATGGCACGGCCACCGACGGCAAAATCATAGCCTGCGGTGTTAATAGCAAGCTCTGTGGCGCCCTGTTGCACATTTTGCAACAACTGGGCGCTGTCGGTAAAAATCGGCATCGACCACACGCGCGAGCCATAACCATACAGCGCCGCTTCAAAGCGCTGACGCTCAA  >  1/92914‑93210                                                              |                                                                                                                                                                                                                                            ccgccgaacgcagcgggtAAACAGTGGTTAAAGTGCGCTCATAATTGCCCACGACCGCGCGCTTGAGCTTGATTTTACTCGGTTCAAGGGCGCTGATTTCGACCATTTCTTTGTTGCTTGCGGTCATCAAAATCGCGCGCCCGCCGATGGCAAAATCATAGCCTGCGGTGTTAATGGCAAGCTCTGTCGCCCCTTGTTGCACATTTTGCAACAACTGCGCGCTGTCGGTAAAAATc                                                               >  1:839159/18‑236 (MQ=255)        ccgcagcGGATAAACGGTGTTTAAAGTGCGCTAATATTTCCCCACGCCCGCGCGTTTGAGTTTGTTTTTCCTCGGTTCAGGGGCGTTGTTTTCGACCTTTTTTTTGTTGCTTGCGGTCATAAAAACCGCGCGCCCGCCGATGGCAAAATCATAGCCTGCGGTGTTAATGGCAAGCTCTGTCGCCCCTTGTTGCACATTTTGCAACAACTGCGCGCTGTCGGTAAAAATCGGCATCGACCATACGCGCGAGc                                         <  1:653339/245‑1 (MQ=255)                ggtAAACAGTGGTTAAAGTGCGCTCATAATTGCCCACGACCGCGCGCTTGAGCTTGATTTTACTCGGTTCAAGGGCGCTGATTTCGACCATTTCTTTGTTGCTTGCGGTCATCAAAATCGCGCGCCCGCCGATGGCAAAATCATAGCCTGCGGTGTTAATGGCAAGCTCTGTCGCCCCTTGTTGCACATTTTGCAACAACTGCGCGCTGTCGGTAAAAATCGGCATCGACCATAcg                                                >  1:187236/3‑236 (MQ=255)                ggtAAACAGTGGTTAAAGTGCGCTCATAATTGCCCACGACCGCGCGCTTGAGCTTGATTTTACTCGGTTCAAGGGCGCTGATTTCGACCATTTCTTTGTTGCTTGCGGTCATCAAAATCGCGCGCCCGCCGATGGCAAAATCATAGCCTGCGGTGTTAATGGCAAGCTCTGTCGCCCCTTGTTGCACATTTTGCAACAACTGCGCGCTGTCGGTAAAAATCGGCATCGACCATACGCGCGAGCCATAACCa                                 >  1:886681/3‑251 (MQ=255)                ggtAAACAGTGGTTAAAGTGCGCTCATAATTGCCCACGACCGCGCGCTTGAGCTTGATTTTACTCGGTTCAAGGGCGCTGATTTCGACCATTTCTTTGTTGCTTGCGGTCATCAAAATCGCGCGCCCGCCGATGGCAAAATCATAGCCTGCGGTGTTAATGGCAAGCTCTGTCGCCCCTTGTTGCACATTTTGCAACAACTGCGCGCTGTCGGTAAAAATCGGCATCGACCATACGCGCGAGCCATAACCa                                 >  1:80681/3‑251 (MQ=255)                ggtAAACAGTGGTTAAAGTGCGCTCATAATTGCCCACGACCGCGCGCTTGAGCTTGATTTTACTCGGTTCAAGGGCGCTGATTTCGACCATTTCTTTGTTGCTTGCGGTCATCAAAATCGCGCGCCCGCCGATGGCAAAATCATAGCCTGCGGTGTTAATGGCAAGCTCTGTCGCCCCTTGTTGCACATTTTGCAACAACTGCGCGCTGTCGGTAAAAATCGGCATCGACCATACGCGCGAGCCATAACCa                                 >  1:249910/3‑251 (MQ=255)                ggtAAACAGTGGTTAAAGTGCGCTCATAATTGCCCACGACCGCGCGCTTGAGCTTGATTTTACTCGGTTCAAGGGCGCTGATTTCGACCATTTCTTTGTTGCTTGCGGTCATCAAAATCGCGCGCCCGCCGATGGCAAAATCATAGCCTGCGGTGTTAATGGCAAGCTCTGTCGCCCCTTGTTGCACATTTTGCAACAACTGCGCGCTGTCGGTAAAAATCGGCATCGACCATACGCGCGAGCCATAACCa                                 >  1:543314/3‑251 (MQ=255)                   aaaCAGTGGTTAAAGTGCGCTCATAATTGCCCACGACCGCGCGCTTGAGCTTGATTTTACTCGGTTCAAGGGCGCTGATTTCGACCATTTCTTTGTTGCTTGCGGTCATCAAAATCGCGCGCCCGCCGATGGCAAAATCATAGCCTGCGGTGTTAATGGCAAGCTCTGTCGCCCCTTGTTGCACATTTTGCAACAACTGCGCGCTGTCGGTAAAAATCGGCATCGACCATACGCGCGAGCCATAACCATag                              >  1:862101/1‑250 (MQ=255)                        gtggtTAAAGTGCGCTCATAATTGCCCACGACCGCGCGCTTGAGCTTGATTTTACTCGGTTCAAGGGCGCTGATTTCGACCATTTCTTTGTTGCTTGCGGTCATCAAAATCGCGCGCCCGCCGATGGCAAAATCATAGCCTGCGGTGTTAATGGCAAGCTCTGTCGCCCCTTGTTGCACATTTTGCAACAACTGCGCGCTGTCGGTAAAAATCGGCATCGACCATACGCGCGAGCCATAACCATAGAgcg                          >  1:14596/1‑250 (MQ=255)                         tggtTAAAGTGCGCTCATAATTGCCCACGACCGCGCGCTTGAGCTTGATTTTACTCGGTTCAAGGGCGCTGATTTCGACCATTTCTTTGTTGCTTGCGGTCATCAAAATCGCGCGCCCGCCGATGGCAAAATCATAGCCTGCGGTGTTAATGGCAAGCTCTGTCGCCCCTTGTTGCACATTTTGCAACAACTGCGCGCTGTCGGTaaa                                                                   >  1:462896/1‑208 (MQ=255)                          ggtTAAAGTGCGCTCATAATTGCCCACGACCGCGCGCTTGAGCTTGATTTTACTCGGTTCAAGGGCGCTGATTTCGACCATTTCTTTGTTGCTTGCGGTCATCAAAATCGCGCGCCCGCCGATGGCAAAATCATAGCCTGCGGTGTTAATGGCAAGCTCTGTCGCCCCTTGTTGCACATTTTGCAACAACTGCGCGCTGTCGGTAAAAATCGGCATCGACCATACGCGCGAGCCATAAc                                   >  1:28546/1‑239 (MQ=255)                          ggtTAAAGTGCGCTCATAATTGCCCACGACCGCGCGCTTGAGCTTGATTTTACTCGGTTCAAGGGCGCTGATTTCGACCATTTCTTTGTTGCTTGCGGTCATCAAAATCGCGCGCCCGCCGATGGCAAAATCATAGCCTGCGGTGTTAATGGCAAGCTCTGTCGCCCCTTGTTGCACATTTTGCAACAACTGCGCGCTGTCGGTAAAAATCGGCATCGACCATACGCGCGAGCCATAACCATAGAgc                           >  1:491700/1‑247 (MQ=255)                                 gTGCGCTCATAATTGCCCACGACCGCGCGCTTGAGCTTGATTTTACTCGGTTCAAGGGCGCTGATTTCGACCATTTCTTTGTTGCTTGCGGTCATCAAAATCGCGCGCCCGCCGATGGCAAAATCATAGCCTGCGGTGTTAATGGCAAGCTCTGTCGCCCCTTGTTGCACATTTTGCAACAACTGCGCGCTGTCGGTAAAAATCGGCATCGAcc                                                     >  1:104724/1‑214 (MQ=255)                                     gctCATAATTGCCCACGACCGCGCGCTTGAGCTTGATTTTACTCGGTTCAAGGGCGCTGATTTCGACCATTTCTTTGTTGCTTGCGGTCATCAAAATCGCGCGCCCGCCGATGGCAAAATCATAGCCTGCGGTGTTAATGGCAAGCTCTGTCGCCCCTTGTTGCACATTTTGCAACAACTGCGCGCTGTCGGTAAAAATCGGCATCGACCATACGCGCGAGCCATAACCATAGAGCGCCGCTTCAAAgcg             >  1:752436/3‑250 (MQ=255)                                        cATAATTGCCCACGACCGCGCGCTTGAGCTTGATTTTACTCGGTTCAAGGGCGCTGATTTCGACCATTTCTTTGTTGCTTGCGGTCATCAAAATCGCGCGCCCGCCGATGGCAAAATCATAGCCTGCGGTGTTAATGGCAAGCTCTGTCGCCCCTTGTTGCACATTTTGCAACAACTGCGCGCTGTCGGTAAAAATCGGCATCGACCATACGCGCGAGCCATAACCATAGAGCGCCGCTTCAAAGCGCTGa         >  1:555352/1‑251 (MQ=255)                                        cATAATTGCCCACGACCGCGCGCTTGAGCTTGATTTTACTCGGTTCAAGGGCGCTGATTTCGACCATTTCTTTGTTGCTTGCGGTCATCAAAATCGCGCGCCCGCCGATGGCAAAATCATAGCCTGCGGTGTTAATGGCAAGCTCTGTCGCCCCTTGTTGCACATTTTGCAACAACTGCGCGCTGTCGGTAAAAATCGGCATCGACCATACGCGCGAGCCATAACCATAGAGCGCCGCTTCAAAGCGCTGa         >  1:788551/1‑251 (MQ=255)                                        cATAATTGCCCACGACCGCGCGCTTGAGCTTGATTTTACTCGGTTCAAGGGCGCTGATTTCGACCATTTCTTTGTTGCTTGCGGTCATCAAAATCGCGCGCCCGCCGATGGCAAAATCATAGCCTGCGGTGTTAATGGCAAGCTCTGTCGCCCCTTGTTGCACATTTTGCAACAACTGCGCGCTGTCGGTAAAAATCGGCATCGACCATACGCGCGAGCCATAACCATAGAGCGCCGCTTCAAAGCGCTGa         >  1:1960/1‑251 (MQ=255)                                               gCCCACGACCGCGCGCTTGAGCTTGATTTTACTCGGTTCAAGGGCGCTGATTTCGACCATTTCTTTGTTGCTTGCGGTCATCAAAATCGCGCGCCCGCCGATGGCAAAATCATAGCCTGCGGTGTTAATGGCAAGCTCTGTCGCCCCTTGTTGCACATTTTGCAACAACTGCGCGCTGTCGGTAAAAATCGGCATCGACCATACGCGCGAGCCATAAc                                   >  1:518439/1‑218 (MQ=255)                                               gCCCACGACCGCGCGCTTGAGCTTGATTTTACTCGGTTCAAGGGCGCTGATTTCGACCATTTCTTTGTTGCTTGCGGTCATCAAAATCGCGCGCCCGCCGATGGCAAAATCATAGCCTGCGGTGTTAATGGCAAGCTCTGTCGCCCCTTGTTGCACATTTTGCAACAACTGCGCGCTGTCGGTAAAAATCGGCATCGACCATACGCGCGAGCCATAACCATAGAGCGCCGCTTCAAAGCGCTGACGCTCaa  >  1:13718/1‑251 (MQ=255)                                                              |                                                                                                                                                                                                                                            TTCCGGCTGCAATCGGATAAACGGTGGTTAAAGTGCGGTCATAATTGCCCACGACCGCGCGTTTGAGCTTGATTTTGTTGGATTCAAGCGCACTAATCTCGACCATTTCTTTGTTTTTTCCTGTAATCAAAATGGCACGGCCACCGACGGCAAAATCATAGCCTGCGGTGTTAATAGCAAGCTCTGTGGCGCCCTGTTGCACATTTTGCAACAACTGGGCGCTGTCGGTAAAAATCGGCATCGACCACACGCGCGAGCCATAACCATACAGCGCCGCTTCAAAGCGCTGACGCTCAA  >  1/92914‑93210

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