Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,760,381	T→C	32.2%	S290G (AGT→GGT)	selA ←	selenocysteine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,760,381	0	T	C	32.2%	73.0 / 8.7	60	S290G (AGT→GGT)	selA	selenocysteine synthase
Reads supporting (aligned to +/- strand):  ref base T (25/13);  new base C (2/18);  total (28/32)
Fisher's exact test for biased strand distribution p-value = 6.30e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAA‑GACGC‑AACGTGGCTTCCAGCGCCGC‑‑GAGGG‑TCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCA‑GGCGGGCG‑ATCATCTCTTTT‑TT‑A‑CCAAC‑AATAATTCCTGCCTGCG‑GC‑CCGCCTAACAAC‑TTGTCGC‑CGGAGAAAC‑TCACCAGACTGACGC‑CCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTTTCATTCACCGC  >  NC_000913/3760144‑3760609                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                      gtgCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAA‑GACGC‑AACGTGGCTTCCAGCGCCGC‑‑GAGGG‑TCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGGGGCTTTGCA‑GGCGGGCG‑ATCATCTCTTTT‑TT‑A‑CCAAC‑AATAATTCCTGCCTGCG‑GC‑CCGCCTTACAAC‑TTTTCGC‑CGGGGAAAC‑Tc                                                                                                                                                                                                                                      >  2:289769/1‑239 (MQ=255)  ttccTGGATTTGAATGCCCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAA‑GACGC‑AACGTGGCTTCCAGCGCCGC‑‑GAGGG‑TCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCA‑GGCGGGCG‑ATCATCTCTTTT‑TT‑A‑CCAAC‑AATAATTCCTGCCTGCG‑GC‑CCGCCTAACAAC‑TTGTCGC‑CGGAGAAAC‑TCa                                                                                                                                                                                                                                     <  2:22498/237‑1 (MQ=255)   gCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCCGGGTGTAAATAAA‑GACGC‑AACGTGGCTTCCAGCGCCGC‑‑GAGGG‑TCATTTTATCCGCGCGTAATGCACGCTTCAGCGGCTGGCTTTGCC‑GGCGGGCG‑ATCATCTCTTTT‑TT‑A‑CCAAC‑AATAATTCCTGCCTCCG‑GC‑ACGCCACACACC‑TTTTCCC‑CCGGGAATC‑TCa                                                                                                                                                                                                                                     >  2:252867/1‑238 (MQ=255)   gCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAA‑GACGC‑AACGTGGCTTCCAGCGCCGC‑‑GAGGG‑TCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCA‑GGCGGGCG‑ATCATCTCTTTT‑TT‑A‑CCAAC‑AATAATTCCTGCCTGCG‑GC‑CCGCCTAACAAC‑TTGTGGC‑CGGGGAAAA‑TCa                                                                                                                                                                                                                                     >  2:679944/1‑238 (MQ=255)        gATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAA‑GACGC‑AACGTGGCTTCCAGCGCCGC‑‑GAGGG‑TCATTTTATCCGCGCGGAATGCACGCTTCAGCGGGTGGCTTTGCA‑GGCGGGCG‑AGCAACTCTTTT‑TT‑T‑ACAAA‑AACAATTCCTGCCTGCG‑GC‑CCCCCCAAAAAA‑CTTTTGC‑CGGGGAAAA‑ACACCacaa                                                                                                                                                                                                                               >  1:200311/1‑236 (MQ=255)                              cgGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAA‑GACGC‑AACGTGGCTTCCAGCGCCGC‑‑GAGGG‑TCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCA‑GGCGGGCG‑GTCATCCCTTTT‑TT‑A‑CCAAC‑AAAAATTCCTGCCTTCG‑GC‑CCGCCTAAACAC‑TTGTCCC‑CGGAGAAAC‑TCACCAACCATACGC‑CCGCCCCaactaaac                                                                                                                                                                                                         >  2:263677/1‑232 (MQ=255)                              cgGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAA‑GACGC‑AACGTGGCTTCCAGCGCCGC‑‑GAGGG‑TCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCA‑GGCGGGCG‑ATCATCTCTTTT‑TT‑A‑CCAAC‑AATAATTCCTGCCTGCG‑GC‑CCGCCTAACAAC‑TTTTCGC‑CGGAGAAAC‑TCAACAGACTGACGC‑CCGCCGCACTCAAcc                                                                                                                                                                                                         >  2:363342/1‑237 (MQ=255)                              cgGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAA‑GACGC‑AACGTGGCTTCCAGCGCCGC‑‑GAGGG‑TCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCA‑GGCGGGCG‑ATACTCTCTTTT‑TT‑A‑CCAAC‑AAAAATTCCTGCCTGCG‑GC‑CCGCCTAACAAC‑CTGTCCC‑CGCAGAAAA‑CCACCCAACTTAccc‑ccccccacaacacac                                                                                                                                                                                                         >  1:874836/1‑228 (MQ=255)                                          cgcAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAA‑GACGC‑AACGTGGCTTCCAGCGCCGC‑‑GAGGG‑TCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCA‑GGCGGGCG‑ATCATCTCTTTT‑TT‑A‑CCAAC‑AATAATTCCTGCCTCCG‑GC‑CCGCCTAACCAA‑TTTTCGC‑CCGAGAAAA‑TCACCCAACTGACGC‑CCCCCCCAAACAACCCCTGCGGGACTgg                                                                                                                                                                                            >  1:877231/1‑239 (MQ=255)                                                           ttCACTCAGAGCTTCAGGGTGTAAATAAA‑GACGC‑AACGTGGCTTCCAGCGCCGC‑‑GAGGG‑TCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCA‑GGCGGGCG‑ATCATCTCTTTT‑TT‑A‑CCAAC‑AATAATTCCTGCCTGCG‑GC‑CCGCCTAACAAC‑TTGTCGC‑CGGAGAAAC‑TCACCAGACTGACGC‑CCGCCGCAATCAACTCCTGCGGCATTTGCTCTTTCGGCCAACCGt                                                                                                                                                                           >  2:745271/1‑239 (MQ=255)                                                           ttCACTCAGAGCTTCAGGGTGTAAATAAA‑GACGC‑AACGTGGCTTCCAGCGCCGC‑‑GAGGG‑TCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCA‑GGCGGGCG‑ATCATCTCTTTT‑TT‑A‑CCAAC‑AATAATTCCTGCCTGCG‑GC‑CCGCCTAACAAC‑TTGTCGC‑CGGAGAAAC‑TCACCAGACTGACGC‑CCGCCGCAATCAACCCCCGCGTCATTGGCTCTTTCGGCaacacg                                                                                                                                                                            >  2:458107/1‑234 (MQ=255)                                                            tcccacagggtgtgaGGGTGTAAATAAA‑GACGG‑CACGGGGCTGCCGTGGCGGG‑‑GGGGG‑GCATTTTCTCGGCGCTTAGTCACCTCTCCGGCGGGTGGCTTTGCG‑GGGGTGCG‑ACTATCTTTTTT‑TT‑C‑CCAAC‑AATAATCCCTGCCTGCG‑AC‑CCGCCTAACACT‑TTGTCGC‑CGGAGAACC‑TCACCAGACAGACGC‑CCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTa                                                                                                                                                                          <  1:661577/225‑1 (MQ=255)                                                            cccccaggagCTTGGGGGTG‑AAAAAAA‑GAGCCAAACGTGGCTCCCA‑CGCCCC‑‑GGGGGTTTTTTTTATCGGCGGTTAAGGCCCCCTTCAGGGGGTGGCTTGGCG‑GGCGGGCT‑TTCATCTTTTTT‑CT‑A‑CCAAC‑AATAATCCCTGCCTGCG‑GC‑CCGCCTAACAAC‑TTGTCGC‑GGGAGAAAC‑TCACCAGCCTGACGC‑CCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTa                                                                                                                                                                          <  1:639443/232‑1 (MQ=255)                                                            agttgtgtataggagaaagtaaaaaaaa‑GCCGC‑AGCGTGGCTCCCAGCGCGGG‑‑GAGGT‑TCATTGTATCCGTGCTTAATGCACGTTTCGGCGGGTGCTTTGGCG‑GGCGGGCG‑ATCATCTTTTTT‑TT‑A‑CCAAA‑AATAATCCCTGCCGGCG‑GC‑CCGCCTAACAAC‑TTGTCGC‑CGGAGAAAC‑CCCCCAGCCTGACGC‑CCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTa                                                                                                                                                                          <  1:672882/218‑1 (MQ=255)                                                             tcctcagggtttgtgggTGTAAAAAAA‑GAAGA‑AGCG‑GGGTTCCCGGGCCGCGGGAGGG‑TCTTTTTCCCCGCGTGAAAGCCCCGCTCCAGGGGGTGTCTTCGGG‑GGGGGGAA‑ATTACTTCTTTT‑TT‑C‑CAAAC‑AATATTCCCGGCCTGGG‑GC‑CCGCTTAACACC‑TTTTCCC‑CGGAGAAAC‑GCCCCACACGGACCC‑CCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTa                                                                                                                                                                          <  1:679944/225‑1 (MQ=255)                                                             ctcccagagtgttaGGGTGTAAAAAAACGACGC‑ACCGGGGCTCCCGGCGCGGC‑‑GGGGG‑TCTTTTGAGCCGCGCTTAATCCCCCTTCCCGGGGGTGGTTTTGC‑‑GGCGGGGGTATCATTCCTTTT‑TT‑C‑CCAAC‑ACAATTTCCTGCCTGCG‑GC‑CCTCTTACCAAC‑TGTTCGC‑CGGAAAACC‑CCCCCAGACTGACCC‑CCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTa                                                                                                                                                                          <  1:391815/226‑1 (MQ=255)                                                                   gagCTTCAGGGTGTAAATAAA‑GACGC‑AACGTGGCTTCCAGCGCCGC‑‑GAGGG‑TCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCA‑GGCGGGCG‑ATCCTCTCTTTT‑TT‑C‑CCAAC‑AATAATTCCTGCCTGCG‑GC‑CCGCCTCACAAC‑TTTTCGC‑CGGGGAAAC‑TACCCCTACTGACGC‑CCCCCCCAATAAACTCCCCCGGGATTTGCCCTTTCGGCAAAccacaacggcca                                                                                                                                                                   >  1:711341/1‑229 (MQ=255)                                                                      cTTCAGGGTGTAAATAAA‑GACGC‑AACGTGGCTTCCAGCGCCGC‑‑GAGGG‑TCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCA‑GGCGGGCG‑ATCATCTCTTTT‑TT‑A‑CCAAC‑AATAATTCCTGCCTGCG‑GC‑CCGCCTAACAAC‑TTGTCGC‑CGGAGAAAC‑TCAcc                                                                                                                                                                                                                                   <  2:118193/173‑1 (MQ=255)                                                                              tgtAAATAAA‑GACGC‑AACGTGGCTTCCAGCGCCGC‑‑GAGGG‑TCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCA‑GGCGGGCG‑ATCATCTCTTTT‑TT‑T‑CCCAC‑AAAAATTCCTTCCTGCG‑GC‑CCGCCCAAACAC‑TTGTCGC‑CGGAGAAAC‑TCACCAGACGTACGC‑CCGCCGCACTCAACTCCTGCCGCCTTGGACCTTTCGGTCAAccacaccagaccaaaaacaaacg                                                                                                                                                        >  1:657945/1‑218 (MQ=255)                                                                              tgtAAATAAA‑GACGC‑AACGTGGCTTCCAGCGCCGC‑‑GAGGG‑TCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCA‑GGCGGGCG‑ATCATCTCTTTT‑TT‑T‑CCAAC‑AATAATTCCTGCCTGCG‑GC‑CCGCCTAACAAC‑TTGTCGC‑CGGAGAAAA‑TCACCAGACTGACGC‑CCGCCGCCATCAACCCCCGCGGCATTGGCTCTTTCGACAAACCGAACTGGCTaaaa                                                                                                                                                                >  2:305657/1‑229 (MQ=255)                                                                               ttaaaaaaa‑GACCC‑AAGGTGGGCTCCGGAGCGGG‑‑GGGGG‑TTCTTTTCTCCCCGCGTAATGGACCGTCTAGCGGGGGGCTTTGGAGGGCGGGC‑‑AACATTTCTTTT‑TTCA‑CAAAA‑AATAATTCCCGCCTCGG‑CC‑CCCCCTAAAACC‑TTGCCGC‑CGGAAAACC‑CCCCCAGACTGCCCC‑CCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTGAGATCGACCAg                                                                                                                                                        <  1:654286/238‑1 (MQ=255)                                                                                       aacggcccatcGTTGCCTCCGGCGCCCG‑‑GGGGG‑TTTTTTTACCCCAGCGTAGTGCACGTTTGGGGGGGTGTTTTTGCG‑GGCGGGAT‑A‑CCTCTTTTTTGTC‑A‑CCCAA‑AAAAATGCCTCCCGGCG‑CC‑CCTCCTAACAAC‑TTGTCGC‑CGGAACAAC‑TCACCAGCCTGACCC‑CCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCa                                                                                                                                                 <  1:714484/229‑1 (MQ=255)                                                                                           cGC‑AACGTGGCTTCCAGCGCCGC‑‑GAGGG‑TCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCA‑GGCGGGCG‑ATCATCTCTTTT‑TT‑A‑CCAAC‑AATAATTCCTGCCTGCG‑GC‑CCGCCTAACAAC‑TTTTCGC‑CGGGGAAAC‑CCACCCGACTGACCC‑CCCCCGCAATCAAATCCTTCGGCATTGGCCCTTTCGGCAAACCGTACCGGCCAAGaactgtcccttttacccaat                                                                                                                                             >  2:672882/1‑219 (MQ=255)                                                                                                                  gcgggtggttttttttCTCCGCGGATATTCCCCGCTTCGGGGGGTTGCTTGGGG‑GGCGGCCA‑TACCTTTTTTTT‑TT‑A‑CCAAC‑AAAAACTCCCGCCGGGG‑CC‑CCCCCTAAAACC‑TTGCCGC‑GGGAAAAAC‑CCACCAGCCGGCCCC‑CCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGa                                                                                                                         <  1:289769/227‑1 (MQ=255)                                                                                                                  ccggggggtttttttttcccGCGCTTAATCCACGCTGCAGCGGGGGGTGTTGCA‑GGGGGGCG‑ATCCTCTTTTTT‑CT‑A‑CCAAA‑AATAATCCCTGCGTGGG‑GC‑CCGCCTAACAAC‑TTGTCGC‑GGGAAAAAC‑TCCCAAGCCGAACGC‑CCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGa                                                                                                                         <  1:719915/221‑1 (MQ=255)                                                                                                                       ggg‑TCTTTTTACCCGCGCTTAATGCACGCTTCGGGGGGTGGCTTTGCA‑GGCGGGCG‑ATCATCTCTTTT‑TT‑A‑CCAAC‑AATAATCCCTGCCTGCG‑GC‑CCGCCTAACAAC‑TTGTCGC‑CGGAAAAAC‑TCCCCAGACTGACGC‑CCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATc                                                                                                                    <  2:513707/238‑1 (MQ=255)                                                                                                                                                              tttcttggga‑GGCGGCCG‑TTCATCTTTTTT‑TC‑C‑CCAAA‑AATACTTCCGCCCTGCG‑CC‑CCGCTTAACACC‑TTCTCGC‑CGGAAAAAC‑CCACCAGACTGACGC‑CCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTa                                                                              <  2:200311/229‑1 (MQ=255)                                                                                                                                                               ggCTTTGCA‑GGCGGGCG‑ATCATCTCTTTT‑TT‑A‑CCAAC‑AATAATTCCTGCCTGCG‑GC‑CCGCCTAACAAC‑TTGTCGC‑CGGAGAAAC‑TCACCAGACTGACGC‑CCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCa                                                                                                                                                                                 <  2:567603/138‑1 (MQ=255)                                                                                                                                                               ggCTTTGCA‑GGCGGGCG‑ATCATCTCTTTT‑TT‑A‑CCAAC‑AATAATTCCTGCCTGCG‑GC‑CCGCCTAACAAC‑TTGTCGC‑CGGAGAAAC‑TCACCAGACTGACGC‑CCGCCGCAATCAACTCCTGCGGCATTGGCCCTTTCGGCa                                                                                                                                                                                 >  1:567603/1‑138 (MQ=255)                                                                                                                                                                                 tatcttcTTTTTT‑CT‑A‑CCAAA‑AATACCCCCCGCCGGCC‑CC‑CCCCCACACACC‑TTGCCGC‑GGGAAAACC‑TCCCCAACCGGACCC‑CCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAACCCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCt                                                            <  1:641800/238‑1 (MQ=255)                                                                                                                                                                                  aacattttcttt‑tt‑t‑CCAACAAAAAACTCCTGCCGGCG‑GC‑CCCCCTAACAAC‑TTGCCCC‑CGGAAAACC‑CCCCAAGACTGCCCC‑CCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCt                                                            <  2:877231/230‑1 (MQ=255)                                                                                                                                                                                   tacttcctTTT‑TT‑C‑CAAAA‑AAAAAATCCCGCC‑GCGCGC‑CCCCCTAAAAACTTTGCCGC‑GGGAAAACC‑CCCCAAGACTGCCGC‑CCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCt                                                            <  2:874836/232‑1 (MQ=255)                                                                                                                                                                                    cattcctTTT‑TT‑A‑CCAAC‑AATAATTCCTGCCTGCG‑GC‑CCGCTTAACAAC‑TTGCCGC‑CGGAAAAAC‑CCCCCAGCCTGCCGC‑CCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGaa                                                         <  2:639820/234‑1 (MQ=255)                                                                                                                                                                                     cattttttt‑tt‑ACCCAAC‑AATAATGCCCGCCGGGG‑GC‑CCCCCTAACAAC‑TTGCCGC‑CGGAAAAAC‑CCCCCAGGCTGACGC‑CCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGaa                                                         <  2:443414/234‑1 (MQ=255)                                                                                                                                                                                        tttttt‑tt‑C‑CCAAA‑AAAAATCCTTGCCTGCG‑GC‑CCGCCAACCAAC‑TGGCCGC‑GGGAGAACC‑CCACCAGCCTGCCGC‑CCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCt                                                     <  2:33428/237‑1 (MQ=255)                                                                                                                                                                                           ttttttgcacaaac‑aaTAATCCCCGCCGTCG‑CC‑CCGCCTAACCAC‑TCGCCGC‑CGGAGAAAC‑CCACCAGACTGACGC‑CCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCt                                                     <  2:101221/228‑1 (MQ=255)                                                                                                                                                                                            tt‑ta‑A‑CCAAC‑AATAATCCCGGCCTGCG‑GC‑CCGCCTAACAAC‑TTGTCGC‑CGGAGAACC‑CCACCAGACTGACGC‑CCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAg                                                 <  2:242498/239‑1 (MQ=255)                                                                                                                                                                                                t‑A‑CCAAC‑AATAATTCCTGCCTGCG‑GC‑CCGCCTAACAAC‑TTGTCGC‑CGGAGAAAC‑TCACCAGACTGACGC‑CCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGtt                                               >  2:181573/1‑238 (MQ=255)                                                                                                                                                                                                  acccAAA‑AAAAACTCCCGCC‑GCG‑GCGCCCCCAAACAAC‑TGGCCGC‑CGAAAAACC‑CCACAGTACTGCCCC‑CCGCCGCAATCAACTCCTGCGGCATTGGATCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAGCCCCTGAATGCTGTAGTTa                                              <  2:813069/237‑1 (MQ=255)                                                                                                                                                                                                   accaaaaaataacTCCTGCTTGCG‑GC‑CCGCACAACATC‑TTG‑CGCGCGGACAAAC‑TCCCCAGACTGCCGC‑CCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTa                                              <  2:498117/232‑1 (MQ=255)                                                                                                                                                                                                       ac‑aaTAATTCCTGCCTGCG‑GC‑CCGCCTAACAAC‑TTGTCGC‑CGGAGAAAC‑TCACCAGCCTGACGC‑CCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTa                                                                                                                                                                          <  2:405237/110‑1 (MQ=255)                                                                                                                                                                                                       ac‑aaTAATTCCTGCCTGCG‑GC‑CCGCCTAACAAC‑TTGTCGC‑CGGAGAAAC‑TCACCAGACTGACGC‑CCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTa                                                                                                                                                                          >  1:405237/1‑110 (MQ=255)                                                                                                                                                                                                       ac‑aaTAATTCCTGCCTGCG‑GC‑CCGCCTAACAAC‑TTGTCGC‑CGGAGAAAC‑TCACCAGACTGACGC‑CCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCCCTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCCAGCGCCCCCAGTTCCGCTTCCTCTTTCGCTTTGGTTAACCCCCAAATGCTTTTGTTTCGGGt                                         >  1:448582/1‑239 (MQ=255)                                                                                                                                                                                                          aataatTCCTGCCTGCG‑GC‑CCGCCTAACAAC‑TTGTCGC‑CGGAGAAAC‑TCACCAGACTGACGC‑CCGCCGCAATCACCTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGATCGCCACCAGTTCCGCTTCATCTATCGCCTTGGTGAACCCCTGAATGCTGTAGTTACTGGAAt                                       >  2:61396/1‑239 (MQ=255)                                                                                                                                                                                                             cttcccgccctgcG‑CC‑CCCCCAAAAATT‑TTGCCCC‑GGGAAAACC‑TCACCGAACTGCCGC‑CCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTa                                    <  1:204300/231‑1 (MQ=255)                                                                                                                                                                                                             atccccggcctgcG‑GC‑CCCCCAACAAAC‑TGGTCGC‑CGGAGACAC‑CCCCCAGACTGCCGC‑CCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTa                                    <  2:711341/232‑1 (MQ=255)                                                                                                                                                                                                              attcccgccctggcgC‑CCCCCTAACAAC‑TTGTCGC‑CGGAAAAAC‑CCACCAGCCTGCCCC‑CCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTa                                    <  2:206754/225‑1 (MQ=255)                                                                                                                                                                                                               tcccggcctccG‑CC‑CCCCATAACACC‑TTGCCGC‑CGGGAAAAC‑CCCCAAGACGGCCGC‑CCTCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACt                                  <  2:830425/234‑1 (MQ=255)                                                                                                                                                                                                               cccctgccggCG‑CC‑CCGCCTAACAAC‑TTTCCGC‑TGGAAAACC‑CCCCAAGACTGCCCC‑CCGCCGCAATCAACCCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACt                                  <  2:764682/237‑1 (MQ=255)                                                                                                                                                                                                                    gcctgcG‑GC‑CCGCCTAACAAC‑TTGTCGC‑CGGAGAAAC‑TCACCAGACTGACGC‑CCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATTGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGGTACTGGTATGTACTTtca                              >  2:661463/1‑238 (MQ=255)                                                                                                                                                                                                                    gcctgcG‑GC‑CCGCCTAACAAC‑TTGTCGC‑CGGAGAAAC‑TCACCAGACTGACGC‑CCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCCCCAGTTCCGCTTCTTCTACCGCTTTGGGTAACCCCTGAAAGTGGTAGTTACCGGTAAGAACTTtaat                             >  1:455828/1‑239 (MQ=255)                                                                                                                                                                                                                    gcctgcG‑GC‑CCGCCTAACAAC‑TTGTCGC‑CGGAGAAAC‑TCACCAGACTGACGC‑CCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCCCTTTGCCGAGCGCCACCAGTTCCGCTTCCTCTATCGCTTTGGTAAACCCCTGAATGCTGTAGTTACTGGGATGAACTTtcat                             >  1:453033/1‑239 (MQ=255)                                                                                                                                                                                                                       tgcG‑GC‑CCGCCTAACAAC‑TTGTCGC‑CGGAGAAAC‑TCACCAGACTGACGC‑CCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCGACGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTTCCGAGGGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAAAGCTGCAGTTACTGGTATGTACTTtcatac                           >  1:586827/1‑236 (MQ=255)                                                                                                                                                                                                                       tgcG‑GC‑CCGCCTAACAAC‑TTGTCGC‑CGGAGAAAC‑TCACCAGACTGACGC‑CCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTtcatca                           >  2:204300/1‑238 (MQ=255)                                                                                                                                                                                                                       tgcG‑GC‑CCGCCTAACAAC‑TTGTCGC‑CGGAGAAAC‑TCACCAGACTGACGC‑CCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCaa                          >  2:77887/1‑239 (MQ=255)                                                                                                                                                                                                                         gggcc‑ccGCCAAAC‑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)                                                                                                                                                                                                                                                    cGGAAAAAC‑CCCCCGGACTCACCC‑CCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTCCTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTTTCATTCACCGc  <  2:619304/238‑1 (MQ=255)                                                                                                                                                                                                                                                       agaAAC‑TCACCAGACTGACGC‑CCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTttcattca      >  2:790170/1‑231 (MQ=255)                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                      GTGCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAA‑GACGC‑AACGTGGCTTCCAGCGCCGC‑‑GAGGG‑TCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCA‑GGCGGGCG‑ATCATCTCTTTT‑TT‑A‑CCAAC‑AATAATTCCTGCCTGCG‑GC‑CCGCCTAACAAC‑TTGTCGC‑CGGAGAAAC‑TCACCAGACTGACGC‑CCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTTTCATTCACCGC  >  NC_000913/3760144‑3760609

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: ‑