Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,760,381	T→C	19.0%	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	19.0%	110.6 / 3.4	60	S290G (AGT→GGT)	selA	selenocysteine synthase
Reads supporting (aligned to +/- strand):  ref base T (33/14);  new base C (4/9);  total (37/23)
Fisher's exact test for biased strand distribution p-value = 2.13e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGT‑AAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGT‑AATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACA‑ATAATTCCTGCCTGCGGCCCGCCTAACAACTTGT‑CGCCGGAG‑AAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGC‑CACTGCCTAAATCAGTCACTA‑CGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATC‑TATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTTTCATTCAC  >  NC_000913/3760146‑3760606                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                     gCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGT‑AAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGT‑AATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGGTCATCTCTTTTTTACCAACA‑ATAATTCCTCCCTGCGGCCCGCCTAAAACCTTGT‑CGCCGGGG‑AAACTcc                                                                                                                                                                                                                                    >  2:1230866/1‑237 (MQ=255) gCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGT‑AAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGT‑AATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTGTTTTTTACCCACA‑ATAATTCCTGCCTACGGCCCGCCTAACAACTTGT‑CGCCCGCG‑AAAATCa                                                                                                                                                                                                                                    >  2:687135/1‑238 (MQ=255)                aCCTCTGCGCTGCGGGTAAGCAGGCGCCGGGTCGGTAATTTTTCACTCAGAGCCTCAGGGTGT‑AAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCCCGT‑AATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACA‑ATAATTCCTGCCTCCGGCCCGCCTAACAACTTGT‑CGCCGGCG‑AAACCCCCCCAACCGAcgcccg                                                                                                                                                                                                                     >  2:355019/1‑238 (MQ=255)                aCCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGT‑AAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGT‑AATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCCCTTTTTTACCCACA‑ATAATTCCTGCCTGCGGCCCGCCTAACAACTTTG‑CGCCGGGG‑AAACTCACCCACCTGAcccccc                                                                                                                                                                                                                     >  2:1143056/1‑237 (MQ=255)                     tgcgctgcgGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGT‑AAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGT‑AATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACA‑ATAATTCCTGCCTGCGGCCCGCCTAACAACTTTT‑CGCCGGGG‑AAACTCACCCGACTGAACCCCGCCCCa                                                                                                                                                                                                                >  2:867865/1‑238 (MQ=255)                                  aGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGT‑AAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGT‑AATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGGTCATCTCTTTTTTACCAACC‑ATAATTCCTGCCTGCGGCCCGCCCTACCACTTTT‑CCCCGGGG‑AAACCCCCCCGCCCGACCCCCCCCCCCATCAAccccctggg                                                                                                                                                                                                  >  2:587246/1‑234 (MQ=255)                                    cAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGT‑AAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGT‑AATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGGTCATCTCTTTTTTACCAACA‑ATAATTCCTGCCTGCGGCCCGCCTAACAACTTGT‑CGCCGGAG‑AAACTCACCAGACCGACGCCCGCCGCAATCCACCCCCGCGGCa                                                                                                                                                                                                >  2:1079695/1‑239 (MQ=255)                                        cgcAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGT‑AAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGT‑AATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGGCCATCTCTTTTTTACCAACA‑ATAAATCCTGCCTTCGGCCCGCCCACAAACTTGT‑TGCCGGCGCAAAACCCCCACCTTGAC‑CCCGCCCCCAACAACTCCTGCggacatgt                                                                                                                                                                                            >  1:825006/1‑233 (MQ=255)                                                    tttttttgcaccaggtgttgaGGGTGTAAAAAAAAGACGCAGCGTGTGTGCCGGGGCGGGGGTGGTCATTTTTCCCGCGGTT‑AATACGCGCTTAGGGGGGGGGTTTGGCAGGGGGGCGACCTTCCTTTTTTTCCAAAAA‑ATAATCCCGCCCTGGGCCCCCCCTACCAACTGGC‑GGCGGGAA‑AAACTCCCCAGACTACCGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCa                                                                                                                                                                                 <  1:988026/219‑1 (MQ=255)                                                         ttCACTCAGAGCTTCAGGGTGT‑AAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGT‑AATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGGACATCTCTTTTTTACCAAAA‑AAAATTCCTGCCCGGTGCCCGCCTAAAACCTTTT‑TGCCGGAG‑AAACTCACCCACATCACGCCCGCACCCACCAACCCCTCCGGGAGTGGCGCTTTCGGCAAACCgc                                                                                                                                                                           >  1:1044934/1‑238 (MQ=255)                                                         ttCACTCAGAGCTTCAGGGTGT‑AAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGT‑AATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACA‑ATAATTCCTGCCTGCGGCCCGCCTAACAACTTGT‑CGCCGGAG‑AAACTCACCCGACTGACGCCCCCCCCAATCAACTCCTTCGGCATTTGCCCTTTTGGCCAACCg                                                                                                                                                                            >  2:844745/1‑238 (MQ=255)                                                          tCACTCAGAGCTTCAGGGTGA‑AAAAAAAGGCGCAACGTGGCTTCCAGCGCCGGGGGGGTCATTTCATCCGGGCGT‑AATGCACGCTTCAGCGGGTGGCTTGGCAGGCGGGCGATCATCTTTTTTTTACCAACA‑ATAATCCCTGCCTGCGGCCCGCCTAACAACTTGT‑CGCCGGAG‑AAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTa                                                                                                                                                                          <  2:300985/239‑1 (MQ=255)                                                                   gtttCAGGGTGT‑AAAAAAAGCCGCAACGGGGCTTCCAGCCCCGCGGGGGTTTTTTTACCCGCGCGT‑AATGCCCGTTTCAGGGGGTGGCTTTGGAGGGGGGCGATCATCTTTTTTTTACCAACA‑ATAATTCCTGCCTGCGGCCCGCCTAACAACTTGT‑CGCCGGAG‑AACCTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTa                                                                                                                                                                   <  2:3121/235‑1 (MQ=255)                                                                          ggTGT‑AAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGT‑AATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGGTCATCTCTTTTTTACCAACC‑ATAATTCCTGCCTGCGGCCCGCCTAACAACTTGT‑CGCCGGAG‑AAACTCACCAGACTCACGCCCCCCGCAACCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGCCCTGGcgcacgactacc                                                                                                                                                          >  2:464348/1‑228 (MQ=255)                                                                              t‑aaataaaGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGT‑AATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCCCTTTTTTAACAACA‑ATAATTCCCGCCGCCGGCCCGCCTACCAAATTGT‑TCCCCGAG‑AAAACCACCAGCACGACACCCGCCGCAACCACCTCCTTCGGCACTGCCCCTCTCGCGCAACCGCACtccccaaaaacgacacac                                                                                                                                                       >  1:975123/1‑221 (MQ=255)                                                                                 aataaaGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGT‑AATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACA‑ATAATTCCTGCCTGCGGCCCGCCTAACAACTTGT‑CGCCGGAG‑AAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTCCGGCATTGGCCCTTTCGGCAAACCCGACTGGCTAAAATCGGCCAGCGa                                                                                                                                                     >  2:345643/1‑238 (MQ=255)                                                                                           cAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGT‑AATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACA‑ATAATTCCTGCCTGCGGCCCGCCTAACAACTTTT‑CGCCGGAG‑AAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCCCCTTCGGCAAACCGTACTGGCTAAAATCGACCACCGAGC‑CACTGCCtc                                                                                                                                         >  2:190514/1‑238 (MQ=255)                                                                                                   ctccgaGCGCCGGGAGGGTATTTTTAGCCGC‑CGTAAATGCCCGTTTCGGGGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACA‑ATAATTCCTGCCTGCGGCCCGCTTAACAACTTGC‑CGCCGGAG‑AAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGC‑CACTGCCTAAATCAGt                                                                                                                                  <  2:876444/233‑1 (MQ=255)                                                                                                   cTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGT‑AATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACA‑ATAATTCCTGCCTGCGGCCCGCCTAACAACTTTT‑CGCCGGAG‑AAACTCACCAGACTGACGCCCGCCGCAATCAACTCCCTCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAAATCGACCAGGCAGC‑CCCTGCCTaaaaaagg                                                                                                                                  >  2:378091/1‑233 (MQ=255)                                                                                                                       aTTTTATCCGCGCGT‑AATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACA‑ATAATTCCTGCCTGCGGCCCGCCTAACAACTTGT‑CGCCGGAG‑AAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCCCTTTCGGCCAACCGTACTGGCTAAAATCGACCAGCGAGC‑CACTGCCTAAATCAGTCACTA‑CGGGAAAATCCAGc                                                                                                              >  2:255606/1‑237 (MQ=255)                                                                                                                                   cgT‑AATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACA‑ATAATTCCTGCCTGCGGCCCGCCTAACAACTTGT‑CGCCGGAG‑AAACTCACCAGACTGACGCCCGCCGCCATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGC‑CACTGCCTAAATCAGTCACTA‑CGGGAACATCCCCCCCTTTTCCCAgcg                                                                                                 >  2:793566/1‑238 (MQ=255)                                                                                                                                                 tggggggggggTTTGGCGGGCGGGCGACCACTTTTTTTTTACCAACA‑ATAATCCCTGCCTGCGGCCCGCCAAACAATTTGC‑CGCCGGAG‑AAACCCCCCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGC‑CACTGCCTAAATCAGTCACTA‑CGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCt                                                                                   <  2:376264/234‑1 (MQ=255)                                                                                                                                                   aGCGTGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACA‑ATAATTCCTGCCTGCGGCCCGCCTAACAACTTGT‑CGCCGGAG‑AAACTAACCCGAATGGCGCCCACCGCAATATACTCCTTCGGCATTGGCCCTTTCGGCAAAACTGACTGGGGAA‑ATCGACCAGCGGAC‑CACCGCCCAAAACAGACACCACCGGGAAAAACAAGTCCTTTGCCCACCGCC‑CCCGTTCCCCTTc                                                                                 >  1:916634/1‑238 (MQ=255)                                                                                                                                                         tGGTTTTGCGGGCGGGCGAACATCTCTTTTTACCCAAAA‑ATAATTCCTGCCTGGGGCCCGCCTAACAACTGGT‑CGCGGGAG‑AACCTCCCCAGACTGCCGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGC‑CACTGCCTAAATCAGTCACTA‑CGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATC‑Ta                                                                           <  2:390367/238‑1 (MQ=255)                                                                                                                                                          ggCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACA‑ATAATTCCTGCCTGCGGCCCGCCTAACAACTTGT‑CGCCGGAG‑AAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCCAACCGTACTGGCTAAGATCGACCAGCGAGC‑CACTGCCTAAATCAGGCACTA‑CGGGAACATCCAGCTCTTTGCCGAACGCCCCCAGTTCCGCTTCATC‑TAt                                                                          >  2:269219/1‑238 (MQ=255)                                                                                                                                                          ggCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACA‑ATAATTCCTGCCTGCGGCCCGCCTAACAACTTGT‑CGCCGGAG‑AAACTCACCAGACTGACGCCCGCCCCAAACAACTCCTGCGGCAATGGCTCCTTCGGCAAACCGTAATGGCTAAAATCGACCCACGAGC‑CCCATCCCAAAACAACCACCA‑CGGGGACAACCAGCCCTCTTCCACGACCCACCACATCCGCCTCacacac                                                                           >  1:525163/1‑233 (MQ=255)                                                                                                                                                          ggCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACA‑AATATTCCTGCCTGCGGCCCGCCTAACAACTTTT‑CGCAGGGG‑AAAATCACCAGACTCACGCCCGCCCCGATCAACTCCCCCGGCATTGGCTATTTCGTCAAAACGCACTCCATAAGGTCGCCCA‑CGAGCACCACACCTATATCCAACACAA‑CGGGAAACTACCGCTCTTTTCCCAGCTCCCCCAATTACGCTTCCAC‑TAt                                                                          >  1:822992/1‑238 (MQ=255)                                                                                                                                                                    gggggcgATCATCTCTTTTTTACCAAAA‑ATAATTCCTGCCTGCGGCCCGCCTAACAACTTGC‑CGCGGGAA‑AAACTCACCAGACGGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGC‑CACTGCCTAAATCAGTCACTA‑CGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATC‑TATCGCTTTGGTg                                                                <  2:975123/236‑1 (MQ=255)                                                                                                                                                                    ggagctccatcatttttttttCCCAAAATATAATCCCGGCCCGCGGCCCGCATAACATCTTGA‑CTCAGAAA‑AACCCCACCAACCTCACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTACGGCAAACCGTACTGGCTAAGATCGACCAGCGAGC‑CACTGCCTAAATCAGTCACTA‑CGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATC‑TATCGCTTTGGTg                                                                <  1:404042/224‑1 (MQ=255)                                                                                                                                                                      gggcgATCATCTCTTTTTTACCAACA‑ATAATTCCTGCCTGCGGCCCGCCTAACCACTTGT‑CGCCGGAG‑AAAATCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGC‑CACTGCCTAAATCAGTCACTA‑CGGGAACATCCAGCTCTTTGCCGGGCGCCACCAGTTCCGCTTCATC‑TATCGCTTTGGTGaa                                                              >  2:1157421/1‑238 (MQ=255)                                                                                                                                                                       ggcgGTCATCTCTTTTTTACCAACA‑ATAATTCCTGCCTGCGGCCCGCCTAACAACTTTT‑CGCCGGAG‑AAACTCACCAGACTGACGCCCCCCGCAAACCACTCCTGCGGCATTGGCTCTTTCGGCCAACCGGACTGGCTAAAATCGACCAACGAGC‑CCACGCCCAAATCAGTCACTA‑CCGGCAACCCCAGCTCTTTTCCCAGCGCCCCCCGTTCCCCCCCCTC‑TATCCCTTTGGTGAAc                                                             >  1:853902/1‑238 (MQ=255)                                                                                                                                                                       gcgcgtttctttttttttACCACAA‑ATATTTCCCGCCGGGGGCCCCACTAACACTTTGC‑CGCGGGAG‑AACCTCACCAGCCCGCCGCCCGCCGAAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGC‑CACTGCCTAAATCAGTCACTA‑CGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATC‑TATCGCTTTGGTGAAcc                                                            <  1:1143056/227‑1 (MQ=255)                                                                                                                                                                            tcatcTCTTTTTTACAAAAA‑ATATTCCCCGCCTGCGCCCCCCCTAACACCTTGC‑CGCCGGAA‑AACCCCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGC‑CACTGCCTAAATCAGTCACTA‑CGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATC‑TATCGCTTTGGTGAACCCCTGa                                                       <  2:1044934/239‑1 (MQ=255)                                                                                                                                                                                      tttACCAACA‑ATAATTCCTGCCTGCGGCCCGCCTAACAACTTGT‑CTCCGGAG‑AAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGCCTCTTTCGGCAAACCGTACTGGCTAAGATCGAACAGCGAGC‑CACGGCCCAACTCAGTCCCTA‑CGGGAACATCCAGCTCCTTGCCGAGCGCCACCAGTTCCCCTTCATC‑TATCGCTTTGCTGGAACCCCGCATGCTgcaa                                              >  1:190431/1‑235 (MQ=255)                                                                                                                                                                                      tttACCAACA‑ATAATCCCTCCCTGCGGCCCCCCTAACAATTTGC‑CGCCGGAA‑AAACCCCCCAGACGGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGC‑CACTGCCTAAATCAGTCACTA‑CGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATC‑TATCGCTTTGGTGAACCCCTGAATGCTGTAGt                                             <  2:191951/239‑1 (MQ=255)                                                                                                                                                                                        tcccAACA‑ATAACCCCGGCCGGCGGCCCGCTTAACAACTTGT‑CGCCGGAG‑AACCTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGC‑CACTGCCTAAATCAGTCACTA‑CGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATC‑TATCGCTTTGGTGAACCCCTGAATGCTGTAGTTa                                           <  2:825006/237‑1 (MQ=255)                                                                                                                                                                                        tACCAACA‑ATAATTCCTGCCTGCGGCCCGCCTAACAACTTGT‑CGCCGGAG‑AAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTg                                                                                                                                                                                                     >  1:1124879/1‑88 (MQ=255)                                                                                                                                                                                        tACCAACA‑ATAATTCCTGCCTGCGGCCCGCCTAACAACTTGT‑CGCCGGAG‑AAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTg                                                                                                                                                                                                     <  2:1124879/88‑1 (MQ=255)                                                                                                                                                                                        ccccAACA‑AAAATCCCCGCCCGGCGCCCCCCAAACAACTTGC‑CGCCGGAA‑AAACTCACCAGACGGCCGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGC‑CACTGCCTAAATCAGTCACTA‑CGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATC‑TATCGCTTTGGTGAACCCCTGAATGCTGTAGTTa                                           <  2:1223912/237‑1 (MQ=255)                                                                                                                                                                                             aca‑aTAATTCCTGCCTGCGGCCCGCCTAACAACTTGT‑CGCCGGAG‑AAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGC‑CACTGCCTAAATCAGTCCCCA‑CGGGAACATCCCGCCCTTTGCCGAGCGCCCCCAGTTCCGCTTCATC‑TATCCCTTTGGTGAACCCCTGAATGCCGGATTTACTggg                                      >  1:1223912/1‑238 (MQ=255)                                                                                                                                                                                               a‑ataatCCCGGCCTGCGGCCCCCCTAACACTTTGT‑CGCCGGAG‑AACCTCCCCAGACTGCCGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGC‑CACTGCCTAAATCAGTCACTA‑CGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATC‑TATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTAt                                    <  2:711271/239‑1 (MQ=255)                                                                                                                                                                                                   tatccctcccgggggCCCGCCAAACAACTTGT‑CGCCGGAA‑AAACCCCCCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGC‑CACTGCCTAAATCAGTCACTA‑CGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATC‑TATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTa                                 <  2:853902/226‑1 (MQ=255)                                                                                                                                                                                                   agcccgggcacgcGGCCCGCTTAACACTTTGG‑CGCGAGAG‑AACCCCCACAGACGGGCGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGC‑CACTGCCTAAATCAGTCACTA‑CGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATC‑TATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTa                                 <  2:332376/228‑1 (MQ=255)                                                                                                                                                                                                   actccctgcctgCGGCCCGCCTACCAACTTGT‑CGCCGGAA‑AACCCCCCCAGACGGCCGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGC‑CACTGCCTAAATCAGTCACTA‑CGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATC‑TATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTa                                 <  2:349300/235‑1 (MQ=255)                                                                                                                                                                                                   actcccggcctgcGGCCCGCCTAACAACTGGC‑CGCCGGAA‑AACCTCCCCAACCTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGC‑CACTGCCTAAATCAGTCACTA‑CGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATC‑TATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTa                                 <  2:759039/232‑1 (MQ=255)                                                                                                                                                                                                   aactCCGGCCTGCGCCCCGCCAAAAAACTTGG‑CGCGGGAA‑AAGCTCACCAGCCCGACGCCCGCCGCAATCAACTTCTGCGGCATTGGCTTTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGC‑CACTGCCTAAATCAGTCACTA‑CGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATC‑TATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTa                                 <  2:810825/236‑1 (MQ=255)                                                                                                                                                                                                   aacccccgccggcGCCCCGTCAAACAAGTTGC‑CGCCAGAG‑AACCCAACCAGCCAGACCCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGC‑CACTGCCTAAATCAGTCACTA‑CGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATC‑TATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTa                                 <  2:46052/232‑1 (MQ=255)                                                                                                                                                                                                     ttCCCGCCGGCGGCCCCCCTAA‑AACCTGTGCCGGGGAG‑AAACCCACCAGCATGCCGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGC‑CACTGCCTAAATCAGTCACTA‑CGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATC‑TATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACt                               <  2:854896/239‑1 (MQ=255)                                                                                                                                                                                                       cctgcctgCGGCCCGCCTAACAACTTGT‑CGCCGGAG‑AAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGC‑CACTGCCTAAATCAGTCACTA‑CGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATC‑TATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACtt                              >  2:1065417/1‑238 (MQ=255)                                                                                                                                                                                                         tgcctgcGGCCCGCCTAACAACTTGT‑CGCCGGAG‑AAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGC‑CACTGCCTAAATCAGTCACTA‑CGGGAACATCCAGCTCTTTGCCGGACGCCACCAGTTCCGCTTCCTC‑TATCGCTTTGTGGAACCCCCGAATGCTGTATTTACGGGGAttaccttt                             >  1:554810/1‑230 (MQ=255)                                                                                                                                                                                                          gcctgcGGCCCGCCTAACAACTTGT‑CGCCGGAG‑AAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGC‑CACTGCCTAAATCAGTCACTA‑CGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATC‑TATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTtca                           >  2:3758/1‑238 (MQ=255)                                                                                                                                                                                                          gcctgcGGCCCGCCTAACAACTTGT‑CGCCGGAG‑AAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGC‑CACTGCCTAAATCAGTCACTA‑CGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATC‑TATCGCTTTCGTGAACCCCTGACTGCTGTAGTTCCTTGTATGTACTTtca                           >  2:554768/1‑238 (MQ=255)                                                                                                                                                                                                          gcctgcGGCCCGCCTAACAACTTGT‑CGCCGGAG‑AAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGC‑CACTGCCTAAATCAGTCACTA‑CGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATC‑TACCGCTTTGGTGAACCCCTGGAGGCTTTAGTTACTGGTATGTACTTtcat                          >  1:124896/1‑239 (MQ=255)                                                                                                                                                                                                          gcctgcGGCCCGCCTAACAACTTGT‑CGCCGGAG‑AAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGC‑CACTGCCTAAATCAGTCACTA‑CGGGAACATCCAGCTCTTTGCCAATCGCCACCAGTTCCGCTTCATC‑TATAGCTTTGGTGAACCCCCGAACGCTGTCATTACTGGAATGGACTTtcat                          >  1:345223/1‑239 (MQ=255)                                                                                                                                                                                                          gcctgcGGCCCGCCTAACAACTTGT‑CGCCGGAG‑AAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGC‑CACTGCCTAAATCAGTCACTA‑CGGGAACATCCACCCCTTTGCCCGCCGCCACCCGTTCCGCTTCACC‑TTTCCCTTGGGTGAACCCCTGAATGCGTCATTTACTGGGAttgcctttact                          >  1:903131/1‑229 (MQ=255)                                                                                                                                                                                                          gcctgcGGCCCGCCTAACAACTTGT‑CGCCGGAG‑AAACTAACCCGACTGACGCCCGCCGCAATAAACTCCTGCGGCATTGGCTCTTTCGGCAAAACGTACTGGCTAAGATCGACCAGCCAGC‑CACTGCCTAAATCAGTCACTC‑CGGGAACATCCAGCTCTTTGCCGAGCGCCACCCGCTCCCCTT‑ATCTTATCGCCTAGGTGAACCCCCGACTACTTAAGTTCCTGGTCTGTAAttttca                          >  1:876646/1‑236 (MQ=255)                                                                                                                                                                                                                ggCCCGCCTAACAACTTGT‑CGCCGGAG‑AAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGC‑CACTGCCTAAATCAGTCACTA‑CGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATC‑TATCGCTTTGGGGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCaa                       >  2:650883/1‑236 (MQ=255)                                                                                                                                                                                                                   ccGCCTAACAACTTCT‑CGGCGGAA‑AACCTCACGAGACTGACGCCCGACGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGC‑CACTGCCTAAATCAGTCACTA‑CGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATC‑TATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGc                 <  2:942544/239‑1 (MQ=255)                                                                                                                                                                                                                                 gT‑CGCCGGAG‑AAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGC‑CACTGCCTAAATCAGTCACTA‑CGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATC‑TATCGCTTTGGTGAAACCCTGAATGCTTTCGTTCATGGAATGTACTTTCTTAAACAGGGCGGGTTTttcattc    >  2:449394/1‑238 (MQ=255)                                                                                                                                                                                                                                    cGCCGGAG‑AAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGC‑CACTGCCTAAATCAGTCACTA‑CGGGAACATCCAGCTCTTTGCCCAGCGCCACCAGTTCCGCTTCATC‑TATCGCTTTGGTGAACCCCTGAATCCCGTAGTTACTGGTATGTACTTTCATCAACCGTGCGGCATTTTCATTCAc  >  2:892733/1‑238 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                     GCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGT‑AAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGT‑AATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACA‑ATAATTCCTGCCTGCGGCCCGCCTAACAACTTGT‑CGCCGGAG‑AAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGC‑CACTGCCTAAATCAGTCACTA‑CGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATC‑TATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTTTCATTCAC  >  NC_000913/3760146‑3760606

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