Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,956,003	G→A	D225N (GAT→AAT)	ilvA →	threonine deaminase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,956,003	0	G	A	95.0%	227.3 / NA	82	D225N (GAT→AAT)	ilvA	threonine deaminase
Reads supporting (aligned to +/- strand):  ref base G (4/0);  new base A (27/49);  total (32/50)
Fisher's exact test for biased strand distribution p-value = 1.99e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCGATCGAACTGTCACAGCAGCAGGGGTTCACCTGGGTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTG‑CCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGA‑TGCCGCAAATCAAAGTG‑ATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGGA‑TGCGGG‑TCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACG‑AAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATGAAAAAATATATCGCCCTGCACAACATTCGC  >  NC_000913/3955753‑3956248                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                         aGCGATCGAACTGTCACAGCAGCAGGGGTTCACCTGGGTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTG‑CCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGGTCCAACCACTGA‑TGCCGCAAATCAAAGTG‑ATCGCCGTAGAAGAGGGAGACTCCCCCTTCCTGAAAACCGCGCTgg                                                                                                                                                                                                                                                          >  1:536292/1‑251 (MQ=255)          aCTGTCACAGCAGCAGGGGTTCACCTGGGTGCCGCCGTTCGACCATCCGAGGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTG‑CCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGA‑TGCCGCAAATCAAAGTG‑ATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAA‑TGCGGG‑Tc                                                                                                                                                                                                                                               <  1:237051/251‑1 (MQ=255)                       cgggggTTCACCGGGGTCCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTG‑CCAGTCGGCGGCGGCGGTCTGGCTGCGGGCGTGGCGGTGCTGATCAAACAACTGA‑TGCCGCAAATCAAAGTG‑ATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAA‑TGCGGG‑TCATCCGGTTGATCt                                                                                                                                                                                                                                  <  2:223089/249‑1 (MQ=255)                         ggggttccacgtggggggtccgTTCGGACATCCGAGGGTGATTGCCGGGCAAGGCGCGCTGGCGCGGGAACTGCGCCAGCAGGACGCCCATCGCGACCGCGTATTTGTG‑CCAGACGGAGGCGGCGGTCTTGCTGGTGGCGTGGCGGTGCTGATCAAACAACTGA‑AGTCGCAAAGCAAAGTG‑ATCGCCGGAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAA‑TTCGGG‑GCATCCGGTTGATCTGc                                                                                                                                                                                                                                <  2:335912/232‑1 (MQ=255)                               cacggggggtccgccgTGCGACCAGCCGATGGGGATTGCCGGGCAAGGCTCGCTGGCTCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTG‑CCAGTCGGCGGCGGCGTTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGA‑TGCCGCAAATCAAAGTG‑ATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAA‑TGCGGG‑TCATCCGGTTGATCTGCcgcgcg                                                                                                                                                                                                                          <  2:137394/241‑1 (MQ=255)                                          cgccgTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTG‑CCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGA‑TGCCGCAAATCAAAGTG‑ATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAA‑TGCGGG‑TCATCCGGTGGATCTGCCGCGCGTAGGGCTAttt                                                                                                                                                                                                               <  1:489633/251‑1 (MQ=255)                                          cgccgTGCGTCCACCCGAGGGTGATTGTCGGGCAAGGCTCGCTGGCGATGGAACTGCTCCCGCAGGACGCCCATCACGCCCGCGTATTTGTG‑CCAGTCGGAGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGCTCAAACAACTGA‑TGCCGCAAATAAAAGTG‑ATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAA‑TGCGGG‑TCATCCGGTTGATCTGCCGCGCGTAGGGCTAttt                                                                                                                                                                                                               <  2:426675/251‑1 (MQ=255)                                              gTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTG‑CCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGGTCCAACAACTGA‑TGCCGCCAATCAAAGGGGAT‑GCCGTAGAAGCGGAAAACTCCCCCTCCCTGAAAGCCGCGCTGAA‑TGCGGG‑TTCTCCGGGTGGTCTGCcgcgcgggggggttttttctt                                                                                                                                                                                                           >  1:279401/1‑236 (MQ=255)                                                   aCCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTG‑CCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGA‑TGCCGCAAATCAAAGTG‑ATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAA‑TGCGGG‑TCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGc                                                                                                                                                                                                      <  1:18627/251‑1 (MQ=255)                                                      agccgaggtggATTGCAGGGCAAGGCCCGTGGGCGCTGGAACGGCGCCAGCAGGACGCCCATCTCGACCGCGTATGTGTG‑CCAGTCGGCGGCGGCGGTCTGGCGGCTGGCGTGGCGGTGCTGATCAAACAACGGA‑TGCCGCAAATCAAAGTG‑ATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAA‑TGCGGG‑TCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCGGAAGGCGTa                                                                                                                                                                                                   <  2:160647/241‑1 (MQ=255)                                                          gATGTTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTG‑CCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGA‑TGCCGCAAATCAAAGTG‑ATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAA‑TGCGGG‑TCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCgg                                                                                                                                                                                               <  1:496640/251‑1 (MQ=255)                                                                    gCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTG‑CCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGA‑TGCCGCAAATCAAAGTG‑ATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAA‑TGCGGG‑TCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCa                                                                                                                                                                                      <  1:23973/250‑1 (MQ=255)                                                                     ccGGGCAAGGCCCGCTGGCGCTGGAACTGCCCCAGCAGGACGCCCACCGCGACCGCGTATTGGGG‑CCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGA‑TGCCGCAAATCAAAGTG‑ATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAA‑TGCGGG‑TCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATc                                                                                                                                                                                    <  2:326785/251‑1 (MQ=255)                                                                     ccGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTG‑CCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGA‑TGCCGCAAATCAAAGTG‑ATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAA‑TGCGGG‑TCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATc                                                                                                                                                                                    <  1:127616/251‑1 (MQ=255)                                                                       gggCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTG‑CCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTTCTGATCAAACAACTGA‑TGCCGCAAATCCAAGTG‑ATCGCCGTAGAAGCGGGAGACTCCGCCTGCCTGAAAGCAACGGTGAA‑TTCGGG‑TCATCCGGTTGGTCTGGCGCCCGGAGGGCACATTGCGGAAGGGGTTGCGGTCAAACGCAACgg                                                                                                                                                                                  >  2:87925/1‑251 (MQ=255)                                                                          cAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCCTCTAGACCGCGTATTTGTG‑CCAGTCGGCGGAGGAGGTCTGGCTGCTGGCGTGGCGGTGGAGGTAAAACACATGA‑AGCCGCAAAACCAAGGG‑CACGCCGGAAAATCGGGAAACTCCCCCTGCCTGAGAGAAAGGGTAGA‑AGCGGG‑GCAACCGGGTGATCTACCGCGAGTAGGGATATTTGCGGAAggaggggcgggaaacagaaacggggg                                                                                                                                                                               >  2:339289/1‑227 (MQ=255)                                                                           aaGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTG‑CCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGA‑TGCCGCAAATCAAAGTG‑ATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAA‑TGCGGG‑TCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGAc                                                                                                                                                                              >  1:462025/1‑251 (MQ=255)                                                                             ggCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTG‑CCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGGTCAAACAACTGG‑AGCCGCAAATCAAAGGG‑ATCGCCGTAGAAGCGGAAAACTCCGCCTGCCTGAAAGCAGCGC‑GTACAGCGGG‑GCAACCGGTTGATCTGCCGCGCGGAGGGCTAATTGTTGAAAGCGTAGCGGTAAAAAGCATCGGTGACg                                                                                                                                                                             >  1:372090/1‑250 (MQ=255)                                                                                  gctggcgctggAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTG‑CCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGGTCAAACAACTGA‑TGCCGCAAATCAAAGTG‑ATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGGAAGCAGCGCTGAA‑TGCGGG‑TCATCCGGTTGATCTGCCGCGCGGTGGGCTATTTGCTGAAGGCGTAGCGGGAAAAAGCCACGGTGAAG‑aaaacc                                                                                                                                                                      >  1:566319/1‑248 (MQ=255)                                                                                     ggcgctggAACTGCCCCAGCAGGACGCCCATCGCGACCGCGTATTTGGG‑CCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGA‑TGCCGCAAATCAAAGTG‑ATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAA‑TGCGGG‑TCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACG‑AAACCTTc                                                                                                                                                                    <  2:152185/250‑1 (MQ=255)                                                                                       cgctggAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTG‑CCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGCGCTGCTCAAACAACTGA‑TGCCGCAAATCAAAGTG‑ATCGCCGTAGAAGAGGAAGACTCCGCCCGCCTGCAAGCAGCGCCGAA‑TGCGGG‑TCATCCGGGTGATCTGCCGCGCGGAGGGCTATTTGCTGAAGGCGTAGCGGTAAAAAGCATCGGGGACG‑AAACCTTCCGt                                                                                                                                                                 >  1:253002/1‑251 (MQ=255)                                                                                        gctggAACTGCCCCAGCGGGACGCCCAGCTCGACTGCGGTTTGGCG‑CCAGCCGGCGTCGGCGGTCTGGCTGATTGCGTGGCGGCGAGGAGTAAACAACTGA‑TGCCGCAAACCAAAGTG‑ATCGCCGTATAAGCGGAAGACTCCGCCTGCCGTAAAGCAGCGTTGAA‑TGCGGG‑GCATCCGGTGGATCTGCCGTGCGTAGGGCTATTTGCTGAAGGCGGAGCGGTAAAACGCATCGGTGACG‑AAACCTTCCGtt                                                                                                                                                                <  2:332670/251‑1 (MQ=255)                                                                                              aCTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTG‑CCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGTTCACACAACTGG‑TGCCGCAAATCAAAGTG‑ATCGCCCGAGGAGGGGAAGGCCCCACCTTCCTGGAAGCAGCGCTGAG‑TGCGGG‑TCATCCGGTTTATCTGCCGCGCGGAGGGCTATTTTCTGAAAGCGTAGCGGTAAAACGCATAGGCGGCG‑AAACCtttcttgtttgcc                                                                                                                                                          >  2:538290/1‑240 (MQ=255)                                                                                                tGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTG‑CCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTTCTTGTCAAACAACTGT‑TACCGCAAATCAAAGGG‑AACGCCGTAGAAGCGGAAGACACCACCCGCCCGAAAACGGCGCTGAA‑TGCGGG‑TCATCCGGGTGATCTGCCCCGCGGAGGGGTATTTGCTGGAGGCGAGGCGGTAAAAAGCATCGGTGAGG‑AAACCTTCCGTTTATGcc                                                                                                                                                          >  2:431702/1‑249 (MQ=255)                                                                                                    ccAGCAGGACGCCCATCTCGACCGCGTATTTGTG‑CCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGA‑TGCCGCAAATCAAAGTG‑ATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAA‑TGCGGG‑TCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACG‑AAACCTTCCGTTTATGCCAGGa                                                                                                                                                      <  1:439740/249‑1 (MQ=255)                                                                                                     cagcagGACGCCCATCTCGACCGCGTATTTGTG‑CCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGGGCTGATCAAACAACTGA‑AGCCGCAAATCAAAGTG‑GTCGCCGTAGAAGCGGGAGACTCCGCCTGCCTGAAAGCAGCGCTGAA‑TGCGGG‑TCATCCGGGTGATCTTCCGCGCGGAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGGGACG‑AAACCTTCCGTTTATGGCAGGAGt                                                                                                                                                    >  1:444874/1‑250 (MQ=255)                                                                                                        cagGACGCCCATCTCGACCGCGTATTTGTG‑CCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGGCCTGATCCAACCACTTAGTGC‑GAAAATAAAAGAG‑ATCGCCGTAAAAACGGGAGACACCGCCTGCCTGAAAGCAACGCTGGA‑TGCGGG‑GCATCCCGTTGATATGGCGCGCGTCGGGCGATTTTCTGCAGGCGTTGCGGTAAAACAGATCGGGGACG‑AAACCTTCCCTTTATGCCAgggggatct                                                                                                                                                >  1:185896/1‑244 (MQ=255)                                                                                                          gtgatgcccttctggagctctttttGTGGCCAGTCGGCGGCGGCTGGCTGGCGGCGGGCGTGGCGGTGCTGATCAAACAACTGA‑TGCCGCAAATCAAAGTG‑ATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAA‑TGCGGG‑TCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTATCGTTAAACCGCATCGGTGACG‑AAACCTTCCGTTTATGCCTGGAGTAtctc                                                                                                                                               <  2:66325/229‑1 (MQ=255)                                                                                                          ggACGCCCATCTCGACCGCGTATTTGTG‑CCAGGCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGA‑TGCCGCAAATCAAAGTG‑ATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAA‑TGCGGG‑TCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACG‑AAACCTTCCGTTTATGCCAGGAGTATCTcg                                                                                                                                              <  1:293421/251‑1 (MQ=255)                                                                                                            gcGCCCATCTCGACCGCGTATTTGTG‑CCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGA‑TGCCGCAAATCAAAGTG‑ATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAA‑TGCGGG‑TCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACG‑AAACCTTCCGTTTATGCCAGGAGTATCTcgac                                                                                                                                            <  1:225183/250‑1 (MQ=255)                                                                                                             cGCCCATCTCGACCGCGTATTTGTG‑CCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGA‑TGCCGCAAATCAAAGTG‑ATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAA‑TGCGGG‑TCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACG‑AAACCTTCCGTTTATGCCAGGAGTATCTcgacg                                                                                                                                           <  1:100055/251‑1 (MQ=255)                                                                                                                    ctcGACCGCGTATTTGTG‑CCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGA‑TGGCGCAAATCAAAGTG‑ATCGCCGGAGAAGCGGGAGACTCCGCCTGCCCGAAAGCAGCGCAGGA‑TGCGGG‑TTATCCGGTTTTTCTGCCGCGCGGGGGGCGATTTGCGGGAGGCGTAGCGGTAAAAAGCATCGGTGGCG‑GAACCTTCTTCTTATGCCAGGGGTATTTCGACcacagaat                                                                                                                                    >  2:152069/1‑247 (MQ=255)                                                                                                                     tcGACCGCGTATTTGTG‑CCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGGTCAAACAACTGA‑TGCCGCAAATCAAAGTG‑ATCGCCGTAGAAGCGGGAGACTCCGCCTGCCTGAAAGCAGCGCTGAA‑TGCGGG‑GCATCCGGTTGATCTTCCGCGCGGAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGGGGCG‑AAACCTTCCGTTTATGCCAGGGGTATCTCgacgacctcctc                                                                                                                                   >  1:514800/1‑245 (MQ=255)                                                                                                                       gACCGCGTATTTGTG‑CCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGGTCAAACAACTTA‑TGCCGCCAATCAAAGTG‑ATCGCCGTAGAAGCGGGAGACTCCGCCTGCCTGAAAGCAGCGCTGAA‑TGCGGG‑TCATCCGGGTGATCTGCCGCGCGGAGGGCTTTTTTCTGAAGGCGTAGCGGGAAAACCCATCGGGGGAG‑AAACCTTCCGTTTATGCCAGGAGGATATCGCCGGCATCATCAc                                                                                                                                 >  1:350140/1‑251 (MQ=255)                                                                                                                       gACCGCGTATTTGGG‑CCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGA‑TGCCGCAAACCAAAGTG‑ATCGCCGTAGAAGCGGAAGACGCCGCCTGCCTGAAAGCAGCGCTGAA‑TGCGGG‑TCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACG‑AAACCTTCCGTGTATGCCAGGAGTATCTCGACGACAGCATCAc                                                                                                                                 <  2:514800/251‑1 (MQ=255)                                                                                                                                     g‑CCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGA‑TGCCGCAAATCAAAGTG‑ATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAA‑TGCGGG‑TCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACG‑AAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAgcgatg                                                                                                                   <  1:87925/251‑1 (MQ=255)                                                                                                                                       ccAGCCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGA‑TGCCGCAAATCAAAGTG‑ATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAA‑TGCGGG‑TCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACG‑AAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAgcgatgc                                                                                                                  <  2:214299/251‑1 (MQ=255)                                                                                                                                          gTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGA‑TGCCGCAAATCAAAGTG‑ATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAA‑TGCGGG‑TCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACG‑AAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAgcgatgcgat                                                                                                               <  1:36728/251‑1 (MQ=255)                                                                                                                                          gTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAAATGA‑TGCCGCAAATCAAAGGG‑ATCGCCGTAGAAGCGGAAAAATCCGCCTGCCTGAAAACAGCGCTGAA‑AGCGGG‑TCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGGAGCGGTAAAACGCATCGGGGACG‑AAACCTTCCGGTTATGCCAGGAGGATCTCCACGACATCATCAACGTCGATAgcggtgcggt                                                                                                               >  1:395029/1‑249 (MQ=255)                                                                                                                                              gcggcggcggTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGA‑TGCCGCAAATCAAAGTG‑ATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAA‑TGCGGG‑TCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACG‑AAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCtgt                                                                                                           <  1:561030/251‑1 (MQ=255)                                                                                                                                               cggcggcggTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGA‑TGCCGCAAATCAAAGTG‑ATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAA‑TGCGGG‑TCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACG‑AAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCtgtg                                                                                                          <  1:287630/251‑1 (MQ=255)                                                                                                                                                ggcggcggTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGA‑TGCCGCAAATCAAAGTG‑ATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAA‑TGCGGG‑TCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACG‑AAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTgc                                                                                                         <  1:496434/251‑1 (MQ=255)                                                                                                                                                ggcggcggTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGA‑TGCCGCAAATCAAAGTG‑ATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAA‑TGCGGG‑TCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACG‑AAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTgc                                                                                                         <  1:112318/251‑1 (MQ=255)                                                                                                                                                ggcggcggTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGA‑TGCCGCAAATCAAAAGG‑ATCGCCGGAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAA‑TGCGGG‑GCATCCGGGTGATCTGCCGCGCGGAGGGCTATTTGCTGAAGGCGGAGCGGTAAAAAGCATCGGGGACG‑AAACCTTCCGTTTATGCCCGGAGTATCTCGACGGCATCATCAACGTCGATAGCGATGCGATCtgggc                                                                                                         >  1:160693/1‑248 (MQ=255)                                                                                                                                                ggcggcggTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAAATGA‑TGCCGCAAATAAAAAGG‑ATCGCCGGAGAAACGGGAAACTCCGCCCGCCTTAAAACAGCGCGGGA‑AGAGGG‑GCATCCGGTTGATCTTCCGCGCGGAGGGCTTATTGCTTAAAGCGTAGCGGGAAAAAGCATCGGTGACGGAAACTTTCCGTTT‑TGCCAGGAGGATCTCGACGACATCATCACCGTCgcaagcgaagacacatttgg                                                                                                         >  1:486741/1‑232 (MQ=255)                                                                                                                                                           ggCTGCTGGCGTGGGGGGGCTGATCAAACAACTGG‑TGCCGCAAAACAAAGTG‑ATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGCA‑TGCGGG‑TCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACG‑AAACCTTCCGTTTATGCCAGGAGTATCTCGACGACAGCATCACCGTCGATATCGATGCGATCTGTGCGGCGATGa                                                                                                 <  2:408155/248‑1 (MQ=255)                                                                                                                                                           ggCTGCTGGCGTGGCGGTGCTGATCAAACAACTGA‑TGCCGCAAATCAAAGTG‑ATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAA‑TGCGGG‑TCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACG‑AAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGaa                                                                                                <  1:12760/249‑1 (MQ=255)                                                                                                                                                               gctgGCGTGGCGGTGCTGATCAAACAACGGA‑TGCCGCAAATCAAAGTG‑ATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAA‑TGCGGG‑TCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACG‑AAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGAt                                                                                            <  2:138048/249‑1 (MQ=255)                                                                                                                                                                      tggcgGTGCTGATCAAACAACTGA‑TGCCGCAAATCAAAGTG‑ATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAA‑TGCGGG‑TCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACG‑AAACCTTCCGTTTATg                                                                                                                                                            <  1:144754/178‑1 (MQ=255)                                                                                                                                                                      tggcgGTGCTGATCAAACAACTGA‑TGCCGCAAATCAAAGGG‑ATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGGAAGCAGAGCTGAC‑TGCGGG‑TCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGCCG‑AAACCTTCCGGTTagg                                                                                                                                                            >  2:144754/1‑176 (MQ=255)                                                                                                                                                                          ggggCTGATAAAACACCTGA‑GGCCGCAAATCAAAGTG‑AACGCCGTAGAAGCGGAAGCCTCCGCCTGCCTGAAAGCAGCGCTGAA‑TGCGGG‑TCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACG‑AAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGAtg                                                                               <  2:536292/248‑1 (MQ=255)                                                                                                                                                                          ggTGCTGATCAAACAACTGA‑TGCCGCAAATCAAAGTG‑ATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAA‑TGCGGG‑TCATCCGGTTGATCTTCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACG‑AAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGAtg                                                                               <  1:418927/251‑1 (MQ=255)                                                                                                                                                                             gCTGATCAAACAACTGA‑TGCCGCAAATCAAAGTG‑ATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAA‑TGCGGG‑TCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACG‑AAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTgc                                                                            <  1:558721/251‑1 (MQ=255)                                                                                                                                                                                gATCAAACAACTGA‑TGCCGCAAATCAAAGTG‑ATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAA‑TGCGGG‑TCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACG‑AAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTgcg                                                                           <  1:446697/249‑1 (MQ=255)                                                                                                                                                                                 aTCAAACAACTGA‑TGCCGCAAATCAAAGGG‑ATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGTAAGCAGCGCTGAA‑TGCGGG‑TCATCCGGTTGATCTGCCGCGCATAGTGCTATTTGCTGAAGGCGGAGCGGTAAAACGCATCGGTGACA‑AAACCTTCCGTTTATACCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCTAAGGTGTgctcgc                                                                        >  1:434671/1‑251 (MQ=255)                                                                                                                                                                                     accaacTTA‑TGCCGCAAATAAAAGGG‑CTCGCCGTAGACGGTGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAA‑TGCGGG‑TCATCCGGTTGATCTGCCGCTCGTAGGGCTATTTGCTGAAGGCATAGCGGTAAAACGCATCGGTGACG‑AAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATTTGCGCGCGGTg                                                                    <  2:568721/249‑1 (MQ=255)                                                                                                                                                                                     aacaacTGA‑TGCCGCAAATCAAAGTG‑ATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAA‑TGCGGG‑TCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACG‑AAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTg                                                                    >  2:118277/1‑251 (MQ=255)                                                                                                                                                                                     aacaacTGA‑TGCCGCAAATCAAAGTG‑AGTGCCGTAGAAGCGGAAGACGCCGCCTGCCTGAAAGCAGCGCTGAA‑TGCGGG‑TCATCCGGTTGAGCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACG‑AAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAATATGTGCGCGCGGTg                                                                    <  2:187170/251‑1 (MQ=255)                                                                                                                                                                                       caacTGA‑TGCCGCAAATCAAAGTG‑TTCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAA‑TGCGGG‑TCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACG‑AAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGAGTTATTCGAAGATGTGCGCGCGGTGGc                                                                  <  2:462025/251‑1 (MQ=255)                                                                                                                                                                                        caaatgaggCCGCAAATCAAAGTG‑CTCGCCGTCGAAGGGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAA‑TGCGGG‑TCATCCGGTTGCTCAGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACG‑AAACATTCAGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCACTGCGATCGGTGCGCAGATGAAGGATTTATTCGAAGTTGTGAGCGCGGTGGc                                                                  >  2:388683/9‑251 (MQ=255)                                                                                                                                                                                                   gCAAATCAAAGTG‑ATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAA‑TGCGGG‑TCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACG‑AAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTg                                                       >  1:198714/1‑251 (MQ=255)                                                                                                                                                                                                   gCAAATCAAAGTG‑ATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAA‑TGCGGG‑TCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACG‑AAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTg                                                       <  1:119275/251‑1 (MQ=255)                                                                                                                                                                                                            aGTG‑ATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGAGCTGAA‑TGCGGG‑TCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACG‑AAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCgcgctgg                                              <  1:475286/251‑1 (MQ=255)                                                                                                                                                                                                                       gTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAA‑TGCGGG‑TCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACG‑AAACCTTCCGTTTATGCCAGGAGTATCTCGACGAcatca                                                                                                                                     >  1:183019/1‑154 (MQ=255)                                                                                                                                                                                                                       gTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAA‑TGCGGG‑TCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACG‑AAACCTTCCGTTTATGCCAGGAGTATCTCGACGAcatca                                                                                                                                     <  2:183019/154‑1 (MQ=255)                                                                                                                                                                                                                       gTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAA‑TGCGGG‑TCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACG‑AAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCggg                                    <  1:468939/251‑1 (MQ=255)                                                                                                                                                                                                                         aGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAA‑TGCGGG‑TCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACG‑AAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGa                                   >  1:140529/1‑250 (MQ=255)                                                                                                                                                                                                                            aGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAA‑TGCGGG‑TCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACG‑AAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATGa                               <  1:152069/251‑1 (MQ=255)                                                                                                                                                                                                                             gCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAA‑TGCGGG‑TCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACG‑AAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATGa                               >  1:296799/1‑250 (MQ=255)                                                                                                                                                                                                                              cGGAAGACTCCGCCTGTCTGAAAGCAGTGCTGAA‑TGCGGG‑TCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACG‑AAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATGaaa                             <  2:395029/251‑1 (MQ=255)                                                                                                                                                                                                                               ggAAGACTCCGCCTGCCTGAAAGCAGCGCTGAA‑TGCGGG‑TCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACG‑AAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATGaaa                             >  1:512249/1‑250 (MQ=255)                                                                                                                                                                                                                                gAAGACTCCGCCTGCCTGAAAGCAGCGCTGAA‑TGCGGG‑TCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACG‑AAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATGaaaa                            >  1:439004/1‑250 (MQ=255)                                                                                                                                                                                                                                       ccGCCTGCCTGAAAGCAGCGCTGAA‑TGCGGG‑TCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACG‑AAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATGAAAAAATATATc                    <  1:225006/251‑1 (MQ=255)                                                                                                                                                                                                                                       ccGCCTGCCTGAAAGCAGCGCTGAA‑TGCGGG‑TCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACG‑AAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGAGCGGTGGCGGAACCCTCTGGGGCGCTGGCGCTGGCGGGGAAGaaaaaaaaaat                     >  2:562253/1‑245 (MQ=255)                                                                                                                                                                                                                                          cctgcctgAAAGCATCGCTGAA‑TGCGGG‑TCATCCGGTTGATCTGCCGCGCGTAGGGGTATTTGCTGAAGGCGGAGCGGTAAAACGCATCGGTGACG‑AAACCTTCCGTTTATGCCAGGAGTACCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGTCGCTGGCGGGAATGAAAAAATATATCGcc                 <  2:238936/251‑1 (MQ=255)                                                                                                                                                                                                                                            tgcctgAAAGCAGCGCTGAA‑TGCGGG‑TCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACG‑AAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGGATGAAAAAAGATATCccccc               >  2:95973/1‑250 (MQ=255)                                                                                                                                                                                                                                            ggcctgAAAGCAGTGCTGAA‑TGCGGGGTC‑TCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAATGCGTAGCGGTAAATCGCAGCGGTGACG‑AACCCTTCCGTTTATGCCAGGAGTATCTCGAGTACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTGATTCGAAGATGTGCGTGCGGTGGCGGAACCCTCTGGCGCGCTGGCTCTGGTGTGAAAGAAAAAATATATCGCCCt               <  2:535628/250‑1 (MQ=255)                                                                                                                                                                                                                                                 gAAAGCAGCGCTGAA‑TGCGGG‑TCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACG‑AAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATGAAAAAATATATCGCCCTGcac           <  1:212007/250‑1 (MQ=255)                                                                                                                                                                                                                                                 gAAAGCAGCGCTGAA‑TGCGGG‑TCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACG‑AAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATGAAAAAATATATCGCCCTGca            <  1:186137/249‑1 (MQ=255)                                                                                                                                                                                                                                                    agcagcGCTGAA‑TGCGGG‑TCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACG‑AAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATGAAAAAATATATCGCCCTGCacaac        >  1:398156/1‑250 (MQ=255)                                                                                                                                                                                                                                                      cagcGCTGAA‑TGCGGG‑TCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACG‑AAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGACGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATAACAAAAAATATCGCCCTGCACAACAtt     >  1:124985/1‑251 (MQ=255)                                                                                                                                                                                                                                                      cagcGCTGAA‑TGCGGG‑TCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACG‑AAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATGAAAAAATATATCGCCCTGCCCAACAtt     >  1:198206/1‑251 (MQ=255)                                                                                                                                                                                                                                                         cgcTGAA‑TGCGGG‑TCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACG‑AAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATGAAAAAATATATCGCCCTGCACAACATTcgc  <  1:249584/251‑1 (MQ=255)                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                         AGCGATCGAACTGTCACAGCAGCAGGGGTTCACCTGGGTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTG‑CCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGA‑TGCCGCAAATCAAAGTG‑ATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGGA‑TGCGGG‑TCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACG‑AAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATGAAAAAATATATCGCCCTGCACAACATTCGC  >  NC_000913/3955753‑3956248

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