BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,774,222	A→G	5.3%	intergenic (+28/‑202)	mtlA → / → mtlD	mannitol‑specific PTS enzyme: IIA, IIB and IIC components/mannitol‑1‑phosphate dehydrogenase, NAD‑dependent

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,774,222	0	A	G	5.3%	3128.9 / 95.7	928	intergenic (+28/‑202)	mtlA/mtlD	mannitol‑specific PTS enzyme: IIA, IIB and IIC components/mannitol‑1‑phosphate dehydrogenase, NAD‑dependent
Reads supporting (aligned to +/- strand): ref base A (489/386); new base G (31/19); total (523/405)
Fisher's exact test for biased strand distribution p-value = 4.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.29e-03

Only 100 of 950 total aligned reads displayed.

TTGAAGCGAAAGATCGCGTACTGAAAACGGGCGTCGTGTTCTGCCAGTACCCGGAAGGCGTGCGCTTCGGTGAAGAAGAAGATGACATTGCCCGTCTGGTGATTGGTATTGCTGCCCGTAACAACGAGCACATTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTG‑CTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGGAGAGGGTTTGGGTGAGGGAAAAGCCTCACCCCAGCCCTCTCGGGTAAAAACATTGATGAAGGTTAATACTATGAAAGCATTACATTTTGGCGCAGGTAATATCGGTCGTGGCTTTATCGGTAAACTGCTGGCAGACGCGGGTA  >  NC_000913/3773950‑3774496

     gCGAAAGATCGCGTACTGAAAACGGGCGTCGTGTTCTGCCAGTACCCGGAAGGCGTGCGCTTCGGTGAAGAAGAAGATGACATTGCCCGTCTGGTGATTGGTATTGCTGCCCGTAACAACGAGCACATTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTG‑CTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACCTGGAGAAGGTggg                                                                                                                                                                                                                                                                              >  1:198236/1‑274 (MQ=255)

                     tGAAAACGGGCGTCGTGTTCTGCCAGTACCCGGAAGGCGTGCGCTTCGGTGAAGAAGAAGATGACATTGCCCGTCTGGTGATTGTTATTGCTGCCCGTAACAACGAGCACATTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTG‑CTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATgcg                                                                                                                                                                                                                                                             <  2:460198/275‑1 (MQ=255)

                        gaaCGCGAGTGGTGTTGTGCCAGTGCCCGGGAGGTGTGCGGTTGGGGGAAGAAGAGGAGGACATTGCCCGTCTGGTGATTGGTATTGGGACCCGTAACAACGAGGACATTCAGGTTAGCACCAGCCTGACCACTGTACTGGATGATGAGACCGGCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTG‑CTGGAATTGCTGGCAGGGCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGGGGGGTTAATTGCCTGATGCGCTa                                                                                                                                                                                                                                                          <  4:187099/274‑1 (MQ=255)

                           cggggggcggtttcttcCAGTACCCGGAAGGCGTGCGGTTCGGTGAAGAAGAAGATGACATTGCCCGTCTGTTGATTGGTATTGCTGCCCGTAAGAACGAGCACATTCAGGTTATCACCAGCCTGACCAATGCGCTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTG‑CTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTAc                                                                                                                                                                                                                                                         <  4:167701/262‑1 (MQ=255)

                                         gtccagagccgggAAGGCGGGTGTTTCGGTGAAGAAGAAGATGACATTGACCGTCTGGTGATTGGTATTGCTGCCCGTAACAACGAGCACATTCAGGTTATGACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTG‑CTGGAACTGCTGGCAGGTCGTAAGGAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACa                                                                                                                                                                                                                                         <  4:482634/264‑1 (MQ=255)

                                               tACCCGGAAGGCGTGCGCTTCGGTGAAGAAGAAGATGACATTGCCCGTCTGGTGATTGGTATTGCTGCCCGTAACAACGAGCACATTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCNGTCATCGAGCGTCTGGCACCCACCACCAGCGTGGATGAAGTC‑CGGGAACTGCTGGCAGGTCGTAAGTAANCCAATCCCACCCCCCCCACATGGCGAAGGAGGGGTTAATTGCCTCNTNCCCTCACCCAATCAGGCCTACaccg                                                                                                                                                                                                                                      >  3:69804/1‑269 (MQ=255)

                                                                        aagagaaatcggACATAGCCCGTGCTGCGATTGGTTTTGCTGCCCGCAACGACGAACAAATTCAGGTTAGACCCAGCATTCCCAAGGCACTGGATAAGGAGTCCGACATCGGGCGTCTGGCACACACCACCAGCGTGGAGGAAGTG‑CTGGAACTGCGGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATTGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACg                                                                                                                                                                                                           <  4:35600/263‑1 (MQ=255)

                                                                                              tCTGGTGATTGGTATTGCTGCCCGTAACAACGAGCACATTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTG‑CTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAATgcgc                                                                                                                                                                                    >  2:624579/1‑275 (MQ=255)

                                                                                              tCTGGTGATTGGTATTGCTGCCCGTAACAACGAGCACATTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTG‑CTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGgcgc                                                                                                                                                                                    >  2:315770/1‑275 (MQ=255)

                                                                                                                 gcCCGTAACAACCGGACCATTCAGGTTCNNNCCAGCCTGACCAATGCACTGGATGATGGGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTG‑CTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGgcgc                                                                                                                                                                  <  3:198236/274‑1 (MQ=255)

                                                                                                                                cacaTTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTG‑CTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATgcg                                                                                                                                                                                                                                                             <  4:692592/168‑1 (MQ=255)

                                                                                                                                cacaTTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTG‑CTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACg                                                                                                                                                  >  4:237395/1‑275 (MQ=255)

                                                                                                                                 acaTTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTG‑CTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTgg                                                                                                                                                                                                                                                                               <  6:1012171/149‑1 (MQ=255)

                                                                                                                                         gTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTG‑CTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGTGGTTAAtt                                                                                                                                                                                                                                                                       <  8:207059/149‑1 (MQ=255)

                                                                                                                                         gTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTG‑CTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAAtt                                                                                                                                                                                                                                                                       >  7:426154/1‑149 (MQ=255)

                                                                                                                                         gTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTG‑CTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAAtt                                                                                                                                                                                                                                                                       <  8:636646/149‑1 (MQ=255)

                                                                                                                                         gTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTG‑CTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGGAGAGGGtt                                                                                                                                         >  4:571020/1‑275 (MQ=255)

                                                                                                                                             tCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTG‑CTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCt                                                                                                                                                                                                                                                                   <  8:561701/149‑1 (MQ=255)

                                                                                                                                             tCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTG‑CTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGtt                                                                                                                                                                                                         >  4:276067/1‑207 (MQ=255)

                                                                                                                                                ccAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTG‑CTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGAt                                                                                                                                                                                                                                                                >  5:672993/1‑149 (MQ=255)

                                                                                                                                                 cAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTG‑CTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGc                                                                                                                                                                                                <  4:373479/212‑1 (MQ=255)

                                                                                                                                                 cAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTG‑CTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGcccc                                                                                                                                                            <  2:506921/248‑1 (MQ=255)

                                                                                                                                                  aGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGTACACACCACCAGCGTGGATGAAGTG‑CTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATgc                                                                                                                                                                                                                                                              >  7:322337/1‑149 (MQ=255)

                                                                                                                                                    ccTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTG‑CTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGGAGAGGGTTTGGGTGAGGGa                                                                                                                              <  2:364817/275‑1 (MQ=255)

                                                                                                                                                       gACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTG‑CTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTAc                                                                                                                                                                                                                                                         >  5:428023/1‑149 (MQ=255)

                                                                                                                                                       gACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTG‑CTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTAt                                                                                                                                                                                                                                                   <  2:53058/155‑1 (MQ=255)

                                                                                                                                                       gACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTG‑CTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCa                                                                                                                                                                                                             <  4:522616/193‑1 (MQ=255)

                                                                                                                                                       gACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTG‑CTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCa                                                                                                                                                                                                             <  2:737781/193‑1 (MQ=255)

                                                                                                                                                             tGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTG‑CTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCTTCTCCATGTGGAGAGGGTggg                                                                                                                                                                                                                                                                              >  5:393632/1‑122 (MQ=255)

                                                                                                                                                             tGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTG‑CTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGGAGAGGGTTTGGGTGAGGGAAAAGCCTCa                                                                                                                     >  2:172038/1‑275 (MQ=255)

                                                                                                                                                                 cTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTG‑CTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAgg                                                                                                                                                                                                                                               <  8:759839/149‑1 (MQ=255)

                                                                                                                                                                  tGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTG‑CTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAAtt                                                                                                                                                                                                                                                                       >  8:5083/1‑124 (MQ=255)

                                                                                                                                                                  tGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTG‑CTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGa                                                                                                                                                                                                                                                                 >  1:115000/1‑130 (MQ=255)

                                                                                                                                                                  tGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTG‑CTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGc                                                                                                                                                                                                                                              <  6:290184/149‑1 (MQ=255)

                                                                                                                                                                  tGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTG‑CTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGc                                                                                                                                                                                                                                              >  8:291513/1‑149 (MQ=255)

                                                                                                                                                                   ggATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTG‑CTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAg                                                                                                                                                                                                                                                                                   <  5:152940/111‑1 (MQ=255)

                                                                                                                                                                    gatgatGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTG‑CTGGAACTGCTGGCAAGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGGAGAGGGTTTGGGTGAGGGAAAAGCCTCACCCCAGc                                                                                                              >  2:131233/1‑275 (MQ=255)

                                                                                                                                                                      tgatgaGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTG‑CTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATgcg                                                                                                                                                                                                                                                             <  3:115000/130‑1 (MQ=255)

                                                                                                                                                                       gatgaGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTG‑CTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACa                                                                                                                                                                                                                                         >  8:949162/1‑149 (MQ=255)

                                                                                                                                                                         tggggccgcaacggAGCGTCTGGCACACACCACCGGCGTGGATGAAGTG‑GTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGGAGAGGGTTTGGGTGAGGGAAAAGCCTCACCCCAGCCctct                                                                                                         <  4:627274/262‑1 (MQ=255)

                                                                                                                                                                           aGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTG‑CTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGAt                                                                                                                                                                                                                                     <  7:589912/149‑1 (MQ=255)

                                                                                                                                                                           aGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTG‑CTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGAt                                                                                                                                                                                                                                     <  8:877400/149‑1 (MQ=255)

                                                                                                                                                                              ccGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTG‑CTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGGAGAGGGTTTGGGTGAGGGAAAAGCCTCACCCCAGCCCTCTCGGGt                                                                                                    >  2:570122/1‑275 (MQ=255)

                                                                                                                                                                              ccGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTG‑CTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTACCCTCACGGAGAGGGTTTGGGTGAGGGAAAAGCCTCACCCCAGACCTTTCggta                                                                                                    >  2:413894/1‑273 (MQ=255)

                                                                                                                                                                               cGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTG‑CTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGcc                                                                                                                                                                                                                                                                    >  6:431720/1‑114 (MQ=255)

                                                                                                                                                                               cGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTG‑CTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCAt                                                                                                                                                                                                                                 <  6:76161/149‑1 (MQ=255)

                                                                                                                                                                                   aTCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTG‑CTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATcaca                                                                                                                                                                                                                             >  8:419787/1‑149 (MQ=255)

                                                                                                                                                                                       aGCGTCTGGCACACACCACCAGCGTGGATGAAGTG‑CTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAAttt                                                                                                                                                                                                                         <  5:475910/149‑1 (MQ=255)

                                                                                                                                                                                              ggCACACACCACCAGCGTGGATGAAGTG‑CTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAAt                                                                                                                                                                                                                  >  7:704629/1‑149 (MQ=255)

                                                                                                                                                                                               gCACACACCACCAGCGTGGATGAAGTG‑CTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAAtt                                                                                                                                                                                                                 <  6:361965/149‑1 (MQ=255)

                                                                                                                                                                                                    cacagccaGTGTGGATG‑AGTGTGTGGAACTGGTGTCGGGTTGTAAGTAATCCATCCCCACCCGCTCCACAGGGGGANGTTGGGGTGAATTGCCTGATGCTCTACGTTTATCAGGCCTACAGGATGCATCTGAATTTGTTGAATTTGCCCGTTCTTGGAGGTCTGATAAGGCGCTTGCGCCGCATCCGGCGCTGCCCCTCTCCTCACGGAGAGGGTGTGGGTGAGGGAAAAGCATCACCACAGCCCTCTCGGGTAACAACATTGATGATGGTTaa                                                                               <  3:166556/269‑1 (MQ=255)

                                                                                                                                                                                                     accaccAGCGTGGATGAAGTG‑CTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTAt                                                                                                                                                                                                                                                   >  8:457158/1‑109 (MQ=255)

                                                                                                                                                                                                        accaGCGTGGATGAAGTG‑CTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTc                                                                                                                                                                                                        >  6:485145/1‑149 (MQ=255)

                                                                                                                                                                                                               tGGATGAAGTG‑CTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGAAAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGGAGAGGGGTTGGGTGAGGGAAAAGCGTCACACCAGCCCTCTCGGGTAAAAACATTGATGAGGGTGAATACTGTGAAAGc                                                                   >  2:220838/1‑275 (MQ=255)

                                                                                                                                                                                                                 gATGAAGTG‑CTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCt                                                                                                                                                                                                                                                           <  7:227829/89‑1 (MQ=255)

                                                                                                                                                                                                                 gATGAAGTG‑CTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTg                                                                                                                                                                                                               >  1:45326/1‑133 (MQ=255)

                                                                                                                                                                                                                 gATGAAGTG‑CTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGcc                                                                                                                                                                                               >  7:529458/1‑149 (MQ=255)

                                                                                                                                                                                                                 gATGAAGTG‑CTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGGAGAGGGTTTGGGTGAGGGaaaa                                                                                                                           >  1:70405/1‑217 (MQ=255)

                                                                                                                                                                                                                     aaGTG‑CTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTctcct                                                                                                                                                      >  2:105346/1‑186 (MQ=255)

                                                                                                                                                                                                                            tGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGgcg                                                                                                                                                                                     >  8:379301/1‑149 (MQ=255)

                                                                                                                                                                                                                             ggAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGgcgc                                                                                                                                                                                    >  5:225830/1‑149 (MQ=255)

                                                                                                                                                                                                                              gAACTGCAGGGAGGTCGTAAGGAAACCAATCCCACCATCTCCAAATGGGGAAGGTGGGGGTAATTGCCTGATGGGTTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGGAGAGGGTTTGGGTGAGGGAAAAGCCTCACCCCAGCCCTCTCGGGTAAAAACATTGATGAAGGTTAATACTATGAAAGCATTACATTTTGgc                                                      <  4:43746/274‑1 (MQ=255)

                                                                                                                                                                                                                                 ctgctgGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTAc                                                                                                                                                                                <  7:761761/149‑1 (MQ=255)

                                                                                                                                                                                                                                  tgctgGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGACTTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTAcg                                                                                                                                                                               >  6:123413/1‑149 (MQ=255)

                                                                                                                                                                                                                                       gCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCa                                                                                                                                                                          >  8:1001133/1‑149 (MQ=255)

                                                                                                                                                                                                                                         aGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATc                                                                                                                                                                        >  8:706813/1‑149 (MQ=255)

                                                                                                                                                                                                                                         aGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATc                                                                                                                                                                        >  6:1076803/1‑149 (MQ=255)

                                                                                                                                                                                                                                         aGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATc                                                                                                                                                                        <  5:85216/149‑1 (MQ=255)

                                                                                                                                                                                                                                              gtaagtaaTCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCg                                                                                                                                                                      <  5:103821/146‑1 (MQ=255)

                                                                                                                                                                                                                                              gtaagtaaTCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGgcg                                                                                                                                                                   >  7:819621/1‑149 (MQ=255)

                                                                                                                                                                                                                                              gtaagtaaTCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGgcg                                                                                                                                                                   >  6:675835/1‑149 (MQ=255)

                                                                                                                                                                                                                                                  gtaaTCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGc                                                                                                                                                               <  5:115052/149‑1 (MQ=255)

                                                                                                                                                                                                                                                    aatccaatccCACCCTCTCCACAAGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGccc                                                                                                                                                             <  6:87960/149‑1 (MQ=255)

                                                                                                                                                                                                                                                     atccaatccCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAg                                                                                                                                                                                                  <  5:604201/111‑1 (MQ=255)

                                                                                                                                                                                                                                                     aattaaccccACCTTCTCCATGTGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGcccc                                                                                                                                                            >  5:299906/8‑149 (MQ=25)

                                                                                                                                                                                                                                                      tccaatccCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCct                                                                                                                                                           >  8:134409/1‑149 (MQ=255)

                                                                                                                                                                                                                                                      tccaatccCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCa                                                                                                                                                    >  2:661697/1‑156 (MQ=255)

                                                                                                                                                                                                                                                       ccaatccCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCctc                                                                                                                                                          <  7:278607/149‑1 (MQ=255)

                                                                                                                                                                                                                                                          accccACCTTCTCCATGTGGAGAGGGTGGGATTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCAc                                                                                                                                                                                                            <  7:361410/94‑1 (MQ=25)

                                                                                                                                                                                                                                                            cccACCTTCTCCATGTGGAGAGGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTctcctc                                                                                                                                                     >  8:352110/1‑149 (MQ=255)

                                                                                                                                                                                                                                                             ccACCTTCTCCATGTGGAGAGGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCa                                                                                                                                                    >  6:225555/1‑149 (MQ=255)

                                                                                                                                                                                                                                                              cACCTTCTCCATGTGGAGAGGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCAc                                                                                                                                                   >  5:545397/1‑149 (MQ=255)

                                                                                                                                                                                                                                                               aCCTTCTCCATGTGGAGAGGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACg                                                                                                                                                  >  7:478399/1‑149 (MQ=255)

                                                                                                                                                                                                                                                               aCCTTCTCCATGTGGAGAGGGTGGGATTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACg                                                                                                                                                  >  6:827588/1‑149 (MQ=25)

                                                                                                                                                                                                                                                               aCCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCtctc                                                                                                                                                        >  7:136932/1‑143 (MQ=255)

                                                                                                                                                                                                                                                                cccTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACgg                                                                                                                                                 >  5:953558/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                  ctctCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGgag                                                                                                                                               >  7:1012630/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                   tctcCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGgaga                                                                                                                                              <  5:752351/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                   tccccACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGgaga                                                                                                                                              <  5:343141/146‑1 (MQ=255)

                                                                                                                                                                                                                                                                    ctcCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGGAGAGGGTTTGGGTGAGGGAAAAGCCTCACCCCAGCCCTCTCGGGTAAAAACATTGATGAAGGTTaa                                                                               <  4:119478/211‑1 (MQ=255)

                                                                                                                                                                                                                                                                     tcCACATGGAGTAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGGAGAgg                                                                                                                                            >  7:935248/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                     tcCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGGAGAGGGTTTGGGTGAGGGAAAAGCCTCACCCCAGCCCTCTCGGGTAAAAACAtt                                                                                           >  1:132983/1‑198 (MQ=255)

                                                                                                                                                                                                                                                                         caTGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGGAGAGGGttt                                                                                                                                        <  6:179250/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                          aTGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGAt                                                                                                                                                                                                                                     >  6:922807/1‑55 (MQ=255)

                                                                                                                                                                                                                                                                          aTGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGAt                                                                                                                                                                                                                                     <  5:922807/55‑1 (MQ=255)

                                                                                                                                                                                                                                                                           tGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGGAGAGGGTTTgg                                                                                                                                      <  6:573097/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                            ggAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGGAGAGGGTTTGGGTGAGGGAAAAGCCTCACCCCAGcc                                                                                                             >  2:25702/1‑173 (MQ=255)

                                                                                                                                                                                                                                                                               gaAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGGAGAGGGTTTGGGTGAGGGAAAAGCCTCACCCCAGCCCTCTCGGGTAAAAACATTGATGAAGGTTAATACTATGAAAGAATTACATTTTGGCGCAGGt                                                >  4:84998/1‑231 (MQ=255)

                                                                                                                                                                                                                                                                                aaGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGGAGAGGGTTTGGGTGAg                                                                                                                                 <  7:88095/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                                 aGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGGAGAGGGTTTg                                                                                                                                       >  1:171473/1‑142 (MQ=255)

TTGAAGCGAAAGATCGCGTACTGAAAACGGGCGTCGTGTTCTGCCAGTACCCGGAAGGCGTGCGCTTCGGTGAAGAAGAAGATGACATTGCCCGTCTGGTGATTGGTATTGCTGCCCGTAACAACGAGCACATTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTG‑CTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGGAGAGGGTTTGGGTGAGGGAAAAGCCTCACCCCAGCCCTCTCGGGTAAAAACATTGATGAAGGTTAATACTATGAAAGCATTACATTTTGGCGCAGGTAATATCGGTCGTGGCTTTATCGGTAAACTGCTGGCAGACGCGGGTA  >  NC_000913/3773950‑3774496

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