Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,760,384	A→C	26.7%	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	26.7%	125.6 / 17.4	78	V289G (GTG→GGG)	selA	selenocysteine synthase
Reads supporting (aligned to +/- strand):  ref base A (37/19);  new base C (2/20);  total (39/39)
Fisher's exact test for biased strand distribution p-value = 7.56e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTC‑ATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTG‑CCTGCGGCCCGCC‑TAACAA‑CTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTTTCATTC  >  NC_000913/3760147‑3760605                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                              gCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTC‑ATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGTCCATCTCTTTTTTTCTAACACTTCTTTCTG‑ACTGCGGCCCGCC‑CCACAA‑CCTCTCGCCCGAGAAAATCa                                                                                                                                                                                                                               >  2:931746/1‑238 (MQ=255)   cTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTC‑ATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCGGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTG‑CCTGCGGCCCGCC‑TAACAA‑CTTGTCGCCGGAGAAACTCACCa                                                                                                                                                                                                                            <  2:1035558/239‑1 (MQ=255)      gATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTC‑ATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCCTCTCTTTTTTACCCACCATAATTCCTT‑CCCGCGCCCCGCG‑CACCTC‑CTTCTCTCCGCAGAAACTCAcccca                                                                                                                                                                                                                          >  2:7356/1‑235 (MQ=255)       aTTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAA‑TTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTC‑ATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGTCGGTCTGTCATCTCTTTTTTATCCACAATAATTCCCG‑CCTGCGGCCCGCC‑TAAACA‑CTTTTCTCCTCAGAAACTTACCAGACt                                                                                                                                                                                                                        >  1:1226360/1‑238 (MQ=255)                aCCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTC‑ATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGGTCATCTCTTTTTTACCAACAATAATTCCTG‑CCTGCGGCCCGCC‑TAACAA‑CTTGTCGCCGGAGAAAATCACCAGACTGGagcccg                                                                                                                                                                                                                >  2:251531/1‑238 (MQ=255)                aCCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTC‑ATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTTCTG‑CCTGCGGCCCGCC‑TAACAA‑CTTGTCGCCGGAGAAACTCACCAGACTGAcgcccg                                                                                                                                                                                                                >  2:506297/1‑238 (MQ=255)                     tgcgctgcgGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTC‑ATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTG‑CCTGCGGCCCGCC‑TAACAA‑CTTGTCGCCGGAGAAACTCACCAGACTGACGCCCCCCGCaa                                                                                                                                                                                                          >  1:672520/1‑239 (MQ=255)                       cgctgcgGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTC‑ATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGGTCATCTCTTTTTTACCAACCATAATTCCTT‑GCAGCGGGCCCCC‑TAACAA‑ATTTTCGCCGGAGAAACTCACCAGACTGGCCCCCCCCCCaaac                                                                                                                                                                                                        >  2:472268/1‑237 (MQ=255)                            cgGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTC‑ATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTG‑CCTGCGGCCCGCC‑TAACAA‑CTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACt                                                                                                                                                                                                    >  1:427091/1‑238 (MQ=255)                               gTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTC‑ATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGGTCATCTCTTTTTTTCCAACCATAATTCCTG‑CCTGCGGCCCGCC‑TAACAA‑CCTGTCGCCGGAGAAACTCACCCGCCTGACGCCcgccgccacccacccct                                                                                                                                                                                                 >  2:1221413/1‑227 (MQ=255)                                    cAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTC‑ATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTG‑CCTGCGGCCCGCC‑TAACAA‑CTTGTTGCCGGAGAAAATCCCCAGACTGACCCCCGCCGCAATCAACTCCTGCGGCa                                                                                                                                                                                           >  1:781847/1‑239 (MQ=255)                                         gcAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTC‑ATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTG‑CCTGCGGCCCGCC‑TAACAA‑CTTGTCGCCGGAGAAACTCACCAGACTTACGCCCGCCGCAATCAACTCCTGCGGCATTg                                                                                                                                                                                        >  1:274485/1‑237 (MQ=255)                                                     ttttttCACGCAGAGCTTAAGGGTGTAAATAAGGCCGCACCGGGGCTCCCAGCGCGGCG‑GGGTCTTTTTTCTCCGCGCGTAAGGCCCGCTCCGGCGGGTGCTTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATATTTCCTC‑CCTGCGGCCCCCC‑TAACAC‑CTTGCCGCCGGAGAAACTCCCCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCaaa                                                                                                                                                                          <  2:398189/238‑1 (MQ=255)                                                     aTTTTTCACTCAGAGCTTCAGGGTGTAAAAAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTC‑ATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATCCCTG‑CCTGCGGCCCGCC‑TAACAA‑CTTGTCGCCGGAGAAACTCCCCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCaaa                                                                                                                                                                          <  1:428420/239‑1 (MQ=255)                                                       ttttCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTC‑ATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGGTCATCTCTTTTTTACCAACAATAATTCCTG‑CCTGCGGCCCGCC‑TAACAA‑CTTGTCGCCGGAGAAACTCACCCGACTGACGCCCGCCGCCATCCACTCCTGCGGCATTGGCCCTTTCGGCAAAc                                                                                                                                                                         >  2:428420/1‑238 (MQ=255)                                                         ttCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTC‑ATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTG‑CCTGCGGCCCGCC‑TAAAAA‑CTTTTCGCCGGAGAAACTCACCAGACTGACGCCCGCCCCAATCAACTCCTGCGGCCTTGGCTCTTTCGGCAAACCGt                                                                                                                                                                      >  1:884711/1‑239 (MQ=255)                                                          tCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTCCCAGCGCCGCGAGGGTC‑ATTTTACCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTG‑CCTGCGGCCCGCC‑TAACAA‑CTTGTCGCCGGAGAAACTCCCCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTa                                                                                                                                                                     <  1:1184920/239‑1 (MQ=255)                                                          tCACGCAGGGCTTCAGGGTGTAAATAAAGCCGCAGGGTGGCTTCCAGCGCCGCGGGGGTC‑ATTTTACCCGCGCGTAATGCCCGCTTCAGGGGGTGGTTTTGCAGGGGGGCGATCATCTCTTTTTTACCAACAATAATCCCTC‑CCTGCGGCCCGCC‑TAAAAA‑CTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTa                                                                                                                                                                     <  2:164954/239‑1 (MQ=255)                                                          tCACGCAGAGCTCAGGGGTTTAAATAAAGACGCAACGTGGTTTCCAGCGCCGGGAGGGTC‑ATTTTAGCCGGACTTAATGCCCGCTTGAGCGGGTGGCTTTGCGGGCGGGCGATCATTTCTTTTTTACCAACAAAAATCCCGG‑CCTGCGCCCCGCC‑TAACAA‑CTTGTCGCCGGAGAAACTCCCCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTa                                                                                                                                                                     <  2:756627/239‑1 (MQ=255)                                                          ccACTCAGAGCTTCGGGGTAGAAAAAAAGACGCAACTGGGCCCCCAGCGCGGAGGGGGTT‑ATTTTATCAGCACTTCAGGCACGCTACTGCGGGTGGTTTGGCAGGCGTGATAACATCTCTTTTTTACCAAAAAAAATCCCTG‑CCCGCGCCCCACC‑TAACAC‑CTTGTCGCCGGATAACCCCACCAGACGGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTa                                                                                                                                                                     <  2:339980/238‑1 (MQ=255)                                                                   gtttCAGGTTGTAAAAAAAGACCCAACGTGGCTTCCAGCGCCGCGAGGGTC‑ATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTG‑CCTGCGGCCCGCC‑TAACAA‑CTTGCCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTa                                                                                                                                                              <  1:1103395/235‑1 (MQ=255)                                                                        aGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTC‑ATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTG‑CCTGCGGCCCGCC‑TAACAA‑CTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCCAAGATCg                                                                                                                                                        >  1:281930/1‑238 (MQ=255)                                                                            tgtAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTC‑ATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTG‑CCTGCGGGCCGCC‑TAACAA‑ATTGTCGCCCGAGAAACTCCCCAGACTTACGCCCGCCGCCATCAACTCCTGCGGCCTTGGCTCTTTCGGCACCCCGTACTGGCTAAGAACgcaca                                                                                                                                                    >  2:735659/1‑234 (MQ=255)                                                                            tgtAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTC‑ATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTG‑CCTGCGGCCCGCC‑TAACAA‑CTTTTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCa                                                                                                                                                    >  1:412842/1‑238 (MQ=255)                                                                            tgtAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTC‑ATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTG‑CCTGCGGCCCGCC‑TAACAA‑CTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCa                                                                                                                                                    >  1:336286/1‑238 (MQ=255)                                                                            tgtAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTC‑ATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTG‑CCTGCGGCCCGCC‑TAACAA‑CTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCa                                                                                                                                                    >  2:124621/1‑238 (MQ=255)                                                                            tgtAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTC‑ATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTG‑CCTGCGGCCCGCC‑TAACAA‑CTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCa                                                                                                                                                    >  1:300629/1‑238 (MQ=255)                                                                             gaaaataaaGAGCCAACGGCGCTCCCGGCGCGGCGTGTGTC‑TTTTTCCCCCCCCGTAACGCACGATTCAGCGGTTGCCTTGGCGGGCGGGCGATCATCTTTTTTTTACAAAAAATATTCCCTC‑CCTGCGGCCCGCC‑AAACAA‑TGTGTCGCCGGAGAAACTCACCAGCATGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGc                                                                                                                                                  <  2:781847/237‑1 (MQ=255)                                                                             aaaaataaaGCGGCACGGGGGCTCACGGCGCGGCGGGGGTC‑TTTTTCCCCGCGCTAAATGCCCGTTCCGGCGGGGGGTTTGGCGGGCGGGCTATCATCTTTTTTTTACAAAAAATAATCCCTG‑CCGGCGGCCCGCC‑AAACAC‑CTTGTCCCCGAAGAACCCCCCCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGc                                                                                                                                                  <  2:70288/237‑1 (MQ=255)                                                                              taaataaaGACGCAACGTGGCTTCCAGCGCCGCGAGGGTC‑ATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTG‑GCTGCGGCCCACC‑TAACAA‑CTTGTCGCCGGAGAAACTCACCAGACTGAAGCCCGCCGCAATCAACTCCTGCGGCATTGGACCTTTCGGCAAACCGTAATGGCTAAGATCGCCCAGCg                                                                                                                                                 >  1:720619/1‑239 (MQ=255)                                                                                   aaaGACGCAACGTGGCTTCCAGCGCCGGGAGGGTC‑TTTTTACCCGCGCGTAATGCACGCTTCGGCGGGGGGCTTTGCGGGCGGGCTATCATCTCTTTTTTACCAACAATAATTCCTG‑CCTGCGGCCCGCC‑TAACAA‑CTTGTCGCCGGAGAACCTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCa                                                                                                                                            <  1:1192754/239‑1 (MQ=255)                                                                                     aGACGCAACGTGGCTTCCAGCGCCGCGAGGGTC‑ATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTG‑CCTGCGGCCCGCC‑CAACAA‑CTTTTCGCCGGAAAAACTCACCCGACTGACCCCCGCCCCAAACAACTCCTGCGGGCTTGGCCCTTTCCGCAAACCCTACCGGCCAAAATCGccccacgcgcccc                                                                                                                                           >  2:10432/1‑228 (MQ=255)                                                                                         gCAACGGGGCTCCCAGCGCCGCGGGGGTC‑ATTTTACCCGCGCGTAATCCACGCTTCAGCGGGGGGCTTTGCGGGCGGCCGATCATCTCTTTTTACCCAAAAATAACCCCTG‑CCGGCGCCCCGCC‑AAACAC‑CTTGCCGCGGGAGAACCTCCCCAGACGGACCCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCt                                                                                                                                      <  2:427091/239‑1 (MQ=255)                                                                                             cGTGGCTTCCAGCGCCGCGAGGGTC‑ATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGGTCATCTCTTTTTTACCAACAATAATTCCTG‑CCTGCGGCCCGCC‑TAACAA‑CTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTaaa                                                                                                                                   >  1:86682/1‑238 (MQ=255)                                                                                                  ctcccggcgcGGCGGGGGTG‑TTTTTACCCGCGCTTACTCCCCGCTTCAGCGGGTGGTTTGGCGGGGGGGCAATCATCTTTTTTTCCCAAACAAAATTCCCTC‑CCGGCGCCCCCCC‑AAACAT‑CTTGCCCCCGGAGAAACTCACCAGACGGCCGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGt                                                                                                                              <  2:1192601/232‑1 (MQ=255)                                                                                                      cggcgcCGGGGGGTTC‑TTTTTACCCGCGCGAAATGCACTCTTCGGCGGGTGGCTTGCCGGGGGGGCGATCATCTCTTTTTCACAAACAATAATCCCCG‑CCGGCGCCCCGCC‑AAACAC‑CTTTTCGGCGGAGAAACCCACCAGCCTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTa                                                                                                                         <  2:1204801/237‑1 (MQ=255)                                                                                                         cgccgcGAGGGTC‑ATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTG‑CCTGCGGCCCGCC‑TAACAA‑CTTGTCGCCGGAGAAACTCACCAGACTGACGCccgccg                                                                                                                                                                                                             >  2:358871/1‑152 (MQ=255)                                                                                                         cgccgcGAGGGTC‑ATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAAAAATTCCTG‑CCTGCGGCCCGCC‑TAACAA‑CTTGTCGCCGGAGAAACTCACCAGACTGACGCccgccg                                                                                                                                                                                                             <  1:358871/152‑1 (MQ=255)                                                                                                           ccggGAGGTTC‑ATTTTATCCGCGCGTAATGCCCGCTTCGGCGGGTGGCTTGGCGGGGGGGCGATCATCTCTTTTTTACCAACAATAATCCCTG‑CCGGCGGCCCGCC‑TAACAA‑CTTGTCGCCGGAGAACCTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGa                                                                                                                    <  1:1085299/239‑1 (MQ=255)                                                                                                           ccgcGAGGGTC‑ATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTG‑CCTGCGGCCCGCC‑TAACAA‑CTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCCAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGa                                                                                                                    >  1:216650/1‑239 (MQ=255)                                                                                                            cgggAGGGTT‑TTTTTACCCGCGGTTAATGCACTCTTCGGCGGGGGCCTTTGCAGGGGGGGGATCATCTCTTTTTTACAAACAATAATCCCTG‑CCGGCGCCCCGCC‑TAACAA‑CTTGCCGCCGGAGAACCCCCCCAGACTGCCGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGa                                                                                                                    <  1:10432/235‑1 (MQ=255)                                                                                                            cgagggggTC‑TTTTTACCCGCGCGAAATGCACGTTTCCGCGTGGGGCTTTGCAGGCGGCCGCTCATTTGTTTTTTCCAAACAATATTCCCTG‑CCGCCGGCCCGCC‑‑AACAACCTTGGCGCCGGAAAAACTCCCAAAACGAACCCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGa                                                                                                                    <  2:86682/233‑1 (MQ=255)                                                                                                                              ccGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTG‑CCTGCGGCCCGCC‑TAACAA‑CTTGTCGCCGGAGAAACTCCCCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCg                                                                                                  <  1:124621/239‑1 (MQ=255)                                                                                                                                 cgcgTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTG‑CCTGCGGCCCGCC‑TAACAA‑CTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAgc                                                                                               >  1:993410/1‑239 (MQ=255)                                                                                                                                 cgcgTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTG‑CCTGCGGCCCGCC‑TAACAA‑CTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAgc                                                                                               >  1:474911/1‑239 (MQ=255)                                                                                                                                   cgTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTG‑CCTGCGGCCCGCC‑TAACAA‑CTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTTGCTCTTTCGGCAAACCGTACTGGCTTAGATCTACCCGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAgcg                                                                                              >  1:654084/1‑238 (MQ=255)                                                                                                                                     tAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTG‑CCGGCGGCCCGCC‑TAACAA‑CTTGCCGCCGGAGAAACCCCCCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGcca                                                                                           <  1:1221405/239‑1 (MQ=255)                                                                                                                                        tGCCCGCTTCGGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTCCCAACAATAATTCCTG‑CCTGCGGCCCGCC‑TAACAC‑CTTGCCCCCGGAAAACCTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGccacca                                                                                        <  1:506297/239‑1 (MQ=255)                                                                                                                                           gcgtgccggCGGGTGGCTTGGCGGGCGGGCTATCTTCTCTTTTTAACAAACAAGAATCCCTC‑CCGGCGCCCCGCC‑AAAAAC‑CTGGTCGCCGAAAAAACCCCCCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGt                                                                                      <  2:101493/230‑1 (MQ=255)                                                                                                                                                         ggCTTTGCAGGCGGGCGATCATCTCTTTTTTACCCACAATAATTCCTG‑CCTGCGGCCCGCC‑TAACAA‑CTTGTCGCCGGAGAAAATCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGGTAAGATCGACCAGCGGGGCAATGCCTAAATCAGTCACTACCGGGACATCCAGCTCTTTGCCCAGCGCCCCCAGTTCCGCtttctccat                                                                        >  2:527232/1‑231 (MQ=255)                                                                                                                                                         ggCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTG‑CCTGCGGCCCGCC‑TAACAA‑CTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTAt                                                                        >  1:816558/1‑238 (MQ=255)                                                                                                                                                         ggCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTG‑CCTGCGGCCCGCC‑TAACAA‑CTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAAATCGACCAGCCAGCCACCGCCTAAATCAGTCACTACGGGAACATCCCGCCCTTTGCCCGGCGCCCCCACTTCCGCTTCATCTAt                                                                        >  2:1102773/1‑238 (MQ=255)                                                                                                                                                                   gcgggcgTTCATCTC‑TTTTTACCAACAAAAATTCCCGCCCGGGGGCCCCCC‑TAAAAA‑CTTGCCCCCGGAAAAACCCCCCAGACTGCCGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTg                                                              <  2:816558/238‑1 (MQ=255)                                                                                                                                                                     gggcgATCATTTTTTTTTTACCAACAATAATTCCTG‑CCTGCGGCCCGCC‑TAACAA‑CTTGTCGCCGGAGAAACTCACCAGCCTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGaa                                                            <  1:904694/238‑1 (MQ=255)                                                                                                                                                                         catcatcGTTTTTCCACCAACAA‑AACTCCCG‑CCCGCCGCCCCCCATAACAC‑CTTCTGCCCGGAAAACCTCCCAAGACTGACCCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCt                                                       <  2:1097620/238‑1 (MQ=255)                                                                                                                                                                            catcTCTTTTTTACCAACAATAATTCCTG‑CCTGCGGCCCGCT‑TAACAA‑CTTGTCGCCGGAGAACCTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGaa                                                    <  1:1142723/239‑1 (MQ=255)                                                                                                                                                                            aagctctTTTTACCAAAAAATAATACCTC‑CCGGGCGCCCGCA‑AAACAC‑CGGCGCGCGGGAGAACCTCCCCGGACGGCCCCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGaa                                                    <  2:281930/236‑1 (MQ=255)                                                                                                                                                                                  ttttttCCCAACAATAATCCCTG‑CCGGCGGCCCGCC‑TAACAA‑CTTGTCGCCGGAGAAACCCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGt                                              <  1:1221413/239‑1 (MQ=255)                                                                                                                                                                                     ttaaCCAACAATAATTCCTG‑CCGGCGGCCCGCC‑TAACAA‑CTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGt                                           <  1:251531/236‑1 (MQ=255)                                                                                                                                                                                       tcccAACAATAATCCCTC‑CCGGCGGCCCCCC‑TAACAA‑CTTGTCCCCGGAGAAACTCCCCAGACTGCCGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTa                                         <  1:173678/237‑1 (MQ=255)                                                                                                                                                                                       tACCAACAATACTCCCTC‑CCTGCGGCCCGCC‑TAAAAA‑CTTGTCGCCGAAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTa                                         <  1:142117/239‑1 (MQ=255)                                                                                                                                                                                       aaccacaaataatTCCTG‑CCGGCGGCCCGCC‑AAAAAA‑CTTGCCGCGGGAGAACCCCCCCAGACTGACCCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTa                                         <  1:735659/232‑1 (MQ=255)                                                                                                                                                                                         ccAACAATAATTCCTG‑CCTGCGGCCCGCC‑TAACAA‑CTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGtt                                          >  1:854232/1‑236 (MQ=255)                                                                                                                                                                                            acaaTAATTCCTG‑CCTGCGGCCCGCC‑TAACAA‑CTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGt                                    >  1:515742/1‑239 (MQ=255)                                                                                                                                                                                                 attccctc‑ccggCGCCCCCCC‑AACCAC‑CTTGCCGCCGGAGAAACTCACCGGCCTGCCCCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTa                               <  1:1074605/235‑1 (MQ=255)                                                                                                                                                                                                 atcccccg‑ccggcGCCCCCCC‑AAACAT‑CTTGCCGCCGAAAAAACCCCCAAGCGGCACCCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTa                               <  2:474911/232‑1 (MQ=255)                                                                                                                                                                                                 aatCCCTG‑CCGCCGGCCCGCC‑AAACAA‑TTGGTCCCCGGAGAACCTCCCCAGCCTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTa                               <  1:447008/239‑1 (MQ=255)                                                                                                                                                                                                     cctg‑cctgCGGCCCGCC‑TAACAA‑CTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACtt                            >  2:1195983/1‑238 (MQ=255)                                                                                                                                                                                                       tg‑cctgcGGCCCGCC‑TAACAA‑CTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGaa                                                    >  1:439955/1‑212 (MQ=255)                                                                                                                                                                                                        g‑cctgcGGCCCGCC‑TAACAA‑CTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTTTGTACTTtca                         >  2:709488/1‑238 (MQ=255)                                                                                                                                                                                                        g‑cctgcGGCCCGCC‑TAACAA‑CTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTtca                         >  1:908662/1‑238 (MQ=255)                                                                                                                                                                                                        g‑cctgcGGCCCGCC‑TAACAA‑CTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTtca                         >  2:938292/1‑238 (MQ=255)                                                                                                                                                                                                            cgcGGCACGCC‑AAACAA‑TTTGCCCCGGGAAAACCCCCCCAGCGCGCCGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCaa                     <  2:672520/238‑1 (MQ=255)                                                                                                                                                                                                                 cccGCT‑TAACAT‑CTTGTCGCGGTAGAAACTCACCAGCCTGAAGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTg                <  1:775349/239‑1 (MQ=255)                                                                                                                                                                                                                  ccacc‑TAACAA‑CTGGTCGCGGGAAAAACCCACCAGACGGAGCCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTg                <  2:888741/235‑1 (MQ=255)                                                                                                                                                                                                                     cccaaacaa‑tTTGGCGCCGGAGAAACTCCCCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCgg             <  1:414956/235‑1 (MQ=255)                                                                                                                                                                                                                     cc‑TAACAA‑CTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGt            >  2:582547/1‑239 (MQ=255)                                                                                                                                                                                                                                  gTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTttcattc  >  2:783140/1‑238 (MQ=255)                                                                                                                                                                                                                                  gTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTATAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTttcattc  >  2:949435/1‑238 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                              GCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTC‑ATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTG‑CCTGCGGCCCGCC‑TAACAA‑CTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTTTCATTC  >  NC_000913/3760147‑3760605

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