BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,642,361	T→A	6.1%	intergenic (+30/+19)	dtpB → / ← rsmJ	dipeptide and tripeptide permease B/16S rRNA m(2)G1516 methyltransferase, SAM‑dependent

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,642,361	0	T	A	6.1%	752.5 / 27.9	213	intergenic (+30/+19)	dtpB/rsmJ	dipeptide and tripeptide permease B/16S rRNA m(2)G1516 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand): ref base T (97/103); new base A (5/8); total (102/111)
Fisher's exact test for biased strand distribution p-value = 5.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.86e-01

Only 100 of 213 total aligned reads displayed.

AGCATTTGATGGGCTTTATTCTCGGGATGTGGTTCCTGACGCAGGCTGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCGGACAACATTACCGATCCGCTTGAGACGTTGCCCGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGTTTGGTCGCCAGCAAGCGTGCAGGCTCCAGTAATCCATCGGCATCAAGATCCGGTCCCACCAGCGACTGAAAGACACGCATCTCTTTTTTCACTAGCGCGCTTTTCTGCTTATGCGGGAACATCGGGTCGAGATAAACC  >  NC_000913/3642090‑3642632

ggcaTTTGTGTGGCGTTTTTTTCGGGAGGTGGTGCCGTCAGAAGGCTGCCGGGTTTTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCGGACAACAGTACCGATCCGCTTGAGACGTTGCCAGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCg                                                                                                                                                                                                                                                                               <  4:157469/273‑1 (MQ=255)

