Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,879,047	C→T	E242E (GAG→GAA)	cas1 ←	multifunctional nuclease Cas1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,879,047	0	C	T	88.9%	16.7 / NA	9	E242E (GAG→GAA)	cas1	multifunctional nuclease Cas1
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (5/3);  minor base G (1/0);  total (6/3)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTC‑ACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCTTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCCTTTTCCCAGTCTTTCGGATCGTAGC  >  NC_000913/2878816‑2879282                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                             gggggTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGACTGCATCTCCCAGTGCAACAGGAAGCGGAATGGCAACAGGCTGTGCCTCTTCAGGTGGGGCCGGCGGTTGCATTTCTCAAGCCGCCAGGAAGTACTCTATAACCGGAATCCATTTTGGTATTGTTTTAATACTGCGAAAAATATTCCTGCCAGCCAAACGGAATTCCCGGTACGGTTC‑CCCAGGGATACGAcgcg                                                                                                                                                                                                                          >  1:189658/4‑251 (MQ=255)    ggTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTTTCCAGCCGCAAACACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGTTC‑ACCAGGGttttgctatcaa                                                                                                                                                                                                                        >  1:195294/1‑240 (MQ=255)                                    cATTTCAGCTACTCAGATCGCCTGCATCTCACAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATTAACGGGATCAATTTGGCTAATTTTTTAATACTGCGGAAAATATCCCTGCACGCCAAACGGACTTTCCGGTCCGGGTCCACAGGGGTTA‑GACGCGCTATCTCCAAAACTTTTGGTacaaaaa                                                                                                                                                                                               >  1:72434/1‑240 (MQ=255)                                          nnnnnnncnaanncccTNCATCTCCCCGTGAAACAGGAAGCGNAATGGCAACGGGGTGTGCATCTGCAGGGGGGNCNGGCGGTTGTTTTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCNATTTGGGTAANGTTNNACTANTNCNAAAAATAGCCCTGCGCGCCAAACGGACTTNCCGGTCCGGTTC‑ACCAGGGTTACNNCGNNNTGNNTCaaa                                                                                                                                                                                                                <  2:189658/206‑1 (MQ=255)                                                                                                                                                                                         ctactGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGTTA‑ACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCTTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGt                                  >  1:188350/1‑251 (MQ=255)                                                                                                                                                                                                             gCACGCCAAACGGACTTCCCGGTCCGGTTC‑ACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATAGCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATcgc                                                                                           <  1:171434/174‑1 (MQ=255)                                                                                                                                                                                                                cGCCAAACGGACTTCCCGGTCCGGTTC‑ACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTAt                                                                                                                                      >  1:209132/1‑128 (MQ=255)                                                                                                                                                                                                                         gACTTCCCGGTCCGGTTC‑ACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACCACCGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAACTGGTGCATAACCAGCTGCAAGTATCGCCGCTTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCCTTTTCCCAGTCTTTCGGATCGTAGc  >  1:112890/1‑251 (MQ=255)                                                                                                                                                                                                                           cTTCCCGGTCCGGTTC‑ACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAggn                                                                                                                                                <  1:66456/107‑2 (MQ=255)                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                             GGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTC‑ACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCTTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCCTTTTCCCAGTCTTTCGGATCGTAGC  >  NC_000913/2878816‑2879282

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