BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,642,358	T→A	6.1%	intergenic (+27/+22)	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,358	0	T	A	6.1%	742.1 / 28.4	213	intergenic (+27/+22)	dtpB/rsmJ	dipeptide and tripeptide permease B/16S rRNA m(2)G1516 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand): ref base T (96/104); new base A (5/8); total (101/112)
Fisher's exact test for biased strand distribution p-value = 5.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.67e-02

Only 100 of 219 total aligned reads displayed.

TGCCGCAGCATTTGATGGGCTTTATTCTCGGGATGTGGTTCCTGACGCAGGCTGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCGGACAACATTACCGATCCGCTTGAGACGTTGCCCGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGTTTGGTCGCCAGCAAGCGTGCAGGCTCCAGTAATCCATCGGCATCAAGATCCGGTCCCACCAGCGACTGAAAGACACGCATCTCTTTTTTCACTAGCGCGCTTTTCTGCTTATGCGGGAACATCGGGTCGAGATAAACC  >  NC_000913/3642084‑3642632

tGCCGCAGCATTTGATGGGCTTTATTCTCGGGATGTGGTTCCTGACGCAGGCTGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCGGACAACATTACCGATCCGCTTGAGACGTTGCCCGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGt                                                                                                                                                                                                                                                                                    <  4:538945/275‑1 (MQ=255)

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

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

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

                                   tgGTTCCTGACGCAGGCTGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCGGACAACATTACCGATCCGCTTGAGACGTTGCCCGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTAc                                                                                                                                                                                                                                                 <  2:206132/275‑1 (MQ=255)

                                         cTGACGCAGGCTGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCGGACAACATTACCGATCCGCTTGAGACGTTGCCCGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCAt                                                                                                                                                                                                                                           <  2:189472/275‑1 (MQ=255)

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

                                                 ggCTGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCGGACAACATTACCGATCCGCTTGAGACGTTGCCCGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTaca                                                                                                                                                                                                                                                          <  2:373756/252‑1 (MQ=255)

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

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

                                                                                   aCATTTACCGCGGTGCCGGACAACATTACCGATCCGCTTGAGACGTTGCCCGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACt                                                                                                                                                                                                                                                        >  2:425586/1‑220 (MQ=255)

                                                                                               gTGCCGGACAACATTACCGATCCGCTTGAGACGTTGCCCGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCg                                                                                                                                                                                     >  2:538402/1‑275 (MQ=255)

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

                                                                                                       caacaTTACCGGTCCGCTTGGGACGTTGACCGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGt                                                                                                                                                                             <  4:169989/275‑1 (MQ=255)

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

                                                                                                                    ccGCTTGAGACGTTGCCCGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGa                                                                                                                                                                                                                >  4:126030/1‑227 (MQ=255)

                                                                                                                    ccGCTTGAGACGTTGCCCGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCg                                                                                                                                                                >  4:607453/1‑275 (MQ=255)

                                                                                                                        ttGAGACGTTGCCCGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTaca                                                                                                                                                                                                                                                          <  2:349420/181‑1 (MQ=255)

                                                                                                                               gTTGCCCGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGtt                                                                                                                                                                                                                                                                                   <  5:536762/149‑1 (MQ=255)

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

                                                                                                                                  gCCCGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATc                                                                                                                                                                                                                                                                                <  5:192977/149‑1 (MQ=255)

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

                                                                                                                                     cGTNTNTACCAACGTGTTTGGTAAGATTGGTCTGGTCACGNTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGAt                                                                                                                                                                                                                            <  4:697073/198‑1 (MQ=255)

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

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

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

                                                                                                                                        cTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGATAACGCTAcc                                                                                                                                                                                                                                                                          <  6:840819/149‑1 (MQ=255)

                                                                                                                                         tataCCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGATAACGCTAccc                                                                                                                                                                                                                                                                         <  8:401885/149‑1 (MQ=255)

                                                                                                                                          attactaCGTGTTTGGTAAGATTGGTCTGGTTACCGCGGGCGGTGGAGTAATTATGCTGTTGATGGTGCCGTTGCTGAAACGCCTGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCAATAACGCTACCCt                                                                                                                                                                                                                                                                        <  6:590359/143‑1 (MQ=255)

                                                                                                                                               aaCGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTgtttt                                                                                                                                                                                                                                                                   >  7:752286/1‑149 (MQ=255)

                                                                                                                                               aaCGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAACAGGGTAGCGATAACGCTACCCTgcttt                                                                                                                                                                                                                                                                   <  8:173811/149‑1 (MQ=255)

                                                                                                                                                aCGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTgtttta                                                                                                                                                                                                                                                                  <  5:328599/149‑1 (MQ=255)

                                                                                                                                                 cGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGtt                                                                                                                                                                                                                                                                                   >  2:84179/1‑131 (MQ=255)

                                                                                                                                                 cGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTgttttag                                                                                                                                                                                                                                                                 <  6:549247/149‑1 (MQ=255)

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

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

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

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

                                                                                                                                                         gTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTgttttagt                                                                                                                                                                                                                                                                >  7:238747/1‑142 (MQ=255)

                                                                                                                                                         gTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTgttttagt                                                                                                                                                                                                                                                                <  8:238747/142‑1 (MQ=255)

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

                                                                                                                                                           aaGATTGGTCTGGTCACGTTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTg                                                                                                                                                                                                                                                       <  5:301983/149‑1 (MQ=255)

                                                                                                                                                            aGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTgg                                                                                                                                                                                                                                                      >  7:856330/1‑149 (MQ=255)

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

                                                                                                                                                              aTTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCg                                                                                                                                                                                                                                                    >  6:825084/1‑149 (MQ=255)

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

                                                                                                                                                                       gtcACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTTTCGCTACCCTGTTTTAGTTTTACACTGGCGTACCAGCAt                                                                                                                                                                                                                                           >  8:674910/1‑149 (MQ=255)

                                                                                                                                                                        tcACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATa                                                                                                                                                                                                                                          >  6:567688/1‑149 (MQ=255)

                                                                                                                                                                              tGGGCGTTGTAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATc                                                                                                                                                                                                                                    >  7:908366/1‑149 (MQ=255)

                                                                                                                                                                                   gTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGc                                                                                                                                                                                                                               >  5:966402/1‑149 (MQ=255)

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

                                                                                                                                                                                     tGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGa                                                                                                                                                                                                                             <  6:1080436/149‑1 (MQ=255)

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

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

                                                                                                                                                                                           agtGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCt                                                                                                                                                                                                                       <  8:752286/149‑1 (MQ=255)

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

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

                                                                                                                                                                                               aTGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTATCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGcgcg                                                                                                                                                                       <  4:632095/193‑1 (MQ=255)

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

                                                                                                                                                                                                        aTGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAg                                                                                                                                                                                                                   <  8:111832/140‑1 (MQ=255)

                                                                                                                                                                                                          ggTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACcgcg                                                                                                                                                                                                        >  7:399452/1‑149 (MQ=255)

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

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

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

                                                                                                                                                                                                               cGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCCCCACGCGTGTGGTCGCCAGCAAGCGGGCAGGCCCCGGTTACCCATGAGCACCAAGATCCGGTCCCACACGCGACTGcaa                                                                     >  4:388032/1‑275 (MQ=255)

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

                                                                                                                                                                                                                gTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCt                                                                                                                                                                                                                                               >  7:530191/1‑104 (MQ=255)

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

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

                                                                                                                                                                                                                     tGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTtggcgtggc                                                                                                                                                                                             >  6:1091270/1‑149 (MQ=255)

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

                                                                                                                                                                                                                           gCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACAtt                                                                                                                                                                                       >  5:706654/1‑149 (MQ=255)

                                                                                                                                                                                                                             aTGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCg                                                                                                                                                                                     >  5:873028/1‑149 (MQ=255)

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

                                                                                                                                                                                                                                  tGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTAcc                                                                                                                                                                                                                                                                          <  3:147596/59‑1 (MQ=255)

                                                                                                                                                                                                                                    cGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCgg                                                                                                                                                                              >  8:103094/1‑149 (MQ=255)

                                                                                                                                                                                                                                      aCGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTg                                                                                                                                                                            <  7:608632/149‑1 (MQ=255)

                                                                                                                                                                                                                                      aCGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTg                                                                                                                                                                            <  7:819271/149‑1 (MQ=255)

                                                                                                                                                                                                                                             aaaGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTa                                                                                                                                                                     <  5:567688/149‑1 (MQ=255)

                                                                                                                                                                                                                                                gCCATTAATTATTCTTGCAGAACGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTc                                                                                                                                                                  <  7:663834/149‑1 (MQ=255)

                                                                                                                                                                                                                                                  cATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCt                                                                                                                                                                                                                                               >  4:23178/1‑70 (MQ=255)

                                                                                                                                                                                                                                                   attaattaTTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCgg                                                                                                                                                               <  5:788426/149‑1 (MQ=255)

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

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

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

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

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

                                                                                                                                                                                                                                                       attattCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGgcgc                                                                                                                                                           >  5:397088/1‑149 (MQ=255)

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

                                                                                                                                                                                                                                                          attCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCttt                                                                                                                                                        <  6:1007688/149‑1 (MQ=255)

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

                                                                                                                                                                                                                                                             cTTGCAGAAATCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGTTTGGTCGCCAGCAAGCGTGCAGGCTCCAGTAATCCATCGGCATCAAGATCCGGTCCCACCAGCGACTGAAAGACACGCATCTCTTTTTTCACTAGCGCGCTTTTCTGCTTATGCggg                       <  2:412726/275‑1 (MQ=255)

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

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

                                                                                                                                                                                                                                                                      aGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGt                                                                                                                                            >  7:762378/1‑149 (MQ=255)

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

                                                                                                                                                                                                                                                                       gCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGtt                                                                                                                                           >  6:663106/1‑149 (MQ=255)

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

                                                                                                                                                                                                                                                                              aGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGTTTGGTCGc                                                                                                                                    <  7:975573/149‑1 (MQ=255)

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

TGCCGCAGCATTTGATGGGCTTTATTCTCGGGATGTGGTTCCTGACGCAGGCTGCCGCGTTCTTGCTGGGCGGCTATGTGGCAACATTTACCGCGGTGCCGGACAACATTACCGATCCGCTTGAGACGTTGCCCGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGTTTGGTCGCCAGCAAGCGTGCAGGCTCCAGTAATCCATCGGCATCAAGATCCGGTCCCACCAGCGACTGAAAGACACGCATCTCTTTTTTCACTAGCGCGCTTTTCTGCTTATGCGGGAACATCGGGTCGAGATAAACC  >  NC_000913/3642084‑3642632

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