Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,760,367	G→T	24.5%	D294E (GAC→GAA)	selA ←	selenocysteine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,760,367	0	G	T	24.5%	95.4 / 4.3	65	D294E (GAC→GAA)	selA	selenocysteine synthase
Reads supporting (aligned to +/- strand):  ref base G (20/27);  new base T (14/3);  total (34/31)
Fisher's exact test for biased strand distribution p-value = 9.46e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.98e-01

TGTAAACGTTGTGCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAA‑GACGC‑AACGTGGCTTCCAGCGCCGC‑‑GAGGG‑TCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCA‑GGCGGGCG‑ATCATCTCTTTT‑TT‑A‑CCA‑AC‑AATAATTCCTGCCTGCG‑GC‑CCGCCTA‑ACAAC‑TTGTCGC‑CGGAGAAAC‑TCACCAGACTGACGC‑CCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTTTCATTC  >  NC_000913/3760134‑3760604                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                 tGTAAACGTTGTGCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAA‑GACGC‑AACGTGGCTTCCAGCGCCGC‑‑GAGGG‑TCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCA‑GGCGGGCG‑ATCATCTCTTTT‑TT‑A‑CCA‑AC‑CATAATTCCTGCCTGCG‑GG‑CCGCCTC‑ACAAC‑TTTTCGc‑c                                                                                                                                                                                                                                            >  2:654286/1‑239 (MQ=255)   tAAACGTTGTGCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAA‑GACGC‑AACGTGGCTTCCAGCGCCGC‑‑GAGGG‑TCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCA‑GGCGGGCG‑ATCCTCTCTTTT‑TT‑T‑CCC‑AC‑AATAATTCCTGCCTGCG‑GC‑CCGCCTA‑ACAAA‑TTTTGGC‑Cgg                                                                                                                                                                                                                                          >  2:17228/1‑239 (MQ=255)   tAAACGTTGTGCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAA‑GACGC‑AACGTGGCTTCCAGCGCCGC‑‑GAGGG‑TCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCA‑GGCGGGCG‑ATCATTCCTTTT‑TT‑T‑ACA‑AC‑AAAAAATTCCGCCCGCG‑GC‑CCGCCCC‑ACAAA‑TTTTTGc‑ccg                                                                                                                                                                                                                                          >  1:33428/1‑237 (MQ=255)          tgtgCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAA‑GACGC‑AACGTGGCTTCCAGCGCCGC‑‑GAGGG‑TCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCA‑GGCGGGCG‑ATCCTCTCTTTT‑TT‑A‑CCA‑AC‑AATAATTCCTGCCCGCG‑GG‑CCGCCCA‑ACAAA‑TTTTCGC‑CGGAaaaaaa                                                                                                                                                                                                                                   >  1:242498/1‑237 (MQ=255)          tgtgCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAA‑GACGC‑AACGTGGCTTCCAGCGCCGC‑‑GAGGG‑TCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCA‑GGCGGGCG‑ATCATCTCTTTT‑TT‑A‑CCA‑AC‑AATAATTCCTGCCTGCG‑GC‑CCGCCTACACAA‑‑TTTTCGC‑CGGGGAAAc                                                                                                                                                                                                                                    >  2:135750/1‑238 (MQ=255)           gtgCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAA‑GACGC‑AACGTGGCTTCCAGCGCCGC‑‑GAGGG‑TCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGGGGCTTTGCA‑GGCGGGCG‑ATCATCTCTTTT‑TT‑A‑CCA‑AC‑AATAATTCCTGCCTGCG‑GC‑CCGCCTT‑ACAAC‑TTTTCGC‑CGGGGAAAC‑Tc                                                                                                                                                                                                                                 >  2:289769/1‑239 (MQ=255)            ttccTGGATTTGAATGCCCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAA‑GACGC‑AACGTGGCTTCCAGCGCCGC‑‑GAGGG‑TCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCA‑GGCGGGCG‑ATCATCTCTTTT‑TT‑A‑CCA‑AC‑AATAATTCCTGCCTGCG‑GC‑CCGCCTA‑ACAAC‑TTGTCGC‑CGGAGAAAC‑TCa                                                                                                                                                                                                                                <  2:22498/237‑1 (MQ=255)            tgCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAA‑GACGC‑AACGTGGCTTCCAGCGCCGC‑‑GAGGG‑TCATTTTATCCGCCCGGAATGCACGCTTCAGCGGGGGGCTTTGCA‑GGCCGGCG‑CT‑ATCTCCTTT‑TT‑T‑CCACAC‑AACAATGCCCACCCGGG‑GC‑CGCCCCA‑ACAAA‑TTTTCCC‑CGGAGaaagacc                                                                                                                                                                                                                                 >  1:397244/1‑235 (MQ=255)             gCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCCGGGTGTAAATAAA‑GACGC‑AACGTGGCTTCCAGCGCCGC‑‑GAGGG‑TCATTTTATCCGCGCGTAATGCACGCTTCAGCGGCTGGCTTTGCC‑GGCGGGCG‑ATCATCTCTTTT‑TT‑A‑CCA‑AC‑AATAATTCCTGCCTCCG‑GC‑ACGCCAC‑ACACC‑TTTTCCC‑CCGGGAATC‑TCa                                                                                                                                                                                                                                >  2:252867/1‑238 (MQ=255)             gCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAA‑GACGC‑AACGTGGCTTCCAGCGCCGC‑‑GAGGG‑TCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCA‑GGCGGGCG‑ATCATCTCTTTT‑TT‑A‑CCA‑AC‑AATAATTCCTGCCTGCG‑GC‑CCGCCTA‑ACAAC‑TTGTGGC‑CGGGGAAAA‑TCa                                                                                                                                                                                                                                >  2:679944/1‑238 (MQ=255)                  gATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAA‑GACGC‑AACGTGGCTTCCAGCGCCGC‑‑GAGGG‑TCATTTTATCCGCGCGGAATGCACGCTTCAGCGGGTGGCTTTGCA‑GGCGGGCG‑AGCAACTCTTTT‑TT‑T‑ACA‑AA‑AACAATTCCTGCCTGCG‑GC‑CCCCCCA‑AAAAA‑CTTTTGC‑CGGGGAAAA‑ACACCacaa                                                                                                                                                                                                                          >  1:200311/1‑236 (MQ=255)                             ccTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAA‑GACGC‑AACGTGGCTTCCAGCGCCGC‑‑GAGGG‑TCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCA‑GGCGGGCG‑ATCATCCCTTTT‑TT‑T‑CCA‑AC‑AATAATTCCTTCCTCCC‑GC‑CCGCCTA‑AAACA‑TTTTTAC‑CGGAGaaaacccccacacatgccccccc                                                                                                                                                                                                                 >  1:639820/1‑218 (MQ=255)                                        cgGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAA‑GACGC‑AACGTGGCTTCCAGCGCCGC‑‑GAGGG‑TCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCA‑GGCGGGCG‑GTCATCCCTTTT‑TT‑A‑CCA‑AC‑AAAAATTCCTGCCTTCG‑GC‑CCGCCTA‑AACAC‑TTGTCCC‑CGGAGAAAC‑TCACCAACCATACGC‑CCGCCCCaactaaac                                                                                                                                                                                                    >  2:263677/1‑232 (MQ=255)                                        cgGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAA‑GACGC‑AACGTGGCTTCCAGCGCCGC‑‑GAGGG‑TCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCA‑GGCGGGCG‑ATCATCTCTTTT‑TT‑A‑CCA‑AC‑AATAATTCCTGCCTGCG‑GC‑CCGCCTA‑ACAAC‑TTTTCGC‑CGGAGAAAC‑TCAACAGACTGACGC‑CCGCCGCACTCAAcc                                                                                                                                                                                                    >  2:363342/1‑237 (MQ=255)                                        cgGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAA‑GACGC‑AACGTGGCTTCCAGCGCCGC‑‑GAGGG‑TCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCA‑GGCGGGCG‑ATACTCTCTTTT‑TT‑A‑CCA‑AC‑AAAAATTCCTGCCTGCG‑GC‑CCGCCTA‑ACAAC‑CTGTCCC‑CGCAGAAAA‑CCACCCAACTTAccc‑ccccccacaacacac                                                                                                                                                                                                    >  1:874836/1‑228 (MQ=255)                                                    cgcAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAA‑GACGC‑AACGTGGCTTCCAGCGCCGC‑‑GAGGG‑TCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCA‑GGCGGGCG‑ATCATCTCTTTT‑TT‑A‑CCA‑AC‑AATAATTCCTGCCTCCG‑GC‑CCGCCTA‑ACCAA‑TTTTCGC‑CCGAGAAAA‑TCACCCAACTGACGC‑CCCCCCCAAACAACCCCTGCGGGACTgg                                                                                                                                                                                       >  1:877231/1‑239 (MQ=255)                                                                     ttCACTCAGAGCTTCAGGGTGTAAATAAA‑GACGC‑AACGTGGCTTCCAGCGCCGC‑‑GAGGG‑TCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCA‑GGCGGGCG‑ATCATCTCTTTT‑TT‑A‑CCA‑AC‑AATAATTCCTGCCTGCG‑GC‑CCGCCTA‑ACAAC‑TTGTCGC‑CGGAGAAAC‑TCACCAGACTGACGC‑CCGCCGCAATCAACTCCTGCGGCATTTGCTCTTTCGGCCAACCGt                                                                                                                                                                      >  2:745271/1‑239 (MQ=255)                                                                     ttCACTCAGAGCTTCAGGGTGTAAATAAA‑GACGC‑AACGTGGCTTCCAGCGCCGC‑‑GAGGG‑TCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCA‑GGCGGGCG‑ATCATCTCTTTT‑TT‑A‑CCA‑AC‑AATAATTCCTGCCTGCG‑GC‑CCGCCTA‑ACAAC‑TTGTCGC‑CGGAGAAAC‑TCACCAGACTGACGC‑CCGCCGCAATCAACCCCCGCGTCATTGGCTCTTTCGGCaacacg                                                                                                                                                                       >  2:458107/1‑234 (MQ=255)                                                                      tcccacagggtgtgaGGGTGTAAATAAA‑GACGG‑CACGGGGCTGCCGTGGCGGG‑‑GGGGG‑GCATTTTCTCGGCGCTTAGTCACCTCTCCGGCGGGTGGCTTTGCG‑GGGGTGCG‑ACTATCTTTTTT‑TT‑C‑CCA‑AC‑AATAATCCCTGCCTGCG‑AC‑CCGCCTA‑ACACT‑TTGTCGC‑CGGAGAACC‑TCACCAGACAGACGC‑CCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTa                                                                                                                                                                     <  1:661577/225‑1 (MQ=255)                                                                      cccccaggagCTTGGGGGTG‑AAAAAAA‑GAGCCAAACGTGGCTCCCA‑CGCCCC‑‑GGGGGTTTTTTTTATCGGCGGTTAAGGCCCCCTTCAGGGGGTGGCTTGGCG‑GGCGGGCT‑TTCATCTTTTTT‑CT‑A‑CCA‑AC‑AATAATCCCTGCCTGCG‑GC‑CCGCCTA‑ACAAC‑TTGTCGC‑GGGAGAAAC‑TCACCAGCCTGACGC‑CCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTa                                                                                                                                                                     <  1:639443/232‑1 (MQ=255)                                                                      agttgtgtataggagaaagtaaaaaaaa‑GCCGC‑AGCGTGGCTCCCAGCGCGGG‑‑GAGGT‑TCATTGTATCCGTGCTTAATGCACGTTTCGGCGGGTGCTTTGGCG‑GGCGGGCG‑ATCATCTTTTTT‑TT‑A‑CCA‑AA‑AATAATCCCTGCCGGCG‑GC‑CCGCCTA‑ACAAC‑TTGTCGC‑CGGAGAAAC‑CCCCCAGCCTGACGC‑CCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTa                                                                                                                                                                     <  1:672882/218‑1 (MQ=255)                                                                       tcctcagggtttgtgggTGTAAAAAAA‑GAAGA‑AGCG‑GGGTTCCCGGGCCGCGGGAGGG‑TCTTTTTCCCCGCGTGAAAGCCCCGCTCCAGGGGGTGTCTTCGGG‑GGGGGGAA‑ATTACTTCTTTT‑TT‑C‑CAA‑AC‑AATATTCCCGGCCTGGG‑GC‑CCGCTTA‑ACACC‑TTTTCCC‑CGGAGAAAC‑GCCCCACACGGACCC‑CCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTa                                                                                                                                                                     <  1:679944/225‑1 (MQ=255)                                                                       ctcccagagtgttaGGGTGTAAAAAAACGACGC‑ACCGGGGCTCCCGGCGCGGC‑‑GGGGG‑TCTTTTGAGCCGCGCTTAATCCCCCTTCCCGGGGGTGGTTTTGC‑‑GGCGGGGGTATCATTCCTTTT‑TT‑C‑CCA‑AC‑ACAATTTCCTGCCTGCG‑GC‑CCTCTTA‑CCAAC‑TGTTCGC‑CGGAAAACC‑CCCCCAGACTGACCC‑CCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTa                                                                                                                                                                     <  1:391815/226‑1 (MQ=255)                                                                             gagCTTCAGGGTGTAAATAAA‑GACGC‑AACGTGGCTTCCAGCGCCGC‑‑GAGGG‑TCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCA‑GGCGGGCG‑ATCCTCTCTTTT‑TT‑C‑CCA‑AC‑AATAATTCCTGCCTGCG‑GC‑CCGCCTC‑ACAAC‑TTTTCGC‑CGGGGAAAC‑TACCCCTACTGACGC‑CCCCCCCAATAAACTCCCCCGGGATTTGCCCTTTCGGCAAAccacaacggcca                                                                                                                                                              >  1:711341/1‑229 (MQ=255)                                                                                cTTCAGGGTGTAAATAAA‑GACGC‑AACGTGGCTTCCAGCGCCGC‑‑GAGGG‑TCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCA‑GGCGGGCG‑ATCATCTCTTTT‑TT‑A‑CCA‑AC‑AATAATTCCTGCCTGCG‑GC‑CCGCCTA‑ACAAC‑TTGTCGC‑CGGAGAAAC‑TCAcc                                                                                                                                                                                                                              <  2:118193/173‑1 (MQ=255)                                                                                        tgtAAATAAA‑GACGC‑AACGTGGCTTCCAGCGCCGC‑‑GAGGG‑TCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCA‑GGCGGGCG‑ATCATCTCTTTT‑TT‑T‑CCC‑AC‑AAAAATTCCTTCCTGCG‑GC‑CCGCCCA‑AACAC‑TTGTCGC‑CGGAGAAAC‑TCACCAGACGTACGC‑CCGCCGCACTCAACTCCTGCCGCCTTGGACCTTTCGGTCAAccacaccagaccaaaaacaaacg                                                                                                                                                   >  1:657945/1‑218 (MQ=255)                                                                                        tgtAAATAAA‑GACGC‑AACGTGGCTTCCAGCGCCGC‑‑GAGGG‑TCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCA‑GGCGGGCG‑ATCATCTCTTTT‑TT‑T‑CCA‑AC‑AATAATTCCTGCCTGCG‑GC‑CCGCCTA‑ACAAC‑TTGTCGC‑CGGAGAAAA‑TCACCAGACTGACGC‑CCGCCGCCATCAACCCCCGCGGCATTGGCTCTTTCGACAAACCGAACTGGCTaaaa                                                                                                                                                           >  2:305657/1‑229 (MQ=255)                                                                                         ttaaaaaaa‑GACCC‑AAGGTGGGCTCCGGAGCGGG‑‑GGGGG‑TTCTTTTCTCCCCGCGTAATGGACCGTCTAGCGGGGGGCTTTGGAGGGCGGGC‑‑AACATTTCTTTT‑TTCA‑CAA‑AA‑AATAATTCCCGCCTCGG‑CC‑CCCCCTA‑AAACC‑TTGCCGC‑CGGAAAACC‑CCCCCAGACTGCCCC‑CCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTGAGATCGACCAg                                                                                                                                                   <  1:654286/238‑1 (MQ=255)                                                                                                 aacggcccatcGTTGCCTCCGGCGCCCG‑‑GGGGG‑TTTTTTTACCCCAGCGTAGTGCACGTTTGGGGGGGTGTTTTTGCG‑GGCGGGAT‑A‑CCTCTTTTTTGTC‑A‑CCC‑AA‑AAAAATGCCTCCCGGCG‑CC‑CCTCCTA‑ACAAC‑TTGTCGC‑CGGAACAAC‑TCACCAGCCTGACCC‑CCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCa                                                                                                                                            <  1:714484/229‑1 (MQ=255)                                                                                                     cGC‑AACGTGGCTTCCAGCGCCGC‑‑GAGGG‑TCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCA‑GGCGGGCG‑ATCATCTCTTTT‑TT‑A‑CCA‑AC‑AATAATTCCTGCCTGCG‑GC‑CCGCCTA‑ACAAC‑TTTTCGC‑CGGGGAAAC‑CCACCCGACTGACCC‑CCCCCGCAATCAAATCCTTCGGCATTGGCCCTTTCGGCAAACCGTACCGGCCAAGaactgtcccttttacccaat                                                                                                                                        >  2:672882/1‑219 (MQ=255)                                                                                                                            gcgggtggttttttttCTCCGCGGATATTCCCCGCTTCGGGGGGTTGCTTGGGG‑GGCGGCCA‑TACCTTTTTTTT‑TT‑A‑CCA‑AC‑AAAAACTCCCGCCGGGG‑CC‑CCCCCTA‑AAACC‑TTGCCGC‑GGGAAAAAC‑CCACCAGCCGGCCCC‑CCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGa                                                                                                                    <  1:289769/227‑1 (MQ=255)                                                                                                                            ccggggggtttttttttcccGCGCTTAATCCACGCTGCAGCGGGGGGTGTTGCA‑GGGGGGCG‑ATCCTCTTTTTT‑CT‑A‑CCA‑AA‑AATAATCCCTGCGTGGG‑GC‑CCGCCTA‑ACAAC‑TTGTCGC‑GGGAAAAAC‑TCCCAAGCCGAACGC‑CCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGa                                                                                                                    <  1:719915/221‑1 (MQ=255)                                                                                                                                 ggg‑TCTTTTTACCCGCGCTTAATGCACGCTTCGGGGGGTGGCTTTGCA‑GGCGGGCG‑ATCATCTCTTTT‑TT‑A‑CCA‑AC‑AATAATCCCTGCCTGCG‑GC‑CCGCCTA‑ACAAC‑TTGTCGC‑CGGAAAAAC‑TCCCCAGACTGACGC‑CCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATc                                                                                                               <  2:513707/238‑1 (MQ=255)                                                                                                                                                                        tttcttggga‑GGCGGCCG‑TTCATCTTTTTT‑TC‑C‑CCA‑AA‑AATACTTCCGCCCTGCG‑CC‑CCGCTTA‑ACACC‑TTCTCGC‑CGGAAAAAC‑CCACCAGACTGACGC‑CCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTa                                                                         <  2:200311/229‑1 (MQ=255)                                                                                                                                                                         ggCTTTGCA‑GGCGGGCG‑ATCATCTCTTTT‑TT‑A‑CCA‑AC‑AATAATTCCTGCCTGCG‑GC‑CCGCCTA‑ACAAC‑TTGTCGC‑CGGAGAAAC‑TCACCAGACTGACGC‑CCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCa                                                                                                                                                                            <  2:567603/138‑1 (MQ=255)                                                                                                                                                                         ggCTTTGCA‑GGCGGGCG‑ATCATCTCTTTT‑TT‑A‑CCA‑AC‑AATAATTCCTGCCTGCG‑GC‑CCGCCTA‑ACAAC‑TTGTCGC‑CGGAGAAAC‑TCACCAGACTGACGC‑CCGCCGCAATCAACTCCTGCGGCATTGGCCCTTTCGGCa                                                                                                                                                                            >  1:567603/1‑138 (MQ=255)                                                                                                                                                                                           tatcttcTTTTTT‑CT‑A‑CCA‑AA‑AATACCCCCCGCCGGCC‑CC‑CCCCCAC‑ACACC‑TTGCCGC‑GGGAAAACC‑TCCCCAACCGGACCC‑CCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAACCCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCt                                                       <  1:641800/238‑1 (MQ=255)                                                                                                                                                                                            aacattttcttt‑tt‑t‑CCA‑ACAAAAAACTCCTGCCGGCG‑GC‑CCCCCTA‑ACAAC‑TTGCCCC‑CGGAAAACC‑CCCCAAGACTGCCCC‑CCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCt                                                       <  2:877231/230‑1 (MQ=255)                                                                                                                                                                                             tacttcctTTT‑TT‑C‑CAA‑AA‑AAAAAATCCCGCC‑GCGCGC‑CCCCCTA‑AAAACTTTGCCGC‑GGGAAAACC‑CCCCAAGACTGCCGC‑CCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCt                                                       <  2:874836/232‑1 (MQ=255)                                                                                                                                                                                              cattcctTTT‑TT‑A‑CCA‑AC‑AATAATTCCTGCCTGCG‑GC‑CCGCTTA‑ACAAC‑TTGCCGC‑CGGAAAAAC‑CCCCCAGCCTGCCGC‑CCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGaa                                                    <  2:639820/234‑1 (MQ=255)                                                                                                                                                                                               cattttttt‑tt‑ACCCA‑AC‑AATAATGCCCGCCGGGG‑GC‑CCCCCTA‑ACAAC‑TTGCCGC‑CGGAAAAAC‑CCCCCAGGCTGACGC‑CCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGaa                                                    <  2:443414/234‑1 (MQ=255)                                                                                                                                                                                                  tttttt‑tt‑C‑CCA‑AA‑AAAAATCCTTGCCTGCG‑GC‑CCGCCAA‑CCAAC‑TGGCCGC‑GGGAGAACC‑CCACCAGCCTGCCGC‑CCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCt                                                <  2:33428/237‑1 (MQ=255)                                                                                                                                                                                                     ttttttgcacaa‑ac‑aaTAATCCCCGCCGTCG‑CC‑CCGCCTA‑ACCAC‑TCGCCGC‑CGGAGAAAC‑CCACCAGACTGACGC‑CCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCt                                                <  2:101221/228‑1 (MQ=255)                                                                                                                                                                                                      tt‑ta‑A‑CCA‑AC‑AATAATCCCGGCCTGCG‑GC‑CCGCCTA‑ACAAC‑TTGTCGC‑CGGAGAACC‑CCACCAGACTGACGC‑CCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAg                                            <  2:242498/239‑1 (MQ=255)                                                                                                                                                                                                          t‑A‑CCA‑AC‑AATAATTCCTGCCTGCG‑GC‑CCGCCTA‑ACAAC‑TTGTCGC‑CGGAGAAAC‑TCACCAGACTGACGC‑CCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGtt                                          >  2:181573/1‑238 (MQ=255)                                                                                                                                                                                                            acccA‑AA‑AAAAACTCCCGCC‑GCG‑GCGCCCCCAA‑ACAAC‑TGGCCGC‑CGAAAAACC‑CCACAGTACTGCCCC‑CCGCCGCAATCAACTCCTGCGGCATTGGATCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAGCCCCTGAATGCTGTAGTTa                                         <  2:813069/237‑1 (MQ=255)                                                                                                                                                                                                              accaaaaaataacTCCTGCTTGCG‑GC‑CCGCACA‑ACATC‑TTG‑CGCGCGGACAAAC‑TCCCCAGACTGCCGC‑CCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTa                                         <  2:498117/232‑1 (MQ=255)                                                                                                                                                                                                                  ac‑aaTAATTCCTGCCTGCG‑GC‑CCGCCTA‑ACAAC‑TTGTCGC‑CGGAGAAAC‑TCACCAGCCTGACGC‑CCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTa                                                                                                                                                                     <  2:405237/110‑1 (MQ=255)                                                                                                                                                                                                                  ac‑aaTAATTCCTGCCTGCG‑GC‑CCGCCTA‑ACAAC‑TTGTCGC‑CGGAGAAAC‑TCACCAGACTGACGC‑CCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTa                                                                                                                                                                     >  1:405237/1‑110 (MQ=255)                                                                                                                                                                                                                  ac‑aaTAATTCCTGCCTGCG‑GC‑CCGCCTA‑ACAAC‑TTGTCGC‑CGGAGAAAC‑TCACCAGACTGACGC‑CCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCCCTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCCAGCGCCCCCAGTTCCGCTTCCTCTTTCGCTTTGGTTAACCCCCAAATGCTTTTGTTTCGGGt                                    >  1:448582/1‑239 (MQ=255)                                                                                                                                                                                                                     aataatTCCTGCCTGCG‑GC‑CCGCCTA‑ACAAC‑TTGTCGC‑CGGAGAAAC‑TCACCAGACTGACGC‑CCGCCGCAATCACCTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGATCGCCACCAGTTCCGCTTCATCTATCGCCTTGGTGAACCCCTGAATGCTGTAGTTACTGGAAt                                  >  2:61396/1‑239 (MQ=255)                                                                                                                                                                                                                        cttcccgccctgcG‑CC‑CCCCCAA‑AAATT‑TTGCCCC‑GGGAAAACC‑TCACCGAACTGCCGC‑CCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTa                               <  1:204300/231‑1 (MQ=255)                                                                                                                                                                                                                        atccccggcctgcG‑GC‑CCCCCAA‑CAAAC‑TGGTCGC‑CGGAGACAC‑CCCCCAGACTGCCGC‑CCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTa                               <  2:711341/232‑1 (MQ=255)                                                                                                                                                                                                                         attcccgccctggcgC‑CCCCCTA‑ACAAC‑TTGTCGC‑CGGAAAAAC‑CCACCAGCCTGCCCC‑CCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTa                               <  2:206754/225‑1 (MQ=255)                                                                                                                                                                                                                          tcccggcctccG‑CC‑CCCCATA‑ACACC‑TTGCCGC‑CGGGAAAAC‑CCCCAAGACGGCCGC‑CCTCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACt                             <  2:830425/234‑1 (MQ=255)                                                                                                                                                                                                                          cccctgccggCG‑CC‑CCGCCTA‑ACAAC‑TTTCCGC‑TGGAAAACC‑CCCCAAGACTGCCCC‑CCGCCGCAATCAACCCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACt                             <  2:764682/237‑1 (MQ=255)                                                                                                                                                                                                                               gcctgcG‑GC‑CCGCCTA‑ACAAC‑TTGTCGC‑CGGAGAAAC‑TCACCAGACTGACGC‑CCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATTGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGGTACTGGTATGTACTTtca                         >  2:661463/1‑238 (MQ=255)                                                                                                                                                                                                                               gcctgcG‑GC‑CCGCCTA‑ACAAC‑TTGTCGC‑CGGAGAAAC‑TCACCAGACTGACGC‑CCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCCCCAGTTCCGCTTCTTCTACCGCTTTGGGTAACCCCTGAAAGTGGTAGTTACCGGTAAGAACTTtaat                        >  1:455828/1‑239 (MQ=255)                                                                                                                                                                                                                               gcctgcG‑GC‑CCGCCTA‑ACAAC‑TTGTCGC‑CGGAGAAAC‑TCACCAGACTGACGC‑CCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCCCTTTGCCGAGCGCCACCAGTTCCGCTTCCTCTATCGCTTTGGTAAACCCCTGAATGCTGTAGTTACTGGGATGAACTTtcat                        >  1:453033/1‑239 (MQ=255)                                                                                                                                                                                                                                  tgcG‑GC‑CCGCCTA‑ACAAC‑TTGTCGC‑CGGAGAAAC‑TCACCAGACTGACGC‑CCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCGACGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTTCCGAGGGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAAAGCTGCAGTTACTGGTATGTACTTtcatac                      >  1:586827/1‑236 (MQ=255)                                                                                                                                                                                                                                  tgcG‑GC‑CCGCCTA‑ACAAC‑TTGTCGC‑CGGAGAAAC‑TCACCAGACTGACGC‑CCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTtcatca                      >  2:204300/1‑238 (MQ=255)                                                                                                                                                                                                                                  tgcG‑GC‑CCGCCTA‑ACAAC‑TTGTCGC‑CGGAGAAAC‑TCACCAGACTGACGC‑CCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCaa                     >  2:77887/1‑239 (MQ=255)                                                                                                                                                                                                                                    gggcc‑ccGCCAA‑AC‑AC‑TTGGTCC‑CGGGGAAACTTCACC‑GACTGGCGCGCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGAGCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCaa                     <  2:456413/234‑1 (MQ=255)                                                                                                                                                                                                                                                         tGTCGC‑CGGAGAAAC‑TCACCAGACTGACGC‑CCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGGGGGATTttcatt   >  1:811380/1‑238 (MQ=255)                                                                                                                                                                                                                                                         tGTCGC‑CGGAGAAAC‑TCACCAGACTGACGC‑CCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGAATGTACTTTCATCAACAGTGCGGTATTttcattc  >  2:789713/1‑239 (MQ=255)                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                 TGTAAACGTTGTGCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAA‑GACGC‑AACGTGGCTTCCAGCGCCGC‑‑GAGGG‑TCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCA‑GGCGGGCG‑ATCATCTCTTTT‑TT‑A‑CCA‑AC‑AATAATTCCTGCCTGCG‑GC‑CCGCCTA‑ACAAC‑TTGTCGC‑CGGAGAAAC‑TCACCAGACTGACGC‑CCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTTTCATTC  >  NC_000913/3760134‑3760604

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