BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,891,504	T→A	8.9%	intergenic (+21/‑111)	tnaB → / → mdtL	tryptophan transporter of low affinity/multidrug efflux system protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,891,504	0	T	A	8.9%	798.3 / 55.7	246	intergenic (+21/‑111)	tnaB/mdtL	tryptophan transporter of low affinity/multidrug efflux system protein
Reads supporting (aligned to +/- strand): ref base T (114/110); new base A (8/14); total (122/124)
Fisher's exact test for biased strand distribution p-value = 2.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.31e-05

Only 100 of 250 total aligned reads displayed.

CGTTGTTGTATCTGATCTTCCCGAACGGCTTTATTTACGGGATCGGCGGTGCCGGGCTGTGCGCCACCATCTGGGCGGTCATTATTCCCGCAGTGCTTGCAATCAAAGCTCGCAAGAAGTTTCCCAATCAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCTGGTTTTACTTTATCCCGCCGGGATTGATATGTACCTCGTTGGTTTACCGCGCATCGCCGCCGATCTCAATGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGGCAGCTGCG  >  NC_000913/3891239‑3891776

cGTTGTTGTATCTGATCTTCCCGAACGGCTTTATTTACGGGATCGGCGGTGCCGGGCTGTGCGCCACCATCTGGGCGGTCATTATTCCCGCAGTGCTTGCAATCAAAGCTCGCAAGAAGTTTCCCAATCAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGct                                                                                                                                                                                                                                                                         >  2:151607/1‑275 (MQ=255)

          tCTGATCTTCCCGAACGGCTTTATTTACGGGATCGGCGGTGCCGGGCTGTGCGCCCCCATCTGGGCGGTCATTATTCCCGCAGTGCTTGCAATCAAAGCTCGCAACAAGTTTCCCAATCAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTTTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTTTTTGACGTGTTACCTACATTTGGCTAAATCCTTCAAGAAGCCCGCCATTCGCTGGCTTCTTGCcttt                                                                                                                                                                                                                                                               >  2:485796/1‑273 (MQ=255)

                                           cggcggTGCCGGGCTGTGCGCCACCATCTGGGCGGTCATTATTCCCGCAGTGCTTGCAATCAAAGCTCGCAAGAAGTTTCCCAATCAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCaaa                                                                                                                                                                                                                                         >  4:695288/1‑264 (MQ=255)

                                                        ctggcggcCACAATCGGGGAGGTCATTATTCCCGAAGTTCTTGCAATAAAAGAGAGAAAGAAGTTTCCAAATCAGATGGTCACGGTCGGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATGTTGTGCTGGTTCGGCAACGGCTTTAACGTGTTACCTAAATTTGGCTAAATACTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTc                                                                                                                                                                                                                 <  1:136348/269‑1 (MQ=255)

                                                                    aTCTGGGCGGGCATTATTCCTGCAGTGCTGGGAATGAAAGCCTGCAAAAAGTTTGCCAAGAAGAGGTGAACGGTATGGGGCGGCAATCTTATTCCGGCGATTGTAAGTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTGACGTGTTACCTAAATTTTGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTCCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGc                                                                                                                                                                                                      <  4:485794/274‑1 (MQ=255)

                                                                       tgtgccgacATTATTCCCGCAGTGATTGCAATCAAAGCTCGCAAGAAGTTTCCCAATACGATGTTCCCGGTTGGGGGCTGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGTTTTTGTGCTGGTTCTGCAACGTCTTTAACGTGTTACCTAAATTTTGCTAAATCCTTCAAGAAGCCAGCCATTCGCGGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGTTCCTCCTTCACCACGTATGCTTTg                                                                                                                                                                                                  <  4:713045/267‑1 (MQ=255)

                                                                          gCGGTCATTATTCCCGTGTTGCTTGCAATCGAAGCTCCCATGAAGTTTCCCAATCAGATGTTCACGGTCTGGGGCGGCAATCTGTTTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGTTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCg                                                                                                                                                                                                <  4:322160/274‑1 (MQ=255)

                                                                             gTCATTATTCCCGCAGTGCTTGCAATCAAAGCTCGCAAGAAGTTTCCCAATCAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCa                                                                                                                                                                                             >  2:209016/1‑274 (MQ=255)

                                                                                        cGCAGTGCTTGCAATCAAAGCTCGCAAGAAGTTTCCCAATCAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTCCGTCACCTTACTATCAg                                                                                                                                                                                 >  2:696355/1‑275 (MQ=255)

                                                                                                ttGCCATCCAAGGTTGGAAGGAGTTTCCCAGGCAGATGTTCCCAGTCCTTGGCTGCAAGCTTATTCCGGCGATTGTCATTCTCTTTTGTATAACGGTGATTTTGTGCTGGTTAGGCAACGTCGTTAACGTGTTACCTAAATGTGGCTAAATCCTTCCAGAAGCCAGCCATTCGCTGGCTTCTTACCTCTCAGGAAATCACTTATGTCCAAATGGCGACTCGCCTGATCCTCCTTCACCGCGTATGCTTTGCGTCTCCTTACTATCAGGACGCttt                                                                                                                                                                         <  4:199389/275‑1 (MQ=255)

                                                                                                                         tCCCAATCAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCAtt                                                                                                                                                                                                                                                                                 >  7:767269/1‑146 (MQ=255)

                                                                                                                            cAATCAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCGAATGGCTgg                                                                                                                                                                                                                                                                           <  8:494278/149‑1 (MQ=255)

                                                                                                                             aaTCAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTgg                                                                                                                                                                                                                                                                           >  5:11748/1‑148 (MQ=255)

                                                                                                                              aTCAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGc                                                                                                                                                                                                                                                                          >  6:768480/1‑148 (MQ=255)

                                                                                                                               tCAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCGAATggctggct                                                                                                                                                                                                                                                                         <  5:481788/148‑1 (MQ=255)

                                                                                                                               tCAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGctt                                                                                                                                                                                                                                                                        >  8:894204/1‑149 (MQ=255)

                                                                                                                               tCAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTGCAAGAAGCCAGCGAATGGCTGGctt                                                                                                                                                                                                                                                                        <  7:474402/149‑1 (MQ=255)

                                                                                                                                 aGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCGAATGGCTGGcttc                                                                                                                                                                                                                                                                       <  6:563229/148‑1 (MQ=255)

                                                                                                                                  gATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGcttct                                                                                                                                                                                                                                                                      <  5:768480/148‑1 (MQ=255)

                                                                                                                                         acgtgcTGGGGCGGAAAACTTATTCCGGCGAGTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGGTTTACCTAAATTTTGCTAAACCCTTCCAGAAGCCAGCGATTCGCTGTCTTCTTTCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCTGGtt                                                                                                                                <  4:184623/270‑1 (MQ=255)

                                                                                                                                         aCGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCTTCACCTTCCTATCAGGACGCTTTAGCCCATTTCCCGCTTTTTGATTTTTTGTTTTTCCCTGGt                                                                                                                                 >  4:374839/1‑274 (MQ=255)

                                                                                                                                            gTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAg                                                                                                                                                                                                                                                            >  5:246134/1‑148 (MQ=255)

                                                                                                                                               tGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGaa                                                                                                                                                                                                                                                         <  5:738120/148‑1 (MQ=255)

                                                                                                                                                ggggCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGaaa                                                                                                                                                                                                                                                        <  6:426953/148‑1 (MQ=255)

                                                                                                                                                 gggCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCATGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCt                                                                                                                                                                                                     <  4:149033/198‑1 (MQ=255)

                                                                                                                                                  ggcggcAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGTTTCGGCAACGTCTTTAACGTGTTACCTAAATTTTGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAAta                                                                                                                                                                                                                                                      >  8:555784/1‑147 (MQ=255)

                                                                                                                                                       cAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCctc                                                                                                                                                                                                                                                                >  2:761896/1‑133 (MQ=255)

                                                                                                                                                       cAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCctc                                                                                                                                                                                                                                                                >  3:183532/1‑133 (MQ=255)

                                                                                                                                                       cAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCTGGTTTTACTTTATCccgc                                                                                                                  >  4:425979/1‑275 (MQ=255)

                                                                                                                                                       cAATCTTATTCCGGCGATTGTCATTCTCTTTAGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCTGGTTTTACTTTATCccgc                                                                                                                  <  2:50269/275‑1 (MQ=255)

                                                                                                                                                       cAAGCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCGGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCTGGTTTTACTTTATCccgc                                                                                                                  <  2:294871/275‑1 (MQ=255)

                                                                                                                                                         aTCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCGAATGGCTGGcttct                                                                                                                                                                                                                                                                      >  5:606153/1‑125 (MQ=255)

                                                                                                                                                         aTCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCATGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTGACTAt                                                                                                                                                                                    <  1:101320/207‑1 (MQ=255)

                                                                                                                                                          tCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGt                                                                                                                                                                                                                                              >  6:893001/1‑148 (MQ=255)

                                                                                                                                                           cTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCTGGTTTTACTTTATCccgc                                                                                                                  <  2:692485/271‑1 (MQ=255)

                                                                                                                                                               ttCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCaaa                                                                                                                                                                                                                                         <  6:606264/148‑1 (MQ=255)

                                                                                                                                                               ttCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCaaa                                                                                                                                                                                                                                         <  5:893001/148‑1 (MQ=255)

                                                                                                                                                                    gCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGc                                                                                                                                                                                                                              >  4:562856/1‑154 (MQ=255)

                                                                                                                                                                       aTTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTc                                                                                                                                                                                                                                >  7:612305/1‑149 (MQ=255)

                                                                                                                                                                        ttGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGa                                                                                                                                                                                                                          <  2:562858/154‑1 (MQ=255)

                                                                                                                                                                        ttGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCCAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGa                                                                                                                                                                               <  2:42540/197‑1 (MQ=255)

                                                                                                                                                                            cATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTg                                                                                                                                                                                                                           >  7:226941/1‑149 (MQ=255)

                                                                                                                                                                            cATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTg                                                                                                                                                                                                                           >  8:371875/1‑149 (MQ=255)

                                                                                                                                                                              ttCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGAt                                                                                                                                                                                                                         <  8:256713/149‑1 (MQ=255)

                                                                                                                                                                              ttCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGAt                                                                                                                                                                                                                         >  8:692925/1‑149 (MQ=255)

                                                                                                                                                                               tctcTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGAt                                                                                                                                                                                                                         >  5:434280/1‑148 (MQ=255)

                                                                                                                                                                                  ctTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTg                                                                                                                                                                                                                                                                    >  5:552546/1‑102 (MQ=255)

                                                                                                                                                                                  ctTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGAtcct                                                                                                                                                                                                                      >  6:233096/1‑148 (MQ=255)

                                                                                                                                                                                   tttGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGAtcctc                                                                                                                                                                                                                     <  6:293833/148‑1 (MQ=255)

                                                                                                                                                                                    ttGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGAtcctcc                                                                                                                                                                                                                    >  6:125524/1‑148 (MQ=255)

                                                                                                                                                                                       gTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGAtc                                                                                                                                                                                                                        <  4:800531/141‑1 (MQ=255)

                                                                                                                                                                                        tataACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGa                                                                                                                                                                                                                          <  5:829575/138‑1 (MQ=255)

                                                                                                                                                                                          taACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGaa                                                                                                                                                                                                                                                         <  5:156202/105‑1 (MQ=255)

                                                                                                                                                                                           aaCCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGcttctt                                                                                                                                                                                                                                                                     >  2:75963/1‑92 (MQ=255)

                                                                                                                                                                                           aaCCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGc                                                                                                                                                                                                                              >  4:85229/1‑131 (MQ=255)

                                                                                                                                                                                               gTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGaa                                                                                                                                                                                                                                                         <  8:144680/100‑1 (MQ=255)

                                                                                                                                                                                                tGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACg                                                                                                                                                                                                           >  6:205844/1‑145 (MQ=255)

                                                                                                                                                                                                 gATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATCGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCTGGTTTTACTTTATTCCGCCTGGATTGATGTGTACCTCGTTGGTTTACCGCGCATcgccc                                                                         >  1:179608/1‑273 (MQ=255)

                                                                                                                                                                                                  aTTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTNTGTCCAAATGGCAACTCGCCTGANNCTCCTTCANCACGTATGCtt                                                                                                                                                                                                    >  3:55176/1‑150 (MQ=255)

                                                                                                                                                                                                  aTTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTcac                                                                                                                                                                                                               >  2:406511/1‑139 (MQ=255)

                                                                                                                                                                                                    tttGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCttt                                                                                                                                                                                                   <  7:894204/149‑1 (MQ=255)

                                                                                                                                                                                                      tgtgCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCg                                                                                                                                                                                                <  1:55177/150‑1 (MQ=255)

                                                                                                                                                                                                          cTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTc                                                                                                                                                                                              >  5:652847/1‑148 (MQ=255)

                                                                                                                                                                                                             gTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGc                                                                                                                                                                                                 >  7:554552/1‑142 (MQ=255)

                                                                                                                                                                                                             gTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCt                                                                                                                                                                                          >  8:210154/1‑149 (MQ=255)

                                                                                                                                                                                                                      cGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCAc                                                                                                                                                                                                                                                    <  6:606140/82‑1 (MQ=255)

                                                                                                                                                                                                                      agTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCTGGTTTTACTTTATCCCGCCGGGATTGata                                                                                                       <  4:144558/222‑1 (MQ=255)

                                                                                                                                                                                                                           ttAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACg                                                                                                                                                                             >  5:297753/1‑148 (MQ=255)

                                                                                                                                                                                                                             aaCGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCTGGTTTTACTTTATCCCGCCGGGATTGATATGTACCTCGTTGGTTTACCGCGCATCGCCGCCGATCTCAATGCCAGCGAAGCGCAGtt                                             <  2:80598/274‑1 (MQ=255)

                                                                                                                                                                                                                                 tgtTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTg                                                                                                                                                                                                  >  2:100830/1‑121 (MQ=255)

                                                                                                                                                                                                                                  gtTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGc                                                                                                                                                                      <  6:834988/148‑1 (MQ=255)

                                                                                                                                                                                                                                   ttACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCTGGTTTTACTTTATCCCGCCGGGATTGATATGTACCTCGTTGGTTTACCGCGCATCGCCGCCGATCCCAATGCCAGCGAAGCGCAGTTGCATAtt                                      >  2:603991/1‑275 (MQ=255)

                                                                                                                                                                                                                                   ttACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCTGGTTTTACTTTATCCCGCCGGGATTGATATGTACCTCGTTGGTTTACCGCCCATCGCCGCCGGTCCCAATGCCAGCGAAGCGCAGTTGCATAtt                                      >  4:626278/1‑275 (MQ=255)

                                                                                                                                                                                                                                   ttACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCTGGTTTTACTTTATCCCGCCGGGATTGATATGTACCTCGTGGGTTTACCGCCCCTCGCCTCCCATCTCAATTCCCGCGAACCGCCGTGTCatat                                       >  4:687426/1‑274 (MQ=255)

                                                                                                                                                                                                                                   ttACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCGGGTTTTACTTTATCCCGCCGGGATTGATATGTACCTCGTTGGTTTACCGCGCATCGCCGCCGATCTCAATGCCAGCGCAGCGCAGTTGCATAtt                                      >  2:437264/1‑275 (MQ=255)

                                                                                                                                                                                                                                    tACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCTGGTTTTACTTTATCccgcc                                                                                                                 >  4:407346/1‑199 (MQ=255)

                                                                                                                                                                                                                                     aCCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTc                                                                                                                                                                                              <  4:100829/121‑1 (MQ=255)

                                                                                                                                                                                                                                       ntNAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCNCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCTGGTTTTACTTTATCccgc                                                                                                                  <  4:761445/194‑1 (MQ=255)

                                                                                                                                                                                                                                        tAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGAc                                                                                                                                                                              >  7:190212/1‑134 (MQ=255)

                                                                                                                                                                                                                                        tAAATTTGGCTAAATCCATCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCTGGTTTTACTTTATCCCGCCGGGATTGATATGTACCTCGTTGGTTTACCGCGCATCGCCGCCGATCTCAATGCCAGCGAAGCGCAGTTGCATATTGCGtt                                 <  2:534909/275‑1 (MQ=255)

                                                                                                                                                                                                                                            tttGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCTGGTTTTACTTTATCccgcc                                                                                                                 >  1:184509/1‑191 (MQ=255)

                                                                                                                                                                                                                                            tttGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCTGGTTTTACTTTATCCCGCCGGGATTGa                                                                                                         >  2:271792/1‑199 (MQ=255)

                                                                                                                                                                                                                                              tGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCt                                                                                                                                                          <  6:112838/148‑1 (MQ=255)

                                                                                                                                                                                                                                                gCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCtttt                                                                                                                                                       >  8:100615/1‑149 (MQ=255)

                                                                                                                                                                                                                                                gCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCTGGTTTTACTTTATCCCGCCGGGATTGATATg                                                                                                     <  4:271791/199‑1 (MQ=255)

                                                                                                                                                                                                                                                  tAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTa                                                                                                                                                                                        >  7:668/1‑114 (MQ=255)

                                                                                                                                                                                                                                                  tAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTg                                                                                                                                                     >  8:413149/1‑149 (MQ=255)

                                                                                                                                                                                                                                                   aaaTCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGa                                                                                                                                                    <  7:100615/149‑1 (MQ=255)

                                                                                                                                                                                                                                                       ccTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATTTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGAttt                                                                                                                                                 >  5:221918/1‑148 (MQ=255)

                                                                                                                                                                                                                                                         ttCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCt                                                                                                                                                                                          >  7:753248/1‑105 (MQ=255)

                                                                                                                                                                                                                                                         ttCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCt                                                                                                                                                                                          <  8:753248/105‑1 (MQ=255)

                                                                                                                                                                                                                                                          tCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGAtcctc                                                                                                                                                                                                                     >  6:737027/1‑77 (MQ=255)

                                                                                                                                                                                                                                                          tCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCTGGTTTTACTTTATc                                                                                                                      >  4:500806/1‑172 (MQ=255)

                                                                                                                                                                                                                                                            aagaagCCAGCGAATGGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCtt                                                                                                                                                                          <  8:124967/118‑1 (MQ=255)

                                                                                                                                                                                                                                                            aagaagCCAGCGAATGGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCtt                                                                                                                                                                          >  7:124967/1‑118 (MQ=255)

                                                                                                                                                                                                                                                            aagaagCCAGCGAATGGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGt                                                                                                                                            >  6:544214/1‑148 (MQ=255)

                                                                                                                                                                                                                                                                  ccagccaTTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGccc                                                                                                                                     >  8:839809/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                   cagccaTTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCt                                                                                                                                    >  7:34868/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                      ccaTTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGAtt                                                                                                                                                  <  7:385301/132‑1 (MQ=255)

                                                                                                                                                                                                                                                                      ccaTTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCTGGTTTTACTTTATCccgcc                                                                                                                 >  4:238335/1‑165 (MQ=255)

CGTTGTTGTATCTGATCTTCCCGAACGGCTTTATTTACGGGATCGGCGGTGCCGGGCTGTGCGCCACCATCTGGGCGGTCATTATTCCCGCAGTGCTTGCAATCAAAGCTCGCAAGAAGTTTCCCAATCAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCTGGTTTTACTTTATCCCGCCGGGATTGATATGTACCTCGTTGGTTTACCGCGCATCGCCGCCGATCTCAATGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGGCAGCTGCG  >  NC_000913/3891239‑3891776

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