Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,512,501	C→T	A303A (GCG→GCA)	fecC ←	ferric citrate ABC transporter membrane subunit FecC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,512,501	0	C	T	96.0%	195.6 / NA	76	A303A (GCG→GCA)	fecC	ferric citrate ABC transporter membrane subunit FecC
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (35/38);  minor base G (3/0);  total (38/38)
Fisher's exact test for biased strand distribution p-value = 2.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.77e-01

GCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAAT‑CAGCGCCAGCACTGC‑GCCTGCGGGCAGAT‑CT‑CCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACC‑GGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGC  >  NC_000913/4512258‑4512739                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                gCAGTCGGTAATCCATCAATACATACTAATGCTCGCGTCAGGCCAGCCAGTCGGACAGCAGAGCGCGCCACGGCACGGGGATCACTCCAATATGGAGTGATACCAGCGCACAGCCCGCCAGGGCAAGGGTGATGAACATAACCAGCGCAATTTTCATCCTAGCCTCCTCACCATCCCGCAAAAGCAAGGGCTGCCAAT‑CAGCGCCATCACTGC‑GCCTGCGGGCAGAT‑CT‑CCGGGGCAGTACAGGGCGCGTg                                                                                                                                                                                                                                          >  2:199165/1‑251 (MQ=255)      cGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAAT‑CAGCGCCAGCACTGC‑GCCTGCGGGCAGAT‑CT‑CCGGGGAAGGCCAGTGCGCGTGCGAGt                                                                                                                                                                                                                                     <  1:496885/251‑1 (MQ=255)         tACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGCTGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAAT‑CAGCGCCAGCTCTGC‑GCCTGCGGGCAGAT‑CT‑CCAGGGAAGGCCAGTGCGCGTGCGAGTAc                                                                                                                                                                                                                                   >  2:349777/1‑250 (MQ=255)         tACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAAT‑CAGCGCCAGCACTGC‑GCCTGCGGGCAGAT‑CT‑CCGGGGAAGGCCAGGGCGCGTGCGAGTAc                                                                                                                                                                                                                                   >  1:275679/1‑250 (MQ=255)             ccATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGCGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAAT‑CAGCGCCAGCACTGC‑GCCTGCGGGCAGAT‑CT‑CCGGGGAAGGCCAGTGCGCGTGCGAGTACATCt                                                                                                                                                                                                                               <  1:355242/250‑1 (MQ=255)                       ataataATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGCGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAAT‑CAGCGCCAGCACTGC‑GCCTGCGGGCAGAT‑CT‑CCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAgcagcac                                                                                                                                                                                                                    >  2:144380/1‑250 (MQ=255)                       ataataAGGCCCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAAT‑CAGCGCCAGCACTGC‑GCCTGCGGGCAGAT‑CT‑CCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCAt                                                                                                                                                                                                                    <  2:215487/251‑1 (MQ=255)                         aataatGCTCGCGTCCGGCCTGCCAGGCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAAT‑CAGCGCCAGCACTGC‑GCCTGCGGGCAGAT‑CT‑CCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATc                                                                                                                                                                                                                   <  1:70726/250‑1 (MQ=255)                               gctcgcttCCGGCCTGCCAGTCGGCCTGCAGCGCGCGCCACGGCCCGGGGATCACTCCCATATGGCGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGAGGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAAC‑CAGCGCCAGCACTGC‑GCCGGCGGGCAGAT‑CA‑CCTGGGACGGCCAGTGCGCGTGCGAGAACATCTGcccgacgcctcagcgcgg                                                                                                                                                                                                            >  2:142288/8‑235 (MQ=255)                               gCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAAT‑CAGCGCCAGCACTGC‑GCCTGCGGGCAGAT‑CT‑CCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTgg                                                                                                                                                                                                            >  1:164347/1‑251 (MQ=255)                                       ccGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAAT‑CAGCGCCAGCACTGC‑GCCTGCGGGCAGAT‑CT‑CCGGGGACGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCcagc                                                                                                                                                                                                    >  2:42055/1‑251 (MQ=255)                                         ggCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAAT‑CAGCGCCAGCACTGC‑GCCTGCGGGCAGAT‑CT‑CCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCcagcag                                                                                                                                                                                                  <  1:396826/251‑1 (MQ=255)                                            cTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAAT‑CAGCGCCAGCACTGC‑GCCTGCGGGCAGAT‑CT‑CCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCAt                                                                                                                                                                                               >  2:288981/1‑251 (MQ=255)                                                 aGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAAT‑CAGCGCCAGCACTGC‑GCCTGCGGGCAGAT‑CT‑CCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTc                                                                                                                                                                                           >  2:83857/1‑250 (MQ=255)                                                       tCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAAT‑CAGCGCCAGCACTGC‑GCCTGCGGGCAGAT‑CT‑CCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGc                                                                                                                                                                                     >  2:95116/1‑250 (MQ=255)                                                        cagcagCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAAT‑CAGCGCCAGCACTGC‑GCCTGCGGGCAGAT‑CT‑CCGGGGAAGGCCAGTGCGCGTGCGGGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCCTGCTCCCTGGc                                                                                                                                                                                     >  2:113622/1‑249 (MQ=255)                                                            agcGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAAT‑CAGCGCCAGCACTGC‑GCCTGCGGGCAGAT‑CT‑CCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCCGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACg                                                                                                                                                                               >  1:74926/1‑251 (MQ=255)                                                                    ccACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAAT‑CAGCGCCAGCACTGC‑GCCTGCGGGCAGAT‑CT‑CCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTg                                                                                                                                                                       <  1:220343/251‑1 (MQ=255)                                                                       cggcacggGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGAGTGATGAAAATAACCAGGGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCGGCAAAT‑AAGCGCCTGCACTGC‑GCCTCCGAGACGAT‑CT‑CCGGGGAAGGCCAGTGCGCGTGCGATTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCAATAGCGGTCCGTTGCCCTGagc                                                                                                                                                                    >  2:36433/1‑249 (MQ=255)                                                                       cggcacggGGATCACGCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAAT‑CAGCGCCAGCACTGC‑GCCTGCGGGCAGAT‑CT‑CCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATc                                                                                                                                                                    <  2:157431/251‑1 (MQ=255)                                                                       cggcacggGGATAACTCCCATATGGAGTGATAACAGCGCGCAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCACGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAAT‑CAGCGCCAGCACTGC‑GCCTGCGGGCAGAT‑CG‑CCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATc                                                                                                                                                                    <  1:275921/251‑1 (MQ=255)                                                                        ggcacggGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGCGTGATGAAAATAACCAGAGCAATTTTCATCCTCGCCTCCTCACAAGCAAGCCAAAGCAAGGGCTGCCAAT‑CAGCGCCAGCACTGA‑GCCTGCGGGCAGAT‑CT‑CCGGGGAAGGCCAGTGCGCGAGCGAGTAAATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCACGCTCACCGGCAGTACggtggtcggatcg                                                                                                                                                                   >  2:360095/1‑239 (MQ=255)                                                                        ggcacggGGATCACTCCCATATGGAGTGATAACAGCGCAAAGCCCGCCAGGGCAAGGGTGATGGAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAAT‑CAGGGCCAGCACTGC‑GCCTGCGGACAGAT‑CT‑CCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGAATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCg                                                                                                                                                                   >  1:189733/1‑251 (MQ=255)                                                                        ggcacggGGATCACGCGCATATGGAGTTATAACAGCGCACAGCCCGCCATGGCAAGGGTGATGAAAAAAACCAGCGCAAGTTTCATCCTCGCCTCCTCTCAAGCCAGACAAAGCAAGTGCTGCAAAT‑CAGCGCCACCACTGC‑GCCTGCGGGCAGAT‑CG‑CCGGGGAAGGCCAGTGCGCGTGCTAGGACATCTTCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCg                                                                                                                                                                   <  2:320362/251‑1 (MQ=255)                                                                           acggGGATCGCTCCCATATGGAGTGATAACTGCGCACAGTCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAAT‑CAGCGCCAGCACTGC‑GCCTGCGGGCAGAT‑CT‑CCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCAGCAGCGTGGACCCCAGCAGCATGCTCACGGGCAGTACGTTGCGCTGATCGaat                                                                                                                                                                <  2:78294/251‑2 (MQ=255)                                                                           acggGGATCACTACCATATGGAGTGATAAAAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACACAGCAAGGGCTGCCAAT‑CAGCGCCAGCACTGC‑GCCTGCAGGCAGAT‑CT‑CCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGGGGCCCCCAGCAGCCTGCTCACTGGCAGTACGTTGCGCTGATCGa                                                                                                                                                                  >  1:424082/1‑249 (MQ=255)                                                                              gggATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGGGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAAT‑CAGCGCCAGCACTGC‑GCCTGCGGGCAGAT‑CT‑CCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCt                                                                                                                                                             >  1:386581/1‑251 (MQ=255)                                                                                ggtCACTCCCATATGGAGTGAAAACAGCGCACAGCCAGCCAGGGCCATTGTGCTGTAAATAACCAGCGCAATTTTCATACTCGCCTCCTCACAGGCCAGACAAAGCAAGGGCTGCCAAT‑CAGCGCCAGCACTGC‑GCCTGCGGGCAGAT‑CC‑CCGGGGAAGGCCAGGGCGCGTGCGCGTACATCTGCCAGCAGCATCAGCGTGCCCCCCAGCAGCATGCCCACTGGCAGTACGTTGAGCTGATCGCAGCGTGc                                                                                                                                                           >  2:4342/3‑251 (MQ=255)                                                                                gATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGACAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCAAAAGCCAGACAAAGCAAGGGCTGCCAAT‑CAGCGCCAGCACTGC‑GCCTGCGGGCAGAT‑CT‑CCGGGGAAGGCCAGTGAGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGc                                                                                                                                                           >  2:438773/1‑251 (MQ=255)                                                                                    ccTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAAT‑CAGCGCCAGCACTGC‑GCCTGCGGGCAGAT‑CT‑CCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAg                                                                                                                                                       <  1:9763/250‑1 (MQ=255)                                                                                       cccATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAAT‑CAGCGCCAGCACTGC‑GCCTGCGGGCAGAT‑CT‑CCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGaa                                                                                                                                                     <  1:126669/250‑1 (MQ=255)                                                                                       cccATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAAT‑CAGCGCCAGCACTGC‑GCCTGCGGGCAGAT‑CT‑CCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAg                                                                                                                                                    <  1:326316/251‑1 (MQ=255)                                                                                           tatGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAAT‑CAGCGCCAGCACTGC‑GCCTGCGGGCAGAT‑CT‑CCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAgcgcg                                                                                                                                                >  2:22600/1‑251 (MQ=255)                                                                                            atGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAAT‑CAGCGCCAGCACTGC‑GCCTGCGGGCAGAT‑CT‑CCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAgcgcgc                                                                                                                                               >  1:431672/1‑251 (MQ=255)                                                                                                            gcACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAAT‑CAGCGCCAGCACTGC‑GCCTGCGGGCAGAT‑CT‑CCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACcagc                                                                                                                               >  1:125470/1‑251 (MQ=255)                                                                                                            gcACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAAT‑CAGCGCCAGCACTGC‑GCCTGCGGGCAGAT‑CT‑CCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACcagc                                                                                                                               >  1:417939/1‑251 (MQ=255)                                                                                                              acaGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAAT‑CAGCGCCAGCACTGC‑GCCTGCGGGCAGAT‑CT‑CCGGGGAAGGCCAGTGCGCCTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGACGCGCGCCAGATGCGGCACcagc                                                                                                                               >  2:446754/1‑249 (MQ=255)                                                                                                                aGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAAT‑CAGCGCCAGCACTGC‑GCCTGCGGGCAGAT‑CT‑CCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGAc                                                                                                                           >  1:321762/1‑251 (MQ=255)                                                                                                                         gggCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCGCGCCTCCCCACAAGCCAGACAAAGCAAGGGCTGCCAAT‑CAGCGCCAGCACTGC‑GCCTGCGGGCAGAT‑CT‑CCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACg                                                                                                                  <  1:117490/251‑1 (MQ=255)                                                                                                                              aGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAAT‑CAGCGCCAGCACTGC‑GCCTGCGGGCAGAT‑CT‑CCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATTCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCAcc                                                                                                             <  1:214093/251‑1 (MQ=255)                                                                                                                                  ggatgaAAATAACCAGCGTAATTTTCATCCCCGCCTCCTCACCAGCCAGACAAAGCAAGGGCTGCCAAC‑CAGCGCCAGCACTGC‑GCCTGCGGGCAGAT‑CT‑CCGGGGAAGGCCAGTGCGCGGGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACC‑GGAc                                                                                                        <  2:45930/250‑1 (MQ=255)                                                                                                                                   gatgaAAAAACCCAGCGATATTCTAATCCTCTCCGCCTCACAAGCCAGCCAAATCAAGGGTTGCCAAC‑AAGCGCCAGCACTGC‑GCCTGCGGTCAGAT‑CT‑CCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCAGCAGTGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTTATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCCCCAGCAGACCGATAAACGCCTCC‑GGAc                                                                                                        <  2:270921/250‑1 (MQ=255)                                                                                                                                      gaAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAAT‑CAGCGCCAGCACTGC‑GCCTGCGGGCAGAT‑CT‑CCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACC‑GGACCTg                                                                                                     <  1:227934/250‑1 (MQ=255)                                                                                                                                            aaCCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAAT‑CAGCGCCAGCACTGC‑GCCTGCGGGCAGAT‑CT‑CCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGAc                                                                                                                           >  1:378272/1‑223 (MQ=255)                                                                                                                                            aaCCAGCGCAATTTTCATCCCCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAAT‑CAGCGCCAGCACTGC‑GCCTGCGGGCAGAT‑CT‑CCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGAc                                                                                                                           <  2:378272/223‑1 (MQ=255)                                                                                                                                                     aaTTTTCATCCTCGCCTCCTCACAAGCCAGACAAGGCAAGGGCTGCCAAT‑CAGCGCCAGCACTGC‑GCCTGCGGGCAGCT‑CT‑CCGGGGAAGTCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGCCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCGGCCCAGAAGCGCGCCCGATGTGGCACCAGCAGACCGCTCAAC‑CCACCCGGACCTGCCACGCtaacccacac                                                                                     >  2:78870/1‑242 (MQ=255)                                                                                                                                                        tttCATCCGCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAAT‑CAGCGCCAGCACTGC‑GCCTGCGGGCAGAT‑CT‑CCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACC‑GGACCTGCCACGCTGACgcacgcac                                                                                   <  1:83857/250‑1 (MQ=255)                                                                                                                                                             tCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAAT‑CAGCGCCAGCACTGC‑GCCTGCGGGCAGAT‑CT‑CCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCca                                                                                                                                                                                                      <  1:419517/131‑1 (MQ=255)                                                                                                                                                             tCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAAT‑CAGCGCCAGCACTGC‑GCCTGCGGGCAGAT‑CT‑CCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCca                                                                                                                                                                                                      >  2:419517/1‑131 (MQ=255)                                                                                                                                                             cccccccctcccCACACGCCAGAAAAAGCAAGGGC‑GGCAATACAGCGCCAGCCCTGC‑GCCTGCGGGCAGAT‑CT‑CCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACC‑GGACCTGCCACGCTGACGCACGCaccaacca                                                                             <  2:406016/245‑1 (MQ=255)                                                                                                                                                              ccTCGCCTCCTCCCAAGCCAGACAAAGCAAGGGCTGCCAAT‑CAGCGCCAGCACTGC‑GCCGGCGGGCAGAT‑CT‑CCGGGAAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCAGCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACC‑GGACCTGCCACGCTGACGCACGCACCAACCAg                                                                            <  2:28929/251‑1 (MQ=255)                                                                                                                                                                 cGCCTCCGCACAAGCCAGACAAAGCAAGGGCTGCCAAT‑CAGCGCCAGCACTGC‑GCCTGCGGGCAGAT‑CT‑CCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCAGGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAACCGCCACC‑GGGCCTGCCACGCTGACGCACGCACCAACCAGAag                                                                         <  1:20986/251‑1 (MQ=255)                                                                                                                                                                   cctcctCACAAGCCAGACAAAGCAAGGGCTGCCAAT‑CAGCGCCAGCACTGC‑GCCTGCGGGCAGAT‑CT‑CCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACC‑GGACCTGCCACGCTGACGCACGCACCAACCAGAagc                                                                        >  1:127626/1‑250 (MQ=255)                                                                                                                                                                      cctcACAAGCCAGACAAAGCAAGGGCTGCCAAT‑CAGCGCCAGCCCGGC‑GCCGGCGGGCAGAT‑CT‑CCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACC‑GGCCCTGCCACGCTGACGCACGCACCAACCAGAAgcagca                                                                    <  2:483227/251‑1 (MQ=255)                                                                                                                                                                      cccaacaAGCAAGAAAAAGCAAGGGCTGCCAAT‑CAGGGCCAGCACTGC‑GCCTGCGGGCAGAC‑CT‑CCGGGGAAGGCCAGTGCGCGGGCGAGTACATCTGCCAGCAGCATAAGCGTGGCCCCCAGCAGCATGCTCACTGGCCGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGAAGACCGATAAACGCCACG‑GGACCAGCCACGCTGACGCACGCACCAACCAGAAgcagca                                                                    <  2:224703/247‑1 (MQ=255)                                                                                                                                                                                     aaaGCAAGGGCTGCCAAT‑CAGCGCCAGCACTGC‑GCCTGCGGGCAGAT‑CT‑CCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGACGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGCTAAACGCCCCC‑GGACCTGCCCCGCTGACGCACGCACCACCCAGAAGCCACACTAACCTATtgatc                                                      >  2:400198/1‑249 (MQ=255)                                                                                                                                                                                     aaaGCAAGGGCTGCCAAC‑AAGCGCCAGCACTGC‑GCCTGCGGGCAGAT‑CT‑CCGGGGAAGGCCAGTGCGCGGGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACC‑GGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATtgat                                                       <  1:400198/249‑1 (MQ=255)                                                                                                                                                                                      aaGCAAGGGCTGCCAAT‑CAGCGCCAGCACTGC‑GCCTGCGGGCAGAT‑CT‑CCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACC‑GGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGAc                                                    <  1:288981/251‑1 (MQ=255)                                                                                                                                                                                           aGGGCTGCCAAT‑CAGCGCCAGCACTGC‑GCCTGCGGGCAGAT‑CT‑CCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTTATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACC‑GGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAAc                                               <  1:167653/251‑1 (MQ=255)                                                                                                                                                                                             gccggCCAAT‑CACCGCCAGCAC‑GCAGCCTGCGGGC‑GATGCC‑CCGGGGAAGGCCAGTGCGCGTGCGAGCACATCTGCCACCAGCAGCGGCTTGGCCCCCAGCAGAATGCCCACTGGCAGTAGTTTGCGCTGTCCGAAGCCTGCCCAGAAGCGCGCCAGAGGTGGCACCAGCAGACCGATAAACGCCACC‑GGAACTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGt                                             <  2:305247/247‑1 (MQ=255)                                                                                                                                                                                             gccggCCAAC‑CACCGCCAGCACTGC‑GCCGGCGGGCAGA‑‑CTGCCGGGGAAGGCCAGTGCGCGGGCGAGTCCAGCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACC‑GGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGt                                             <  2:417939/247‑1 (MQ=255)                                                                                                                                                                                                  ccAAT‑CAGCGCCAGCACTGC‑GCCTGCGGGCAGAT‑CT‑CCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACC‑GGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGccc                                        >  2:417566/1‑250 (MQ=255)                                                                                                                                                                                                    aaT‑CAGCGCCAGCACTGG‑CCCTGCGGGCAGAT‑CT‑CCGGGGAAGGCCAGTGCGCGTGCGGGTACATCTGCCAGCAGCAGCAGCGTGGCCCCCAGCAGCAGGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCTATAAACGCCTCC‑GGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCg                                      <  2:480266/251‑1 (MQ=255)                                                                                                                                                                                                     aT‑CAGCGCCAGCACTGC‑GCCTGCGGGCAGAT‑CT‑CCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACC‑GGCCCTGCCACGCTGACGCCCGCCACAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGt                                     >  2:152084/1‑251 (MQ=255)                                                                                                                                                                                                        cagcgccagcgcTGC‑GCCTGCGGGCAGAT‑CT‑CCGGGGAAGACCAGTGCGCGTGCGAGTACATCTGCCAACAGCATCAGATTGGCCCCCAGCAGCATGCTCACTGGCAGTCCGTTGCGCTGATCGAAGCCTGCCAAGAAGCGCGCAAGATGTGGCACCAGCCGACCGATAAACGCCACC‑GGACCTCCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACACATTGCTGACCAGACGCACCCTCgcc                                    >  2:352871/1‑248 (MQ=255)                                                                                                                                                                                                              catcACTGC‑GCCTGCGGGCAGAT‑CT‑CCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCACAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACC‑GTACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCa                             <  1:170926/248‑1 (MQ=255)                                                                                                                                                                                                                 cACTGC‑GCCTGCGGGCAGAT‑CT‑CCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACC‑GGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTc                          >  2:91177/1‑251 (MQ=255)                                                                                                                                                                                                                  aCTGC‑GCCTGCGGGCAGAT‑CT‑CCGGGGAAGGCCAGTGCGCGTGCGAGTACCTCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACC‑GGACCTGCCACGCTGACGCCCGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTcc                         >  2:60617/1‑251 (MQ=255)                                                                                                                                                                                                                     gc‑gcCTGCGGGCAGAT‑CT‑CCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACC‑GGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTccc                        <  1:349777/249‑1 (MQ=255)                                                                                                                                                                                                                            gCGGGCAGAC‑CT‑CCGGGGAAGGCCAGTGCGCGTGCGAGGACCTCTGCCCGAAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCTGTCCGTTGCGCTGATCGAAGCCTGCCCAGAGGCGCGCCAGATGTGGCTCCAGCAGACCGATAAACGCCACC‑GGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCCCTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTa                  <  2:275679/249‑1 (MQ=255)                                                                                                                                                                                                                             ggtgccagtcct‑cCGGGGCAGGCCCGTGCGCGTGCGAGTACATCTGCCAACAGCCTCAGCGTGCCCCCCAGCAGCTTGCTCACTGGCAGTCCGTTGCGCTGACCGAAGCCTGCCCAGAAGCGCGCCAGAGGTGGCACCAGCAGACCGATAACCGCCACC‑GGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATg                <  2:386581/241‑1 (MQ=255)                                                                                                                                                                                                                                          ccGGGGAGGGCCCGTGCGCGTGCGAGTACATCCGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCACAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACC‑GGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTg      <  2:268333/249‑1 (MQ=255)                                                                                                                                                                                                                                          ccGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACC‑GGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTc    >  1:465644/1‑251 (MQ=255)                                                                                                                                                                                                                                          ccGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACC‑GGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTc    <  1:505107/251‑1 (MQ=255)                                                                                                                                                                                                                                           cGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCCTGCTCCCTGGCAGTACGTTGCACTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACTCCACC‑GGCCCTGCCCCGCTGACGCCCGCACCACCCAGAAGCAGCACTAACATATTGAGGACACAACGAAACCTCGTCAGGTTCCCTCCCCGCGTATGGGAGGTgaggtcc   >  2:164414/1‑245 (MQ=255)                                                                                                                                                                                                                                            ggggAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACC‑GGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGc  >  1:366981/1‑251 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                GCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAAT‑CAGCGCCAGCACTGC‑GCCTGCGGGCAGAT‑CT‑CCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACC‑GGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGC  >  NC_000913/4512258‑4512739

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