BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,116,937	T→G	81.7%	K96Q (AAG→CAG)	glpK ←	glycerol kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,116,937	0	T	G	81.7%	696.9 / 133.6	251	K96Q (AAG→CAG)	glpK	glycerol kinase
Reads supporting (aligned to +/- strand): ref base T (23/23); new base G (115/90); total (138/113)
Fisher's exact test for biased strand distribution p-value = 5.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.36e-01

Only 100 of 252 total aligned reads displayed.

GGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAG‑CTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCG  >  NC_000913/4116668‑4117197

           tCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTc                                                                                                                                                                                                                                                        <  2:90108/274‑1 (MQ=255)

           gaacaTGTACAGGGCCCTGAGTCATTTTCCAGATAAGCCAAGTATCAACCGTAGCAAACAGCAATTCACCAGGACGTGCACGCTCGCGGGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTAAGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTcc                                                                                                                                                                                                                                                       <  4:146510/273‑1 (MQ=255)

                tGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGa                                                                                                                                                                                                                                                   <  4:57302/274‑1 (MQ=255)

                                                               aCCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAt                                                                                                                                                                                                   <  2:51365/275‑1 (MQ=255)

                                                                        cAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGCTAGGCTGGCCGGGTTCTTTTTCCCACACAAATGGTGTTTCAAGCGGGTTCGTAATCCCGATATCTACACTTTGCTCGGAACtt                                                                                                                                                                                           >  4:630017/1‑273 (MQ=255)

                                                                           ttCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGGACGGGTCAATAACCAGACCGGTATTGCCGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGATTCTTTTTCACAGACAATAGTGGTTTCACGCTGGTTAGTAATACAGATAGCTGCCATTTGATAGGAACTGATCTCGGc                                                                                                                                                                                   >  1:96372/1‑279 (MQ=255)

                                                                                  cgacgTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATGGCTGAAATTTGACCGGAAATGATATCGGCtt                                                                                                                                                                                 >  4:480967/1‑274 (MQ=255)

                                                                                  cgacgTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCtt                                                                                                                                                                                 >  4:121999/1‑274 (MQ=255)

                                                                                     cgTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTACATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATACTTNNNATTTGTTAGGAACTTATATCTGCTTTCGCCAg                                                                                                                                                                          >  3:221792/1‑278 (MQ=255)

                                                                                       tGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCa                                                                                                                                                                           <  2:480967/275‑1 (MQ=255)

                                                                                         cACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGc                                                                                                                                                                         >  2:69791/1‑275 (MQ=255)

                                                                                         cACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCGGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACGGAAATCGGCTTTCCCCAGc                                                                                                                                                                         >  4:377912/1‑275 (MQ=255)

                                                                                                                 ggTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCCGCACTTTCACCAGCGTGGAA‑CCTTg                                                                                                                                                 >  2:616959/1‑274 (MQ=255)

                                                                                                                      aGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGACCTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAG‑CTTTGGGTGGCCCAg                                                                                                                                       >  1:186742/1‑279 (MQ=255)

                                                                                                                            cacttcacttTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAG‑CTTTGGGTGGCCCAGAt                                                                                                                                     <  2:317155/275‑1 (MQ=255)

                                                                                                                             acttcacttTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAG‑CTTTGGGTGGCCCAGAtt                                                                                                                                    >  2:435386/1‑275 (MQ=255)

                                                                                                                             acttcacttTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGG‑‑ATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATANGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCTATAGCTGCAATTTTATCGGAACTGNAATCGGCTTTCGCCCGCAATTTTATNAGCTGGGGG‑TTTTGGGTGGGCCAGAttttcttt                                                                                                                              >  3:134290/1‑274 (MQ=255)

                                                                                                                             acttcacttTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGG‑‑ATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCTATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAG‑TTTTGGGTGGCCCAGAt                                                                                                                                     <  1:134290/272‑1 (MQ=255)

                                                                                                                                 cacttTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCa                                                                                                                                                                           <  4:465778/233‑1 (MQ=255)

                                                                                                                                 cacttTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAG‑CTTTGGGTGGCCCAGATTTCCAtt                                                                                                                              <  1:111997/277‑1 (MQ=255)

                                                                                                                                    tttGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCtt                                                                                                                                                                                                                                                            >  5:107295/1‑149 (MQ=255)

                                                                                                                                     ttGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCtt                                                                                                                                                                                                                                                            >  7:610083/1‑148 (MQ=255)

                                                                                                                                      tGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTACACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGCCGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAG‑CTTTGGGTGGCCCAGATTTCCATTggg                                                                                                                           >  2:660767/1‑275 (MQ=255)

                                                                                                                                        gTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTc                                                                                                                                                                                                                                                        <  6:1006388/149‑1 (MQ=255)

                                                                                                                                        gTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAG‑CTTTGGGTGGCCCAGATTTCCATTGGGTc                                                                                                                         <  2:377912/275‑1 (MQ=255)

                                                                                                                                             aGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAg                                                                                                                                                                                                                                                    <  7:223554/148‑1 (MQ=255)

                                                                                                                                             aGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAg                                                                                                                                                                                                                                                    >  7:987237/1‑148 (MQ=255)

                                                                                                                                             aGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGa                                                                                                                                                                                                                                                   >  6:450360/1‑149 (MQ=255)

                                                                                                                                             aGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGa                                                                                                                                                                                                                                                   >  6:472660/1‑149 (MQ=255)

                                                                                                                                                        ggTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGt                                                                                                                                                                                                                                         >  8:255043/1‑148 (MQ=255)

                                                                                                                                                            aaTCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGa                                                                                                                                                                                                                                                   >  1:212100/1‑134 (MQ=255)

                                                                                                                                                                 ccAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTgg                                                                                                                                                                                                                               >  6:896682/1‑149 (MQ=255)

                                                                                                                                                                  cAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGt                                                                                                                                                                                                                              >  6:296093/1‑149 (MQ=255)

                                                                                                                                                                    ggccGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACCATGGCGTTATAGATAGGCAGGCCGGTTTCTTTTTCCCCGGCAATACTGGTTTNAGGCTGGTTCGTAGTACACATAGCTGCACTTTGATCGGAACTGCTATCGGCTTTAGCCCGATATTCTACCAGAGTGGGGCCTTCGAGCGGCCC‑GACTTCAGTTGGGNCnngtgc                                                                                                                   >  3:157568/3‑251 (MQ=255)

                                                                                                                                                                         gTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTg                                                                                                                                                                                                                                           >  5:379422/1‑129 (MQ=255)

                                                                                                                                                                           aTTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATAc                                                                                                                                                                                                                     >  5:31905/1‑149 (MQ=255)

                                                                                                                                                                              gctgcGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAG‑CTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTGGATTTGCTCaa                                                                                   >  2:153359/1‑275 (MQ=255)

                                                                                                                                                                                   ggATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCt                                                                                                                                                                                                             >  6:649456/1‑149 (MQ=255)

                                                                                                                                                                                      tataATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCa                                                                                                                                                                                                          >  5:87148/1‑149 (MQ=255)

                                                                                                                                                                                      tataATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCa                                                                                                                                                                                                          >  5:173287/1‑149 (MQ=255)

                                                                                                                                                                                        taATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAAt                                                                                                                                                                                                        >  6:760473/1‑149 (MQ=255)

                                                                                                                                                                                          aTCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAg                                                                                                                                                                                                                                             <  8:564017/110‑1 (MQ=255)

                                                                                                                                                                                          aTCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTgg                                                                                                                                                                                                                               <  8:228431/124‑1 (MQ=255)

                                                                                                                                                                                           tcttctAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAAttt                                                                                                                                                                                                      >  7:22456/1‑148 (MQ=255)

                                                                                                                                                                                            cttctAAACCGTTACATTTTAAATGCTCCACGGTTTCGGCGGTACGCCGGCACTGCCAGACAATGGCCTTATGGATAGGCTGGCCGGTTGCGTATTCCCAGACGATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAG‑CTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGac                                                                     <  4:310251/275‑1 (MQ=255)

                                                                                                                                                                                            cttctAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTg                                                                                                                                                                                                     <  8:5217/148‑1 (MQ=255)

                                                                                                                                                                                             ttctAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAt                                                                                                                                                                                                   >  6:291290/1‑149 (MQ=255)

                                                                                                                                                                                             ttctAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGAt                                                                                                                                                                                                   <  5:63946/149‑1 (MQ=255)

                                                                                                                                                                                              tctAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATc                                                                                                                                                                                                  <  6:932299/149‑1 (MQ=255)

                                                                                                                                                                                              tctAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATc                                                                                                                                                                                                  <  5:496396/149‑1 (MQ=255)

                                                                                                                                                                                              tctAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATc                                                                                                                                                                                                  <  6:496668/149‑1 (MQ=255)

                                                                                                                                                                                               ctAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCCCTTCTACCAGCGTGGAG‑CTTTGGGTGGCCCAGATTTCAATTGGGTCGTGTTCTACCCAACCTGGTTGTGGGTAGCTTTGCTCAAATTCGCGCttgactctt                                                                  >  2:577041/1‑267 (MQ=255)

                                                                                                                                                                                               ctAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAG‑CTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCt                                                                          >  2:224837/1‑267 (MQ=255)

                                                                                                                                                                                                   aCCGTCACGTTTTAAATGCTCGCAGATTTTTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAG‑CTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGa                                                                      <  4:224837/267‑1 (MQ=255)

                                                                                                                                                                                                     cGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACt                                                                                                                                                                                            >  7:32835/1‑148 (MQ=255)

                                                                                                                                                                                                     cGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTg                                                                                                                                                                                           >  5:1023233/1‑149 (MQ=255)

                                                                                                                                                                                                      gTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGa                                                                                                                                                                                          >  5:656382/1‑149 (MQ=255)

                                                                                                                                                                                                      gTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGa                                                                                                                                                                                          <  6:357898/149‑1 (MQ=255)

                                                                                                                                                                                                       tCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAG‑CTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGa                                                           >  2:301127/1‑274 (MQ=255)

                                                                                                                                                                                                       tCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAG‑CTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTCCGACCCGCTTAtta                                                           >  4:566/1‑272 (MQ=255)

                                                                                                                                                                                                        cACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGat                                                                                                                                                                                         >  7:977399/1‑148 (MQ=255)

                                                                                                                                                                                                           gTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCg                                                                                                                                                                                     <  6:811867/149‑1 (MQ=255)

                                                                                                                                                                                                            ttttAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCgg                                                                                                                                                                                    >  5:908745/1‑149 (MQ=255)

                                                                                                                                                                                                            ttttAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAG‑CTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTgg                                                     <  4:291624/275‑1 (MQ=255)

                                                                                                                                                                                                             tttAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTgg                                                                                                                                                                                                                                          <  7:715974/94‑1 (MQ=255)

                                                                                                                                                                                                             tttAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCgg                                                                                                                                                                                    >  7:415100/1‑148 (MQ=255)

                                                                                                                                                                                                               tAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAAc                                                                                                                                                                                             <  8:393391/137‑1 (MQ=255)

                                                                                                                                                                                                                 aaTGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTc                                                                                                                                                                               >  5:269630/1‑149 (MQ=255)

                                                                                                                                                                                                                       cGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAG‑CTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACAAGCTAATGATATTGGCATCGTGATcc                                          >  2:282884/1‑275 (MQ=255)

                                                                                                                                                                                                                        gCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGc                                                                                                                                                                         <  8:977399/148‑1 (MQ=255)

                                                                                                                                                                                                                         cAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCa                                                                                                                                                                        >  7:142434/1‑148 (MQ=255)

                                                                                                                                                                                                                                          acgGCACTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCa                                                                                                                                                                           >  8:141432/1‑128 (MQ=255)

                                                                                                                                                                                                                                           cgGCACTGCCAGACCATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAG‑CTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGc                       <  2:566/274‑1 (MQ=255)

                                                                                                                                                                                                                                            ggCACTGCCGGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAG‑Ctt                                                                                                                                                   >  5:923683/1‑149 (MQ=255)

                                                                                                                                                                                                                                               aCTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAG‑CTTTgg                                                                                                                                                <  6:551685/149‑1 (MQ=255)

                                                                                                                                                                                                                                               aCTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAG‑CTTTgg                                                                                                                                                <  6:1049064/149‑1 (MQ=255)

                                                                                                                                                                                                                                               aCTGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAG‑CTTTg                                                                                                                                                 >  8:1044363/1‑148 (MQ=255)

                                                                                                                                                                                                                                                cTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCttttt                                                                                                                                                                                                                                                         >  8:365834/1‑44 (MQ=255)

                                                                                                                                                                                                                                                cTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCttttt                                                                                                                                                                                                                                                         <  7:365834/44‑1 (MQ=255)

                                                                                                                                                                                                                                                cTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAG‑ttt                                                                                                                                                   >  4:407734/1‑145 (MQ=255)

                                                                                                                                                                                                                                                 tGCCAGACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAG‑CTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCtggtgg                 <  4:69791/274‑1 (MQ=255)

                                                                                                                                                                                                                                                    cAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAG‑TTTTggg                                                                                                                                               <  2:407734/145‑1 (MQ=255)

                                                                                                                                                                                                                                                     aGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAG‑CTTTGGGTGGcc                                                                                                                                          >  6:709814/1‑149 (MQ=255)

                                                                                                                                                                                                                                                     aGACAATGGCGTTATAGATAGGCTTGCCGGTTACTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAG‑TTTTGGGTGGc                                                                                                                                           <  8:862367/148‑1 (MQ=255)

                                                                                                                                                                                                                                                      gACAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAG‑CTTTGGGTGGcc                                                                                                                                          >  7:570286/1‑148 (MQ=255)

                                                                                                                                                                                                                                                        cAATGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAG‑CTTTGGGTGGCCCa                                                                                                                                        >  7:883538/1‑148 (MQ=255)

                                                                                                                                                                                                                                                          aTGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAG‑TTTTGGGTGGCCCAGAt                                                                                                                                     <  5:291290/149‑1 (MQ=255)

                                                                                                                                                                                                                                                          aTGGCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAG‑CTTTGGGTGGCCCAGa                                                                                                                                      >  7:292546/1‑148 (MQ=255)

                                                                                                                                                                                                                                                             gCGTTATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAG‑CTTTGGGTGGCCCAGAttc                                                                                                                                   >  8:736359/1‑147 (MQ=255)

                                                                                                                                                                                                                                                                ttATAGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAG‑CTTTGGGTGGCCCAGATTTCCAt                                                                                                                               >  6:212592/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                 tataGATAGGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAG‑CTTTGGGTGGCCCAGATTTCCAtt                                                                                                                              <  6:386611/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                  atagatagGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACtt                                                                                                                                                                     >  2:750825/1‑110 (MQ=255)

                                                                                                                                                                                                                                                                      atagGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAG‑TTTTGGGTGGCCCAGATTTCCATTGGGTc                                                                                                                         <  6:22367/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                      atagGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTAc                                                                                                                                                                 <  4:750825/110‑1 (MQ=255)

                                                                                                                                                                                                                                                                      atagGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAG‑Ctt                                                                                                                                                   >  7:1026395/1‑123 (MQ=255)

                                                                                                                                                                                                                                                                      atagGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAG‑CTTTGGGTGGCCCAGATTTCCATTGGGTc                                                                                                                         <  5:212592/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                      atagGCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAG‑CTTTGGGTGGCCCAGATTTCCATTGGGTc                                                                                                                         <  5:748198/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                          gCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAG‑CTTTGGGTGGc                                                                                                                                           >  4:152652/1‑127 (MQ=255)

                                                                                                                                                                                                                                                                          gCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAG‑CTTTGGGTGGCCCAGATTTCCATTGGGTCgtgt                                                                                                                     <  6:727228/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                          gCTGGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAG‑CTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGa                         >  4:613473/1‑241 (MQ=255)

GGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCGGATATAATCTTCTAAACCGTCACGTTTTAAATGCTCGCAGATTTCTGCGGTACGACGGCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAG‑CTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCG  >  NC_000913/4116668‑4117197

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 11 ≤` `ATCG/ATCG` `< 25 ≤` `ATCG/ATCG` `< 37 ≤` `ATCG/ATCG` `< 40 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`