BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,800,515	C→A	5.6%	intergenic (+40/‑208)	ygaM → / → nrdH	putative membrane‑anchored DUF883 family ribosome‑binding protein/hydrogen donor for NrdEF electron transport system; glutaredoxin‑like protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,800,515	0	C	A	5.6%	857.4 / 28.2	250	intergenic (+40/‑208)	ygaM/nrdH	putative membrane‑anchored DUF883 family ribosome‑binding protein/hydrogen donor for NrdEF electron transport system; glutaredoxin‑like protein
Reads supporting (aligned to +/- strand): ref base C (112/124); new base A (8/6); total (120/130)
Fisher's exact test for biased strand distribution p-value = 5.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.01e-02

Only 100 of 253 total aligned reads displayed.

CTGGGGCAGCGACGCCAAAGGGGAAGCTGAAGCCGCACGCAGTAAAGCTCAGGCATTGCTGAAAGAAACCAGGGCGCGAATGCATGGTCGTACTCGCGTCCAGCAAGCCGCGCGCGATGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTGTATCAACAGAGAGACAACCCGACGCGTATCATCGCGCCGTATCTTCATTTTAAACGGAAATACGAATCATGCGCATTACTATTTACACTCGTAACGATTGCGTTCAGTGCCACGCCACCAAACGGG  >  NC_000913/2800247‑2800780

