Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	78,831	T→A	5.4%	intergenic (+32/+17)	setA → / ← leuD	broad specificity sugar efflux system/3‑isopropylmalate dehydratase small subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	78,831	0	T	A	5.4%	318.6 / 5.2	90	intergenic (+32/+17)	setA/leuD	broad specificity sugar efflux system/3‑isopropylmalate dehydratase small subunit
Reads supporting (aligned to +/- strand):  ref base T (40/45);  new base A (3/2);  total (43/47)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.37e-01

TGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCCACGTCGAAATGGATCCCCGGATTAGCTTTCACCAGCGCAAACAGTTCGTCCACTTCTGCATCGCTTAATTTCACCGGCAG  >  NC_000913/78565‑79096                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                          tGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTCCTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCCTTATTTTTAACCGCAAAGGTTAAAGTCGTTTGATGCCGTGGACGATTGCGGAAAGCCCGGGCATTTGCCCggg                                                                                                                                                                                                                                                                   >  4:367113/1‑275 (MQ=255) tGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCggg                                                                                                                                                                                                                                                                   >  2:228427/1‑275 (MQ=255)   tATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCa                                                                                                                                                                                                                                                                 >  4:384808/1‑275 (MQ=255)     tttATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTCAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAg                                                                                                                                                                                                                                                               <  2:94670/275‑1 (MQ=255)     tttATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTAAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAg                                                                                                                                                                                                                                                               <  2:427770/275‑1 (MQ=255)     tttATCGGCATTGTTGCGGGTATTGGGATGCTATAGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAg                                                                                                                                                                                                                                                               <  2:376249/275‑1 (MQ=255)                   tgCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATa                                                                                                                                                                                                                                                 <  2:478868/275‑1 (MQ=255)                            tGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGgt                                                                                                                                                                                                                                        >  4:447469/1‑275 (MQ=255)                                aTGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGa                                                                                                                                                                                                                                                                        <  4:176716/239‑1 (MQ=255)                                 tGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGttt                                                                                                                                                                                                                                   <  2:347854/275‑1 (MQ=255)                                   cTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAATGGTTAACGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCCAGGGGATTAATTCATAAACGCAGGTTGTTTTg                                                                                                                                                                                                                                 >  2:122563/1‑275 (MQ=255)                                     atcgTTTCAGGATCTAATGTCTGGAAGAGCGAGGGCAGCCACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATCCTGGCTGGCGCTATTCCGGGAGCAATTGCACAAAGTTGGTGGCACTTTGCTGTCTACTGGTTAATTTCGGTTATTGCAGTTGTCGCATTTTTTTTAACCGCATCGGTTACCGACGTTTTATGATGTGGACGATAGCGGACAGCCCGGTCCGTTGACCGGGCAAGGGGATTATTGCATATACGACGGTTNGTTTgat                                                                                                                                                                                                                               >  4:512177/4‑273 (MQ=255)                                                 aggggacggttgAAGAGGTGGGGTCGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGGAATTCTGGCCGGCGTTATTACGGGAGCAATTGCGCAAAGTGTGGAGCACTTTGCGGTCTACTGGGTAATTGCGGTTATTTCGGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGc                                                                                                                                                                                                                    <  4:311707/263‑1 (MQ=255)                                                                    ggggCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTTCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCATGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAg                                                                                                                                                                                                                        >  2:311708/1‑251 (MQ=255)                                                                        cAGCTACCACCTTATTTCATAACATGATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGt                                                                                                                                                                                            <  4:233992/275‑1 (MQ=255)                                                                        cAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGt                                                                                                                                                                                            <  4:228426/275‑1 (MQ=255)                                                                        cAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGt                                                                                                                                                                                            <  2:447470/275‑1 (MQ=255)                                                                        cAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGt                                                                                                                                                                                            <  2:367114/275‑1 (MQ=255)                                                                        cAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGt                                                                                                                                                                                            <  2:343667/275‑1 (MQ=255)                                                                        cAGCTAACTCCTTTTTTGCTAACATAATTTCTTCCGGGGGAATTCTGGCTGGGGTGTTTCAGGGAGAAATTGCAAAAAGTTGTGGGAAATTTGCTGTGTAATGGGTAATGGCGGTTATTTCGGTTGTCGAATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGt                                                                                                                                                                                            <  4:483246/275‑1 (MQ=255)                                                                         atgtACCACCTTTTTTAATAAAAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGt                                                                                                                                                                                            <  4:446065/271‑1 (MQ=255)                                                                                   ttATTTACGAACATTATTTCTACCGGGGTAATTCTTGGTGGGGATATTAAGGGGGTAAATGAAAAAAGTAGGGGCAACTTTGCTGTCTACAGGGTAATTGCGGTTATTTCTGTTTTCGCATTATTTTTAACCGCAAAGTTTAAAGACGTTTGATGACGTTGACGATAGCGGAAGGCCCGGTCATTTGACCGGGCAAGGAGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATTGCGTCGTCGTGCTGCAAGGTAAGCCCAATa                                                                                                                                                                                  <  4:507599/274‑1 (MQ=255)                                                                                   gtATTTACTATCACTATTTTTGCGCGGGTAATTTTTGCTGGCGTTTTTCAGGGCGCAATTGCACAAAGTTGGGGGCATGTTGGTTTTTCCTACGGAATTGCGGTTATTGCTGTTGTCGCATTATGTTTAGCCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATa                                                                                                                                                                                  <  4:122562/273‑1 (MQ=255)                                                                                        tACTAACATTAGGTCTACCGGGGTAAGTCTGGCTGTCGTTATTCAGGTAGCAATTGAACAAAGTTGGGGGCAAGTTGCTGTCTACTGGGTAATTGCTGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTTATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAATCGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTcc                                                                                                                                                                            <  2:96556/275‑1 (MQ=255)                                                                                           tAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGTTGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCGTAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCa                                                                                                                                                                                                >  2:248909/1‑252 (MQ=255)                                                                                           tAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTGTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGa                                                                                                                                                                         <  4:339863/275‑1 (MQ=255)                                                                                           tAACAGTAGTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGa                                                                                                                                                                         <  2:368856/275‑1 (MQ=255)                                                                                            aaCAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGAc                                                                                                                                                                        >  2:455571/1‑275 (MQ=255)                                                                                               aGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGTTGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCGTAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGt                                                                                                                                                                                            <  4:248908/252‑1 (MQ=255)                                                                                                         ccGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGAATTATTTTTAACCGCAAAGGTTAAAGAAGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGAACGGGAAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGc                                                                                                                                                            >  1:47640/1‑274 (MQ=255)                                                                                                             ggTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTgg                                                                                                                                                       <  2:384809/275‑1 (MQ=255)                                                                                                                          ggcGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTg                                                                                                                                                                                                                                                                         <  7:5041/148‑1 (MQ=255)                                                                                                                                 ttCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGt                                                                                                                                   <  4:455570/275‑1 (MQ=255)                                                                                                                                      ggAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCgg                                                                                                                                                                                                                                                                    >  4:539441/1‑141 (MQ=255)                                                                                                                                          cAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCaa                                                                                                                                                                                                                                                                <  2:539442/141‑1 (MQ=255)                                                                                                                                             ttGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAAt                                                                                                                                                                                                                                                      >  6:170521/1‑148 (MQ=255)                                                                                                                                                     aGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCtt                                                                                                                                                                                                                              >  2:257774/1‑164 (MQ=255)                                                                                                                                                         gggggCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATa                                                                                                                                                                                                                          <  4:257773/164‑1 (MQ=255)                                                                                                                                                          ggggCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCCAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGCCGCCGTGCTGCAATGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCCGCCCTGGTAACCCCGTAGGTTTtctcccacgt                                                                                                           >  4:337484/1‑268 (MQ=255)                                                                                                                                                               aCTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATTACGTGGACGATAGCGGAAAGCCCGGTCATTTTACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCCATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCATGGGCCGCGGCATGACACGATGTTAAAGCGATAGGTTTTCTCCCCCCCTTTCAc                                                                                                     >  4:171276/1‑275 (MQ=255)                                                                                                                                                               aCTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTGCTCTCCCGCTTTCAc                                                                                                     >  1:21500/1‑275 (MQ=255)                                                                                                                                                                  ttGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTg                                                                                                                                                                                                                                                                         >  6:300387/1‑108 (MQ=255)                                                                                                                                                                  ttGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTg                                                                                                                                                                                                                                                                         <  5:300387/108‑1 (MQ=255)                                                                                                                                                                           aCTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGAc                                                                                                                                                                                                                                                                       >  5:922467/1‑101 (MQ=255)                                                                                                                                                                           aCTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGAc                                                                                                                                                                                                                                                                       <  6:922467/101‑1 (MQ=255)                                                                                                                                                                           aCTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCCTAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCa                                                                                                                                                                           >  2:268093/1‑193 (MQ=255)                                                                                                                                                                            cggggTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGTGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTc                                                                                         <  2:263514/272‑1 (MQ=255)                                                                                                                                                                               ggTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGAcc                                                                                                                                                                       <  4:268092/193‑1 (MQ=255)                                                                                                                                                                                  aaTTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATc                                                                                   >  4:386976/1‑274 (MQ=255)                                                                                                                                                                                    ttGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATCGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGGGATAGGTTTTCTCTACCGCTTTCACCTCTTGCGCTTCCAGATCCa                                                                                 >  2:89587/1‑274 (MQ=255)                                                                                                                                                                                    ttGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATgg                                                                                                                                                                                                               <  5:170521/148‑1 (MQ=255)                                                                                                                                                                                    ttGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCa                                                                                                                                                                                                >  1:76801/1‑163 (MQ=255)                                                                                                                                                                                    ttGCGGTTATTTCTGTTGTCGCACTATTTTTAAACGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGCTTAATTCATAAACGTAGGTTGTTTTGCTTCATAAGTGGTATTGGCGTCGTCGTGCTGCTAGGTAAGCCACATTCTGTCCAGACCGGTCATCATGCCGTTGCGGCGGAAGGCATCACAGGTAAACGGATAGGTTATCACCCCCGCATTCCACTCTTGCGCCTcaggaaccc                                                                                 >  4:27390/1‑266 (MQ=255)                                                                                                                                                                                        ggTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGt                                                                                                                                                                                            <  3:76801/163‑1 (MQ=255)                                                                                                                                                                                               tCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCCACGTCGAAATggg                                                                     >  2:169774/1‑274 (MQ=255)                                                                                                                                                                                                    tgtCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGAAAGGCATCGATGGTAAAGAGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCTACGTCGCAATGAATccc                                                                 >  1:147590/1‑274 (MQ=255)                                                                                                                                                                                                       cGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCAtt                                                                                                                                                                                                                                                                           >  6:272988/1‑69 (MQ=255)                                                                                                                                                                                                       cGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCAtt                                                                                                                                                                                                                                                                           <  5:272988/69‑1 (MQ=255)                                                                                                                                                                                                        gCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGAt                                                                                                                                      >  2:211636/1‑201 (MQ=255)                                                                                                                                                                                                          attattTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCa                                                                                                                                                                                                                                           >  6:971239/1‑98 (MQ=255)                                                                                                                                                                                                          attattTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCa                                                                                                                                                                                                                                           <  5:971239/98‑1 (MQ=255)                                                                                                                                                                                                            tattTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTa                                                                                                                                  <  4:211635/201‑1 (MQ=255)                                                                                                                                                                                                            tattTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCAAATGACCgg                                                                                                                                                                                                                                                                    <  6:618390/71‑1 (MQ=255)                                                                                                                                                                                                            tattTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCAAATGACCgg                                                                                                                                                                                                                                                                    >  5:618390/1‑71 (MQ=255)                                                                                                                                                                                                                        cAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATaa                                                                                                                                                                                                                                                <  8:884902/79‑1 (MQ=255)                                                                                                                                                                                                                        cAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATaa                                                                                                                                                                                                                                                <  8:695662/79‑1 (MQ=255)                                                                                                                                                                                                                        cAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATaa                                                                                                                                                                                                                                                >  7:884902/1‑79 (MQ=255)                                                                                                                                                                                                                        cAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATaa                                                                                                                                                                                                                                                >  7:695662/1‑79 (MQ=255)                                                                                                                                                                                                                          aaGGTTAAAGACGTTTGATGACTTGGACGATAGAGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCAAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCCACGTCGAAATGGATCCCCGGATTAGCTTTCACCAGCGCaa                                          <  2:474285/275‑1 (MQ=255)                                                                                                                                                                                                                               tAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCCACGTCGAAATGGATCCCCGGATTAGCTTTCACCAGCGCAAACAGt                                     <  2:337485/275‑1 (MQ=255)                                                                                                                                                                                                                                    aCGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCCACGTCGAAATGGATCCCCGGATTAGCTTTCACCAGCGCAAACAGTTCGTc                                <  2:386977/275‑1 (MQ=255)                                                                                                                                                                                                                                     cGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTCTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCCACGTCGAAATGGATCCCCGGATTAGCTTTCACCAGCGCaaa                                         >  2:291645/1‑265 (MQ=255)                                                                                                                                                                                                                                       tttGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTa                                                                                                                                                                                           >  4:370490/1‑117 (MQ=255)                                                                                                                                                                                                                                       tttGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCCCCTCTTGCGCTTCCAGATCCAAGTCGAAATGGATCCCCGGATTAGCTTTCACCAGCTCAAACAGTTCGTCCAc                             >  1:113639/1‑275 (MQ=255)                                                                                                                                                                                                                                       tttGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCCACGTCGAAATGGATCCCCGGATTAGCTTTCACCAGCGCAAACAGTTCGTCCAc                             >  1:33419/1‑275 (MQ=255)                                                                                                                                                                                                                                       tttGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCCACGTCGAAATGGATCCCCGGATTAGCTTTCACAAGCGCAAACAGTTCGTCCAc                             >  2:319287/1‑275 (MQ=255)                                                                                                                                                                                                                                       gttGATGTCGTGGATGATAGCTGAAGGCCCGGTCATTTGACCGGGCCAGGGGATTAATTCATAAACGCAGGTTTGTTTGCTTCATAAGAGGCAATGGCTTCGTCGGGCTGCGAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCCACGTCGAAATGGATCCCCGGATTAGCTTTCACCAGCGCAAACAGTTCGTCCAc                             <  4:89586/274‑1 (MQ=255)                                                                                                                                                                                                                                           atgaCGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGcc                                                                                                                                                                                       <  2:370491/117‑1 (MQ=255)                                                                                                                                                                                                                                           atgaCGTGGACGATAGCGGAAAGCCCGGTCAAATGACCg                                                                                                                                                                                                                                                                     >  7:1128162/1‑39 (MQ=37)                                                                                                                                                                                                                                           atgaCGTGGACGATAGCGGAAAGCCCGGTCAAATGACCg                                                                                                                                                                                                                                                                     <  8:1128162/39‑1 (MQ=37)                                                                                                                                                                                                                                             gaCGTGGACGTTAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCAGAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCCACGTCGAAATGGATCCCCGGATTAGCTTTCACCAGCGCAAACAGTTCGTCCACTTCTGc                       <  2:512178/275‑1 (MQ=255)                                                                                                                                                                                                                                                   gACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCCACGTCGAAATGGATCCCCGGATTAGCTTTCACCAGCGCAAACAGTTCGTCCACTTCTGCATCg                   <  2:155097/273‑1 (MQ=255)                                                                                                                                                                                                                                                                aGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCCACGTCGAAATGGATCCCCGGATTAGCTTTCACCAGCGCAAACAGTTCGTCCACTTCTGCATCGCTTAATTTCACCGgc    >  2:150498/1‑275 (MQ=255)                                                                                                                                                                                                                                                                aGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGCTGTTTTGCTTCATAAGCGCTAATGGCGTCGTAGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGTCATCGCTGGTAAAGCGATAGATTTTCTCCCGCGCTTTCACTTCTTCCGCTTGCAGACCCAAGTCGAAATGTATCCCTGGTTACGCTTGCACCACCGCAAGCGTTTCTCCCAATACTGCACCGCTTAATTTCCCCgta    >  4:136568/1‑273 (MQ=255)                                                                                                                                                                                                                                                                  cccGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCCACGTCGAAATGGATCCCCGGATTAGCTTTCACCAGCGCAAACAGTTCGTCCACTTCTGCATCGCTTAATTTCACCGgcag  >  2:355335/1‑275 (MQ=255)                                                                                                                                                                                                                                                                    cGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCt                                                                                                                                                                                                                               >  6:792201/1‑52 (MQ=255)                                                                                                                                                                                                                                                                    cGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCt                                                                                                                                                                                                                               <  5:792201/52‑1 (MQ=255)                                                                                                                                                                                                                                                                    cGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAGGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGc                                                                                                                                                                                        <  7:854996/91‑1 (MQ=255)                                                                                                                                                                                                                                                                    cGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAGGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGc                                                                                                                                                                                        >  8:854996/1‑91 (MQ=255)                                                                                                                                                                                                                                                                    cGGTCAAATGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAg                                                                                                                               >  7:527078/1‑148 (MQ=255)                                                                                                                                                                                                                                                                      gTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATgg                                                                                                                                    >  6:38876/1‑141 (MQ=255)                                                                                                                                                                                                                                                                      gTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATgg                                                                                                                                    <  5:38876/141‑1 (MQ=255)                                                                                                                                                                                                                                                                           ttGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGcgtcgtcg                                                                                                                                                                                                       >  7:481922/1‑69 (MQ=255)                                                                                                                                                                                                                                                                           ttGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGcgtcgtcg                                                                                                                                                                                                       <  8:481922/69‑1 (MQ=255)                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                          TGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTTGCGCTTCCAGATCCACGTCGAAATGGATCCCCGGATTAGCTTTCACCAGCGCAAACAGTTCGTCCACTTCTGCATCGCTTAATTTCACCGGCAG  >  NC_000913/78565‑79096

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