Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,760,384	A→C	18.2%	V289G (GTG→GGG)	selA ←	selenocysteine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,760,384	0	A	C	18.2%	145.4 / 3.5	73	V289G (GTG→GGG)	selA	selenocysteine synthase
Reads supporting (aligned to +/- strand):  ref base A (42/16);  new base C (2/13);  total (44/29)
Fisher's exact test for biased strand distribution p-value = 4.65e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

AATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGG‑CTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTA‑ATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCA‑TCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAAC‑TTGTCGCC‑GGAGAAACTCACCAGACTGA‑CGCCCGCCGCAATCAACTCCTGC‑GGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTTTCATT  >  NC_000913/3760158‑3760604                                                                                                                                                                                                                                        |                                                                                                                                                                                                                               aaggACCTCTGCGTAGCGGATAAGCAGGCGCAGGGTCGGTAATTGTTCACTCGGAGCTTCAGGGTGTAAATAAAGACGCAACGTGG‑CTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTA‑ATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCA‑TCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAAC‑TTGTCGCC‑GGAGAAACTCAc                                                                                                                                                                                                                               <  2:1988809/225‑1 (MQ=255)     aCCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGG‑CTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTA‑ATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCA‑TCCCTTTTTTATCCACAATAATTCCTTCCTGCGGCCCGCCTACCAAA‑TTTTCGCC‑GTGGAAAATCACCaa                                                                                                                                                                                                                            >  1:1993293/1‑226 (MQ=255)                 cgGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGG‑CTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTA‑ATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGTCGGTCA‑TCCCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCCCCTAACAAC‑CTTTTGCC‑GGAGAAAAGCACCAAAcctagtccgcc                                                                                                                                                                                                                >  2:1497417/1‑217 (MQ=255)                 cgGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGG‑CTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTA‑ATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCA‑TCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTTACAAC‑TTGTTCCC‑CGAGAAAATCCCCAGACTGA‑ccccccc                                                                                                                                                                                                               >  1:890865/1‑225 (MQ=255)                          aGGCACAGGTTCAGTTATTTTCCACAGAGAGATCGGGGGTATAAAAAAGCAGGCAACTGGG‑CTCCCAGCGCCGCGGGGGTCTTTTTCACCGCGCTTA‑ATGCCCCCTACGGCGGGTGGCTTGGCAGGCGGGCGTCCA‑TCTCTTTTTTACCAACAATATTCCCTCCCGGCGCCCCGCCTAACAAC‑TTGGCGCG‑GGAGAACCTCCCCAGACTGA‑cgcccgcctgcaaagcc                                                                                                                                                                                                     <  2:1480845/228‑10 (MQ=255)                                                     ggacttgcAGTGTGAAAAGAGACACGCAGCGGGG‑CTCCCAGCCGCGCGGGTATCATTTTACCCGCGTGAA‑ATGCACGCTTAGGCGGGTGGTTTGGCGGGCGGGCGTTCC‑TCTTTTTTTTCCAAACAATAATCCCTGCCTGCGCCCCGCCTACAAAC‑TTGTCCCG‑GGAAAAACCCACCAGACGGA‑CGCCCGCCGCAATCAACTCCTGC‑GGCATTGGCTCTTTCGGCaaa                                                                                                                                                                         <  2:1776529/221‑1 (MQ=255)                                                      gagCTTCAGGGTGTAAATAAAGACGCAACGTGG‑CTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTA‑ATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCA‑TCTCTTTTTTTCCAACAATAATTCCTGCCTGCGGCCCGCCTAAAAAC‑TTGTCGCC‑GGGGAAAATCACCAGAATTA‑AGCCCGCCGCCATCAAATCCTGC‑GGCAATGGCGCCTTTGGCaaaa                                                                                                                                                                        >  2:530272/1‑227 (MQ=255)                                                       agCTTCAGGGTGTAAATAAAGACGCAACGTGG‑CTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTA‑ATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCA‑TCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAAC‑TTGTCGCC‑GGAGAAATTCCCCAAACTGA‑cgcccccc                                                                                                                                                                                                              >  1:613142/1‑190 (MQ=255)                                                          ttaaggttaaaaataaaGACGACACGTGC‑CTCCCCGCGCCGCGAGGGTCATTTTACCCGCCTGTA‑ATCCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCA‑TCTTTTTTTTACCAAAAAAAATCCCTGCCGGCGCCCCGCCAAACAAC‑TTGCCGCG‑GGAAAAACCCCCCAGACTGC‑CGCCCGCCGCAATCAACTCCTGC‑GGCATTGGCTCTTTCGGCAAACCGTa                                                                                                                                                                    <  2:996620/218‑1 (MQ=255)                                                          ttaaGGGTAAAAATAAAGCCGCAACGTGG‑CTTCCGGCGCGGGGGGGTTCTTTTTACCCGCGCTTA‑AGGCCCGCTTCGGGGGGTGCCTTTGCGGGGGGGCGATCA‑TCTTTTTTTACCCAACAATAATTCCTGCCGGCGCCCCGCCTAACAAC‑TTGCCGCG‑GGAAAACCCCCCCAGCCTGA‑CCCCCGCCGCAATCAACTCCTGC‑GGCATTGGCTCTTTCGGCAAACCGTa                                                                                                                                                                    <  2:1970115/225‑1 (MQ=255)                                                          tacAGGGTGTAAAGAAAGACGAAGCGTCG‑CTCCCAGCGCCGCGGGGGTCTTTTTACCCGCGCGTA‑ATGCACGCTTCAGGGGGTGGCTTGGCGGGGGGGCGATCA‑TCTCTTTTTCACAAAAAAAAATCCCTCCCGGCGGCCCGCCTAACAAC‑TTGTCGCG‑GGAAAAACCCACCAGCCGGA‑CGCCCGCCGCAATCAACTCCTGC‑GGCATTGGCTCTTTCGGCAAACCGTa                                                                                                                                                                    <  2:889312/226‑1 (MQ=255)                                                           tCAGGGTGTAAATAAAGACGCAACGTGG‑CTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTA‑ATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCA‑TCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAAC‑TTGTCGCC‑GGAGAAACTCa                                                                                                                                                                                                                                >  1:670271/1‑169 (MQ=255)                                                            cAGGGTGTAAATAAAGACGCAACGTGG‑CTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTA‑ATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCA‑TCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAAC‑TTGTCGCC‑GGAGAAACTCAc                                                                                                                                                                                                                               <  2:669154/169‑1 (MQ=255)                                                             aGGGTGTAAATAAAGACGCAACGTGG‑CTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTA‑ATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCA‑TCTCTTTTTTACCCACAATAATTCCTGCCTGCGGCCCGCCTAACAAC‑TTGTTGCC‑GGAGAAAATCACCAGACTGA‑CGCCCGCCGCAATTAACTCCTGC‑GGCAATGGGTCCTTTGGCCAACCggtct                                                                                                                                                                  >  2:904502/1‑223 (MQ=255)                                                                   taaataaaGACGCAACGTGG‑CTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTA‑ATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCA‑TCTCTTTTTTTCCAACAATAATTCCTGCCCGCGGGCCGCCTAACAAC‑TTGTTGCCTGGAG‑AACTCACCAGACGGA‑CGCCCCCCGCAATCAACTCCCCT‑CGGATTGGCCCATTCGGGCAAACGTACTgggttag                                                                                                                                                           >  2:1021888/1‑223 (MQ=255)                                                                   taaataaaGACGCAACGTGG‑CTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTA‑ATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCA‑TCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAAC‑CTGTCGCC‑GGAGAAACTCACCAGACTGA‑CGCCCGCCGCAATCCACTCCTGC‑GGCATTGGCCCTTTCGGCAAAACGGACTGGCCAAg                                                                                                                                                           >  1:1325405/1‑228 (MQ=255)                                                                         aaGACGCAACGTGG‑CTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTA‑ATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGGTCA‑TCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAAC‑TTTTTGCC‑GGAGAAAATCACCAGACTGA‑CCCCCGCCCCCATCAACTCCTGCGGGCTTGGGCCCTTTCGCCAAACCG‑ACTGGGCAAGATCGa                                                                                                                                                      >  2:2027560/1‑227 (MQ=255)                                                                         aaGACGCAACGTGG‑CTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTA‑ATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCA‑TCTCTTTTTTTCCAACAATAATTCCTGCCTGCGGGCCGCCTAACAAC‑TTGTCGCC‑CGTGAAAATCACCAGACTGA‑CGCCCGCCCCCATCAACTCCCTC‑GGCATTGGGCCCTTCGGCCAACCGt                                                                                                                                                                     >  1:172274/1‑212 (MQ=255)                                                                          aGACGCAACGTGG‑CTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTA‑ATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCA‑TCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAAC‑TTTTCGCC‑GGAGAAACTCACCAGACTGA‑CGCCCGCCGCAATCAACTCCCTC‑GGCAATGGCTCTTTCGGCAAACCGCACTGGCTCAGATCGAc                                                                                                                                                     >  1:1064843/1‑227 (MQ=255)                                                                          aGACGC‑ACGTGGCCTTCCAGCGCCGGGAGGTTCTTTTTACCCGCGCGTA‑ATCCCCGCTCCAGCGGGTGGCTTGGCAGGGGGGCGTTCA‑TCTCTTTTTTCCCAAAAATAATCCCGCCCTGCGGCCCGCCTAACAAC‑TTGTCGCC‑GGAAAAACTCCCCAGACTGC‑CGCCCGCCGCAATCAACTCCTGC‑GGCATTGGCTCTTTCGGCAAACCGTa                                                                                                                                                                    <  2:171976/212‑1 (MQ=255)                                                                           gACGCAACGTGG‑CTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTA‑ATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCA‑TCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAAC‑TTGTCGCC‑GGAGAAACTCACCAGACTGA‑CGCCCGCCGCCATCAACTCCTGC‑GGCATTGGCTCTTTCGGCCAACCGTACTGGCTAAGATCg                                                                                                                                                       >  1:279889/1‑224 (MQ=255)                                                                           gACGCAACGTGG‑CTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTA‑ATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCA‑TCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAAC‑TTGTCGCC‑GGAGAAACTCACCAGACTGA‑CGCCCGCCGCAATCAACTCCTGC‑GGCATTGGCTCTTTCGGCAAACCga                                                                                                                                                                     >  1:1546799/1‑209 (MQ=255)                                                                                      g‑CTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTA‑ATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCA‑TCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAAC‑TTGTCGCC‑GGAGAAACTCACCAGACCGA‑CGCCCCCCCCAATCCAATCCTGC‑GGCATTTGCTCTTATCGCAAACCCTACTGGCCTAGAACGACCCACGAgg                                                                                                                                             >  2:1369615/1‑222 (MQ=255)                                                                                                   gcgAGGTTCATTTGATCGGCG‑GTATAGGCCCGCTTCAGCGGGTGGCTGTGCGGGCGGCCGATCA‑TCTCTTTTTCCCCAACAATAACTCCGCCCCGCGCCCCGCCTAACATC‑TTGTCGCC‑GGAGAACCCCACCAGCCGGA‑CGCCCGCCGCAATCAACTCCTGC‑GGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGt                                                                                                                             <  2:1631154/227‑1 (MQ=255)                                                                                                                      cgcgTA‑ATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCA‑TCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAAC‑TTGTCGCC‑GGAGAAACTCACCAGACTGA‑CGCCCGCCGCAATCAACCCCTGC‑GGCAt                                                                                                                                                                                         >  1:448543/1‑148 (MQ=255)                                                                                                                       gggTA‑AGGCCCGCTTAGGGGGGTGCCTTGGCGGGGGGCCGCTAC‑TCTCTTTTTTCCAAACAAAAATCCCTCCCGGCGCCCCGCCTAACACC‑TTGTCGCG‑GGAAAACCCCCCCAGACGGC‑CGCCCGCCGCAATCAACTCCTGC‑GGCAtt                                                                                                                                                                                        <  2:447808/146‑1 (MQ=255)                                                                                                                             aTGCACCGTTGAGCGGATGGCTTGGGCGGCGCGCG‑TCACTCTTTTTTGCACCAACAAAACCTCCTGCCGCCGCCCCGCCAAACAAC‑TTGCCGCC‑GGAAAACCCCCCCAGACGGC‑GGCCCGCCGCAATCAACTCCTGC‑GGCATTGGCTCTTTCGGCAACCCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGc                                                                                                   <  2:653983/228‑1 (MQ=255)                                                                                                                                cACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCA‑TCTCTTTTTTACCCACAATAATTCCTGCCTGCGGCCCGCCTAACAAC‑TTGTCGCC‑GGAGAAACTCACCAGACTGA‑CGCCCGCCGCAATCAACTCCTGC‑GGCATTGGCTCTTTCGGCAAAACGCACTGGCTAAGATCGACCAGCGAGCCCCTGCCTAAATAAGTCACTACCGGTACATCCAGCTCTTTgtct                                                                                                 >  2:1195787/1‑224 (MQ=255)                                                                                                                                    cTTCAGCGGGTGGCTTTGCAGGCGGGCGATCA‑TCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAAC‑TTGTCGCC‑GGAGAAAATCACCAGACTTA‑CGCCCGCCGCAATCAACTCCTGC‑GGCCTTGGCTCTTTCGGCAAACCGt                                                                                                                                                                     >  1:1833582/1‑155 (MQ=255)                                                                                                                                    cTTCAGCGGGTGGCTTTGCAGGCGGGCGATCA‑TCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAAC‑TTGTCGCC‑GGAGAAAATCACCAGACTGA‑CGCCCGCCGCAATCAACTCCTTC‑GGCATTGGCTCTTTCTGCAAAACCTACTGGCTAAAATCGACCAGCGAGCCACTGCCTAAATCCGTACCCACCGGTACATCCCGACCTTTTCCGAggg                                                                                             >  2:1803607/1‑225 (MQ=255)                                                                                                                                     ttCAGCGGGTGGCTTTGCAGGCGGGCGATCT‑TCTTTTTTTCACCAACAATATTTCCTGCCTGCGCCCCCCCTAACACC‑TTGCCGCG‑GGAAAAACTCACCAGACTGA‑CGCCCGCCGCAATCAACTCCTGC‑GGCATTGGCTCTTTCGGCAAACCGTa                                                                                                                                                                    <  2:1830337/155‑1 (MQ=255)                                                                                                                                          cGGGTGGCTTTGCAGGCGGGCGATCA‑TCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAAC‑TTGTCGCC‑GGAGAAACTCACCAGACTGA‑CGCCCGCCGCAATCAACTCCTGC‑GGCATTGGCTCTTTCGGCAAACCGt                                                                                                                                                                     >  1:68077/1‑149 (MQ=255)                                                                                                                                               ggCTTTGCAGGCGGGCGATCA‑TCTCTTTTTTTCCAACAATAATTCCTGCCTTCGGCCCCCCCCACAAA‑CTTTCGCC‑GGAGAAAATAAACCGACTGA‑CACCCCCCGCAAATAACTCCTCC‑GGGCTTGGCTCTTTTgggc                                                                                                                                                                           >  2:427988/1‑136 (MQ=255)                                                                                                                                               ggCTTTGCAGGCGGGCGATCA‑TCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAAC‑TTTTCGCC‑GGAGAAACTCACCAGACTGA‑CGCCCGCCGCAATCAACTCCTGC‑GGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCCCTACGGGAACATCCAGCCCCTTGCCGAGCGCCACCCGTTca                                                                                  >  2:92797/1‑226 (MQ=255)                                                                                                                                               ggCTTTGCAGGCGGGCGATCA‑TCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAAC‑TTTTCGCC‑GGAGAAACTCACCAGACTGA‑CGCCCGCCGCAATCAACTCCTGC‑GGCATTGGCTCCTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACAACCAGCCCTTTTCCGAGAGCCACCAGtttc                                                                                  >  2:1717752/1‑225 (MQ=255)                                                                                                                                               ggCTTTGCAGGCGGGCGATCA‑TCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCCAACAAC‑TTTTCGCC‑GGAGAAACTCACCAGACTGA‑CGCCCGCCGCAATCAACTCCTGC‑GGCATTGGGCCTTTCGGCAAACCGTACTGGCTAAGAACGACCAGCGAGCCACTGCCTAAATCAGTCACTACCGGAACATCCCGCTCCTTTCCCAGCGCCACCCGGTcc                                                                                  >  2:1077707/1‑227 (MQ=255)                                                                                                                                               ggCTTTGCAGGCGGGCGATCA‑TCTCTTTTTTACCAAAAATAATTCCTGCCTGCGGCCCGCCTAAAAAC‑TTGTCGCC‑GGAGAAACTCACCAGACTGA‑CGCCCGCCGCAATCAACTCCTGG‑GGCATTGGCTCTTTTGGCAAAACGGACTGGCTAAGATCGACCCGCGAGCCACTGCCTAAATCAGTCACTAAGGGGAAAACCAGCTCTTTGACGAAGGCCACCAgattcg                                                                                 >  2:1100169/1‑223 (MQ=255)                                                                                                                                                ccTTTGCAGGCGGGCGATCA‑TCTTTTTTTTCCAAAAAAAAATTCCTCCCTGCGGCCCGCCTAACAAC‑TTGTCGCC‑GAAGAACCTCACCAGACTGA‑CGCCCGCCGCAATCAACTCCTGC‑GGCATTGGCTCTTTCGGCaa                                                                                                                                                                          <  1:428702/137‑1 (MQ=255)                                                                                                                                                                   a‑tcTCTTTTTAACAAACAATCTTCCGTGCCTGCGGCCCGCCTAACAAC‑TTGTCGCA‑GGAAAACCTCACCAGCCTGA‑CGCCCGCCGCAATCAACTCCTGC‑GGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTg                                                             <  1:92957/228‑1 (MQ=255)                                                                                                                                                                     tcttTTTTTCACCAACAAATACTCCCGCCTGCGCCCCGCCTAACACC‑TTCCCGCC‑GGAGAACCCCCCCCATCTGA‑CCCCCGCCGCAATCAACTCCTGC‑GGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTg                                                             <  1:1102146/227‑1 (MQ=255)                                                                                                                                                                     tctcTTTTTCAACAACCATAATCCCTGCCTGCCCCCCTCCTAACATC‑GTGCCGCC‑GGAAAACCCCACCAGACCGA‑CCCCCGCCGCAATCAACCCCTGC‑GGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTg                                                             <  2:55549/227‑1 (MQ=255)                                                                                                                                                                      ctctTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAAC‑TTGTCGCC‑GGAGAAACTCACCAGACTGA‑CGCCCGCCGCAATCAACTCCTGC‑GGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTAACGCTTTGGTg                                                             >  1:1627840/1‑226 (MQ=255)                                                                                                                                                                          tttttACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAAC‑TTGTCGCC‑GGAGAAACTCACCAGACTGA‑CGCCCGCCGCAATCAACTCCTGC‑GGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGTTGAAcccc                                                       >  1:479855/1‑228 (MQ=255)                                                                                                                                                                           ttctACCAACAATAATTCCTGCCTCCGCCCCCCCTAACAAC‑GTCCCGCC‑GGACAACCCCCCCAGCCTGC‑CGCCCGCCGCAATCAACTCCTGC‑GGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCt                                                      <  2:559387/225‑1 (MQ=255)                                                                                                                                                                            tttACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAAC‑TTGTCGCC‑GGAGAAACTCACCAGACTGA‑CGCCCGCCGCAATCAACTCCTGC‑GGCATTGGCTCTTTCGGCAAACCGt                                                                                                                                                                     >  2:1072858/1‑116 (MQ=255)                                                                                                                                                                             ttACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAAC‑TTGTCGCC‑GGAGAAACTCACCAGACTGA‑CGCCCGCCGCAATCAACTCCTGC‑GGCATTGGCTCTTTCGGCAAACCGTa                                                                                                                                                                    <  1:1074781/116‑1 (MQ=255)                                                                                                                                                                              tACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAAC‑TTGTCGCC‑GGAGAAACTCACCAGACTGA‑CGCCCGCCGCAATCAACTCCTGC‑GGCATTGGCTCTTTCGGCAAACCGTACTGGCCAAGCTCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCCTCTATCGCTTTGGTGAACCCCTgc                                                    >  1:579358/1‑226 (MQ=255)                                                                                                                                                                                  aaaaataatCCCTGCCTGCGCCCCG‑CAAACAACTTTGCCCCC‑GGAAAACCTCCCCAGCCGGA‑CGCCCGCCGCAATCAACTCCTGC‑GGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCt                                               <  1:1079639/225‑1 (MQ=255)                                                                                                                                                                                   acaaTAATTCCTGCCTGCGGCCCGCCTAACAAC‑TTGTCGCC‑GGAGAAACTCa                                                                                                                                                                                                                                >  2:1644045/1‑52 (MQ=255)                                                                                                                                                                                   acaaTAATTCCTGCCTGCGGCCCGCCTAACAAC‑TTGTCGCC‑GGAGAAACTCACCAGACTGA‑CGCCCGCCGCAATCAACTCCTGC‑GGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTg                                              >  2:1403086/1‑228 (MQ=255)                                                                                                                                                                                    caaTAATTCCTGCCTGCGGCCCGCCTAACAAC‑TTGTCGCC‑GGAGAAACTCAc                                                                                                                                                                                                                               <  1:1646958/52‑1 (MQ=255)                                                                                                                                                                                     actactccccgcctgcCGCCCGCCAAACAAT‑TGGGCGCC‑GGACAACCCCCCCAGCCGCA‑CGCCCGCCGCAATCAACTCCTGC‑GGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTa                                            <  2:479079/218‑1 (MQ=255)                                                                                                                                                                                               gcctgcGGCCCGCCTAACAAC‑TTGTCGCC‑GGAGAAACTCACCAGACTGA‑CGCCCGCCGCAATCAACTCCTGC‑GGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTg                                              >  1:869191/1‑216 (MQ=255)                                                                                                                                                                                               gcctgcGGCCCGCCTAACAAC‑TTGTCGCC‑GGAGAAACTCACCAGACTGA‑CGCCCGCCGCAATCAACTCCTGC‑GGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGt                                   >  1:1561662/1‑227 (MQ=255)                                                                                                                                                                                               gcctgcGGCCCGCCTAACAAC‑TTGTCGCC‑GGAGAAACTCACCAGACTGA‑CGCCCGCCGCAATCAACTCCTGC‑GGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTAATggac                                  >  2:439920/1‑226 (MQ=255)                                                                                                                                                                                               gcctgcGGCCCGCCTAACAAC‑TTGTCGCC‑GGAGAAACTCACCAGACTGA‑CGCCCGCCGCAATCAACTCCTGC‑GGCATTGGCTCTTTCGGCAAACCGTACCGGCTAAGATCGACCAGCGAGCCACTTCCTAAATCAGTCACTACcgg                                                                                                                    >  1:1829214/1‑146 (MQ=255)                                                                                                                                                                                                   gcgcCCCGCAAAACACC‑TTGCCGCC‑GGAAAACCTCACCAGCCTGA‑CCCCCGCCGCAATCAACTCCTGC‑GGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTa                              <  1:1197933/228‑1 (MQ=255)                                                                                                                                                                                                          gCCTACCAAC‑TTCGCGCG‑GAAGAACCTCCCAAGCCGCA‑CCCCCGCCGCAATCAACTCCTGC‑GGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTtcat                       <  1:1372085/228‑1 (MQ=255)                                                                                                                                                                                                              taacgtgtcggcGCC‑GGAGAAACTCACCGTACGCA‑CCCCCGCCGCAATCAACTCCGGC‑GGCATTGGCTCTTTCGGCAAGCCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCaa                    <  2:1728002/217‑1 (MQ=255)                                                                                                                                                                                                                   c‑TTGTCGCC‑GGAGAAACTCACCAGACTGA‑CGCCCGCCGCAATCAACTCCTGC‑GGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGc              >  2:670581/1‑228 (MQ=255)                                                                                                                                                                                                                      tGTCGCC‑GGAGAAACTCACCAGACTGA‑CGCCCGCCGCAATCAACTCCTGC‑GGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAGTCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAAcccc                                                       >  1:1154897/1‑185 (MQ=255)                                                                                                                                                                                                                      tGTCGCC‑GGAGAAACTCACCAGACTGA‑CGCCCGCCGCAATCAACTCCTGC‑GGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAg                                                                                                                              >  1:1243910/1‑114 (MQ=255)                                                                                                                                                                                                                       gTCGCG‑GGAAAACCCCCCCAGCCGGA‑CGCCCGCCGCAATCAACTCCTGC‑GGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGt                                                                                                                             <  2:1241688/114‑1 (MQ=255)                                                                                                                                                                                                                       gTCGCC‑GGAGAAACTCACCAGACTGA‑CGCCCGCCGCAATCAACTCCTGC‑GGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAAcccc                                                       >  1:1171329/1‑184 (MQ=255)                                                                                                                                                                                                                       gTCGCC‑GGAGAAACTCACCAGACTGA‑CGCCCGCCGCAATCAACTCCTGC‑GGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTgg                                    >  1:1073807/1‑203 (MQ=255)                                                                                                                                                                                                                       gTCGCC‑GGAGAAACTCACCAGACTGA‑CGCCCGCCGCAATCAACTCCTGC‑GGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGCATGTACTTTCATCAACAGTGCgg            >  1:1445825/1‑227 (MQ=255)                                                                                                                                                                                                                       gTCGCC‑GGAGAAACTCACCAGACTGA‑CGCCCGCCGCAATCAACTCCTGC‑GGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCCCTGGTGAACCCCCGAAGGCTGTAGTTACTGGAATGt                               >  2:1491117/1‑208 (MQ=255)                                                                                                                                                                                                                        ccccc‑GGAGAAACACACCAGACGCACCCCCCCCCGCAATCAACTCCTGC‑GGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGt                                   <  2:1071886/201‑1 (MQ=255)                                                                                                                                                                                                                         gcccgggAGAAGCCCAGCAGCCTGC‑CCCCCGCCGCAATCAACTCCTTC‑GGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTCTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTa                              <  1:1493786/203‑1 (MQ=255)                                                                                                                                                                                                                         cGCC‑GGAGAAACTCACCAGACTGA‑CGCCCGCCGCAATCAACTCCTGC‑GGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTa          >  1:264883/1‑227 (MQ=255)                                                                                                                                                                                                                            c‑GGAGAAACTCACCAGACTGA‑CGCCCGCCGCAATCAACTCCTGC‑GGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACAt                                                                                                               >  2:207131/1‑123 (MQ=255)                                                                                                                                                                                                                              ggAGAAACTCACCAGACTGA‑CGCCCGCCGCAATCAACTCCTGC‑GGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATc                                                                                                              <  1:207487/123‑1 (MQ=255)                                                                                                                                                                                                                                 aaaaCTCCACAGACTGA‑CGCCCGCCGCAATCAACTCCTGC‑GGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTttcatt  <  1:178217/227‑1 (MQ=255)                                                                                                                                                                                                                                       cACCAGACTGA‑CGCCCGCCGCAATCAACTCCTGC‑GGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGc                                                                                >  2:1646834/1‑144 (MQ=255)                                                                                                                                                                                                                                       cACCAGACTGA‑CGCCCGCCGCAATCAACTCCTGC‑GGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTttcat   >  2:177912/1‑221 (MQ=255)                                                                                                                                                                                                                                        aCCAGACTGA‑CGCCCGCCGCAATCAACTCCTGC‑GGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCt                                                                               <  1:1649752/144‑1 (MQ=255)                                                                                                                                                                                                                                        |                                                                                                                                                                                                                               AATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGG‑CTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTA‑ATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCA‑TCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAAC‑TTGTCGCC‑GGAGAAACTCACCAGACTGA‑CGCCCGCCGCAATCAACTCCTGC‑GGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTTTCATT  >  NC_000913/3760158‑3760604

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