Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,815,502	G→A	H89H (CAC→CAT)	pyrE ←	orotate phosphoribosyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,815,502	0	G	A	97.4%	275.0 / NA	79	H89H (CAC→CAT)	pyrE	orotate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (9/67);  minor base C (2/0);  total (12/67)
Fisher's exact test for biased strand distribution p-value = 1.83e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCGCGTTCCTGACGATCGAGCGAAATCAACACGCCAGCAAGCGTCGCGCCATTGGCCTGAATAATCTCCATCGACTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCTTCTTTGCGGTTAAAGCAGTACGGCAGGTCCAGGTCGTGATGCTCCGCCAGTGCCACAG‑CGGTTGTGGTGGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCGGAATCCACCAACGCTTCAGCGTAAAAACGGCCTAAC‑AGTGCC‑AGATCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGA‑TTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAA  >  NC_000913/3815255‑3815743                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       gCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGa                                                                                                                                                                                                                                                     <  1:549590/251‑1 (MQ=255)             gAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCa                                                                                                                                                                                                                                           <  1:158692/249‑1 (MQ=255)             gAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAg                                                                                                                                                                                                                                          <  1:64267/250‑1 (MQ=255)              attcttcGAAATCGCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAAGCAATACGCCAGCAATCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCTCGAATCGCCGTTCCGGCGGTTATCACATCAGCGACCAGCATTACGCGTCCTTTTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAgcg                                                                                                                                                                                                                                        <  2:705323/245‑3 (MQ=255)                  gCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGGCAAGGTCATGATGCTCCCCCAGCGCcac                                                                                                                                                                                                                                    >  1:355911/1‑251 (MQ=255)                            cgccgcGCCCCCGTTCCTGACGAGCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGt                                                                                                                                                                                                                             <  2:566388/247‑1 (MQ=255)                               cgcgCCCCCGTTCCTGACGATCAAGCGAAATCAATCCGCCAGCAAGCTGCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAgtggtc                                                                                                                                                                                                                      <  2:515153/251‑2 (MQ=255)                               cgcgCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAgtggtc                                                                                                                                                                                                                      <  1:52170/251‑2 (MQ=255)                                gcgcCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAgtggtcg                                                                                                                                                                                                                     <  2:696029/251‑3 (MQ=255)                                 cgctcccgTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTGGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAgtggtcgc                                                                                                                                                                                                                    <  1:729566/245‑4 (MQ=255)                                   cccccgTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGGTTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCaa                                                                                                                                                                                                                  <  1:351316/247‑1 (MQ=255)                                   cccccgTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGGCCTGATGCCCCGCCAGCGCCAcgg‑cgggagtgggcccc                                                                                                                                                                                                                   >  1:515487/5‑239 (MQ=255)                                    ccccgTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAAt                                                                                                                                                                                                                 <  2:240874/248‑1 (MQ=255)                                          cccTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACAGCATCTACCAGCATTTCGCGTCCTTGTAACGCGCTACCAACCAGATTTCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑At                                                                                                                                                                                                          <  2:175822/250‑1 (MQ=255)                                           ccTGACGATCAAGCGAAATCAATCCGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTGGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCTCGG‑CGGTAGTGGTCGCAATCGGG‑ATc                                                                                                                                                                                                         <  1:671277/251‑1 (MQ=255)                                             ggACGATCAAGCGAAATCAATACGCCAGCAANCGTCGCCCCATTCGCCTGAATAATCGCCATCGATTCGCGAATCGCCGTTCCGGCGGTGCTCACATCATTTTCCAGCATTACGCGACCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATccc                                                                                                                                                                                                       <  2:368143/250‑1 (MQ=255)                                                cGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTGGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCttg                                                                                                                                                                                                    <  2:487258/251‑2 (MQ=255)                                                  atcaaGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGt                                                                                                                                                                                                  <  1:345796/247‑1 (MQ=255)                                                        cGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAg                                                                                                                                                                                            <  1:221125/251‑1 (MQ=255)                                                        cGAAATCAAAACGCCAGCAAGCGGCGCGCCATTAGCCTGAATAATCTCCAGCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTACCGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGAGGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAg                                                                                                                                                                                            <  2:475513/251‑1 (MQ=255)                                                          aaaTCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAggg                                                                                                                                                                                          <  1:267607/251‑1 (MQ=255)                                                           aaacaatacGCCAGCAAGCGTCGCGCCATTANCCTGAATAATCTCCATCGATGCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTCCCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCGCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGc                                                                                                                                                                                         <  2:350629/244‑1 (MQ=255)                                                            aTCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGcc                                                                                                                                                                                        <  1:145037/251‑1 (MQ=255)                                                                     ccAGCAAGCGTCGCGCCAGTAGCCTTAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTCCCAGCATGACGCGTCCTTGTAACGCGCCACCAACCAGATTGCCGCCTTCACCGTTGTCTTTTGCCTCTTTGCGGGTAAAGCAGCACGGCAGTTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAAcagcag                                                                                                                                                                               <  2:244346/251‑1 (MQ=255)                                                                         cAAGCGTCGCGCCATTAGCCTTAATAATCTCCATCGATTCGCGAAGCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCg                                                                                                                                                                           <  2:577760/251‑1 (MQ=255)                                                                           aGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGGCAAGGTCATGATTCTCCGCCCGCGCCACGG‑CGGGAGTGGTCGCAAACGGG‑GTCCCCTTGTAAGCAGGGCCaaaaagacgaacga                                                                                                                                                                          >  1:84872/1‑239 (MQ=255)                                                                                 gcgcCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCGTTTGCCTCTTTGGGGTTAAAGCAGTACGTCAGGTCAAGGGCATGATGCTCCGCCAGCGCGTCGG‑CGGTAGTTGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAAt                                                                                                                                                                   <  1:435270/251‑1 (MQ=255)                                                                                   gcCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCAGTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGc                                                                                                                                                                 <  1:387123/251‑1 (MQ=255)                                                                                    ccATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCTTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGcc                                                                                                                                                                <  1:708632/251‑1 (MQ=255)                                                                                         agCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCCTGATGCTCCGCCCGCGCC‑CAGGCGGGTGTGGTCGCAAACGGGGATCC‑TTTGTAAGCCGGGGCAAAcacccggtcggaatcaatgcccaaa                                                                                                                                                            >  1:50363/2‑227 (MQ=255)                                                                                                ataatCTCCACCGATGCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCCTTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCGGCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGTCCAATCAGCAGATCGAACTCAATGCCAGAATCCTCCaa                                                                                                                                                    <  2:212675/251‑1 (MQ=255)                                                                                                       ccATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCt                                                                                                                                                <  1:517035/248‑1 (MQ=255)                                                                                                           cgattCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGt                                                                                                                                         <  1:516382/247‑1 (MQ=255)                                                                                                                       cGCCTTTCCGGCGGGGATCACAGCATCTACCAGCATTACGCGTCCTTTTAACGCGCTACCAACCAGATTGCCGCCTTCTCCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCCGCGCCGCGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCt                                                                                                                              <  2:256392/250‑1 (MQ=255)                                                                                                                       cGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTa                                                                                                                             <  1:172190/251‑1 (MQ=255)                                                                                                                        gCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTaa                                                                                                                            <  1:569429/251‑1 (MQ=255)                                                                                                                              ccGGCTGTGATCCCATCATCTACCAGCATTACGCGGCCTTGTAACGCGCGACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCTCCAATGCTTCAGCGTAAAAACGGCCTAAC‑AGTGc                                                                                                                     <  2:654657/251‑1 (MQ=255)                                                                                                                              ccGGCGGTGATCACATCATCTACCAGCATTANGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAAC‑AGTGc                                                                                                                     <  2:355911/251‑1 (MQ=255)                                                                                                                                    gTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGGAGTGGTCGCCATCCGG‑ATCCCTTTGTAAGCCGGGCCCAACAGCCGATCGAACCCCATGCCAGAAACCACCAATGCTTCAGCGTAAAAACGGCCCAAC‑AGTGcccaggc                                                                                                               >  1:377878/1‑246 (MQ=255)                                                                                                                                      gATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGcc                                                                                                                                                                <  1:151135/201‑1 (MQ=255)                                                                                                                                         cccatcatCTAACAGCATTACGCGTCCTTGTAACGCGCTACTAACCAGAGTGCCGCCTTCACCGTGGTCTTTCGCCGCTTGGCGGTTAAAGAAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCGCCAATGCTTCAGCGTAAAAACGGCCTAAC‑ATTGCC‑AGGTCGCGcc                                                                                                         <  2:331598/249‑1 (MQ=255)                                                                                                                                         cacaTCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAAC‑AGTGCC‑AGGTCgcgc                                                                                                          <  1:461766/250‑1 (MQ=255)                                                                                                                                          ccatcatCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAAC‑AGTGCC‑AGGTCGCGccc                                                                                                        <  2:536933/250‑1 (MQ=255)                                                                                                                                                      aGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCTTGATGCTCCGCCAGCGCCACGG‑CGGTAGGGGGCGCCATCGGG‑GTCCCTTTGTTAGCAAGGGCAAACCGCAGATCGCACTCCATGCCAGAATCCCCCAATGCTTCCGCGTAAAAACCGCCTAACCAG‑GCC‑AAGTCGCGCCCCGGATTAAAc                                                                                              >  1:706601/1‑249 (MQ=255)                                                                                                                                                       gCATTTCGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAAC‑AGTGCC‑AGGTCGCGCCCGGTATTAAACAGc                                                                                           <  2:84872/251‑1 (MQ=255)                                                                                                                                                           tACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAAC‑AGTGCC‑AGGTCGCGCCCGGTATTAAACAGCCCg                                                                                        <  1:126447/250‑1 (MQ=255)                                                                                                                                                               cgTCCTTGTACCGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCAGGATGCGCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAAC‑AGTGCC‑AGGTCGCGCCCGGTATTAAACAGCCCGGCGtt                                                                                   <  2:515487/251‑1 (MQ=255)                                                                                                                                                               cgTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCCGCGCCACGG‑CGGGAGTGGTCGCCATCCGG‑AACCCTTTGTAAGCCGGGGCCAACAGCAGATCGAACTCCATGCCAGAATCCACCAAAGCTTCCGCGGAAAAAACGCCTAAA‑AATGCC‑CGGTAGCGCCCGGCATTAAACAGCCCGGCGtt                                                                                   >  2:66617/1‑251 (MQ=255)                                                                                                                                                               cgTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAAC‑AGTGCC‑AGGTCGCGCCCGGTATTAAACAGCCCGGCGtt                                                                                   <  2:3697/251‑1 (MQ=255)                                                                                                                                                                ggccTTGTAACGCGCTACCAACCAGATTGCCNCCTTCACCGGGGTCTTTTGCCTCTGTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCAACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGATCTCAATGCCAGAATCCACCAATGCGTCAGCGTAAAAACGGCCTAAC‑AGTGCC‑AGGTCGCGCCCTGTATTAAACAGCCCGGCGTTg                                                                                  <  2:342186/249‑1 (MQ=255)                                                                                                                                                                  cctcgTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAAC‑AGTGCC‑AGGTCGCGCCCGGTATTAAACAGCCCGGCGTTgaa                                                                                <  1:683579/247‑1 (MQ=255)                                                                                                                                                                             gcTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAAC‑AGTGCC‑AGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGc                                                                     <  1:348279/251‑1 (MQ=255)                                                                                                                                                                                   aaccaGATGGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAAC‑AGTGCC‑AGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTgc                                                               <  2:307266/251‑1 (MQ=255)                                                                                                                                                                                      caGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGGGGTCGCAATCGGG‑ATCCCTTTTTAAGCCGGGCCAAACCGCAGAACGAACTCCATGCCAGAATCCACCCATGCTTCAGCGTAAAAACGGCCCAAC‑AATGCC‑CGGGCGCGCCCGGTATTAAACAGCCCCGCGGTGAAGAAATAAGGGCTTTTgcggc                                                            >  1:654874/1‑249 (MQ=255)                                                                                                                                                                                       aGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAAC‑AGTGCC‑AGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGccc                                                           <  2:763100/251‑1 (MQ=255)                                                                                                                                                                                            gccgccTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAAC‑AGTGCC‑AGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGA‑tt                                                     <  1:409945/251‑1 (MQ=255)                                                                                                                                                                                            gccgccTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCCAAC‑AGTGCC‑AGGTCGCGCCCGGTATTAACCAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCAGGA‑tt                                                     <  1:449889/251‑1 (MQ=255)                                                                                                                                                                                              cgccTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGGGGGCGCCATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCAACAATGCTTCAGCGTAAAAACGGCCTAAC‑AGTGCC‑AGGGCGCGCCCGGGATTAAAAAGCCCGGCGTTGGAGAAAAAAGGGCTTTTGCGCCCGGG‑tttt                                                   >  1:532983/1‑250 (MQ=255)                                                                                                                                                                                                   tCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAAC‑AGTGCC‑AGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGA‑TTTCAGc                                                <  1:12760/249‑1 (MQ=255)                                                                                                                                                                                                     aCCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGAAGTGGTCGCCATACGG‑ATCCCTTTGTAAGAAGGGGCAAACAGCAGAACGAACTCAAAGCCAGAAACCACCCATGCTTCAGCGTAAAAACGGGCTAAC‑AATGCCCAG‑GCGCGCCCGGTATTAAACAGACCGACGTGTAAGAAATAGGGGCTTTTTCGCCCGGAGTTTCacgta                                              >  1:555646/1‑246 (MQ=255)                                                                                                                                                                                                      ccGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAAC‑AGTGCC‑AGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGA‑TTTCAGCGTaaa                                           <  1:254614/251‑1 (MQ=255)                                                                                                                                                                                                        gTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCATAATCCACCAATTCTTCAGCGTAAAAACGGCCTAAC‑AGTGCC‑AGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGA‑TTTCAGCGTAAACt                                         <  1:287802/251‑1 (MQ=255)                                                                                                                                                                                                        gTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGTCCTAAC‑AGTGCC‑AGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGA‑TTTCAGCGTAAACt                                         <  1:111336/251‑1 (MQ=255)                                                                                                                                                                                                          ggTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAAC‑AGTGCC‑AGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGA‑TTTCAGCGTAAACTc                                        <  1:380240/250‑1 (MQ=255)                                                                                                                                                                                                          ggTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAAC‑AGTGCC‑AGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGA‑TTTCAGCGTAAACTc                                        <  1:496505/250‑1 (MQ=255)                                                                                                                                                                                                          ggTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAAC‑AGTGCC‑AGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGA‑TTTCAGCGTAAACTCg                                       <  1:326281/251‑1 (MQ=255)                                                                                                                                                                                                          ggTCTTTTGCCGCTTTGCGGTTAAAGCAGGNCGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAAC‑AGTGCC‑AGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGA‑TTTCAGCGTAAACTc                                        <  2:377878/250‑1 (MQ=255)                                                                                                                                                                                                             cctgtgactttttGCGGTTAACGCTGTACGNGAGGTTGAGGTCATGATGCTCCGCCAGCGCTACGG‑AGGTAGTCGTCGCAATCGGG‑ATCCCTTGGTAAGCTGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCGCCAATGCTTCAGCGTAAAAACGGCCTAAC‑AGTGCC‑AGGTCGCGCCCGGTATTATGCAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGA‑TTTAAGCGTAAAATCGcc                                     <  2:343210/240‑1 (MQ=255)                                                                                                                                                                                                              ttttGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAAC‑AGTGCC‑AGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGA‑TTTCAGCGTAAACTCGCCa                                    <  2:189182/250‑1 (MQ=255)                                                                                                                                                                                                               tttGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAAC‑AGTGCC‑AGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGA‑TTTCAGCGTAAACTCGCCa                                    <  1:713375/249‑1 (MQ=255)                                                                                                                                                                                                                   cctcttTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAAc                                                                                                                           <  2:70603/159‑1 (MQ=255)                                                                                                                                                                                                                        ttGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAAC‑AGTGCC‑AGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGA‑TTTCAGCGTAAACTCGCCAAACTTTaa                            <  1:430913/248‑1 (MQ=255)                                                                                                                                                                                                                            ggTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAAC‑AGTGCC‑AGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGA‑TTTCAGCGTAAACTCGCCAAACTTTAACACCtgc                     <  1:188804/251‑1 (MQ=255)                                                                                                                                                                                                                             gTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCCCGG‑CGGTAGTGGTCGCTATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGACCTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAAC‑AGTGCC‑AGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGA‑TTTCAGCGTAAACTCGCCAAACTTTAACACCtgct                    <  1:383456/251‑1 (MQ=255)                                                                                                                                                                                                                                 aaGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAAC‑AGTGCC‑AGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGA‑TTTCAGCGTAAACTCGCCAAACTTTAACACCtgcttgct                <  1:665828/251‑1 (MQ=255)                                                                                                                                                                                                                                          ggCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGAACGAACTCAATGCCAGAATCCAACAATGCTTCAGCGTAAAAACGGCCTAAC‑AGTGCC‑AGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGA‑TTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCaa        >  1:273190/1‑250 (MQ=255)                                                                                                                                                                                                                                           gCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAAC‑AGTGCC‑AGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGA‑TTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAAt      <  1:535397/251‑1 (MQ=255)                                                                                                                                                                                                                                           gCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAAC‑AGTGCC‑AGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGA‑TTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAAt      <  1:59480/251‑1 (MQ=255)                                                                                                                                                                                                                                             aggtcaaggtcaTGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAAC‑AGTGCC‑AGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGA‑TTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTc    <  1:368407/251‑1 (MQ=255)                                                                                                                                                                                                                                                tcaaggtcATGATGCTCCGCCAGCGCCACGG‑CGGTAGTGGTCGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTAAATGCCAGAATCCACCAATGCTTCAGCGTAAAAACGGCCTAAC‑AGTGCC‑AGGTCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGA‑TTTCAGCGTAAAATCCCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCaa  >  1:719131/4‑250 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       GCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCGCGTTCCTGACGATCGAGCGAAATCAACACGCCAGCAAGCGTCGCGCCATTGGCCTGAATAATCTCCATCGACTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCTTCTTTGCGGTTAAAGCAGTACGGCAGGTCCAGGTCGTGATGCTCCGCCAGTGCCACAG‑CGGTTGTGGTGGCAATCGGG‑ATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCGGAATCCACCAACGCTTCAGCGTAAAAACGGCCTAAC‑AGTGCC‑AGATCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGA‑TTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAA  >  NC_000913/3815255‑3815743

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