BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,258,781	C→T	4.4%	V172I (GTC→ATC)	psuK ←	pseudouridine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,258,781	0	C	T	4.4%	313.2 / 2.5	90	V172I (GTC→ATC)	psuK	pseudouridine kinase
Reads supporting (aligned to +/- strand): ref base C (48/38); new base T (2/2); total (50/40)
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.05e-01

ACACGAAGCAAGTCCCGCCATCATGGCATCGCCCGCTCCGGTAACATTAATAACATTGGTTTTGATCGGCGCAGACCAGCCACTTTCACCGCTGATATCGCTGTAATAAACGCCGTCGCCGCCCATGCTCAATACCAGTCGGTTCAGGCCATGTTGATGGAACCAGGCAGCAACTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTATCCAGAATCCATGCCAGTGCCTCTTCACTGATATTACAGTCCGCGACAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCGAGATATTCAGCTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCAACCAGCATTTCACCGGTATTATCGAGTAATGATAAATAACTCGACGTATTTTCTCCCGGCACAATCAGGCATTTATCGACATAAACGCCAGAT  >  NC_000913/2258500‑2259061

acacGAAGCAAGTCCCGCCATCATGGCATCGCCCGCTCCGGTAACATTAATAACATTGGTTTTGATCGGCGCAGACCAGCCACTTTCACCGCTGATATCGCTGTAATAAACGCCGTCGCCGCCCATGCTCAATACCAGTCGGTTCAGGCCATGTTGATGGAACCAGGCAGCAACTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTc                                                                                                                                                                                                                                                                                   <  1:616063/289‑1 (MQ=255)

      aGCAAGTCCCGCCATCATGGCATCGCCCGCTCCGGTAACATTAATAACATTGGTTTTGATCGGCGCAGACCAGCCACTTTCACCGCTGATATCGCTGTAATAAACGCCGTCGCCGCCCATGCTCAATACCAGTCGGTTCAGGCCATGTTGATGGAACCAGGCAGCAACTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCCGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGcc                                                                                                                                                                                                                                                                             >  2:418117/1‑289 (MQ=255)

      aGCAAGTCCCGCCATCATGGCATCGCCCGCTCCGGTAACATTAATAACATTGGTTTTGATCGGCGCAGACCAGCCACTTTCACCGCTGATATCGCTGTAATAAACGCCGTCGCCGCCCATGCTCAATACCAGTCGGTTCAGGCCATGTTGATGGAACCAGGCAGCAACTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGcc                                                                                                                                                                                                                                                                             >  2:787848/1‑289 (MQ=255)

      aGCAAGTCCCGCCATCATGGCATCGCCCGCTCCGGTAACATTAATAACATTGGTTTTGATCGGCGCAGACCAGCCACTTTCACCGCTGATATCGCTGTAATAAACGCCGTCGCCGCCCATGCTCAATACCAGTCGGTTCAGGCCATGTTGATGGAACCAGGCAGCAACTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACCCCGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGCGTGTGGCTCTGATTTAGACGGCCCTGCACTTTGACACATTTCCATGc                                                                                                                                                                                                                                                                              >  1:315794/1‑288 (MQ=255)

               cGCCATCATGGCATCGCCCGCTCCGGTAACATTAATAACATTGGTTTTGATCGGCGCAGACCAGCCACTTTCACCGCTGATATCGCTGTAATAAACGCCGTCGCCGCCCATGCTCAATACCAGTCGGTTCAGGCCATGTTGATGGAACCAGGCAGCAACTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGa                                                                                                                                                                                                                                                                    <  1:383387/289‑1 (MQ=255)

                gCCATCATGGCATCGCCCGCTCCGGTAACATTAATAACATTGGTTTTGATCGGCGCAGACCAGCCACTTTCACCGCTGATATCGCTGTAATAAACGCCGTCGCCGCCCATGCTCAATACCAGTCGGTTCAGGCCATGTTGATGGAACCAGGCAGCAACTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGa                                                                                                                                                                                                                                                                    >  2:727840/1‑288 (MQ=255)

                     catGGCATCGCCCGCTCCGGTAACATTAATAACATTGGTTTTGATCGGCGCAGACCAGCCACTTTCACCGCTGATATCGCTGTAATAAACGCCGTCGCCGCCCATGCTCAATACCAGTCGGTTCAGGCCATGTTGATGGAACCAGGCAGCAACTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACa                                                                                                                                                                                                                                                              <  2:230699/289‑1 (MQ=255)

                               cccGCTCCGGTAACATTAATAACATTGGTTTTGATCGGCGCAGACCAGCCACTTTCACCGCTGATATCGCTGTAATAAACGCCGTCGCCGCCCATGCTCAATACCAGTCGGTTCAGGCCATGTTGATGGAACCAGGCAGCAACTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGGTCTGATTTAGACGGTCGCGCACTTTGACACATTCCCACGCCGAAACCGGCTCGACAAATACGGGt                                                                                                                                                                                                                                                     >  2:239308/1‑288 (MQ=255)

                                    tCCGGTAACATTAATAACATTGGTTTTGATCGGCGCAGACCAGCCACTTTCACCGCTGATATCGCTGTAATAAACGCCGTCGCCGCCCATGCTCAATACCAGTCGGTTCAGGCCATGTTGATGGAACCAGGCAGCAACTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTCGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTg                                                                                                                                                                                                                                               >  2:361357/1‑289 (MQ=255)

                                     ccGGTAACATTAATAACATTGGTTTTGATCGGCGCAGACCAGCCACTTTCACCGCTGATATCGCTGTAATAAACGCCGTCGCCGCCCATGCTCAATACCAGTCGGTTCAGGCCATGTTGATGGAACCAGGCAGCAACTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTgg                                                                                                                                                                                                                                              <  1:37746/289‑1 (MQ=255)

                                        gTAACATTAATAACATTGGTTTTGATCGGCGCAGACCAGCCACTTTCACCGCTGATATCGCTGTAATAAACGCCGTCGCCGCCCATGCTCAATACCAGTCGGTTCAGGCCATGTTGATGGAACCAGGCAGCAACTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGATACATTTCCATGCCGAAACCGGCTCGACAAATACGGGTACGTTggcg                                                                                                                                                                                                                                            >  1:524880/1‑288 (MQ=255)

                                          aaCATTAATAACATTGGTTTTGATCGGCGCAGACCAGCCACTTTCACCGCTGATATCGCTGTAATAAACGCCGTCGCCGCCCATGCTCAATACCAGTCGGTTCAGGCCATGTTGATGGAACCAGGCAGCAACTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCa                                                                                                                                                                                                                                         <  2:139870/289‑1 (MQ=255)

                                          aaCATTAATAACATTAGTGTGGGTCGTGGCCGACCAGCCCCTGTACCCGCTGATATCGCTGTAATAAACGCCGTCGCCGCCGATGCTCAATACCAGTCGGTTCAGGCCATGTTGATGGAACCAGGCAGCAACTTTTGCCACGTCTTCACGCCCTGACAGCGCAATCCCACGCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCa                                                                                                                                                                                                                                         <  1:457747/289‑1 (MQ=255)

                                              ttAATAACATTGGTTTTGATCGGCGCAGACCAGCCACTTTCACCGCTGATATCGCTGTAATAAACGCCGTCGCCGCCCATGCTCAATACCAGTCGGTTCAGGCCATGTTGATGGAACCAGGCAGCAACTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTa                                                                                                                                                                                                                                      >  1:116914/1‑288 (MQ=255)

                                              ttAATAACATTGGTTTTGATCGGCGCAGACCAGCCACTTTCACCGCTGATATCGCTGTAATAAACGCCGTCGCCGCCCATGCTCAATACCAGTCGGTTCAGGCCATGTTGATGGAACCAGGCAGCAACTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTa                                                                                                                                                                                                                                      >  2:154030/1‑288 (MQ=255)

                                              ttAATAACATTGGTTTTGATCGGCGCAGACCAGCCACTTTCACCGCTGATATCGCTGTAATAAACGCCGTCGCCGCCCATGCTCAATACCAGTCGGTTCAGGCCATGTTGATGGAACCAGGCAGCAACTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGACCCACAAATACGGGTACGTTGGCGGCAttta                                                                                                                                                                                                                                     >  1:11407/1‑287 (MQ=255)

                                                       ttGGTTTTGATCGGCGCAGACCAGCCACTTTCACCGCTGATATCGCTGTAATAAACGCCGTCGCCGCCCATGCTCAATACCAGTCGGTTCAGGCCATGTTGATGGAACCAGGCAGCAACTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTATCCAGAATcc                                                                                                                                                                                                                            >  1:321765/1‑289 (MQ=255)

                                                       ttGGTTTTGATCGGCGCAGACCAGCCACTTTCACCGCTGATATCGCTGTAATAAACGCCGTCGCCGCCCATGCTCAATACCAGTCGGTTCAGGCCATGTTGATGGAACCAGGCAGCAACTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTATCCAGAATcc                                                                                                                                                                                                                            >  2:245189/1‑289 (MQ=255)

                                                       ttGGTTTTGATCGGCGCAGACCAGCCACTTTCACCGCTGATATCGCTGTAATAAACGCCGTCGCCGCCCATGCTCAATACCAGTCGGTTCAGGCCATGTTGATGGAACCAGGCAGCAACTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTATCCAGAATc                                                                                                                                                                                                                             >  1:27996/1‑288 (MQ=255)

                                                                        aGACCAGCCACTTTCACCGCTGATATCGCTGTAATAAACGCCGTCGCCGCCCATGCTCAATACCAGTCGGTTCAGGCCATGTTGATGGAACCAGGCAGCAACTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTATCCAGAATCCATGCCAGTGCCTCTTc                                                                                                                                                                                                            <  2:153887/288‑1 (MQ=255)

                                                                           ccAGCCACTTTCACCGCTGATATCGCTGTAATAAACGCCGTCGCCGCCCATGCTCAATACCAGTCGGTTCAGGCCATGTTGATGGAACCAGGCAGCAACTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCTCTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTATCCAGAATCCATGCCAGTGCCTCTTCAc                                                                                                                                                                                                          <  2:309409/287‑1 (MQ=255)

                                                                               ccACTTTCACCGCTGATATCGCTGTAATAAACGCCGTCGCCGCCCATGCTCAATACCAGTCGGTTCAGGCCATGTTGATGGAACCAGGCAGCAACTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTATCCAGAATCCATGCCAGTGCCTCTTCAc                                                                                                                                                                                                          >  2:166925/1‑283 (MQ=255)

                                                                               ccACTTTCACCGCTGATATCGCTGTAATAAACGCCGTCGCCGCCCATGCTCAATACCAGTCGGTTCAGGCCATGTTGATGGAACCAGGCAGCAACTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTATCCAGAATCCATGCCAGTGCCTCTTCAc                                                                                                                                                                                                          <  1:166925/283‑1 (MQ=255)

                                                                               ccACTTTCACCGCTGATATCGCTGTAATAAACGCCGTCGCCGCCCATGCTCAATACCAGTCGGTTCAGGCCATGTTGATGGAACCAGGCAGCAACTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTATCCAGAATCCATGCCAGTGCCTCTTCACTGatat                                                                                                                                                                                                    >  1:450093/1‑289 (MQ=255)

                                                                               ccACTTTCACCGCTGATATCGCTGTAATAAACGCCGTCGCCGCCCATGCTCAATACCAGTCGGTTCAGGCCATGTTGATGGAACCAGGCAGCAACTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTATCCAGAATCCATGCCAGTGCCTCTTCACTGatat                                                                                                                                                                                                    >  2:404355/1‑289 (MQ=255)

                                                                                 aCTTTCACCGCTGATATCGCTGTAATAAACGCCGTCGCCGCCCATGCTCAATACCAGTCGGTTCAGGCCATGTTGATGGAACCAGGCAGCAACTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTATCCAGAATCCATGCCAGTGCCTCTTCACTGATATTa                                                                                                                                                                                                  <  2:207225/289‑1 (MQ=255)

                                                                                    ttCACCGCTGATATCGCTGTAATAAACGCCGTCGCCGCCCATGCTCAATACCAGTCGGTTCAGGCCATGTTGATGGAACCAGGCAGCAACTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGATACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTATCCAGAATCCATGCCAGTGCCTCTTCACTGATATTACAg                                                                                                                                                                                               >  2:769013/1‑289 (MQ=255)

                                                                                        ccGCTGATTATGCTTTAATGAACGCCATAGACGCCCATGCTAAATACCAGTCGGTTCAGGCCATGTTGATGGAACCAGGCAGCAACTTTTGCCACATCTTCACGCCCTGACATCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCCGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTATCCAGAATCCATGCCAGTGCCTCTTCACTGATATTACAGTCcg                                                                                                                                                                                           <  1:289084/289‑1 (MQ=255)

                                                                                        ccGCTGAGGTCGCTGTAATAAACGCCGTCGACGCCCATGCTCAATACCAGTCGGTTCAGGCCATGTTGATGGAACCAGGCAGCAACTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTATCCAGAATCCATGCCAGTGCCTCTTCACTGATATTACAGTCcg                                                                                                                                                                                           <  1:154030/289‑1 (MQ=255)

                                                                                            tGATATCGCTGTAATAAACGCCGTCGCCGCCCATGCTCAATACCAGTCGGTTCAGGCCATGTTGATGGAACCAGGCAGCAACTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCCCTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAACTACGGGTACGTTGGCGGCATTATCCAGCATTCATCTCTGTGCCTCTCTACCTATACTACAGTctgaga                                                                                                                                                                                        >  1:23888/1‑283 (MQ=255)

                                                                                                   gCTGTAATAAACGCCGTCGCCGCCCATGCTCAATACCAGTCGGTTCAGGCCATGTTGATGGAACCAGGCAGCAACTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTATCCAGAATCCATGCCAGTGCCTCTTCACTGATATTACAGTCCGCGACAATGa                                                                                                                                                                                  >  1:862962/1‑287 (MQ=255)

                                                                                                     tGTAATAAACGCCGTCGCCGCCCATGCTCAATACCAGTCGGTTCAGGCCATGTTGATGGAACCAGGCAGCAACTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCCAGGCGGTTCGGCTTGAGCGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGCTCGACAAATACGGGTACGTTGGCGGCATTATCCAGAATCCATGCCAGTGCCTCTTCACTGATATTACAGTCCGCGAGCATGACttt                                                                                                                                                                              >  1:667588/1‑289 (MQ=255)

                                                                                                        aacaaccGCCGTCGCCCCCGATGCTAAATACCAGTCGGTTCAGGCCATGTTGATGGAACCAGGCAGCAACTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTATCCAGAATCCATGCCAGTGCCTCTTCACTGATATTACAGTCCGCGACAATGACCTTTg                                                                                                                                                                            <  1:626807/282‑1 (MQ=255)

                                                                                                               gCCGTCGCCGCCCATGCTCAATACCAGTCGGTTCAGGCCATGTTGATGGAACCAGGCAGCAACTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTATCCAGAATCCATGCCAGTGCCTCTTCACTGATATTACAGTCCGCGACAATGACCTTTGCCCTCTGa                                                                                                                                                                    >  1:158111/1‑289 (MQ=255)

                                                                                                                        gccCATGCTCAATACCAGTCGGTTCAGGCCATGTTGATGGAACCAGGCAGCAACTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGATACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTATCCAGAATCCATGCCAGTGCCTCTTCACTGATATTACAGTCCGCGACAATGACCTTTGCCCTCTGAATAAATTCa                                                                                                                                                           <  2:524880/289‑1 (MQ=255)

                                                                                                                           cATGCTCAATACCAGTCGGTTCAGGCCATGTTGATGGAACCAGGCAGCAACTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTATCCAGAATCCATGCCAGTGCCTCTTCACTGATATTACAGTCCGCGACAATGACCTTTGCCCTCTGAATAAATTCACCg                                                                                                                                                        >  1:269668/1‑289 (MQ=255)

                                                                                                                             tGCTCAATACCAGTCGGTTCAGGCCATGTTGATGGAACCAGGCAGCAACTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTGGACACATTTCCATGCCGAAACCGGGTCGACAAATACGGGTACGTTGGCGGCATTATCCAGAATCCATGCCAGTGCCTCTTCACTGATATTACAGTCCGCGACAATGACCTTTGCCCTCTGAATAAATTCACCGTg                                                                                                                                                      <  1:557505/289‑1 (MQ=255)

                                                                                                                                 cAATACCAGTCGGTTCAGGCCATGTTGATGGAACCAGGCAGCAACTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTATCCAGAATCCATGCCAGTGCCTCTTCACTGATATTACAGTCCGCGACAATGACCTTTGCCCTCTGAATAAATTCACCGTGCtgt                                                                                                                                                  <  2:321765/289‑1 (MQ=255)

                                                                                                                                   aTACCAGTCGGTTCAGGCCATGTTGATGGAACCAGGCAGCAACTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTATCCAGAATCCATGCCAGTGCCTCTTCACTGATATTACAGTCCGCGACAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGc                                                                                                                                                >  1:784364/1‑289 (MQ=255)

                                                                                                                                    tACCAGTCGGTTCAGGCCATGTTGATGGAACCAGGCAGCAACTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGGTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTATCCAGAATCCATGCCAGTGCCTCTTCACTGATATTACAGTCCGCGACAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCg                                                                                                                                               <  1:404355/289‑1 (MQ=255)

                                                                                                                                         gTCGGTTCAGGCCATGTTGATGGAACCAGGCAGCAACTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTATCCAGAATCCATGCCAGTGCCTCTTCACTGATATTACAGTCCGCGACAATGACCTTTTCCCTCTGAATAAATTCACCGTGCTGTGCGAGata                                                                                                                                          >  2:27681/1‑289 (MQ=255)

                                                                                                                                         gTCGGTTCAGGCCATGTTGATGGAACCAGGCAGCAACTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTATCCAGAATCCATGCCAGTGCCTCTTCACTGATATTACAGTCCGCGACAATGACCTTTGCCCTCTGAATAAATTCACCGTg                                                                                                                                                      >  2:887036/1‑277 (MQ=255)

                                                                                                                                         gTCGGTTCAGGCCATGTTGATGGAACCAGGCAGCAACTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTATCCAGAATCCATGCCAGTGCCTCTTCACTGATATTACAGTCCGCGACAATGACCTTTGCCCTCTGAATAAATTCACCGTg                                                                                                                                                      <  1:887036/277‑1 (MQ=255)

                                                                                                                                         gTCGGTTCAGGCCATGTTGATGGAACCAGGCAGCAACTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTATCCAGAATCCATGCCAGTGCCTCTTCACTGATATTACAGTCCGCGACAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCGAGata                                                                                                                                          >  2:142528/1‑289 (MQ=255)

                                                                                                                                         gTCGGTTCAGGCCATGTTGATGGAACCAGGCAGCAACTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTATCCAGAATCCATGCCAGTGCCTCTTCACTGATATTACAGTCCGCGACAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCGAGata                                                                                                                                          >  2:74318/1‑289 (MQ=255)

                                                                                                                                            ggTTCAGGCCATGTTGATGGAACCAGGCAGCAACTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTATCCAGAATCCATGCCAGTGCCTCTTCACTGATATTACAGTCCGCGACAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCGAGATATTc                                                                                                                                       <  1:74318/289‑1 (MQ=255)

                                                                                                                                               tCAGGCCATGTTGATGGAACCAGGCAGCAACTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTATCCAGAATCCATGCCAGTGCCTCTTCACTGATATTACAGTCCGCGACAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCGAGATATTCAg                                                                                                                                     <  2:116914/288‑1 (MQ=255)

                                                                                                                                                    ccATGTTGATGGAACCAGGCAGCAACTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTATCCAGAATCCATGCCAGTGCCTCTTCACTGATATTACAGTCCGCGACAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCGAGATATTCAGCTGTa                                                                                                                                <  2:15045/288‑1 (MQ=255)

                                                                                                                                                        gTTGATGGAACCAGGCAGCAACTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTATCCAGAATCCATGCCAGTGCCTCTTCCCTGATATTACAGTCCGCGACAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCGAGatat                                                                                                                                         >  1:407750/1‑275 (MQ=255)

                                                                                                                                                        gTTGATGGAACCAGGCAGCAACTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTATCCAGAATCCATGCCAGTGCCTCTTCACTGATATTACAGTCCGCGACAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCGAGatat                                                                                                                                         <  2:407750/275‑1 (MQ=255)

                                                                                                                                                        gTTGATGGAACCAGGCAGCAACTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTATCCAGAATCCATGCCAGTGCCTCTTCACTGATATTACAGTCCGCGACAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCGAGATATTCAGCTGTAATAGc                                                                                                                           >  2:514546/1‑289 (MQ=255)

                                                                                                                                                        gTTGATGGAACCAGGCAGCAACTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTATCCAGAATCCATGCCAGTGCCTCTTCACTGATATTACAGTCCGCGACAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCGAGATATTCAGCTGTAATAGc                                                                                                                           >  2:244173/1‑289 (MQ=255)

                                                                                                                                                        gTTGATGGAACCAGGCAGCAACTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTATCCAGAATCCATGCCAGTGCCTCTTCACTGATATTACAGTCCCCGACAATGACCTTTGCCCTCTGAATAAATTCAACGTCCTGTGCGCGATATTCAGCTGTAATAGc                                                                                                                           >  1:658264/1‑289 (MQ=255)

                                                                                                                                                             tGGAACCAGGCAGCAACTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTATCCAGAATCCATGCCAGTGCCTCTTCACTGATATTACAGTCCGCGACAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCGAGATATTCAGCTGTAATAGCGTTGc                                                                                                                      <  2:862962/289‑1 (MQ=255)

                                                                                                                                                              ggAACCAGGCAGCAACTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTATCCAGAATCCATGCCAGTGCCTCTTCACTGATATTACAGTCCGCGACAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCGAGATATTCAGCTGTAATAGCGTTGCt                                                                                                                     <  2:784364/289‑1 (MQ=255)

                                                                                                                                                                   cAGGCAGCAACTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTATCCAGAATCCATGCCAGTGCCTCTTCACTGATATTACAGTCCGCGACAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCGAGATATTCAGCTGTAATAGCGTTGCTAatat                                                                                                                >  1:157824/1‑289 (MQ=255)

                                                                                                                                                                      gcagcaACTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTATCCAGAATCCATGCCAGTGCCTCTTCACTGATATTACAGTCCGCGACAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCGAGATATTCAGCTGTAATAGCGTTGCTAATATTc                                                                                                              <  2:315794/288‑1 (MQ=255)

                                                                                                                                                                           aaCTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTATCCAGAATCCATGCCAGTGCCTCTTCACTGATATTACAGTCCGCGACAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCGAGATATTCAGCTGTAATAGCGTTGCTAATATTCATg                                                                                                           <  2:27996/286‑1 (MQ=255)

                                                                                                                                                                                    cacaTCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGATACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTATCCAGAATCCATGCCAGTGCCTCTTCACTGATATTACAGTCCGCGACAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCGAGATATTCAGCTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCa                                                                                               <  1:769013/289‑1 (MQ=255)

                                                                                                                                                                                       aTCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTATCCAGAATCCATGCCAGTGCCTCTTCACTGATATTACAGTCCGCGACAATGACCTTTGCCCTCTGACTAAATTCACCGTGCTGTGCGAGATATTCAGCTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCAAcc                                                                                            >  2:715819/1‑289 (MQ=255)

                                                                                                                                                                                                 ccTGACAGGGTAAGCCCACTCAGGGTTTGCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGAGCTGATTTAGACCGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACAGGTACGTTGGCGGCATTATCCAGAATCCATGCCAGTGCCTCTTCACTGATATTACTGTCCGCGACAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCGAGATATTCAGCTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCAACCAGCATTTCAc                                                                                  <  1:61710/289‑1 (MQ=255)

                                                                                                                                                                                                 ccTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTATCCAGAATCCATGCCAGTGCCTCTTCACTGATATTACAGTCCGCGACAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCGAGATATTCAGCTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCAACCAGCATTTCAc                                                                                  <  2:158111/289‑1 (MQ=255)

                                                                                                                                                                                                         cgcAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTATCCAGAATCCATGCCAGTGCCTCTTCACTGATATTACAGTCcg                                                                                                                                                                                           >  2:613006/1‑176 (MQ=255)

                                                                                                                                                                                                         cgcAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTATCCAGAATCCATGCCAGTGCCTCTTCACTGATATTACAGTCcg                                                                                                                                                                                           <  1:613006/176‑1 (MQ=255)

                                                                                                                                                                                                         cgcAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTATCCAGAATCCATGCCAGTGCCTCTTCACTGATATTACAGTCCGCGACAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCgag                                                                                                                                             <  2:281971/222‑1 (MQ=255)

                                                                                                                                                                                                         cgcAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTATCCAGAATCCATGCCAGTGCCTCTTCACTGATATTACAGTCCGCGACAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCgag                                                                                                                                             >  1:281971/1‑222 (MQ=255)

                                                                                                                                                                                                         cgcAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTATCCAGAATCCATGCCAGTGCCTCTTCACTGATATTACAGTCCGCGACAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCGAGATATTCAGCTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCAACCAGCATTTCACCGGtatta                                                                          >  2:46626/1‑289 (MQ=255)

                                                                                                                                                                                                          gcAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTATCCAGAATCCATGCCAGTGCCTCTTCACTGATATTACAGTCCGCGACAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCGAGATATTCAGCTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCAACCAGCATTTCACCGGta                                                                             <  1:418117/285‑1 (MQ=255)

                                                                                                                                                                                                          gcAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTATCCAGAATCCATGCCAGTGCCTCTTCACTGATATTACAGTCCGCGACAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCGAGATATTCAGCTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCAACCAGCATTTCACCGGta                                                                             <  1:787848/285‑1 (MQ=255)

                                                                                                                                                                                                           cAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTATCCAGAATCCATGCCAGTGCCTCTTCACTGATATTACAGTCCGCGACAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCGAGATATTCAGTTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCAACCAGCATTTCACCGGTATTATc                                                                        >  2:864961/1‑289 (MQ=255)

                                                                                                                                                                                                             aTCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTATCCAGAATCCATGCCAGTGCCTCTTCACTGATATTACAGTCCGCGACAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCGAGATATTCAGCTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCAACCAGCATTTCACCGGTATTATCGa                                                                      >  1:281318/1‑289 (MQ=255)

                                                                                                                                                                                                              tCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTATCCAGAATCCATGCCAGTGCCTCTTCACTGATATTACAGTCCGCGACAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCGAGATATTCAGCTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCAACCAGCATTTCACCGGTATTATCGAg                                                                     <  2:269668/289‑1 (MQ=255)

                                                                                                                                                                                                                   cTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTATCCAGAATCCATGCCAGTGCCTCTTCACTGATATTACAGTCCGCGACAATGACCTTTGCCCTCTGAATAAATTCACCGTTCTGTTCGAGATATTCAGCTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCAAcc                                                                                            >  1:356425/1‑261 (MQ=255)

                                                                                                                                                                                                                   cTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTATCCAGAATCCATGCCAGTGCCTCTTCACTGATATTACAGTCCGCGACAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCGAGATATTCAGCTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCAAcc                                                                                            <  2:356425/261‑1 (MQ=255)

                                                                                                                                                                                                                   cTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTATCCAGAATCCATGCCAGTGCCTCTTCACTGATATTACAGTCCGCGACAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCGAGATATTCAGCTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCAACCAGCATTTCACCGGTATTATCGAGTAATg                                                                >  1:635220/1‑289 (MQ=255)

                                                                                                                                                                                                                      aGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTATCCAGAATCCATGCCAGTGCCTCTTCACTGATATTACAGTCCGCGACAATGACCTTTGCCCTCTGAATAAATTCACCGTTCTGTGCGAGATATTCCGCTGTAATCGCGTTGCTAATATTCATGTCATTTATAGCAACCAGCATTTCACCGGTATTATCGAGTAATGata                                                             >  1:75917/1‑289 (MQ=255)

                                                                                                                                                                                                                        ggTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTATCCAGAATCCATGCCAGTGCCTCTTCACTGATATTACAGTCCGCGACAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCGAGATATTCAGCTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCAACCAGCATTTCACCGGTATTATCGAGTAATGataaa                                                           >  2:203472/1‑289 (MQ=255)

                                                                                                                                                                                                                          tttCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTATCCAGAATCCATGCCAGTGCCTCTTCACTGATATTACAGTCCGCGACAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCGAGCTATTCAGCTGTAATAGCTTTGCTAATATTCATGTCATTTATAGCAACCAGAATTTCACCGGTATTACCGATTAATGataaata                                                         >  1:783357/1‑289 (MQ=255)

                                                                                                                                                                                                                                  tCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTATCCAGAATCCATGCCAGTGCCTCTTCACTGATATTACAGTCCGCGACAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCGAGATATTCAGCTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCAACCAGCATTTCAACGGTATTATCGAGTAATGATAAATAACTCGACg                                                 >  2:418145/1‑289 (MQ=255)

                                                                                                                                                                                                                                           tttGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTATCCAGAATCCATGCCAGTGCCTCTTCACTGATATTACAGTCCGCGACAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCGAGATATTCAGCTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCAACCAGCATTTCACCGGTATTATCGAGTAATGATAAATAACTCGACGTATTTtct                                         >  2:648945/1‑288 (MQ=255)

                                                                                                                                                                                                                                           tttGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTATCCAGAATCCATGCCAGTGCCTCTTCACTGATATTACAGTCCGCGACAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCGAGATATTCAGCTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCAACCAGCATTTCACCGGTATTATCGAGTAATGATAAATAACTCGACGTATTTtct                                         >  1:537925/1‑288 (MQ=255)

                                                                                                                                                                                                                                                 ttGAGATGGTGGATCTGATTTGACCGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTATCCAGAATCCATGCCAGTGCCTCTTCACTGATATTACAGTCCGCGACAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCGAGATATTCAGCTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCAACCAGCATTTCACCGGTATTATCGAGTAATGATAAATAACTCGACGTATTTTCCCCCGGca                                  <  1:244173/289‑1 (MQ=255)

                                                                                                                                                                                                                                                      agTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTATCCAGAATCCATGCCAGTGCCTCTTCACTGATATTACAGTCCGCGACAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCGAGATATTCAGCTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCAACCAGCATTTCACCGGTATTATCGAGTAATGATAAATAACTCGACGTATTTTCTCCCGGCACAAt                              <  2:450093/288‑1 (MQ=255)

                                                                                                                                                                                                                                                      agTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTATCCAGAATCCATGCCAGTGCCTCTTCACTGATATTACAGTCCGCGACAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCGAGATATTCAGCTGTAATAGCGTTGCTAATATGCATGTCATTTATAGCAACCAGCATTTCACCGGTATTATCGAGTAATGATAAATAACTCGACGTATTTTCTCCCGGCACAAt                              <  2:658264/288‑1 (MQ=255)

                                                                                                                                                                                                                                                                  aTTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTATCCAGAATCCATGCCAGTGCCTCTTCACTGATATTACAGTCCGCGACAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCGAGATATTCAGCTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCAACCAGCATTTCACCGGTATTATCGAGTAATGATAAATAACTCGACGTATTTTCTCCCGGCACAATCAGGCATTTATCg                 >  1:723801/1‑289 (MQ=255)

                                                                                                                                                                                                                                                                            tCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTATCCAGAATCCATGCCAGTGCCTCTTCACTGATATTACAGTCCGCGACAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCGAGATATTCAGCTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCAACCAGCATTTCACCGGTATTATCGAGTAATGATAAATAACTCGACGTATTTTCTCCCGGCACAATCAGGCATTTATCGACATAAACg        >  2:745430/1‑288 (MQ=255)

                                                                                                                                                                                                                                                                            tCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTATCCAGAATCCATGCCAGTGCCTCTTCACTGATATTACAGTCCGCGACAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCGAGATATTCAGCTGTAATAGCGTTGCTAATATTCATGCCATTTATAGCAACCAGCATTTCACCGGTATTATCGAGTAATGATAAATAACTCGACGTATTTTCTCCCGGCACAATCAGGCATTTATCGACATAAACg        >  1:766775/1‑288 (MQ=255)

                                                                                                                                                                                                                                                                             cgcgCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTATCCAGAATCCATGCCAGTGCCTCTTCACTGATATTACAGTCCGCGACAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCGAGATATTCAGCTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCAACCAGCATTTCACCGGTATTATCGAGTAATGATAAATAACTCGACGTATTTTCTCCCGGCACAATCAGGCATTTATCGACATAAACGcc      <  2:635220/289‑1 (MQ=255)

                                                                                                                                                                                                                                                                             cgcgCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTATCCAGAATCCATGCCAGTGCCTCTTCACTGATATTACAGTCCGCGACAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCGAGATATTCAGCTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCAACCAGCATTTCACCGGTATTATCGAGTAATGATAAATAACTCGACGTATTTTCTCCCGGCACAATCAGGCATTTATCGACATAAACGcc      <  1:514546/289‑1 (MQ=255)

                                                                                                                                                                                                                                                                                 cACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTATCCAGAATCCATGCCAGTGCCTCTTCACTGATATTACAGTCCGCGACAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCGAGATATTCAGCTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCAACCAGCATTTCACCGGTATTATCGAGTAATGATAAATAACTCGACGTATTTTCTCCCGGCACAATCAGGCATTTATCGACATAAACGCCAgat  >  2:442227/1‑289 (MQ=255)

ACACGAAGCAAGTCCCGCCATCATGGCATCGCCCGCTCCGGTAACATTAATAACATTGGTTTTGATCGGCGCAGACCAGCCACTTTCACCGCTGATATCGCTGTAATAAACGCCGTCGCCGCCCATGCTCAATACCAGTCGGTTCAGGCCATGTTGATGGAACCAGGCAGCAACTTTTGCCACATCTTCACGCCCTGACAGCGCAATCCCACTCAGGGTTTCCGCTTCAAGGCGGTTTGGCTTGAGAGTGTGGATCTGATTTAGACGGTCGCGCACTTTGACACATTTCCATGCCGAAACCGGATCGACAAATACGGGTACGTTGGCGGCATTATCCAGAATCCATGCCAGTGCCTCTTCACTGATATTACAGTCCGCGACAATGACCTTTGCCCTCTGAATAAATTCACCGTGCTGTGCGAGATATTCAGCTGTAATAGCGTTGCTAATATTCATGTCATTTATAGCAACCAGCATTTCACCGGTATTATCGAGTAATGATAAATAACTCGACGTATTTTCTCCCGGCACAATCAGGCATTTATCGACATAAACGCCAGAT  >  NC_000913/2258500‑2259061

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