agcgTTTGATGGGCTTTATTCTCGGGATGTGGTGCCGGACGCAGGCTGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCGGACAACATTACCGATCCGCTTGAGACGTTGCCCGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCg                                                                                                                                                                                                                                                                               <  4:100653/274‑1 (MQ=255)

                           tgtgGTTCCTGACGCAGGCTGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCGGACAACATTACCGATCCGCTTGAGACGTTGCCCGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCg                                                                                                                                                                                                                                                    <  2:254609/274‑1 (MQ=255)

                           tgtgGTTACTGACTCAGGCTGCCGCGGTTTTGCTGGGCGGCTATGTGGCAACATTCACCGCGGTGCCGGACAACATTACCGATCCGCCTGAGACGTTGCCCGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCAGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCg                                                                                                                                                                                                                                                    <  1:62501/274‑1 (MQ=255)

                                    tGACGCAGGCTGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCGGACAACATTACCGATCCGCTTGAGACGTTGCCCGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGACGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACGCTGTTTTAGTTTTACACTGGCGTACCTGCATa                                                                                                                                                                                                                                          <  1:97696/275‑1 (MQ=255)

                                            gCTGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCGGACAACATTACCGATCCGCTTGAGACGTTGCCCGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCCGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCCTACCTCAATAGATATCaa                                                                                                                                                                                                                                  >  4:284894/1‑275 (MQ=255)

                                                ccGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCGGACAACATTACCGATCCGCTTGAGACGTTGCCCGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTAc                                                                                                                                                                                                                                                 <  2:41835/256‑1 (MQ=255)

                                                     ttCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCGGACAACATTACCGATCCGCTTGAGACGTTGCCCGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGcc                                                                                                                                                                                                                         <  4:352300/275‑1 (MQ=255)

                                                                                         gTGCCGGACAACATTACCGATCCGCTTGAGACGTTGCCCGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTCTTACACTGCCTTACCTGCATAGATATCAAACCGATGCCCTTTAGTGACTCCCGCGTTTGGCTGTGCGACAtttg                                                                                                                                                                                     >  4:412726/1‑273 (MQ=255)

                                                                                             cGGACAACATTACCGATCCGCTTGAGACGTTGCCCGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCt                                                                                                                                                                                                                       <  4:115649/237‑1 (MQ=255)

                                                                                                            aTCCGCTTGAGACGTTGCCCGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAACCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTAAACTGGCGTACCTGCATAGATATCAAAGCGATTCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGt                                                                                                                                                                   >  2:590056/1‑274 (MQ=255)

                                                                                                              ccGCTTGAGACGTTGCCCGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACt                                                                                                                                                                                                                                                        >  4:554599/1‑187 (MQ=255)

                                                                                                                           tGCCCGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAg                                                                                                                                                                                                                                >  2:697073/1‑198 (MQ=255)

                                                                                                                            gCCCGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATc                                                                                                                                                                                                                                                                                <  7:277591/149‑1 (MQ=255)

                                                                                                                              ccGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGc                                                                                                                                                                                                                                                                              >  7:7363/1‑149 (MQ=255)

                                                                                                                               cGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCt                                                                                                                                                                                                                                                                             <  6:872577/149‑1 (MQ=255)

                                                                                                                                gTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGAAGCGTTATCGCTa                                                                                                                                                                                                                                                                            >  6:336697/1‑149 (MQ=255)

                                                                                                                                 tCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTAc                                                                                                                                                                                                                                                                           >  8:841783/1‑149 (MQ=255)

                                                                                                                                 tCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTAc                                                                                                                                                                                                                                                                           <  6:47365/149‑1 (MQ=255)

                                                                                                                                    ataCCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATCGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGATAACGCTACCCt                                                                                                                                                                                                                                                                        <  7:623663/149‑1 (MQ=255)

                                                                                                                                        cAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTgttt                                                                                                                                                                                                                                                                    >  8:231128/1‑149 (MQ=255)

                                                                                                                                         aaCGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTgtttt                                                                                                                                                                                                                                                                   >  5:560192/1‑149 (MQ=255)

                                                                                                                                         aaCGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTgtttt                                                                                                                                                                                                                                                                   >  6:859291/1‑149 (MQ=255)

                                                                                                                                         aaCGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGATAACGCTAccc                                                                                                                                                                                                                                                                         >  7:607471/1‑143 (MQ=255)

                                                                                                                                           cGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTgttttag                                                                                                                                                                                                                                                                 <  5:391551/149‑1 (MQ=255)

                                                                                                                                           cGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCaa                                                                                                                                                                                                                                  >  4:29170/1‑180 (MQ=255)

                                                                                                                                           cGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCaa                                                                                                                                                                                                                                  >  4:670256/1‑180 (MQ=255)

                                                                                                                                            gtgtTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTgttttagt                                                                                                                                                                                                                                                                <  8:490729/149‑1 (MQ=255)

                                                                                                                                             tgtTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTgttttagtt                                                                                                                                                                                                                                                               <  5:293644/149‑1 (MQ=255)

                                                                                                                                              gtTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTgttttagttt                                                                                                                                                                                                                                                              <  8:95968/149‑1 (MQ=255)

                                                                                                                                               tttGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCg                                                                                                                                                                                                                              <  2:29170/180‑1 (MQ=255)

                                                                                                                                                ttGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTgttttagtttta                                                                                                                                                                                                                                                            <  6:777439/149‑1 (MQ=255)

                                                                                                                                                  ggTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTaca                                                                                                                                                                                                                                                          <  6:1017487/149‑1 (MQ=255)

                                                                                                                                                   gTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGGAGCGTTATCGCTACCCTGTTTTAGTTTTacac                                                                                                                                                                                                                                                         >  6:674582/1‑149 (MQ=255)

                                                                                                                                                    tAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTCGCGTTATCGCTACCCTGTTTTAGTTTTACACt                                                                                                                                                                                                                                                        >  6:1003190/1‑149 (MQ=255)

                                                                                                                                                       gATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGc                                                                                                                                                                                                                                                     <  5:3056/149‑1 (MQ=255)

                                                                                                                                                       gATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGTTTGGTCGCCAGCAAGCGTGCa                                                                                                                       <  2:195389/275‑1 (MQ=255)

                                                                                                                                                          tggtctggtcACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGTTTGGTCGCCAGCAAGCGTGCAgg                                                                                                                     <  2:607453/274‑1 (MQ=255)

                                                                                                                                                           ggtctggtcACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTAc                                                                                                                                                                                                                                                 <  7:582814/149‑1 (MQ=255)

                                                                                                                                                            gtctggtcACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTAcc                                                                                                                                                                                                                                                <  6:641682/149‑1 (MQ=255)

                                                                                                                                                             tctggtcACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCaa                                                                                                                                                                                                                                  >  4:224878/1‑162 (MQ=255)

                                                                                                                                                                agtcACTCTGGGCGTTTCAGTAGTGATTCTTTTGATGGTGCCGTGGCTGAAACGCATGATTGCAACGCCGGAAAGCCATTAATTATTCTTTCAGAAAGCAGGGTAGCGTTATCGCTCCCCTGTTTTAGTTTTACACTGTCGTAcc                                                                                                                                                                                                                                                <  4:747905/144‑1 (MQ=255)

                                                                                                                                                                 gtcACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCAt                                                                                                                                                                                                                                           <  5:959808/149‑1 (MQ=255)

                                                                                                                                                                     cGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGat                                                                                                                                                                                                                                       >  5:867332/1‑149 (MQ=255)

                                                                                                                                                                       cTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGatat                                                                                                                                                                                                                                     <  7:509022/149‑1 (MQ=255)

                                                                                                                                                                           gCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGATAACGCTACCCt                                                                                                                                                                                                                                                                        <  8:714848/110‑1 (MQ=255)

                                                                                                                                                                             gTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGc                                                                                                                                                                                                                               <  8:1055728/149‑1 (MQ=255)

                                                                                                                                                                             gTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGc                                                                                                                                                                                                                               >  6:711788/1‑149 (MQ=255)

                                                                                                                                                                             gTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCATACCTGCATAGATATCAAAGc                                                                                                                                                                                                                               >  8:781914/1‑149 (MQ=255)

                                                                                                                                                                                gCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGAt                                                                                                                                                                                                                            <  7:753199/149‑1 (MQ=255)

                                                                                                                                                                                  agtagtAATGCTTTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAACCCATGCCCTTTAGTTAATACCGCGTTTGGCTTGACGCCATACGCCAGCGGTGCCNCGTAATCCGGGCGCTTTACCGCCACCCGCTGGGGCGCCCCCACGCGTAAAGCCACCATTAATCCCTCAGCATaaaa                                                                                             >  3:157492/1‑273 (MQ=255)

                                                                                                                                                                                    tagtGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTCTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTc                                                                                                                                                                  >  2:118018/1‑203 (MQ=255)

                                                                                                                                                                                     agtGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTATCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTgg                                                                                                                                                                           >  2:632095/1‑193 (MQ=255)

                                                                                                                                                                                        gATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTa                                                                                                                                                                                                                    <  7:111289/149‑1 (MQ=255)

                                                                                                                                                                                        gATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTa                                                                                                                                                                                                                    <  8:595503/149‑1 (MQ=255)

                                                                                                                                                                                         aTGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTaca                                                                                                                                                                                                                                                          >  4:618609/1‑110 (MQ=255)

                                                                                                                                                                                         aTGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCaa                                                                                                                                                                                                                                  >  4:51718/1‑134 (MQ=255)

                                                                                                                                                                                         aTGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGTTTGGTCGCCAGCAAGCGTGCAGGCTCCAGTAATCCATCGGCATCAAGATCCGGTc                                                                                     >  2:375085/1‑275 (MQ=255)

                                                                                                                                                                                             tGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCg                                                                                                                                                                                                                              <  2:51718/134‑1 (MQ=255)

                                                                                                                                                                                             tGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGAc                                                                                                                                                                                                               >  6:837408/1‑149 (MQ=255)

                                                                                                                                                                                             tGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGAc                                                                                                                                                                                                               >  5:803289/1‑149 (MQ=255)

                                                                                                                                                                                                   tGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACcgc                                                                                                                                                                                                         >  6:668101/1‑149 (MQ=255)

                                                                                                                                                                                                   tGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACcgc                                                                                                                                                                                                         <  6:560192/149‑1 (MQ=255)

                                                                                                                                                                                                    ggTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTCTTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCCCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGTTTGGCCGCCTGCAAGCGTGCAGGCTCCAGTAATCCATCGGCATCAAGATCTGGTCCCACCAGCGa                                                                           >  2:353327/1‑274 (MQ=255)

                                                                                                                                                                                                      tGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGtt                                                                                                                                                                                                      <  5:1060476/149‑1 (MQ=255)

                                                                                                                                                                                                      tGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGtt                                                                                                                                                                                                      >  7:355362/1‑149 (MQ=255)

                                                                                                                                                                                                       gCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCCCGCGTTTGGTCGCCAGCAAGCGTGCAGGGTCCAGTACTCCATCGGCATCAAGATCCGGTCACACCAGCCACt                                                                         >  1:162387/1‑273 (MQ=255)

                                                                                                                                                                                                         cGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTac                                                                                                                                                      >  4:246140/1‑194 (MQ=255)

                                                                                                                                                                                                          gTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCt                                                                                                                                                                                                                                               <  8:530191/104‑1 (MQ=255)

                                                                                                                                                                                                          gTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTtggc                                                                                                                                                                                                  <  7:161050/149‑1 (MQ=255)

                                                                                                                                                                                                             gCTGAAACGCATGATTGCGACGCAGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTCATCGCCACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGTCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTCGTCCGGGCGCTTTACCACCACGCGTTTGGTCGCCAGCAAGCGTGCAGGCTCCAGTAACCCATCGGCATCAAGATCCGGTCCCACCAGAGACTGAAAGaca                                                                 >  2:167428/1‑275 (MQ=255)

                                                                                                                                                                                                              cTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTtggcgtgg                                                                                                                                                                                              <  8:908366/149‑1 (MQ=255)

                                                                                                                                                                                                               tGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGATAACGCTAcc                                                                                                                                                                                                                                                                          >  8:271743/1‑72 (MQ=39)

                                                                                                                                                                                                                gAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTtggcgtggcg                                                                                                                                                                                            >  6:737092/1‑149 (MQ=255)

                                                                                                                                                                                                                     gCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACAtt                                                                                                                                                                                       >  8:134060/1‑149 (MQ=255)

                                                                                                                                                                                                                        tGATTGCGACGCCGGTAAGCCATTAATGATTCTTGCAGAAAGCAGGGTAGCGTAATCGc                                                                                                                                                                                                                                                                              >  1:147596/1‑59 (MQ=255)

                                                                                                                                                                                                                        tGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGc                                                                                                                                                                                    >  6:580521/1‑149 (MQ=255)

                                                                                                                                                                                                                          aTTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCa                                                                                                                                                                                  <  8:406715/149‑1 (MQ=255)

                                                                                                                                                                                                                               gACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGt                                                                                                                                                                             >  8:1000058/1‑149 (MQ=255)

                                                                                                                                                                                                                               gACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGt                                                                                                                                                                             >  7:203830/1‑149 (MQ=255)

                                                                                                                                                                                                                                           ccATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTc                                                                                                                                                                                      >  7:875742/1‑128 (MQ=255)

                                                                                                                                                                                                                                             attaattaTTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGTTTGGTCGCCAGCAAg                                                                                                                             >  4:85614/1‑183 (MQ=255)

                                                                                                                                                                                                                                              ttaattaTTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCGTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCggg                                                                                                                                                              >  8:930058/1‑149 (MQ=255)

                                                                                                                                                                                                                                              ttaattaTTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCggg                                                                                                                                                              >  7:106555/1‑149 (MQ=255)

                                                                                                                                                                                                                                               taattaTTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGgc                                                                                                                                                             <  7:841783/149‑1 (MQ=255)

                                                                                                                                                                                                                                                aattaTTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCAt                                                                                                                                                                                                                                           <  2:23178/70‑1 (MQ=255)

                                                                                                                                                                                                                                                 attattCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGgcgc                                                                                                                                                           <  7:781914/149‑1 (MQ=255)

                                                                                                                                                                                                                                                 attattCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGTTTGGTCGCCAGCAAGCGTg                                                                                                                         <  2:85614/183‑1 (MQ=255)

                                                                                                                                                                                                                                                   tattCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCtt                                                                                                                                                         >  7:704358/1‑149 (MQ=255)

                                                                                                                                                                                                                                                      tCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGTTTGGTCGCCAGCAAGCGTGCAGGCTCCAGTAATCCATCGGCATCAAGATCCGGTCCCACCAGCGACTGAAAGACACGCATCTCTTTTTTCACTAGCGCGCTTTTCTGCTTATGCgg                        >  1:121390/1‑275 (MQ=255)

                                                                                                                                                                                                                                                        ttGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCttt                                                                                                                                                        <  5:72511/145‑1 (MQ=255)

                                                                                                                                                                                                                                                        tggCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTacca                                                                                                                                                    >  7:925987/3‑149 (MQ=255)

                                                                                                                                                                                                                                                           actaaaacAGGGTAGCGATAACGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTAccacca                                                                                                                                                 >  6:1084796/8‑149 (MQ=255)

                                                                                                                                                                                                                                                            aGAAAGCAGGGTAGCGTTATCGCTACCCTG‑TTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCAcg                                                                                                                                               >  7:342042/1‑149 (MQ=255)

                                                                                                                                                                                                                                                             gAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCAcg                                                                                                                                               >  8:874152/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                aGCAGGGTAGCGTTATCGCTACCCTGTTTTACTTTTACACTGGCTTACCGGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGTTTGGTCGCCAGCAAGCGTGCAGGCTCCAGTAATCCATCGGCATCAAGATCCGGTCCCACCAGCGACTGAAAGACACGCATCTCTTTTTTCACTAGCGCGCTTTTCTGCTTATGCGGGAACATCggg              <  2:253148/275‑1 (MQ=255)

                                                                                                                                                                                                                                                                     ggTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGTTTGGt                                                                                                                                       <  5:43099/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                      ggaGCGTTATCGGTAACCGTGTTAAGGTTGACACGGGCGTACCTGCATAGATTTCAAAGCGATGCCCGTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGTTTGGTCGCCAGCAAGCGTGCAGGCTCCAGTAATCCATCGGCATCAAGATCCGGTCCCACCAGCGACTGAAAGACACGCATCTCTTTTTTCACTAGCGCGCTTTTCTGCTTATGCGGGAACATCGGGTCgaga        <  4:337048/273‑1 (MQ=255)

                                                                                                                                                                                                                                                                      gTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGTTTGGTc                                                                                                                                      <  5:668101/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                            ttATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGTTTGGTCGCCAGc                                                                                                                                >  6:931329/1‑149 (MQ=255)

AGCATTTGATGGGCTTTATTCTCGGGATGTGGTTCCTGACGCAGGCTGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCGGACAACATTACCGATCCGCTTGAGACGTTGCCCGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGTTTGGTCGCCAGCAAGCGTGCAGGCTCCAGTAATCCATCGGCATCAAGATCCGGTCCCACCAGCGACTGAAAGACACGCATCTCTTTTTTCACTAGCGCGCTTTTCTGCTTATGCGGGAACATCGGGTCGAGATAAACC  >  NC_000913/3642090‑3642632

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