Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,320,416	G→A	T83I (ACC→ATC)	phnI ←	carbon‑phosphorus lyase core complex subunit PhnI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,320,416	0	G	A	96.2%	236.4 / ‑4.2	80	T83I (ACC→ATC)	phnI	carbon‑phosphorus lyase core complex subunit PhnI
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base A (44/33);  total (46/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

CGCCAGCCCCTGACGCGCCAGCAGGCTGAAAACGTGCGGCGACGGCTGCTGTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGGCGG‑AGATACGGCGTTC‑GAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGC‑TTCGCCAACGTGGTGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAGCTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGG  >  NC_000913/4320169‑4320653                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                              ccccAGCCCCTGACGCGCCAGCAGGCTGAAAGCGTGCGGCGACGGCTGCTGTTCGCTGTCGGCGGTGGTCAGAGGCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACCACGG‑AGATACGGCTTTC‑GAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGC‑TTCGCCAACGTGATGc                                                                                                                                                                                                                                            <  1:334559/249‑1 (MQ=255)       ccccTGACGCGCCAGCAGGCTGAAAACGTGCGGCGACTGCTGCTGTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTAGTCGTAGTTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGG‑AGATACGGCGTTC‑GAGACGCATCCCGGTGGTGTCTTGCGGCTCGCTTACCGCCAGC‑TTCGCCAACGTGATGCGGTAgg                                                                                                                                                                                                                                      <  1:104918/251‑1 (MQ=255)                   cagcagGCTGAAAACGTGCGGCGACGGCTGCTGTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGG‑AGATACGGCGTTC‑GAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGC‑TTCGCCAACGTGATGCGGTAGGCGCGCAACAgg                                                                                                                                                                                                                           >  1:335974/1‑250 (MQ=255)                       agGCTGAAAACGTGCGGCGACGGCTGCTGTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCCGCAGGGTAAAATCAAGCAGGCGATGGGGGTAGTCGTAGGTTGGGGCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGG‑AGGTACGGCGTTC‑GAGTCGCATCCCGGTGGTGTCGGGCGGCCCGCTTACCGCCCGC‑TTCGCCAACGTGATGCGGTAGGGGCGCAGCAGGaaaa                                                                                                                                                                                                                       >  1:254792/1‑250 (MQ=255)                            gAAAACGTGCGGCGACGGCTGCTGTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTGGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGG‑AGATACGGCGTTC‑GAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGC‑TTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGc                                                                                                                                                                                                                   <  1:478602/249‑1 (MQ=255)                                aCGTGCGGCGACGGCTGCGGTTCGGTGGCGGCGGGGGTCAGGGGCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAAGCGAGCAGGCGATGGGTGTAGGCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGG‑AGATACGGCGTTC‑GAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGC‑TTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAAc                                                                                                                                                                                                             <  2:250112/251‑1 (MQ=255)                                   tGCGGCGCCGGCTGCGGTTCGCTGCCGGCGGGGGTCGGGGTCGGCGCTTCGCCGTTTGCCAGCAGGGAAAAATCGGGCAGGCGATGGGGGTAGTCGTAGGTGGGGCCAGGCAGCGGGCCGCCGGGAATGTCTTTATAAACGACGG‑AGATACGGCGTTC‑GAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGC‑TTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGt                                                                                                                                                                                                           <  2:445778/250‑1 (MQ=255)                                        cgACGGCTGCTGTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTAGGCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGG‑CGGTACGGCGTTC‑GAGACGCATCCCGGTGGTGTCGAGCGGCCCGCTTCCCCCCACCCTCCCCCAAC‑TAGTGCGGTAGGCGCCCAGCAGGAAAAACGCTTCACCGTTATccc                                                                                                                                                                                                     >  2:277779/1‑249 (MQ=255)                                               tgctgTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGG‑AGATACGGCGTTC‑GAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGC‑TTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGc                                                                                                                                                                                              >  1:301844/1‑251 (MQ=255)                                                 ctgTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGG‑AGATACGGCGTTC‑GAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGC‑TTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCt                                                                                                                                                                                            >  1:374277/1‑251 (MQ=255)                                                   gTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGG‑AGATACGGCGTTC‑GAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGC‑TTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGGAAATCGCTTCAACGTTATCGCCGCTGGCCTGt                                                                                                                                                                                          >  1:483489/1‑251 (MQ=255)                                                   gTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGG‑AGATACGGCGTTC‑GAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGC‑TTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGt                                                                                                                                                                                          >  1:209385/1‑251 (MQ=255)                                                      cGCTGTCGGGGGTGGGCGGCGGCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAAACGAGCGGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTAGAAACGACGG‑AGATACGGCGTTC‑GAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGC‑TTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTc                                                                                                                                                                                       <  2:320673/251‑1 (MQ=255)                                                              gcggGGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGGGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGG‑AGATACGGCGTTC‑GAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGC‑TTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTcagcgccag                                                                                                                                                                               <  2:446268/251‑1 (MQ=255)                                                                    gtCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGG‑AGATACGGCGTTC‑GAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGC‑TTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCTCCGCTGGCCTGTTTCAACGCCAGcgccgc                                                                                                                                                                         >  2:137614/1‑251 (MQ=255)                                                                    gtCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGG‑AGATACGGCGTTC‑GAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGC‑TTCCCCAACGTGATGCTGTAGGCGCCCAGCAGGAAAATCGCTTCAACAGTATCGCCGCTGGCCGGTTTTCGCGCCAGcgccgc                                                                                                                                                                         >  2:191077/1‑251 (MQ=255)                                                                                   gcGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAACCGACGG‑AGATACGGCGTTC‑GAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGC‑TTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCgg                                                                                                                                                           <  1:132534/249‑1 (MQ=255)                                                                                      ccGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGG‑CGATACGGCGTTC‑GAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGC‑TTCGCCAACGGGATGCGGTAGGCGCGCCGCAGGACAATCGCTTCACAGTTATCACCGCTGGCCTGTTTCAGACCCCGCGCCGCCAGTTCGCGGTCGGCact                                                                                                                                                       >  2:132589/1‑249 (MQ=255)                                                                                       cGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGG‑AGATACGGCGTTC‑GAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGC‑TTCGCCCACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCCATg                                                                                                                                                      >  2:119279/1‑251 (MQ=255)                                                                                        gtgtgccggcggggTAAAATCGAGCAGGCGATGGGTGTAGTCGTAGGTTGGGCAAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGG‑AGAGACGGCGTTC‑GAGCCGCATCCCGGTGGTGCCGAGCGGCTCGCTTACCGCCAGC‑TTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATgc                                                                                                                                                     <  2:313394/240‑1 (MQ=255)                                                                                        gTTTGCCGGCAGGGTAAACTCGAGCAGGCGCTGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATACACGACGG‑AGCTGCGGCGTTC‑GAGCCGCATCCCGGCGGCGTCGAACGGCTCGCTTACCGCACGC‑CTCGCCAACTCGATGCGGTAGGCGCGCGGCAAGACACTCGCTTCAACGTTCTCGGCGCTGGCCTGGTTCAGCGCCCGCGCCGCCAGTTCCCTGTAGGCAATgc                                                                                                                                                     >  2:45508/1‑251 (MQ=255)                                                                                             ccAGCAGGGTAAAATCGAGCAGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGG‑AGATACGGCGTTC‑GAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGC‑TTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCt                                                                                                                                                >  1:464886/1‑251 (MQ=255)                                                                                              cagcagGGTAAAACCGAGCAGGCGAGGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGCCGG‑AGATACGGCGTTC‑GAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTCCCGCCAGC‑TTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCtt                                                                                                                                               <  1:170550/251‑1 (MQ=255)                                                                                                gcagGGTAAAATCGAGCAGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGG‑AGATACGGCGTCC‑GAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGC‑TTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCg                                                                                                                                             <  1:55902/251‑1 (MQ=255)                                                                                                  agGGTAAAATCGAGAAGGAGGTGGGAGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGG‑AGATACGGCGTTC‑GAGACGAATCCCGGTGGTGTAGAGAGGCTCGCTTACCGCCAGC‑TTCGCCAACGTGATGTGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCCGAGCCAGCGCCGCCAGTTCGCGGTCGGGCATGCCGCCTTCggg                                                                                                                                           >  1:111540/1‑250 (MQ=255)                                                                                                  agGGGAAAATCGAGCAGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGG‑AGATACGGCGTTC‑GAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGC‑TTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGATGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCggg                                                                                                                                           >  2:452897/1‑250 (MQ=255)                                                                                                       aaaaGCGAGCGGGGGAGGGGTGGAGTCGTGGGTGGGGCGAAGCGGCGGGCCGCCGGGACGGTCTTTATAAACTACGG‑AGATACGG‑GTTCCGGGACGCATCCCGGGGGGGTCGCGCGGCTCGCTGACCGCCAGC‑TTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTCCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGtcatca                                                                                                                                      <  2:85012/251‑1 (MQ=255)                                                                                                           ccggggaggggatggggtaagtggTAGGTGGGGCCAAGCAGCGGGCGGTCGGGAAGTTCTTTAAAAACGACGG‑AGATCCGGCGTCC‑GAGCCGCACCCGGGGGGTGTCGGGCGGCCCGCTTACCGCCAGC‑TTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCAcgcg                                                                                                                                  <  2:162606/228‑1 (MQ=255)                                                                                                             gAGCAGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGG‑AGATACGGCGTTC‑GAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGC‑TTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTCCGGTCATCACGCGAt                                                                                                                                >  1:442742/1‑251 (MQ=255)                                                                                                             gAGCAGGCGAGGGGTGTAGTCGTAGGTGGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGG‑AGATACGGCGTTC‑GAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGC‑TTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGAt                                                                                                                                <  1:285431/251‑1 (MQ=255)                                                                                                              ggggggggatggggtgggTCGGAGGGTGGGCCAAGGGGCGGGCCGCCGGGAATGTCTTTATAACCGCGGG‑AGTTACGGGGCTC‑GAGACGCATCCCGGTGGTTTCGAGCGGTTCGCTTACCGCCAGC‑TTCGCCAACGTGAGGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTAGCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATc                                                                                                                               <  2:254792/234‑1 (MQ=255)                                                                                                                cAGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGG‑AGATACGGCGTTC‑GAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGC‑TTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCt                                                                                                                              >  1:373128/1‑250 (MQ=255)                                                                                                                 aGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGG‑AGATACGGCGTTC‑GAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAAC‑TTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAACTAGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCCGTTCGCGGTTGGCAAGTCCGCCCTCGGTCATCCCGCGATCTAc                                                                                                                            >  2:363935/1‑251 (MQ=255)                                                                                                                    cGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGG‑AGATACGGCGTTC‑GAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGC‑TTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACcgc                                                                                                                         >  1:237470/1‑251 (MQ=255)                                                                                                                     gATGGGTGTAGTCGTAGGTTGGGCCAAGCAGATGGCCGCCGGGAATGTCTTTATAAACGACGG‑AGATACGGCGTTC‑GAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGC‑TTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCACCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCCACCTACcgcg                                                                                                                        >  2:498710/1‑251 (MQ=255)                                                                                                                       tGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGG‑AGATACGGCGTTC‑GAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGC‑TTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACcgcg                                                                                                                        >  1:357914/1‑249 (MQ=255)                                                                                                                        gggTGTAGTCGTAGGTTGGGCCAAGCAGCGGGCCGCCGGGAATGTCTTTATAAACGACGG‑AGATACGACGTTC‑GAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGC‑TTCGCCAACGTGATGCGTTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAg                                                                                                                      <  1:52396/250‑1 (MQ=255)                                                                                                                                    aGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGG‑AGATACGGCGTTC‑GAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGC‑TTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAActgctg                                                                                                          <  1:113479/250‑1 (MQ=255)                                                                                                                                        tGGGCCAAGCATCGGGCCGCCGGGAAGTTATTTATAAACGACGG‑AGATACGGCGTTC‑GAGACGCATCCCGGTGGTGGCGAGCGGCTCGCTTACCGCTAGC‑TTCTCCAACGTGATGCGGTAGGCGTGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCTTCACGCGCGCTACCGCGAGGTTAAACTTCTGTTCa                                                                                                      <  1:366598/250‑1 (MQ=255)                                                                                                                                         gggCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGG‑AGATACGGCGTTC‑GAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGC‑TTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAAt                                                                                                    >  1:35867/1‑251 (MQ=255)                                                                                                                                              aaGCAGCTGGCCGCCGGGAATGTCTTTATAAACGCCGG‑AGATACGGCGTTC‑GAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGC‑TTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTggg                                                                                               >  1:392160/1‑251 (MQ=255)                                                                                                                                               agcagcTGGCCGCCGGGAATGTCTTTATAAACGACGG‑AGATACGGCGTTC‑GAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGC‑TTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGGCATCTACCGCGGGGTTAAACAGCTGTTCAATCTgggg                                                                                              >  2:386036/1‑250 (MQ=255)                                                                                                                                                     tGGCCGCCGGGAATGTCTTTATAAACGACGG‑AGATACGGCGTTC‑GAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGC‑TTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCt                                                                                        >  1:343039/1‑251 (MQ=255)                                                                                                                                                      ggCCGCCGGGAATGTCTTTATAAACGACGG‑AGATACGGCGTTC‑GAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGC‑TTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTc                                                                                       >  1:302226/1‑251 (MQ=255)                                                                                                                                                       gccgccGGGAATGTCTTTATAAACGACGG‑AGATACGGCGTTC‑GAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGC‑TTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCCGGCCTGTTTCAGCGACAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCCAGGTTAATCTGCTGTTCAACCTGGGCGACGccc                                                                                       >  2:42310/1‑248 (MQ=255)                                                                                                                                                         cgccgGGAATGTCTTTATAAACGACGG‑AGATACGGCGTTC‑GAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGC‑TTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCCGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCCCGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGt                                                                                    >  2:194211/1‑251 (MQ=255)                                                                                                                                                           ccgGGAATGTCTTTATAAACG‑CGGTGGGTACGGCGTCC‑GTGCCGCATCCCGGTGGTTGCGAGCGGCTCGCTTCCCGCCAGC‑TTCGCCAACGTGATGCGGTAGGCGCGAAGCAGGAAAATCGCTTCAACGTTATCGCCGCGGGCCTGTTTCAGCGCAAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTc                                                                                  <  2:501682/251‑1 (MQ=255)                                                                                                                                                              ggAATGCCTTTATAAACGACGG‑AGATACGGCGTTC‑GAGACGCATCCCGGTGGTGTCGAGCGGCACGCTTACCGCCAGC‑TTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCggg                                                                               <  1:194920/251‑1 (MQ=255)                                                                                                                                                                   tcttttATAACCGACGG‑ACATACGGCTTCC‑GAGACGCATCCCGGTGGTGCCGAGCGGCGCGCCTACCGCCAGC‑TTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCaaa                                                                           <  2:50629/247‑1 (MQ=255)                                                                                                                                                                      tttATAAACGACGG‑AGATACGGCGTTC‑GAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGC‑TTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCaaa                                                                           >  1:50629/1‑247 (MQ=255)                                                                                                                                                                       ttATAAACGGCGG‑AGATACGGCGTTC‑GAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGC‑TTCGCCAACGTGGTGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCGCGCGATCTACCGCGAGGTTAAGCTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGt                                                                      >  1:190061/1‑251 (MQ=255)                                                                                                                                                                              cgacgg‑AGATACGGCGTTC‑GAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGC‑TTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCt                                                               >  1:369583/4‑251 (MQ=255)                                                                                                                                                                              cgacgg‑AGATACGGCGTTC‑GAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGC‑TTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCAAGCGATCTACCGCGAGGTTAAACTGCTGtt                                                                                                        >  1:347868/4‑210 (MQ=255)                                                                                                                                                                              cgacgg‑AGATACGGCGTTC‑GAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGC‑TTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGACAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACAGAGAGGGTAAACTGATGTTCAATCTGGGCGACGCTCAGTTCGGGACAATCGGTATCGCCt                                                               >  1:108009/4‑251 (MQ=255)                                                                                                                                                                                gcgg‑GGATACGGCGTGC‑GGGACGCACCCCGGGGGTGTCGAGCGCCCCGCTACCCGCCAGC‑TTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCGGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCg                                                             <  2:429627/251‑1 (MQ=255)                                                                                                                                                                                 ggg‑AGATACGGCGTTC‑GAGACGCATCCCGGTGGTGTCGAGCGGCCCGCTTACCGCCAGC‑TTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCgg                                                            <  1:378936/250‑1 (MQ=255)                                                                                                                                                                                    gacaTACGGCGTTC‑AGGCCGCATCCGGGGGGTGGCGAGCGGCTCGCTTACCGCCAGC‑TTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCg                                                          <  1:125948/248‑1 (MQ=255)                                                                                                                                                                                     agaTACGGCGTTC‑GAGACGCATCCCGGTGGTGTCGCGCGGCTCGCTTACCGCCAGC‑TTCGCCAAAGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAACTCGGTAACGCCTCGGCGt                                                         >  1:216652/1‑251 (MQ=255)                                                                                                                                                                                     agaTACGGCGTTC‑GAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGC‑TTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCACCGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGGCCTCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGt                                                         >  2:286941/1‑251 (MQ=255)                                                                                                                                                                                     agaTACGGCGTTC‑GAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGC‑TTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGGCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTCAACAGCTGTTCAACATGGGCGACGCTCCGTTCGTGCAAATCGGTATAGCCCCGGCgc                                                         >  2:250247/1‑250 (MQ=255)                                                                                                                                                                                           ggCGTTC‑GAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGC‑TTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGctc                                                   >  1:248474/1‑251 (MQ=255)                                                                                                                                                                                             cGTTC‑GAGACGCATCACGGTGGTGTCGAGCGGCTCGCTTACCGCCAGC‑TTCGCCAACGTGATGCGGGAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGGCATCACGCGATCTACCGCGAGGTTCAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGctctc                                                 >  2:398048/1‑251 (MQ=255)                                                                                                                                                                                                 c‑GAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGC‑TTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTGTCTTGc                                             >  2:23456/1‑251 (MQ=255)                                                                                                                                                                                                    agaCGCCTCCCGGGGGTGTCGAGCGCCTCGCTTACCGCCAGC‑TTCGCCAACGTGATGCGTTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAGACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAg                                           <  1:424177/251‑1 (MQ=255)                                                                                                                                                                                                           ccccGGTGGTGTCGAGCGGCTCGCTTACCGCCAGC‑TTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTgg                                    <  1:386036/250‑1 (MQ=255)                                                                                                                                                                                                            cccGGTGGTGTCGAGCGGCTCGCTTACCGCCAGC‑TTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTggg                                   >  1:121354/1‑251 (MQ=255)                                                                                                                                                                                                              cGGTGGTGTCGAGCGGCTCGCTTACCGCCAGC‑TTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGggcg                                 >  1:51549/1‑251 (MQ=255)                                                                                                                                                                                                                    tgtCGAGCGGCTCGCTTACCGCCAGC‑TTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGAATCGCCCCGGCGTCGGACCTCTTGCAGGGCGTGGGCGGCGTCg                           >  2:154006/1‑251 (MQ=255)                                                                                                                                                                                                                            ggtgtcatgcccGCCAGC‑TGCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGGCGATCGCCtt                   <  2:435694/241‑1 (MQ=255)                                                                                                                                                                                                                                cGCTTACCGCCAGC‑TTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTtctc                <  1:251468/250‑1 (MQ=255)                                                                                                                                                                                                                                 gCTTACCGCCAGC‑TTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGCAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGa                                                                                                                                 >  2:410822/1‑136 (MQ=255)                                                                                                                                                                                                                                 gCTTACCGCCAGC‑TTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGc              >  1:137124/1‑251 (MQ=255)                                                                                                                                                                                                                                 gCTTACCGCCAGC‑TTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCGGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGa                                                                                                                                 <  1:410822/136‑1 (MQ=255)                                                                                                                                                                                                                                   ttACCGCCAGC‑TTCGCCAACGTGATGCGGTAGGCGCGCAGCAGTAAAATCGCTTCAACGTTAACGCCACTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGACCTACCGCGCGGTTAAACTGCTGTTCGATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATAGCCTCGGCGTCGGCTCTCTTGCAGGGCGGGGGCGGCGTCGATCGCCTTCTCGccc            >  2:148168/1‑251 (MQ=255)                                                                                                                                                                                                                                      cccccAGC‑TTCGCCAATGTGATGCGGTGGGCGCGCAGCAGGAAAACCGCTCCAACGTTACCGCCGCTGGCCTGTTGCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCAGCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCt         <  2:56985/248‑1 (MQ=255)                                                                                                                                                                                                                                      ccGCCAGC‑TTCGCAAACGGGGTGAGGGAGGCGCGCAGCAGGAAAATCGCTTCAACGGTATCGTCGCGGGCCTGTTACAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCt         <  2:45755/251‑1 (MQ=255)                                                                                                                                                                                                                                       cGCCAGC‑TTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCtt        >  2:445883/1‑251 (MQ=255)                                                                                                                                                                                                                                        gCCAGC‑TTCGCCAGCGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCACCGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCttt       <  1:479703/251‑1 (MQ=255)                                                                                                                                                                                                                                            gC‑TTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAACCGCTTCACCGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACg   <  1:392656/251‑1 (MQ=255)                                                                                                                                                                                                                                                tCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACgg  <  1:238365/249‑1 (MQ=255)                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                              CGCCAGCCCCTGACGCGCCAGCAGGCTGAAAACGTGCGGCGACGGCTGCTGTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGGCGG‑AGATACGGCGTTC‑GAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGC‑TTCGCCAACGTGGTGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAGCTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCGATCGCCTTCTCGCCCCCTTTCACGG  >  NC_000913/4320169‑4320653

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