cTGGGGCAGCGACGCCAAAGGGGAAGCTGAAGCCGCACGCAGTAAAGCTCAGGCATTGCTGAAAGAAACCAGGGCGCGAATGCATGGTCGTACTCGCGTCCAGCAAGCCGCGCGCGATGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCt                                                                                                                                                                                                                                                                     >  2:599101/1‑275 (MQ=255)

  ggggCAGCGACGCCAAAGGGGAAGCTGAAGCCGCACGCAGTAAAGCTCAGGCATTGCTGAAAGAAACCAGGGCGCGAATGCATGGTCGTACTCGCGTCCAGCAAGCCGCGCGCGATGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAAACGGCGTGCAAAAATGATAATAAATACGAGTCTTTTCACCCGAAGCCTttcttc                                                                                                                                                                                                                                                                   >  4:521994/1‑275 (MQ=255)

                                                  cggccaggcggaaagaaaCCAGGTGGGTAAGTAAGGGTGGGACCCGCGTCCAGCAATCCGCGGGGGATGCCGTTGGCTGAGCGGATTCTTTGGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCa                                                                                                                                                                                                                   <  4:409456/265‑1 (MQ=255)

                                                  aGGCCTTGCTGAAATAACCCAGGGCGCGAATGCATGGTCGTACTCGCGTCCAGCAAGCCGCGCGCGATGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCa                                                                                                                                                                                                                   <  2:440331/275‑1 (MQ=255)

                                                      aTTGCTGAAAGAAACCAGGGCGCGAATGCATGGTCGTACTCGCGTCCAGCAAGCCGCGCGCGATGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGt                                                                                                                                                                                                                                                              <  2:212722/228‑1 (MQ=255)

                                                                    ccAGGGCGCGAATGCATGGTCGTACTCGCGTCCAGCAAGCCGCGCGCGATGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGcc                                                                                                                                                                                                                                                   >  4:128294/1‑225 (MQ=255)

                                                                        ggCGCGAATGCATGGTCGTACTCGCGTCCAGCAAGCCGCGCGCGATGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGt                                                                                                                                                                                                                                               <  2:128296/225‑1 (MQ=255)

                                                                               ggtactGGTCGTACTCGCGTCCAGCCAGGCGCGCTGGATGCCGGTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCATCTGCGGTAGGTATTTGTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTTACCCCGAAGCCTGTCTTCGGTGTTTCTTTTTGCCTGGTTAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACt                                                                                                                                                                                      <  4:256056/270‑1 (MQ=255)

                                                                                 tcgtggtggTACTCGGTTCCAGCAAGCCGCGCGGGAGTCCGTTGGTTGCTCGGATTTTTTTGTTCTGGAAAGATCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGAATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTAt                                                                                                                                                                                    <  4:221519/272‑1 (MQ=255)

                                                                                       gctgaCTCGCGTCCGGCAAGCCGCGCGGGAGGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGAGCCTGGTGTAGCGTGGGTACACCAGCTTCGTTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCt                                                                                                                                                                               <  4:91577/270‑1 (MQ=255)

                                                                                         gTACTCGCGTCCAGCAAGCCGCGCGCGATGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGAc                                                                                                                                                                                          >  2:369778/1‑261 (MQ=255)

                                                                                                 gTCCAGCAAGCCGCGCGCGATGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGAc                                                                                                                                                                                          >  2:256059/1‑253 (MQ=255)

                                                                                                                        cGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGc                                                                                                                                                                                                                                                                           <  5:965119/149‑1 (MQ=255)

                                                                                                                          ttGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCt                                                                                                                                                                                                                                                                         <  7:1619160/149‑1 (MQ=255)

                                                                                                                            ggCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGt                                                                                                                                                                                                                                                                       <  5:1313805/149‑1 (MQ=255)

                                                                                                                               tGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCtt                                                                                                                                                                                                                                                                    >  8:1057427/1‑149 (MQ=255)

                                                                                                                                 cgcgGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCg                                                                                                                                                                                                                                                                  <  6:841278/149‑1 (MQ=255)

                                                                                                                                    ggATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGACAGGCTTCgggg                                                                                                                                                                                                                                                               <  8:827204/149‑1 (MQ=255)

                                                                                                                                    ggATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGACAGGCTTCgggg                                                                                                                                                                                                                                                               <  7:1497355/149‑1 (MQ=255)

                                                                                                                                       ttCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTAt                                                                                                                                                                                    <  1:116295/221‑1 (MQ=255)

                                                                                                                                       ttCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACAACTAt                                                                                                                                                                                    >  3:116295/1‑221 (MQ=255)

                                                                                                                                          ttttGTTGGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGAAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGAGGTTTCtt                                                                                                                                                                                                                                                         >  5:1514080/1‑149 (MQ=255)

                                                                                                                                          ttttGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCtt                                                                                                                                                                                                                                                         <  6:26085/149‑1 (MQ=255)

                                                                                                                                             tGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCttttt                                                                                                                                                                                                                                                      <  7:1319859/149‑1 (MQ=255)

                                                                                                                                               ttCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGc                                                                                                                                                                                                                                                    >  6:1586628/1‑149 (MQ=255)

                                                                                                                                               ttCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGc                                                                                                                                                                                                                                                    >  7:1576914/1‑149 (MQ=255)

                                                                                                                                                  gTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTg                                                                                                                                                                                                                                                 >  6:1177283/1‑149 (MQ=255)

                                                                                                                                                    gAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGt                                                                                                                                                                                                                                               >  8:1603342/1‑149 (MQ=255)

                                                                                                                                                    gAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGt                                                                                                                                                                                                                                               <  8:681656/149‑1 (MQ=255)

                                                                                                                                                    gAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGt                                                                                                                                                                                                                                               <  8:752160/149‑1 (MQ=255)

                                                                                                                                                     aaaGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTg                                                                                                                                                                                                                                              >  8:166649/1‑149 (MQ=255)

                                                                                                                                                     aaaGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTg                                                                                                                                                                                                                                              <  5:41242/149‑1 (MQ=255)

                                                                                                                                                        aaCCCAGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTGTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGt                                                                                                                                                                                                                     <  4:64509/170‑1 (MQ=255)

                                                                                                                                                            cTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATcaca                                                                                                                                                                                                                                       >  6:1312567/1‑149 (MQ=255)

                                                                                                                                                             tGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTGTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACaa                                                                                                                                                                                                                                      <  7:1125253/149‑1 (MQ=255)

                                                                                                                                                                    gCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCTCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAAt                                                                                                                                                            >  1:179302/1‑216 (MQ=255)

                                                                                                                                                                     cGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATcccc                                                                                                                                                                                                                              <  8:180215/149‑1 (MQ=255)

                                                                                                                                                                     cGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATcccc                                                                                                                                                                                                                              <  8:1556822/149‑1 (MQ=255)

                                                                                                                                                                       tGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCa                                                                                                                                                                                                             >  2:760338/1‑164 (MQ=255)

                                                                                                                                                                         ggTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTAc                                                                                                                                                                                                                          >  8:333164/1‑149 (MQ=255)

                                                                                                                                                                         ggTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTAc                                                                                                                                                                                                                          <  6:499560/149‑1 (MQ=255)

                                                                                                                                                                         ggTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTAc                                                                                                                                                                                                                          >  6:1009021/1‑149 (MQ=255)

                                                                                                                                                                         ggTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAg                                                                                                                                        >  4:500699/1‑231 (MQ=255)

                                                                                                                                                                         ggTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTAGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCACCTACCCCGTCACGCTCATATCCAGGGTAATTTCgccaa                                                                                                                                                                                        >  4:158728/1‑179 (MQ=255)

                                                                                                                                                                            aCAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTAcccc                                                                                                                                                                                                                       <  5:1009021/149‑1 (MQ=255)

                                                                                                                                                                            aCAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTAcccc                                                                                                                                                                                                                       >  8:1329585/1‑149 (MQ=255)

                                                                                                                                                                             cagcagCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTAt                                                                                                                                                                                    <  2:158730/183‑1 (MQ=255)

                                                                                                                                                                             cagcagCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTAtc                                                                                                                                    <  2:500703/231‑1 (MQ=255)

                                                                                                                                                                              agcagcTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGt                                                                                                                                                                                                                     >  5:772727/1‑149 (MQ=255)

                                                                                                                                                                                cagcTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTg                                                                                                                                                                                                                                                 <  7:23184/119‑1 (MQ=255)

                                                                                                                                                                                 agcTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCAc                                                                                                                                                                                                                  <  6:1612916/149‑1 (MQ=255)

                                                                                                                                                                                  gcTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATc                                                                                                                                                                                                         <  2:497594/157‑1 (MQ=255)

                                                                                                                                                                                     gCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGNGCAAATNGTANGGCAAAANTGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGAGTAATTTCGACCACTAt                                                                                                                                                                                    <  3:39395/175‑1 (MQ=255)

                                                                                                                                                                                     gCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTAt                                                                                                                                                                                    >  1:39395/1‑175 (MQ=255)

                                                                                                                                                                                      cGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGt                                                                                                                                                                                                                                               <  8:599683/115‑1 (MQ=255)

                                                                                                                                                                                      cGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGt                                                                                                                                                                                                                                               >  7:599683/1‑115 (MQ=255)

                                                                                                                                                                                         taggtaTTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCatat                                                                                                                                                                                                          >  8:1472819/1‑149 (MQ=255)

                                                                                                                                                                                          aggtaTTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATc                                                                                                                                                                                                         <  7:477813/149‑1 (MQ=255)

                                                                                                                                                                                            gtaTTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCa                                                                                                                                                                                                       >  8:406823/1‑149 (MQ=255)

                                                                                                                                                                                            gtaTTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCa                                                                                                                                                                                                       <  7:333164/149‑1 (MQ=255)

                                                                                                                                                                                             taTTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAg                                                                                                                                                                                                      >  5:1167512/1‑149 (MQ=255)

                                                                                                                                                                                             taTTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAg                                                                                                                                                                                                      >  5:365457/1‑149 (MQ=255)

                                                                                                                                                                                              aTTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAgg                                                                                                                                                                                                     <  8:738576/149‑1 (MQ=255)

                                                                                                                                                                                                  ttATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTGTATCAACAGAGAGACAACCCGACGCGTATCATCGCGCCGTATCTTCATTTTAAACGGaaa                                                                    <  2:280703/274‑1 (MQ=255)

                                                                                                                                                                                                       ggCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCg                                                                                                                                                                                            <  6:1570753/149‑1 (MQ=255)

                                                                                                                                                                                                             cTGTTAAGC‑TGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTGTATCAACAGAGAGACAACCCGACGCGTATCATCGCGCCGTATCTTCATTTTAAACGGAAATACGAATCATgcg                                                       <  2:475620/275‑1 (MQ=255)

                                                                                                                                                                                                              tGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGAc                                                                                                                                                                                          <  7:390172/144‑1 (MQ=255)

                                                                                                                                                                                                               gTTAGGCATGCGCAAATCGTAGTGGAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTAt                                                                                                                                                                                    >  5:1459986/1‑149 (MQ=255)

                                                                                                                                                                                                                    gCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCt                                                                                                                                                                               >  7:1205580/1‑149 (MQ=255)

                                                                                                                                                                                                                    gCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCt                                                                                                                                                                               >  8:381444/1‑149 (MQ=255)

                                                                                                                                                                                                                     cATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTGGCta                                                                                                                                                                              >  5:162000/1‑149 (MQ=255)

                                                                                                                                                                                                                        gcgcAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTCTCTATTGTGTGGTTGAATCTTTTTTCAACTACATATAGTATCTCTGTATCAACTGAGAGACAACCCGACGCGTATCATCGCACCGTACCTTAATTTTAAACGGAAATACAAATCactgttccttatgc                                              >  4:478534/1‑261 (MQ=255)

                                                                                                                                                                                                                          gcAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGt                                                                                                                                                                                                                     <  3:205534/105‑1 (MQ=255)

                                                                                                                                                                                                                          gcAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAAtt                                                                                                                                                                                               >  8:1000221/1‑127 (MQ=255)

                                                                                                                                                                                                                               tCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATtgtgt                                                                                                                                                                    <  6:200923/149‑1 (MQ=255)

                                                                                                                                                                                                                                cGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATtgtgtg                                                                                                                                                                   <  7:166649/149‑1 (MQ=255)

                                                                                                                                                                                                                                   aGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTAtt                                                                                                                                                                                   >  2:40605/1‑130 (MQ=255)

                                                                                                                                                                                                                                   aGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTCTGTGGtt                                                                                                                                                                <  6:649773/149‑1 (MQ=255)

                                                                                                                                                                                                                                        aaaaaTGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACTCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGtt                                                                                                                                                                >  5:299361/1‑144 (MQ=255)

                                                                                                                                                                                                                                         aaaaTGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCt                                                                                                                                                          >  6:1212845/1‑149 (MQ=255)

                                                                                                                                                                                                                                           aaTGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCttt                                                                                                                                                        <  5:1642502/149‑1 (MQ=255)

                                                                                                                                                                                                                                           aaTGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCtct                                                                                                                                                        <  7:818898/149‑3 (MQ=255)

                                                                                                                                                                                                                                                 aataaataCGCGCCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCaa                                                                                                                                                  <  7:115740/149‑1 (MQ=255)

                                                                                                                                                                                                                                                   taaataCGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACt                                                                                                                                                >  8:104282/1‑149 (MQ=255)

                                                                                                                                                                                                                                                    aaataCGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATca                                                                                                                                                                                                                                         <  8:1143787/59‑1 (MQ=255)

                                                                                                                                                                                                                                                    aaataCGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGTGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTa                                                                                                                                               <  6:1514080/149‑1 (MQ=255)

                                                                                                                                                                                                                                                       taCGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACAt                                                                                                                                            <  5:21939/149‑1 (MQ=255)

                                                                                                                                                                                                                                                        aCGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATc                                                                                                                                           >  7:669500/1‑149 (MQ=255)

                                                                                                                                                                                                                                                           cgTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAg                                                                                                                                        >  6:982987/1‑149 (MQ=255)

                                                                                                                                                                                                                                                               tttGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACg                                                                                                                                                                                                                 >  1:244397/1‑72 (MQ=255)

                                                                                                                                                                                                                                                               tttGACCCCGAAGCCTGTCTTCGGAGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCttt                                                                                                                                                        <  5:1256568/129‑1 (MQ=255)

                                                                                                                                                                                                                                                                 tGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTatatat                                                                                                                                                                         >  4:676573/1‑110 (MQ=255)

                                                                                                                                                                                                                                                                  gACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATctct                                                                                                                                 <  5:856451/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                    ccccGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTGt                                                                                                                               <  6:390296/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                    ccccGAAGACAGGCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTGt                                                                                                                               >  6:914095/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                     cccGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCtttt                                                                                                                                                       >  8:483832/1‑124 (MQ=255)

                                                                                                                                                                                                                                                                     cccGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCCACTACATCTAGTATCTCTGTa                                                                                                                              <  7:717187/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                           gCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTGTATCAACa                                                                                                                        >  6:761942/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                            ccTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTGTATCAACag                                                                                                                       >  8:1377643/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                            ccTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTGTCAACTACATTTAGTATCTCTGTATCAACAGAGAGACAACCCGACGCGTATCAGCGCGCCGTATCTTCATTTTAAACGGa                                                                      <  2:303848/198‑1 (MQ=255)

CTGGGGCAGCGACGCCAAAGGGGAAGCTGAAGCCGCACGCAGTAAAGCTCAGGCATTGCTGAAAGAAACCAGGGCGCGAATGCATGGTCGTACTCGCGTCCAGCAAGCCGCGCGCGATGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTGTATCAACAGAGAGACAACCCGACGCGTATCATCGCGCCGTATCTTCATTTTAAACGGAAATACGAATCATGCGCATTACTATTTACACTCGTAACGATTGCGTTCAGTGCCACGCCACCAAACGGG  >  NC_000913/2800247‑2800780

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