Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,135,006	G→A	R391H (CGT→CAT)	katG →	catalase/hydroperoxidase HPI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,135,006	0	G	A	98.8%	293.6 / NA	84	R391H (CGT→CAT)	katG	catalase/hydroperoxidase HPI
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (39/43);  minor base C (1/0);  total (40/44)
Fisher's exact test for biased strand distribution p-value = 4.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

GCCATTACCTCTGGTCTGGAAGTAGTCTGGACCCAGACGCCGACCCAGTGGAGCAACTATTTCTTCGAGAACCTGTTCAAGTATGAGTGGGTACAGACCCGCAGCCCGGCTGGCGCAATCCAGTTCGAAGCGGTAGACGCACCGGAAATTATCCCGGATCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTATCGATCTGAAATTCGCGATTGCGGATTCTGGTCTGTCTGTTAGTGAGCTGGTATCGGTGGCCTGGG  >  NC_000913/4134759‑4135251                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      gCCATTACCGCTGGTCTGGAAGTAGTCTGGCCCCAGACGCCGACCCAGTGGAGCAACTATTTCTTCGAGATCCTGTTCAAGTATGAGGGGGTACAGACCCGCAGCCCGGCTGGCGCAATCCAGTTCGAAGCGGTAGACGCACCGGAAATTATCCCGGATCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCtcatcg                                                                                                                                                                                                                                                    <  2:114747/251‑1 (MQ=255)   cATTACCTCTGGTCTGGAAGTAGTCTGGACCCAGACGCCGACCCAGTGGAGCAACTATTTCTTCGAGAACCTGTTCAAGTATGAGTGGGTACAGACCCGCAGCCCGGCTGGCGCAATCCAGTTCGAAGCGGTAGACGCACCGGAAATTATCCCGGATCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCATCGtt                                                                                                                                                                                                                                                  >  1:47368/1‑251 (MQ=255)              gtctggAAGTAGTCTGGACCCAGACGCCGACCCAGTGGAGCAACTATGTCTTCGAGAACCTGTTCAAGTATGAGTGGGTACAGACCCGCAGCCCGGCTGGCGCAATCCAGTTCGAAGCGGTAGACGCACCGGAAATTATCCCGGATCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCATCGTTTCCTCAACGAt                                                                                                                                                                                                                                       <  1:37756/251‑1 (MQ=255)               tctggAAGTAGTCTGGACCCAGACGCCGACCCAGTGGAGCAACTATTTCTTCGAGAACCTGTTCAAGTATGAGTGGGTACAGACCCGCAGCCCGGCTGGCGCAATCCAGTTCGAAGCGGTAGACGCACCGGAAATTATCCCGGATCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCATCGTTTCCTCAACGATc                                                                                                                                                                                                                                      >  1:360096/1‑251 (MQ=255)                 tggAAGTAGTCTGGACCCAGACGCCGACCCAGTGGAGCAACTATTTCTTCGAGAACCTGTTCAAGTATGAGTGGGTACAGACCCGCAGCCCGGCTGGCGCAATCCAGTTCGAAGCGGTAGACGCACCGGAAATTATCCCGGATCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCATCGTTTCCTCAACGATCCg                                                                                                                                                                                                                                    >  2:114282/1‑251 (MQ=255)                    aaGTAGTCCGCACCCAGACGCCGACCCAGTGGAGCAACTATTTCTTAGATAACCTGTTCAAGTATGAGTGGGTACAGACCCGCAGCCCGGCTGGCGCAATCCAGTTCGAAGCGGTAGACGCACCGGAAATTATCCCGGATCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCATCGTTTCCTCAACGATCCGCAg                                                                                                                                                                                                                                 <  1:149852/251‑1 (MQ=255)                              gACCCAGACGCCGACCCAGTGGAGCAACTATTTCTTCGAGAACCTGTTCAAGTATGAGTGGGTACAGACCCGCAGCCCGGCTGGCGCAATCCAGTTCGAAGCGGTAGACGCACCGGAAATTATCCCGGATCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCATCGTTTCCTCAACGATCCGCAGGCGTTCAACg                                                                                                                                                                                                                       <  1:120579/251‑1 (MQ=255)                                 ccAGACGCCGACCCAGTGGAGCAACTATTTCTTCGAGAACCTGTTCAAGTATGAGTGGGTACAGACCCGCAGCCCGGCTGGCGCAATCCAGTTCGAAGCGGTAGACGCACCGGAAATTATCCCGGATCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCATCGTTTCCTCAACGATACGAAGGCGCTCAACGa                                                                                                                                                                                                                      >  2:218490/1‑249 (MQ=255)                                   aGACGCCGACCCAGTGGAGCAACTATTTCTTCGAGAACCTGTTCAAGTATGAGTGGGTACAGACCCGCAGCCCGGCTGGCGCAATCCAGTTCGAAGCGGTAGACGCACCGGAAATTATCCCGGATCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCATCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGcc                                                                                                                                                                                                                  <  1:88249/251‑1 (MQ=255)                                        ccGACCCAGTGGAGCAACTATTTCTTCGAGAACCTGTTCAAGTATGAGTGGGTACAGACCCGCAGCCCGGCTGGCGCAATCCAGTTCGAAGCGGTAGACGCACCGGAAATTATCCCGGATCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCATCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGc                                                                                                                                                                                                             <  1:14039/251‑1 (MQ=255)                                         cggcccAGTGGAGCAACTATTTCTTCGAGACCCTGTTCAAGTATGCGTGGGTACAGACCCGCAGCCCGGCTGGCGCAATCCAGTTCGAAGCGGTAGACGCACCGGAAATTATCCCGGATCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGTCGCTTCGTTTTGATCCTGAGTTCGAGAAGATCTCTCATCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGcc                                                                                                                                                                                                            <  2:238201/248‑1 (MQ=255)                                           aCCCAGTGGAGCAACTATTTCTTCGAGAACCTGTTCAAGTATGAGTGGGTACAGACCCGCAGCCCGGCTGGCGCAATCCAGTTCGAAGCGGTAGACGCACCGGAAATTATCCCGGATCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCATCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCg                                                                                                                                                                                                          >  1:66589/1‑251 (MQ=255)                                            cccAGTGGAGCAACTATTTCTTCGAGAACCTGTTCAAGTATGAGTGGGTACAGACCCGCAGCCCGGCTGGCGCAATCCAGTTCGAAGCGGTAGACGCACCGGAAATTATCCCGGATCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCATCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGt                                                                                                                                                                                                         >  1:34989/1‑251 (MQ=255)                                             ccAGTGGAGCAAATATTTCTTCGAGAACCTGTTCAAGTATGAGTGGGTCCAGCCCCGCAGCCCGGCTGGCGCCATCCAGTTCGAAGCGGTAGACGCACCGGAAATTATCCCGGATCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGAACTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCATCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTg                                                                                                                                                                                                        >  1:263956/1‑251 (MQ=255)                                              cAGTGGAGCAACTATTTCTTCGAGAACCTGTTCAAGTATGAGTGGGTACAGACCCGCAGCCCGGCTGGCGCAATCCAGTTCGAAGCGGTAGACGCACCGGAAATTATCCCGGATCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCATCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGc                                                                                                                                                                                                       <  1:188591/251‑1 (MQ=255)                                                  ggAGCAACTATTTCTTCGAGAACCTGTTCAAGTATGAGTGGGTACAGACCCGCAGCCCGGCTGGCGCAATCCAGTTCGAAGCGGTAGACGCACCGGAAATTATCCCGGATCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCATCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTg                                                                                                                                                                                                    <  2:27864/250‑1 (MQ=255)                                                    aGCAACTATTTCTTCGAGAACCTGTTCAAGTATGAGTGGGTACAGACCCGCAGCCCGGCTGGCGCAATCCAGTTCGAAGCGGTAGACGCACCGGAAATTATCCCGGATCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCATCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGtt                                                                                                                                                                                                 <  1:78774/251‑1 (MQ=255)                                                     gCAACTATTTCTTCGAGAACCTGTTCAAGTATGAGTGGGTACAGACCCGCAGCCCGGCTGGCGCAATCCAGTTCGAAGCGGTAGACGCACCGGAAATTATCCCGGATCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGCAGATCTCTCATCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTc                                                                                                                                                                                                >  2:83696/1‑251 (MQ=255)                                                            tttCTTCGAGAACCTGTTCAAGTATGAGTGGGTACAGACCCGCAGCCCGGCTGGCGCAATCCAGTTCGAAGCGGTAGACGCACCGGAAATTATCCCGGATCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCATCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTg                                                                                                                                                                                          >  2:81419/1‑250 (MQ=255)                                                              tcttcGAGAACCTGTTCAAGTATGAGTGGGTACAGACCCGCAGCCCGGCTGGCGCAATCCAGTTCGAAGCGGTAGACGCACCGGAAATTATCCCGGATCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCATCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGAc                                                                                                                                                                                        <  1:258623/250‑1 (MQ=255)                                                              gcttcGAGACACTGGTCAAGTATGAGGGGGTACAGACCCGCAGCCCGGCTGGCGCAATCTAGTTCGAAGCGGTAGACGCAGCGTAAATGATCCCGGATTCGTTGGATCCGGCGAAGAAATGTAAACCGACAATGCGGGTGACCGACCTGCCGCTGGTTTTTGTTCCTTAGGTCGAGAAGAGCTCGCTTCGGGTCCTCATCGATCCGCAGGTGTTGCACGAAGGCTTTTCCGGTGGGTGGTTCGAACTGAc                                                                                                                                                                                        <  2:270280/249‑1 (MQ=255)                                                                ttcGAGAACCTGTTCAAGTATGAGTGGGTACAGACCCGCAGCCCGGCTGGCGCAATCCAGTTCGAAGCGGTAGACGCACCGGAAATTATCCCGGATCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCATCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAACCTGACGcc                                                                                                                                                                                     >  1:17199/1‑250 (MQ=255)                                                                 tcGAGAACCTGTTCAAGTATGAGTGGGTACAGACCCGCAGCCCGGCTGGCGCAATCCAGTTCGAAGCGGTAGACGCACCGGAAATTATCCCGGATCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCGCATCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACg                                                                                                                                                                                       <  1:378733/248‑1 (MQ=255)                                                                    agaACCTGTTCAAGTATGAGTGGGTACAGACCCGCAGCCCGGCTGGCGCAATCCAGTTCGAAGCGGTAGACGCACCGGAAATTATCCCGGATCCGTTTGATCCGTCGAAGAAACGTAACCCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCATCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAgg                                                                                                                                                                                 >  1:294838/1‑251 (MQ=255)                                                                      aaCCTGTTCAAGTATGAGTGGGTACAGACCCGCAGCCCGGCTGGCGCAATCCAGTTCGAAGCGGTAGACGCACCGGAAATTATCCCGGATCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCATCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAggg                                                                                                                                                                                >  1:189674/1‑250 (MQ=255)                                                                       aCCTGTTCAAGGATGAGTGGGTACAGACCCGCAGCCCGGCTGGCGCAATCCAGTTCGAAGCGGTAGACGCACCGGAAATTATCCCGGATCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCATCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAACCTGACGCACAGGGat                                                                                                                                                                              <  2:88840/251‑1 (MQ=255)                                                                             tCAAGTACGAGTGGGTACAGACCCGCAGCCCGGCTGGCGCAATCCAGTTCGAAGCGGTAGACGCACCGGAAATTATCCCGGATCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCATCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATggg                                                                                                                                                                         <  1:297621/250‑1 (MQ=255)                                                                              cAAGTATGAGTGGGTACAGACCCGCAGCCCGGCTGGCGCAATCCAGTTCGAAGCGGTAGACGCACCGGAAATTATCCCGGATCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCATCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGc                                                                                                                                                                       <  1:80135/251‑1 (MQ=255)                                                                                    taaGTGGGGACAGCCCCGCAGCCCGGCTGGCGCCATCCAGTGCGAAGCGATAGACGCACCGGAAATTATCCCGGATCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCATCGTTTCCTCAACGAACCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAAt                                                                                                                                                                 >  2:284062/3‑251 (MQ=255)                                                                                       gTGGGTACAGACCCGCAGCCCGGCTGGCGCAATCCAGTTCGAAGCGGTAGACGCACCGGAAATTATCCCGGATCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCATCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAAtctc                                                                                                                                                              <  2:326026/251‑1 (MQ=255)                                                                                            tACAGACCCGCAGCCCGGCTGGCGCAATCCAGTTCGAAGCGGTAGACGCACCGGCAATTATCCCGGATCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCATCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATg                                                                                                                                                                           >  2:285193/1‑233 (MQ=255)                                                                                            tACAGACCCGCAGCCCGGCTGGCGCAATCCAGTTCGAAGCGGTAGACGCACCGGAAATTATCCCGGATCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCATCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATg                                                                                                                                                                           <  1:285193/233‑1 (MQ=255)                                                                                                    cGCAGCCCGGCTGGCGCAATCCAGGTCGAAGCGGTAGACGCACCGGAAATTATCCCGGATCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCATCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGcc                                                                                                                                                 <  1:45792/251‑1 (MQ=255)                                                                                                    cGCAGCCCGGCTGGCGCAATCAAGTTCGAAGCGGTAGACGCACCGGACATTATCCCGGATCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCATCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGACTTTGCCCGTGCCTGGTTCAAACTGCCGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGcc                                                                                                                                                 >  2:237585/1‑251 (MQ=255)                                                                                                      cAGCCCGGCTGGCGCAATCCAGTTCGAAGCGGTAGACGCACCGGAAATTATCCCGGATCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCATCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCgg                                                                                                                                               >  1:65010/1‑251 (MQ=255)                                                                                                       aGCCCGGCTGGCGCAATCCAGTTCGAAGCGGTAGACGCACCGGAAATTATCCCGGATCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCATCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCgg                                                                                                                                               >  1:321860/1‑250 (MQ=255)                                                                                                       aGCCCGGCTGGCGCAATCCAGTTCGAAGCGGTAGACGCACCGGAAATTATCCCGGATCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCATCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCgg                                                                                                                                               >  1:196948/1‑250 (MQ=255)                                                                                                         cccGGCTGGCGCAATCCAGTTCGAAGCGGTAGACGCACCGGAAATTATCCCGGATCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCATCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAg                                                                                                                                            <  2:111185/251‑1 (MQ=255)                                                                                                           cGGCTGGCGCAATCCAGTTCGAAGCGGTAGACGCACCGGAAATTATCCCGGATCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCATCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTg                                                                                                                                          >  1:107140/1‑251 (MQ=255)                                                                                                            ggCTGGCGCAATCCAGTTCGAAGCGGTAGACGCACCGGAAATTATCCCGGATCCGTTTGATCCGGCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCATCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGc                                                                                                                                         <  2:363442/251‑1 (MQ=255)                                                                                                                gtcgcAAGCCAGTTCGAAGCGGTAGACGCACCGGAAAGTATCCCGGATCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGGCCTGGCGCTGCGTTTTGATCCTGAGGTCGAGAAGATCTCTCATCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTGTGCCCGTGCCTGGTTCAAACGGACGCGCAGGGATAGGGGGGCGAAATCGCGCTACATCGGGCCGGAAGGGCCGaa                                                                                                                                     <  2:266925/249‑1 (MQ=255)                                                                                                                ggCGCCATCCCGTTCGAAGCGGTAGACGCACCGGTAATTATCCCGGATCCGTTTGATCCGTCGAAGAAACGCAAACCGACAATGCTGGTGACCGACCTTACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCCTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCGAACTGACGACCAGGGAGATGGGACAGAAATCTCGCGACATCGCGCCGGAAGTGCCGa                                                                                                                                      >  2:374099/1‑250 (MQ=255)                                                                                                                   gcAATCCAGTTCGAAGCGGTAGACGCACCGGAAATTATCCCGGATCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCATCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAaaga                                                                                                                                  <  2:147851/251‑1 (MQ=255)                                                                                                                      aTCCAGTTCGAAGCGGTAGACGCACCGGAAATTATCCCGGATCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCATCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAagaaga                                                                                                                               <  1:114282/251‑1 (MQ=255)                                                                                                                                aaGCGGTAGACGCACCGGAAATTATCCCGGATCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCATCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGATCTgg                                                                                                                     <  2:189674/251‑1 (MQ=255)                                                                                                                                 ggCGGTAGACGCACCGGAAATTATCCCGGATCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTGTATCCTGAGTTCGAGAAGATCTCTCATCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGTCCGGAAGTGCCGAAAGAAGATCTGATCTGGc                                                                                                                    <  2:249996/250‑1 (MQ=255)                                                                                                                                       aGACGCACCGGAAATTATCCCGGATCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCATCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGATCTGGCAAGAt                                                                                                               <  1:86496/250‑1 (MQ=255)                                                                                                                                        gACGCACCGGAAATTATCCCGGATCCGTTTGATCCGGCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTGCGAGAAGATCTCTCATCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACGTACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGATCTGGCAAGATcc                                                                                                             <  2:64992/251‑1 (MQ=255)                                                                                                                                             cccGGAAATTATCCCGGATCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCATCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGAGCTGGCAAGTTCCGCTGc                                                                                                        <  2:190644/250‑1 (MQ=255)                                                                                                                                                    aTTATCCCGGATCCGTTTGATCCGGCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCATCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGcc                                                                                                 <  1:120605/251‑1 (MQ=255)                                                                                                                                                     tgaTCCCGGACCCGTTTGTGCCGTCGAAGTAACGTAAACCGACAATGCGTGTGCCCGACCTGACGCTGCGTTTTGATCCTGAGTGCGAGAAGATCTCTCAGCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAATCGGACGGACAGGGATATGGGGTCGATAGCGCGCTTCATCGGGCCGGATGTGCCGACAGAAGAGCTGATCTGGCAAGATCCGCTGCCGCAGcc                                                                                                 <  2:190024/248‑1 (MQ=255)                                                                                                                                                        tCCCGGATCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCATCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATc                                                                                             >  1:255842/1‑251 (MQ=255)                                                                                                                                                                ccGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCATCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAAccc                                                                                     <  1:110724/251‑1 (MQ=255)                                                                                                                                                                  gTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCATCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGCTCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACccg                                                                                    >  2:396129/1‑250 (MQ=255)                                                                                                                                                                         tcGTCGAAGAAACGTAAACCGACAATGCTGGTGTCCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCATCGTTGCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGAGCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAg                                                                              <  1:139084/248‑1 (MQ=255)                                                                                                                                                                         ccGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCATCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCa                                                                            <  2:225603/251‑1 (MQ=255)                                                                                                                                                                             cGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCATCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGAc                                                                        >  1:66442/1‑251 (MQ=255)                                                                                                                                                                             cGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCATCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGATCTGGCAAGATCCGCGGCCGCAGCCGATCTACAACCCGACCGAGCAGGa                                                                         <  1:284475/250‑1 (MQ=255)                                                                                                                                                                              gaagaaACGTAAACCGACAAGGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCATCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACa                                                                       <  2:360096/251‑1 (MQ=255)                                                                                                                                                                                agaaACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCATCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGCTCTGGCAAGATCCGCTGACGCAGCCGATCTACAACCCGACCGAGCAGGACAtt                                                                     >  1:166262/1‑251 (MQ=255)                                                                                                                                                                                   aaCGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCGGAGTTCGAGAAGATCTCTCATCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTATc                                                                  <  2:17199/251‑1 (MQ=255)                                                                                                                                                                                    aCGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCATCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTCCAACCCGACCGAGCAGGACATTATCg                                                                 >  2:393977/1‑251 (MQ=255)                                                                                                                                                                                      gTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCATCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTATCGAt                                                               <  1:224808/251‑1 (MQ=255)                                                                                                                                                                                             gACAATGCTGGTGACCGCCCTGACGCTGCGTTTTGATCCTGAGTTAGAGAAGATCTCTCATCGATTCCTCAACGATCCGCAGGAGTTCAACGAAGCATTTGCCAGTGCCTGGCTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGAGCTGGCAAGATCCGCTGCCGCAGCCGCTCTACAACCCGACCGAGCAGGACAGTAGCGATCTGAAAt                                                        >  2:203791/1‑251 (MQ=255)                                                                                                                                                                                             gACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCATCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTATCGATCTGAAAt                                                        >  1:378820/1‑251 (MQ=255)                                                                                                                                                                                             gACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCATCGTTTCCTCAACGATCCGCAGGCGCTCAACGAAGCCTTTGACCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTATCGATCTGAAAt                                                        >  1:403819/1‑251 (MQ=255)                                                                                                                                                                                              aCAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCATCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGACGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACCTTATCGATCTGAAAtt                                                       >  2:41469/1‑251 (MQ=255)                                                                                                                                                                                               cAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCATCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTATCGATCTGAAATTc                                                      <  1:105843/251‑1 (MQ=255)                                                                                                                                                                                               cAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCATCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTATCGATCTGAAATTc                                                      >  1:173313/1‑251 (MQ=255)                                                                                                                                                                                               cAATGCTGGTGACCGACCTGACGCTGCGTTGTGATCCTGAGTTCGAGAAGATCTCTCATCGTGTCCTCAACGATCCGCAGGCGTGCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATAGGGGGCTGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGATCGGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTATCGATCTGAAATTc                                                      <  1:284062/251‑1 (MQ=255)                                                                                                                                                                                                  tGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCATCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGCAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGAAGTTATCGATCTGAAATTcgcg                                                   >  2:355454/1‑251 (MQ=255)                                                                                                                                                                                                    cTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCATCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTATCGCTCTGAAATTCGCGAt                                                 >  2:120651/1‑251 (MQ=255)                                                                                                                                                                                                        tGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGCAGATCTCTCATCGTTTCCTCAACGATCCGCATGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTAAAACCCGACCGAGCAGGAAATTATCGATCTGAAATTCGCGATTGCg                                             >  1:239623/1‑251 (MQ=255)                                                                                                                                                                                                         gaccgaccTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCATCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTATCGATCTGAAATTCGCGATTGCgg                                            <  1:81419/251‑1 (MQ=255)                                                                                                                                                                                                            cgaccTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCATCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTATCGATCTGAAATTCGCGATTGCGGAtt                                         >  2:177738/1‑251 (MQ=255)                                                                                                                                                                                                             gaccTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCATCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTATCGATCTGAAATTCGCGATTGCGGATTc                                        >  2:155204/1‑251 (MQ=255)                                                                                                                                                                                                                cgtaCGCTGCGTTTTGATCCTGAGGTCGAGAAGATCTCTCATCGTTGCCTCAACGATCCGCAGGCGGTCAACGAAGCCTTTGCCCGTGCCTTGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCGCGCTACATCGGGCCGGAAGGGCCGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTATCGAACTGAAATTCGCGATAGCGGATTCTgg                                     <  2:69183/248‑1 (MQ=255)                                                                                                                                                                                                                cTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCATCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTATCGATCTGAAATTCGCGATTGCGGATTCTgg                                     >  1:399344/1‑251 (MQ=255)                                                                                                                                                                                                                        gCGTTTTGATCCTGAGTTCGAGAAGATCTCTCATCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTATCGATCTGAAATTCGCGATTGCGGATTCTGgtctgtc                               <  1:287155/249‑1 (MQ=255)                                                                                                                                                                                                                        gCGTTTTGATCCTGAGTTCGAGAAGATCTCTCATCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATATACAACCAGACCGAGCAGGACATTATCGATCAGAAATTCGCGCATGCGGATTATGGtctgtctg                             >  2:340119/1‑251 (MQ=255)                                                                                                                                                                                                                                       gTTCGAGAAGATCTCTCATCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTATCGATCTGAAATTCGCGATTGCGGATTCTGGTCTGTCTGTTAGGGAGCTGGTAt              <  2:66589/251‑1 (MQ=255)                                                                                                                                                                                                                                           gagaAGATCTCTCATCGTTTACTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGCGCAGGACATTATCGATCTGAAATTCGCGATTGCGGATTCTGGTCTGTCTGTTAGTGAGCTGGTATCggg          >  2:417986/1‑250 (MQ=255)                                                                                                                                                                                                                                            agaagaTCCCTCATCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTATCGATCTGAAATTCGCGATTGCGGATTCTGGTCTGTCTGTTAGTGAGCGGGTATCGGTg         <  2:399344/251‑1 (MQ=255)                                                                                                                                                                                                                                                   ctctcATCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCCGCCGATCTACAACCCGACCGAGCAGGACATTATCGATCTGAAATTCGCGATTGCGGATTCTGGTCTGTCTGTTAGTGAGCTGGTATCGTTGGCCTggg  >  2:312045/1‑251 (MQ=255)                                                                                                                                                                                                                                                   ctctcATCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGGGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGCAATATCGCTACATCGGGCCGGACGTGCCGAAAGAAGATATGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATCATCGCTCTGAAATTCGAGATTGCGGATTATGGGCTGTCTGGTAGTGAGCTGGCATCGGTGGCCTggg  >  1:217229/1‑251 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      GCCATTACCTCTGGTCTGGAAGTAGTCTGGACCCAGACGCCGACCCAGTGGAGCAACTATTTCTTCGAGAACCTGTTCAAGTATGAGTGGGTACAGACCCGCAGCCCGGCTGGCGCAATCCAGTTCGAAGCGGTAGACGCACCGGAAATTATCCCGGATCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGAAATCTCGCTACATCGGGCCGGAAGTGCCGAAAGAAGATCTGATCTGGCAAGATCCGCTGCCGCAGCCGATCTACAACCCGACCGAGCAGGACATTATCGATCTGAAATTCGCGATTGCGGATTCTGGTCTGTCTGTTAGTGAGCTGGTATCGGTGGCCTGGG  >  NC_000913/4134759‑4135251

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