Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	225,207	0	G	.	54.5%	7.6 / 19.6	11	coding (118/375 nt)	unknown	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (3/2);  new base . (4/2);  total (7/4)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.59e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CCAACTCCTTTAGGTAAGCAATCTGTGCACGTGCTAATGCCAGCTCTCGTTGCAGCTTTTTAAACGCTTTTGGGGAAAAGTCTGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATCTTTGGCTTCACTATTTTAGGGCGTTGAGGGTTTATATAAGGGGATTTTACGCCATTAAGCCCTCTTTCACGAAATGCTTTTAGCCAATAAATGACGAGCGCGTGGGAGATTTTATGAAGTTTCGCCACCTCACGGATACCAAAATCCTGCTCGGTGACCTGTTTGATAGTTTTTAAACGAAATTGATAAGAGTAAGACATAATCTGCACCTCAAATTAGGTGTCCGGGTTTTGGGGTGCAGATCACTTTTTTGTTTCGGAAAGCGGATTATTTAATTACACCACACGCCATACGCGGACCGCCGC  >  1/225006‑225429                                                                                                                                                                                                          |                                                                                                                                                                                                                               ccAACTCCTTTAGGTAAGCAATCTGTGCACGTGCTAATGCCAGCTCTCGTTGCAGCTTTTTAAACGCTTTTGGGGAAAAGTCTGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATCTTTGGCTTCACTATTTTAGGGCGTTGAGGGTTTATATAAGGGGATTTTACGCCATTAAGCCCTCTTTCACGAAATGCTTTTAGCCAATAAATGACGAGCGAATGGGAGATTTTATGAAGTTTCGcc                                                                                                                                                                                     >  1:477617/1‑245 (MQ=255)                                      tGCCAGCTCTCGTTGCAGCTTTTTAAACGCTTTTGGGGAAAAGTCTGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATCTTTGGCTTCACTATTTTAGGGCGTTGAGGGTTTATATAAGGGGATTTTACGCCATTAAGCCCTCTTTCACGAAATGCTTTTAGcc                                                                                                                                                                                                                              >  1:509052/1‑167 (MQ=35)                                      tGCCAGCTCTCGTTGCAGCTTTTTAAACGCTTTTGGGGAAAAGTCTGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATCTTTGGCTTCACTATTTTAGGGCGTTGAGGGTTTATATAAGGGGATTTTACGCCATTAAGCCCTCTTTCACGAAAGGCTTTT‑‑CCAATAAATGACGAGCGAATGGGAGATTTTATg                                                                                                                                                                                               >  1:822717/1‑196 (MQ=11)                                                 gTTGCAGCTTTTTAAACGCTTTTGGGGAAAAGTCTGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATCTTTGGCTTCACTATTTTAGGGCGTTGAGGGTTTATATAAGGGGATTTTACGCCATTAAGCCCTCTTTCACGAAATGCTTTTAGCCaa                                                                                                                                                                                                                            >  1:461201/1‑158 (MQ=37)                                                                   cTTTTGGGGAAAAGTCTGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATCTTTGGCTTCACTATTTTAGGGCGTTGAGGGTTTATATAAGGGGATTTTACGCCATTAAGCCCTCTTTCACGAAATGCTTTTAGCCAATAAATGACGAGCGAATGGGAGATTTTATGAAGTTTCGCCACCTCACGGATACCAAAATCCTGCTCGGTGACCTGTTTGATAGTTTTTAAACGAAATTGGTAAGAGTAAGAc                                                                                                             <  1:897128/251‑1 (MQ=18)                                                                                                                            gCTTCACTATTTTAGGGCGTTGAGGGTTTATATAAGGGGATTTTACGCCATTAAGCCCTCTTTCACGAAATGCTTTTAGCCAATAAATGAc                                                                                                                                                                                                                    <  1:843720/91‑1 (MQ=34)                                                                                                                                                                  gATTTTACGCCATTAAGCCCTCTTTCACGAAAGGCTTTT‑‑CCAATAAATGACTTGAGAACGAGAAATCTGATGAAGTTCAGCCACCTCACGGATACCAAAATCCTGCTCGGTGACCTGTTTGATAATTTCTAAACGAAATTGATAAGAGTAAGACATAATCTGCACCTCAAATTAGGTGTCCAGGTTTTGGGGTGCAGATCACTTTTTTGTTTCGGAAAGCGGATTATTTAATTACACCACACGCCATAcg             <  1:624914/250‑1 (MQ=255)                                                                                                                                                                   aTTTTACGCCATTAAGCCCTCTTTCACGAAAGGCTTTT‑‑CCAATAAATGACTTGAGAACGAGAAATCTGATGAAGTTCAGCCACCTCACGGATACCAAAATCCTGCTCGGTGACCTGTTTGATAATTTCTAAACGAAATTGATAAGAGTAAGACATAATCTGCACCTCAAATTAGGTGTCCAGGTTTTGGGGTGCAGATCACTTTTTTGTTTCGGAAAGCGGATTATTTAATTACACCACACGCCATAcgcg           >  1:661990/1‑251 (MQ=11)                                                                                                                                                                    ttttACGCCATTAAGCCCTCTTTCACGAAAGGCTTTT‑‑CCAATAAATGACTTGAGAACGAGAAATCTGATGAAGTTCAGCCACCTCACGGATACCAAAATCCTGCTCGGTGACCTGTTTGATAATTTCTAAACGAAATTGATAAGAGTAAGACATAATCTGCACCTCAAATTAGGTGTCCAGGTTTTGGGGTGCAGATCACTTTTTTGTTTCGGAAAGCGGATTATTTAATTACACCACACGCCATAcg             <  1:226559/248‑1 (MQ=255)                                                                                                                                                                             aTTAAGCCCTCTTTCACGAAAGGCTTTT‑‑CCAATAAATGACTTGAGAACGAGAAATCTGATGAAGTTCAGCCACCTCACGGATACCAAAATCCTGCTCGGTGACCTGTTTGATAATTTCTAAACGAAATTGATAAGAGTAAGACATAATCTGCACCTCAAATTAGGTGTCCAGGTTTTGGGGTGCAGATCACTTTTTTGTTTCGGAAAGCGGATTATTTAATTACACCACACGCCATAcgc            >  1:303672/1‑240 (MQ=11)                                                                                                                                                                             aTTAAGCCCTCTTTCACGAAAGGCTTTT‑‑CCAATAAATGACTTGAGAACGAGAAATCTGATGAAGTTCAGCCACCTCACGGATACCAAAATCCTGCTCGGTGACCTGTTTGATAATTTCTAAACGAAATTGATAAGAGTAAGACATAATCTGCACCTCAAATTAGGTGTCCAGGTTTTGGGGTGCAGATCACTTTTTTGTTTCGGAAAGCGGATTATTTAATTACACCACACGCCATACGCGGAccgccgc  >  1:675978/1‑250 (MQ=11)                                                                                                                                                                                                          |                                                                                                                                                                                                                               CCAACTCCTTTAGGTAAGCAATCTGTGCACGTGCTAATGCCAGCTCTCGTTGCAGCTTTTTAAACGCTTTTGGGGAAAAGTCTGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATCTTTGGCTTCACTATTTTAGGGCGTTGAGGGTTTATATAAGGGGATTTTACGCCATTAAGCCCTCTTTCACGAAATGCTTTTAGCCAATAAATGACGAGCGCGTGGGAGATTTTATGAAGTTTCGCCACCTCACGGATACCAAAATCCTGCTCGGTGACCTGTTTGATAGTTTTTAAACGAAATTGATAAGAGTAAGACATAATCTGCACCTCAAATTAGGTGTCCGGGTTTTGGGGTGCAGATCACTTTTTTGTTTCGGAAAGCGGATTATTTAATTACACCACACGCCATACGCGGACCGCCGC  >  1/225006‑225429

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