BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,959,515	T→A	6.8%	intergenic (+70/+17)	ilvC → / ← ppiC	ketol‑acid reductoisomerase, NAD(P)‑binding/peptidyl‑prolyl cis‑trans isomerase C (rotamase C)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,959,515	0	T	A	6.8%	350.9 / 10.7	102	intergenic (+70/+17)	ilvC/ppiC	ketol‑acid reductoisomerase, NAD(P)‑binding/peptidyl‑prolyl cis‑trans isomerase C (rotamase C)
Reads supporting (aligned to +/- strand): ref base T (52/43); new base A (4/3); total (56/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.73e-01

Only 100 of 102 total aligned reads displayed.

ACGCTTGTGTGCCGTTGCTGAAACCGTTTATGGCAGAGCTGCAACCGGGCGACCTGGGTAAAGCTATTCCGGAAGGCGCGGTAGATAACGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGAATTCACCTAAATCACCGCCGCGTTTGCCTGATGGGCAAATGGAGTGTTTCTTCGCCAGCTTGCCGAAATCGGCCCCGTTCTTAATCTGCTCGAGAAGATCCAGAGCCAGTTTCTC  >  NC_000913/3959246‑3959771

aCGCTTGTGTGCCGTTGCTGAAACCGTTTATGGCAGAGCTGCAACCGGGCGACCTGGGTAAAGCTATTCCGGAAGGCGCGGTAGATAACGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGgaga                                                                                                                                                                                                                                                             <  2:214178/275‑1 (MQ=255)

                           tgaTGGCAGAGCTGCAACCGGGCGACCTGGGTAAAGCTATTCCGGAAGGCGGGGTAGATAACGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCAAGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAg                                                                                                                                                                                                                                   <  2:215084/272‑1 (MQ=255)

                                    ggctgtgaCCGGGCGATCTGTGTAAAGTTATTTCGGAAGGGGTGGTATATAACGGGGAACTGCGTGTTGTGAACGAAGCTATTCGCAGCCATGCTATTGGGCATGTATGTAAGAAACTGCCAGGCGATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATCAGTTGCGGTACAGCACCTTAat                                                                                                                                                                                                                          <  4:293510/267‑1 (MQ=255)

                                    agCTGCAACCGGGCGACCTGGGTAAAGCTATTCCGGAAGGCGCGGTAGATAACGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTCACGCACCTTAat                                                                                                                                                                                                                          >  2:475418/1‑274 (MQ=255)

                                        gCAACCGGGCGACCTGGGTAAAGCTATTCCGGAAGGCGCGGTAGATAACGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGAtgt                                                                                                                                                                                                                     <  2:356777/275‑1 (MQ=255)

                                                   aCCTGGGTAAAGCTATTCCGGAAGGCGCGGTAGATAACGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCCATATGACAGATATGAAAAGTATTGCTGTTGCGGGTTAAGTGCGGGCTGATGCCCTCACCCCGACCCTCCCCCAAAGGGAGAGGGAGAAAACACCCAAGGCCTTCCTCTGGAAAAGGCCTTGCAAGTAGTTGCGGTACAGCACCTTAAGGATGTGATATGAGAAc                                                                                                                                                                                                          >  4:422997/1‑275 (MQ=255)

                                                                  ttCCGGAAGGCGCGGTAGATAACGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACCCAAGGCCTTCACCTGGAGAAGGCCTTGCTATTAGTTGCAGTACAGCACCTTAATGATCTGATATCCGAACTAGGTGTGCAGCggg                                                                                                                                                                                           >  2:388706/1‑275 (MQ=255)

                                                                  ttCCGGAAGGCGCGGTAGATAACGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCCCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCggg                                                                                                                                                                                           >  2:323880/1‑275 (MQ=255)

                                                                        aaGGCGCGGTAGATAACGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCAAAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGaa                                                                                                                                                                                                                                                            >  2:20601/1‑204 (MQ=255)

                                                                                       aCGGGCAACGGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGa                                                                                                                                                                      <  2:423002/275‑1 (MQ=255)

                                                                                           gCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGt                                                                                                                                                                                                                                       >  2:508931/1‑206 (MQ=255)

                                                                                               cTGNNTGAGGTGAACTAAGCGATTCGCAGCAATGCGATTGAGCAGGTAGGTAANAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTGGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAg                                                                                                                                                                                                                                   <  4:508926/206‑1 (MQ=255)

                                                                                                tgctcgatggtaacgaacgGATTTGTAGCGATGCGATTGTGCAAGTAAGTAAGAAGATGCAGGGCTATATGACAGATATGAAACGTATTTGTGTGGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCACCCACAGGGAGAGGGAGAAAACTCTCGAGGCCTTTGCCTGGAGAAGACCTTGCTTTTGGTTGTGTTACAGCACCTTAATGAAGTGATATCCGAACTAGGTGTGCAGCGGACCGGTCGGCTCCAGTACCGAGCAAGAGaa                                                                                                                                                              <  4:130203/263‑1 (MQ=255)

                                                                                                        gtgAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAat                                                                                                                                                                                                                          <  2:385641/206‑1 (MQ=255)

                                                                                                                     ttCGCAGCCATTCGATTGAGCAGGTCGGTAATAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTCCCGGGTTAAGTGAGCGCTGATGCCCTCTCCCCGACCCTCCCCCACAGGGAGAGGGAGACAACACTCAAGGCCTTCTCCTGGACAATGCCTTTCTGTTAGTTGCGGTACAGCCCCTTAATGATGTGATCACCGCAATGGTTGTTCATCGGGCCGGTCCGCCCCACTTCGGGACCAAAGAAACCCAATTTATCGCACAccagc                                                                                                                                        >  4:493286/1‑272 (MQ=255)

                                                                                                                      tCGCAGCCATGCGTTTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACGCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGaa                                                                                                                                       <  2:13296/275‑1 (MQ=255)

                                                                                                                       tgCAGCCAGGGGAGTGAGAAGGGAGGTAGGAAACTGGGCGGCTATATGACGGATATTAAACGTATTGCTGTGGCGGGTTAAGGGCGGGCTCATGCCCTCAACCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGGGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGaa                                                                                                                                       <  4:475413/273‑1 (MQ=255)

                                                                                                                               tGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGaa                                                                                                                                                                                                                                                            <  5:345500/149‑1 (MQ=255)

                                                                                                                               tGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGaa                                                                                                                                                                                                                                                            <  8:439466/149‑1 (MQ=255)

                                                                                                                                  gATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGc                                                                                                                                                                                                                                                         <  8:258462/149‑1 (MQ=255)

                                                                                                                                   aTTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTg                                                                                                                                                                                                                                                                 >  7:427261/1‑140 (MQ=255)

                                                                                                                                   aTTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTg                                                                                                                                                                                                                                                                 <  8:427261/140‑1 (MQ=255)

                                                                                                                                         cAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTa                                                                                                                                                                                                                                                  <  7:470569/149‑1 (MQ=255)

                                                                                                                                                  aaGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCg                                                                                                                                                                                                                                         >  6:362601/1‑149 (MQ=255)

                                                                                                                                                     aaaCTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTa                                                                                                                                                                                                                                      >  6:202005/1‑149 (MQ=255)

                                                                                                                                                          gcgcGGCGATATGACAGATAGTAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGAATTCACCTAAATCAccg                                                                                                    <  2:118793/274‑1 (MQ=255)

                                                                                                                                                           cgcgGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCCCCTACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCAATTTATCGAACGCCGGAACCAACTGACCCCGGCGGAATTCACCTAAATAAccgcc                                                                                                  >  2:461534/1‑275 (MQ=255)

                                                                                                                                                           cgcgGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCAGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGACGGCCTTGCTATTAGTTGCGGTACAGCGCCTTAATGATGTGATATCCGAACTGGGTGTTCAGCGGGCCGGTCGGCCCCAGTACCGGACAAGAGAACACCACTTTACCGAACGCCGGAACCCACTGAACCTGGCGGACTTCACCACAATCAccgc                                                                                                   >  4:30835/1‑274 (MQ=255)

                                                                                                                                                                cTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGc                                                                                                                                                                                                                                  >  2:67604/1‑142 (MQ=255)

                                                                                                                                                                    atGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCt                                                                                                                                                                                                                              <  4:67603/142‑1 (MQ=255)

                                                                                                                                                                       aCAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGAATTCACCTAAATCACCGCCGCGTTTGCCTGa                                                                                      <  2:30835/275‑1 (MQ=255)

                                                                                                                                                                         aGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTAt                                                                                                                                                                                                                                                 >  6:395849/1‑118 (MQ=255)

                                                                                                                                                                         aGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTAt                                                                                                                                                                                                                                                 <  5:395849/118‑1 (MQ=255)

                                                                                                                                                                          gATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAatgatg                                                                                                                                                                                                                      >  8:100280/1‑144 (MQ=255)

                                                                                                                                                                          gATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAatg                                                                                                                                                                                                                         >  2:112604/1‑141 (MQ=255)

                                                                                                                                                                          gATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGc                                                                                                                                            >  4:421583/1‑218 (MQ=255)

                                                                                                                                                                          gATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGAATTCACCTAAATCACCGCCGCGTTTGCCTGATgg                                                                                   >  2:340389/1‑275 (MQ=255)

                                                                                                                                                                          gATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGAAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAatgatg                                                                                                                                                                                                                      <  7:100280/144‑1 (MQ=255)

                                                                                                                                                                              tGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGAtgt                                                                                                                                                                                                                     <  4:112603/141‑1 (MQ=255)

                                                                                                                                                                              tGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGAtgt                                                                                                                                                                                                                     <  8:238515/141‑1 (MQ=255)

                                                                                                                                                                              tGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGAtgt                                                                                                                                                                                                                     >  7:238515/1‑141 (MQ=255)

                                                                                                                                                                              tGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCg                                                                                                                                                                                                             >  6:395475/1‑149 (MQ=255)

                                                                                                                                                                              tGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGa                                                                                                                                        <  2:421588/218‑1 (MQ=255)

                                                                                                                                                                                aaaCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCAGGAGAAGGCCt                                                                                                                                                                                                                                                       <  7:85274/105‑1 (MQ=255)

                                                                                                                                                                                aaaCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCAGGAGAAGGCCt                                                                                                                                                                                                                                                       >  8:85274/1‑105 (MQ=255)

                                                                                                                                                                                    gTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTgg                                                                                                                                                                                                       <  5:294800/149‑1 (MQ=255)

                                                                                                                                                                                    gTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTgg                                                                                                                                                                                                       >  6:283064/1‑149 (MQ=255)

                                                                                                                                                                                       ttGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGAATTCACCTAAATCACCGCCGCGTTTGCCTGATGGGCAAATGGAgtgt                                                                      >  2:462569/1‑275 (MQ=255)

                                                                                                                                                                                                  ggTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGaa                                                                                                                                                                                                                                                            >  1:78708/1‑82 (MQ=255)

                                                                                                                                                                                                    ttAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGgtgt                                                                                                                                                                                                   >  4:505996/1‑137 (MQ=255)

                                                                                                                                                                                                      aaGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGcc                                                                                                                                                                                                                                                        <  3:78708/82‑1 (MQ=255)

                                                                                                                                                                                                        gTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAg                                                                                                                                                                                               <  2:506001/137‑1 (MQ=255)

                                                                                                                                                                                                             cgcTGATGCCCTCACCCCGACCCTCTCTCACAGGGAGTGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGAGGTACAGCACCTTAATGTTGTGATATTCGAACTGGGTGTGCAGCGGGCCGGTCGGCACCAGTTCTGGACCAGAGAAAACCACTTTATCGAACGCCGGAACCATCTTACCCTGGCGGAATTCACCTGAATCACCGCCCCGTTTGCCCCAGGGCCAAATGGATTGTTTCTTTCACAGCTTTcccgaa                                                 >  4:426483/1‑270 (MQ=255)

                                                                                                                                                                                                              gcTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTc                                                                                                                                                                                     <  8:410927/141‑1 (MQ=255)

                                                                                                                                                                                                              gcTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTc                                                                                                                                                                                     >  7:410927/1‑141 (MQ=255)

                                                                                                                                                                                                                tGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTa                                                                                                                                                                           >  8:297725/1‑149 (MQ=255)

                                                                                                                                                                                                                   tGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGAATTCACCTAAATCACCGCCGCGTTTGCCTGATGGGCAAATGGAGTGTTTCTTCGCCAGCTTGCCGAAATCGGc                                            >  2:475053/1‑273 (MQ=255)

                                                                                                                                                                                                                        tCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACaa                                                                                                                                                                   >  6:183993/1‑149 (MQ=255)

                                                                                                                                                                                                                          aCCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAaga                                                                                                                                                                 <  5:362601/149‑1 (MQ=255)

                                                                                                                                                                                                                          aCCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGACCg                                                                                                                                             >  3:126953/1‑169 (MQ=255)

                                                                                                                                                                                                                            cccGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAgaga                                                                                                                                                               <  7:297725/149‑1 (MQ=255)

                                                                                                                                                                                                                             ccGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGaa                                                                                                                                                              <  5:283064/149‑1 (MQ=255)

                                                                                                                                                                                                                             ccGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGaa                                                                                                                                                              >  6:29522/1‑149 (MQ=255)

                                                                                                                                                                                                                              cGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCgg                                                                                                                                         <  1:126953/169‑1 (MQ=255)

                                                                                                                                                                                                                                    tctcCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACt                                                                                                                                                       >  6:232484/1‑149 (MQ=255)

                                                                                                                                                                                                                                      tcCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGAATTCCCCTAAATCACCGCCGCGTTTCCCTGATGGGCAGATGGACTGTTTCTTCCCCCGCTGTCCACAATCGGCCCCGTTCTTACTctcctccac                       >  4:389752/1‑272 (MQ=255)

                                                                                                                                                                                                                                        ccACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGGGAAGGCCTTGCTATTGGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGAATTCACCTAAATCACCGCCGCGTTTGCCTGATGGGCAAATGGAGTGTTTCTTCGCCAGCTTGCCGAAATCGGCCCCGTTCTTAATCTGCTCGagaa                     <  2:426488/275‑1 (MQ=255)

                                                                                                                                                                                                                                        ccACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCAGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTAt                                                                                                                                                   <  5:395475/149‑1 (MQ=255)

                                                                                                                                                                                                                                         cacaGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATc                                                                                                                                                  >  6:92073/1‑149 (MQ=255)

                                                                                                                                                                                                                                          acaGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGAATTCACCTAAATCACCGCCGCGTTTGCCTGATGGGCAAATGGAGTGTCTCTTCGCCAGCTTGACGAAATCGGCTCCGTTCTTAATCTGCCCGagaaga                   >  2:262565/1‑275 (MQ=255)

                                                                                                                                                                                                                                           caGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGAATTCACCTAAATCACCGCCGCGTTTGCCTGATGGGCAAATGGAGTGTTTCTTCGCCAGCTTGCCGAAATCGGCCCCGTTCTTAATCTGCTCGagaag                    >  2:381132/1‑273 (MQ=255)

                                                                                                                                                                                                                                              ggagagggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGAATTCACCTAAATc                                                                                                        >  4:174476/1‑186 (MQ=255)

                                                                                                                                                                                                                                               gagagggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGAATTCACCTAAATCACCGCCGCGTTTGCCTGATg                                                                                    <  2:226395/205‑1 (MQ=255)

                                                                                                                                                                                                                                                agagggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGACCCATCTGACCCTGGCGGAATTCACCTAAATCACCGCCGCGTTTGCCTGATGGGCAAATGGAGTGTTTCTTACCCCGCTTGCCGAAATCGGCCCCGTTCTTAATCTTCTCgatagaagcca               >  4:458371/1‑264 (MQ=255)

                                                                                                                                                                                                                                                agagggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGAATTCACCTAAATCACCGCCGCGTTTGc                                                                                          >  2:429916/1‑198 (MQ=255)

                                                                                                                                                                                                                                                agagggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGAATTCACCTAAATCACCGCCGCGTTTGCCTGATGGGCAAATGGAGTGTTTCTTCGCCAGCTTGCCGAACTCGGCCCCGTTCTTAATCTGCTCGAGAAGATCCagc             >  2:57717/1‑274 (MQ=255)

                                                                                                                                                                                                                                                 gagggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCg                                                                                                                                          >  5:344999/1‑149 (MQ=255)

                                                                                                                                                                                                                                                  agggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCgg                                                                                                                                         >  8:351050/1‑149 (MQ=255)

                                                                                                                                                                                                                                                  agggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGAATTCACCTAAATCAccg                                                                                                    <  2:174478/186‑1 (MQ=255)

                                                                                                                                                                                                                                                      agaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAAcc                                                                                                                                     <  5:92073/149‑1 (MQ=255)

                                                                                                                                                                                                                                                       gaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCa                                                                                                                                    >  8:472357/1‑149 (MQ=255)

                                                                                                                                                                                                                                                          aaCACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCt                                                                                                                                 >  5:402512/1‑149 (MQ=255)

                                                                                                                                                                                                                                                           aatagcAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGAATTCACCTAAATCACCGCCGCGTTTGCCTGATGGGCAAATGGAGTGTTTCTTCGCCAGCTTGCCGAAATCGGCCCCGTTCTTAATCTGCTCGAGAAGATCCAGAGCCAGTTtctc  >  2:139592/6‑275 (MQ=255)

                                                                                                                                                                                                                                                             acTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGAc                                                                                                                              >  6:72652/1‑149 (MQ=255)

                                                                                                                                                                                                                                                               tCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGAccc                                                                                                                            >  6:156826/1‑149 (MQ=255)

                                                                                                                                                                                                                                                               tCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGAccc                                                                                                                            >  6:47537/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                 aaGGCCTTCTCCAGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTc                                                                                                                                                                                <  6:430336/95‑1 (MQ=255)

                                                                                                                                                                                                                                                                 aaGGCCTTCTCCAGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTc                                                                                                                                                                                >  5:430336/1‑95 (MQ=255)

                                                                                                                                                                                                                                                                  aGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTgg                                                                                                                                                                                                       >  7:88616/1‑71 (MQ=255)

                                                                                                                                                                                                                                                                  aGGCCTTCTCCAGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGc                                                                                                                                                                                          >  5:401997/1‑84 (MQ=255)

                                                                                                                                                                                                                                                                  aGGCCTTCTCCAGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGc                                                                                                                                                                                          <  6:401997/84‑1 (MQ=255)

                                                                                                                                                                                                                                                                    gCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGAtgtg                                                                                                                                                                                                                    >  7:79584/1‑56 (MQ=255)

                                                                                                                                                                                                                                                                    gCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGAtgtg                                                                                                                                                                                                                    <  8:79584/56‑1 (MQ=255)

                                                                                                                                                                                                                                                                      cTTCTCCAGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGa                                                                                                                     >  8:295105/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                       ttCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGaa                                                                                                                    <  5:156826/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                       ttCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGaa                                                                                                                    <  5:47537/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                        tctcCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACa                                                                                                                                                                    >  6:138472/1‑100 (MQ=255)

                                                                                                                                                                                                                                                                        tctcCTGGAGAAGGCCATGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACa                                                                                                                                                                    <  5:138472/100‑1 (MQ=255)

                                                                                                                                                                                                                                                                          tcCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGAATTc                                                                                                                 <  7:351050/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                            cTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGAATTCAc                                                                                                               >  5:390714/1‑149 (MQ=255)

ACGCTTGTGTGCCGTTGCTGAAACCGTTTATGGCAGAGCTGCAACCGGGCGACCTGGGTAAAGCTATTCCGGAAGGCGCGGTAGATAACGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGAATTCACCTAAATCACCGCCGCGTTTGCCTGATGGGCAAATGGAGTGTTTCTTCGCCAGCTTGCCGAAATCGGCCCCGTTCTTAATCTGCTCGAGAAGATCCAGAGCCAGTTTCTC  >  NC_000913/3959246‑3959771

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