BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,477,878	T→A	5.2%	intergenic (+20/+29)	slyX → / ← slyD	phi X174 lysis protein/FKBP‑type peptidyl prolyl cis‑trans isomerase (rotamase)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,477,878	0	T	A	5.2%	345.2 / 6.9	104	intergenic (+20/+29)	slyX/slyD	phi X174 lysis protein/FKBP‑type peptidyl prolyl cis‑trans isomerase (rotamase)
Reads supporting (aligned to +/- strand): ref base T (57/38); new base A (4/2); total (63/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.10e-01

Only 100 of 105 total aligned reads displayed.

GGTGAATCATGCAGGATTTATCATTGGAAGCACGCCTGGCAGAGCTGGAGAGCCGACTGGCTTTTCAGGAAATCACCATTGAAGAACTGAACGTCACGGTGACCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGTGGTGATCGTGCGCGCCGTGAACGTGACCATGAGCCAGTTCTTCTTCAGTCGCTTCGCGAATCGCCACAACTTCAACGTTGAATTTCAGGTTCTGACCGGCCAGCATGTGGTTACCATCAACCACGACG  >  NC_000913/3477632‑3478150

ggttaatcagttaGGTTTTATCATTGGAAGAACGCCTTGAAGCGCTGTAGAAGCGACAGGCTTTTCAAGAAATCACCATTGAAGAACTGAGCGTCACTGTGACCGCGCATGAAATGGAGATGGGGAAACTTCGCGGTCATCGGCGTCTGCGGCCCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACGTTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTa                                                                                                                                                                                                                                                      <  4:172515/263‑1 (MQ=255)

 gTGAATCATGCAGGATTTATCATTGGAAGCACGCCTGGCAGAGCTGGAGAGCCGACTGGCTTTTCAGGAAATCACCATTGAAGAACTGAACGTAACGGTGACCGCTCATGCAATGGAGATGGCACAACTGCGCGATCATCTGCGTCTGCTGACCGAGCAGTTAAACGCCAGCCATCCGTCGAACATCGCGTCGCTGGCTGTCGAAACGCCAACGCCCCATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTTCCTTTTTCGGTAt                                                                                                                                                                                                                                                     >  2:277290/1‑275 (MQ=255)

              gATTTATCATTGGAAGCACGCCTGGCAGAGCTGGAGAGCCGACTGGCTTTTCAGGAAATCACCATTTAAGAACTGAACGTAACGGTGACCGCTCATGACATGGAGATGGAGAAACTTCTCGAACATCTTCGTCTGCTGACCAGGAAGTTAAACGCCAGCCACACGTCAAAACTCGCGTCCCCGGCTGAAGAAACAGCCCCCTCACCTTATTTCGACGCAAAAATAGCGGGCATTCACCCCGTTTCTCTCACTTTTTTGGGTTTAGttcaacacc                                                                                                                                                                                                                                         >  4:596667/1‑266 (MQ=255)

                        tGGAAGCACGCCTGGCAGAGCTGGAGAGCCGACTGGCTTTTCAGGAAATCACCATTGAAGAACTGAACGTCACGGTGACCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCAGTCGAACCTCGCGTCGCAGGCTGAAGAAACGCCACCGCCACACTCTTGAGGCGTAAAAAACACGGGGCTTCCCCGCTTTTTTGTCCCTTTTTCGGTATTAGTGGCACCCGACACCCCCGtg                                                                                                                                                                                                                              >  4:387647/1‑274 (MQ=255)

                             gCACGCCTGGCAGAGCTGGAGAGCCGACTGGCTTTTCAGGAAATCACCATTGAAGAACTGAACGTCACGGTGACCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACAGCCGTTGCCt                                                                                                                                                                                                                          >  2:205012/1‑274 (MQ=255)

                              aaCGCCTGGCAGGGCTGGAGAGCCGACTGGCTTTTCAGGAAATCACCATTGAAGAACTGAACGTCCCGGGGACCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGGGAAGTTAAAAGCCGGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCtt                                                                                                                                                                                                                         <  2:387649/273‑1 (MQ=255)

                                 gCCTGGCAGAGCTGGAGAGCCGACTGGCTTTTCAGGAAATCACCATTGAAGAACTGAACGTCACGGTGACCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTTCCTTTAc                                                                                                                                                                                                                      >  2:209832/1‑274 (MQ=255)

                                                 gagCCGACTGGCTTTTCAGGAAATCACCATTGAAGAACTGAACGTCACGGTGACCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTCGTGGCAACCGCAACCACCGTTGCCTTTACCGCCACAGCAGCCttac                                                                                                                                                                                                     >  2:533802/1‑273 (MQ=255)

                                                      gACTGGCTTTTCAGGAAATCACCATGGAAGAACTGAACGTCACGGTGACCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCAcc                                                                                                                                                                                                <  2:489234/275‑1 (MQ=255)

                                                                                      cTGAACGTCACGGTGACCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGAccgc                                                                                                                                                                >  2:525310/1‑275 (MQ=255)

                                                                                         aaCGTCACGGTGACCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGAGTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCttt                                                                                                                                                                                                                        >  2:242715/1‑216 (MQ=255)

                                                                                         aaCGTCACGGTGACCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACAGCCACATTATTGAGGCGTAAAAACAGCGGGGAGTCCCCGCTTTTTTGTCACTTTTTCGGTATTTGATGCAACCGCAACCCCCATTGCACTTACCGCAACAGCAGACTTCGCCAACGGGTTCATGACCGTGCACGTGCCCATTAccgcccc                                                                                                                                                             >  4:663923/1‑273 (MQ=255)

                                                                                         aaCGTCACGGTGACCGCTCATGAAATGGAGATGGCGAAACTGCGCGACCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCCCTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTTCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTCATCGTGGCCATGAcccacgc                                                                                                                                                             >  4:339599/1‑270 (MQ=255)

                                                                                                   tGACCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGAt                                                                                                                                                                                                                                                                                  >  6:540510/1‑148 (MQ=255)

                                                                                                    gACCGCTCATGAAATGGAGAGGGAGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACcgtcg                                                                                                                                                  <  2:339601/275‑1 (MQ=255)

                                                                                                     aCCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCTTGATCGTGGCCATGACCGCCGCCGCAACcggcgt                                                                                                                                                 >  2:541497/1‑275 (MQ=255)

                                                                                                        gCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCg                                                                                                                                                                                                                                                                            >  8:488699/1‑149 (MQ=255)

                                                                                                           cATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCt                                                                                                                                                                                                                                                                          <  5:154115/148‑1 (MQ=255)

                                                                                                           cATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGAATCCCCGCt                                                                                                                                                                                                                                                                          <  5:314464/148‑1 (MQ=255)

                                                                                                             tGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCtttt                                                                                                                                                                                                                                                                       >  7:71530/1‑149 (MQ=255)

                                                                                                               aaaTGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCtttttt                                                                                                                                                                                                                                                                     <  7:414708/149‑1 (MQ=255)

                                                                                                               aaaTGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGAATCCCCGCttttt                                                                                                                                                                                                                                                                      <  5:437608/148‑1 (MQ=255)

                                                                                                                aaTGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCtttttt                                                                                                                                                                                                                                                                     >  5:443394/1‑148 (MQ=255)

                                                                                                                       aTGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACtt                                                                                                                                                                                                                                                              >  5:168758/1‑148 (MQ=255)

                                                                                                                        tGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACttt                                                                                                                                                                                                                                                             >  5:245026/1‑148 (MQ=255)

                                                                                                                          gCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACAGTATTGAGGCGTAAAAAAAGCGGGGATTACACGCTTTTTTGTCACTTTTTCGGTATTAGGGGCACGCGCAACCGCGCTTGCCTTTACTGCCACCTAAGCAGTCGCCACACTGTTCAGGACCGTGGCCGTGGCCATGAACGGCTGAACAAACGGCGGGAGTGTGATCgtggctgctag                                                                                                                            >  4:76897/1‑268 (MQ=255)

                                                                                                                                cTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTAt                                                                                                                                                                                                                                                     >  6:640170/1‑148 (MQ=255)

                                                                                                                                 tGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTa                                                                                                                                                                                                                                                   >  7:86457/1‑149 (MQ=255)

                                                                                                                                    gcgATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTg                                                                                                                                                                                                                                                <  7:194579/149‑1 (MQ=255)

                                                                                                                                       atcatcTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGgc                                                                                                                                                                                                                                              >  5:289205/1‑148 (MQ=255)

                                                                                                                                        tcatcTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTAtt                                                                                                                                                                                                                                                    <  6:846256/141‑1 (MQ=255)

                                                                                                                                        tcatcTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTAtt                                                                                                                                                                                                                                                    >  5:846256/1‑141 (MQ=255)

                                                                                                                                        tcatcTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGgcaa                                                                                                                                                                                                                                            <  7:488699/149‑1 (MQ=255)

                                                                                                                                                  ctgctgACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGTGGTGATCGTGCGCGCCGTGAACGTGACCAt                                                                                                    >  2:18665/1‑275 (MQ=255)

                                                                                                                                                    gctgACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGt                                                                                                                                                                                  >  2:417293/1‑195 (MQ=255)

                                                                                                                                                        aCCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATc                                                                                                                                                                              <  4:417291/195‑1 (MQ=255)

                                                                                                                                                         ccGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGc                                                                                                                                                                                                                            >  6:582325/1‑148 (MQ=255)

                                                                                                                                                          cGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGc                                                                                                                                                                                                    >  2:499656/1‑171 (MQ=255)

                                                                                                                                                           gagaAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCt                                                                                                                                                                                                                          <  5:608850/148‑1 (MQ=255)

                                                                                                                                                             gaaggaaaaaCCCACCCAGCCGTCGAACGACGCCTCGCGGGGTGACGAACCGCCACGGAACCATTATTGAGGCGTAAAATACGCGGGGATGCCCCGCTTTTTGGTCACTTTTTCGGTATTAGGGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGTGGTGATCGTGCGCGCCGTGAACGTGACCATGAGCCAGttc                                                                                          <  4:525308/268‑1 (MQ=255)

                                                                                                                                                              aaGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCAcc                                                                                                                                                                                                <  4:499654/171‑1 (MQ=255)

                                                                                                                                                                  tAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGc                                                                                                                                                                                                                   >  5:506833/1‑148 (MQ=255)

                                                                                                                                                                   aaaaGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCca                                                                                                                                                                                                                 >  7:545038/1‑149 (MQ=255)

                                                                                                                                                                       gccagccagccGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCAcag                                                                                                                                                                                                              >  5:449970/1‑148 (MQ=255)

                                                                                                                                                                        ccagccagccGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCAcagca                                                                                                                                                                                                            <  7:255855/149‑1 (MQ=255)

                                                                                                                                                                           gccagccGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCt                                                                                                                                                                                                                          >  4:322650/1‑132 (MQ=255)

                                                                                                                                                                           gccagccGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGGGGTGATCGGGCGCGCCGTGAACGTGACCCTGAGCCAGTTCTTCGTCAGTCCCTTc                                                                           >  2:121816/1‑275 (MQ=255)

                                                                                                                                                                           gccagccGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCATTTTTGTCACTTTTTCGGTATTAATTGCAACCGCAACCGCCGTTGCCTTTACCGCCCCAGCAACCTTCGCCACCGTGTTCATGACCGTGATCGTTGCCATGACCGCCGCAGCAACCGTCGTGGTCTTGATCGTGGTGGTGATCGTGCGCGCCCTCAACGTGACCATGAGCCAGTTCTTCTTTAGTccgtt                                                                            >  2:27353/1‑270 (MQ=255)

                                                                                                                                                                              agccGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTc                                                                                                                                                                                                      <  8:71530/149‑1 (MQ=255)

                                                                                                                                                                               gccGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTAc                                                                                                                                                                                                                      <  2:322652/132‑1 (MQ=255)

                                                                                                                                                                                 cggcGAACATGGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGCATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGTGGTGATCGTGCGCGCCGTGAACGTGACCATGAGCCAGTTCTTCTTCAGTCGCTTCGCGAAt                                                                     <  2:424990/272‑1 (MQ=255)

                                                                                                                                                                                   tCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCa                                                                                                                                                                                                  >  6:254305/1‑148 (MQ=255)

                                                                                                                                                                                    cGAACATCGCTTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAACGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAATCGCAACCGCCGTTGCCTTTACCg                                                                                                                                                                                                                    >  5:743237/1‑129 (MQ=255)

                                                                                                                                                                                    cGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCg                                                                                                                                                                                                                    <  6:743237/129‑1 (MQ=255)

                                                                                                                                                                                       aCATCGGGTCGCAGGCTGAAGAAACGCCACCGCCAAATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGTGGTGATCGTGCGCGCCGTGAACGTGACCATGAGCCAGTTCTTCTTCAGTCGCTTCGCGAATCGCcac                                                               <  2:699458/275‑1 (MQ=255)

                                                                                                                                                                                       aCATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCgt                                                                                                                                                                                              >  5:361851/1‑148 (MQ=255)

                                                                                                                                                                                       aCATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCgt                                                                                                                                                                                              >  6:124484/1‑148 (MQ=255)

                                                                                                                                                                                          tCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGtt                                                                                                                                                                                           >  6:308019/1‑148 (MQ=255)

                                                                                                                                                                                          tCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGTGGTGGTCGTGCGCGCCGTGAACGTGACCGTGAGCCAGTTCTTTTTCAGTCGCTTCGCGAATCGCCACAAc                                                            >  2:119986/1‑275 (MQ=255)

                                                                                                                                                                                           cgcgTCGCAGGCTGAAGAAACGCTACCGCCACATTATTGAGGCGTAAACAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCAAAGCAGCCTTCGCCACCGTGTTc                                                                                                                                                                                          >  6:848256/1‑148 (MQ=255)

                                                                                                                                                                                              gTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATg                                                                                                                                                                                       <  6:849789/148‑1 (MQ=255)

                                                                                                                                                                                                   aGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGt                                                                                                                                                                                  <  5:540510/148‑1 (MQ=255)

                                                                                                                                                                                                    ggCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTg                                                                                                                                                                                 <  5:349445/148‑1 (MQ=255)

                                                                                                                                                                                                     gCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGa                                                                                                                                                                                <  6:449970/148‑1 (MQ=255)

                                                                                                                                                                                                     gCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGAt                                                                                                                                                                               >  8:492193/1‑149 (MQ=255)

                                                                                                                                                                                                     gCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGAt                                                                                                                                                                               >  8:49460/1‑149 (MQ=255)

                                                                                                                                                                                                        gaagaaACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGt                                                                                                                                                                            >  7:207983/1‑149 (MQ=255)

                                                                                                                                                                                                           gaaACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTgg                                                                                                                                                                          <  6:813696/148‑1 (MQ=255)

                                                                                                                                                                                                               cgccaccgccacATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCcac                                                                                                                                                                                                                >  2:144021/1‑106 (MQ=255)

                                                                                                                                                                                                                gccaccgccacATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATg                                                                                                                                                                     <  5:124484/148‑1 (MQ=255)

                                                                                                                                                                                                                 ccaccgccacATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGa                                                                                                                                                                    >  6:762530/1‑148 (MQ=255)

                                                                                                                                                                                                                  caccgccacATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTc                                                                                                                                                                                          <  7:887288/125‑1 (MQ=255)

                                                                                                                                                                                                                   accgccacATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCAcagca                                                                                                                                                                                                            <  4:144020/106‑1 (MQ=255)

                                                                                                                                                                                                                   accgccacATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGAccg                                                                                                                                                                 >  7:750485/1‑149 (MQ=255)

                                                                                                                                                                                                                   accgccacATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGTGGTGATCGTGCGCGCCGTGAACGTGACCATGAGCCAGTTCTTCTTCAGTCGCTTCGCGAATCGCCACAACTTCAACGTTGAATTTCAGGTTCTGa                                   >  2:176258/1‑275 (MQ=255)

                                                                                                                                                                                                                     cgccacATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGAccgcc                                                                                                                                                               >  8:79952/1‑149 (MQ=255)

                                                                                                                                                                                                                        cacaTTATTGAGGCGCAAAAAAAGGGGGTCGTCCCCGCTTTTTTGTCCCTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCTCAACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGTGGTGATCGTGCGCGCCGTGAACGTGACCATGAGCCAGTTCTTCTTCAGTCGCTTCGCGAATCGCCACAACTTCAACGTTGAATTTCAGGTTCTGACCgg                               <  2:76898/274‑1 (MQ=255)

                                                                                                                                                                                                                                  aGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACcgtc                                                                                                                                                   >  6:443951/1‑148 (MQ=255)

                                                                                                                                                                                                                                   ggCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACcgtcgt                                                                                                                                                 <  7:492193/149‑1 (MQ=255)

                                                                                                                                                                                                                                   ggCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACcgtcg                                                                                                                                                  <  5:308019/148‑1 (MQ=255)

                                                                                                                                                                                                                                   gcctaaaaaaaaGCGGGGGTTCCCGCCTTTTTTTTCACTTTTACTGTGTTGGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCGTGGCCACCGGGTTCATGACCGTGATCTTGGCCATGACAGCCGCAGCAACCGTCGTGTTCGTGATCGTGGTGGTGATCGTGCGCGCCGTGAACGTGACCATGAGCCAGTTCTTCTTCAGTCGCTTCGCGAATCGCCACAACTTCAACGTTGAATTTCAGGTTCTGACCGGCCAGCATGTGGt                   <  4:176257/270‑1 (MQ=255)

                                                                                                                                                                                                                                    gCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGTGGTGATCGTGCGCGCCGTGAACGTGACCATg                                                                                                   >  4:122255/1‑194 (MQ=255)

                                                                                                                                                                                                                                    gCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGCCGTGGTCGTGATCGTGGTGGTGATCGTTCGCGCCGTGAACGTGACCGTGACCCAGTTCTTCTTCAGTCGCTTAGCGAATCGCTACAACTTCAACGTTGACTTTCCGGTTCTGACCGGCCGCCATGTGGt                   >  4:433248/1‑274 (MQ=255)

                                                                                                                                                                                                                                      gTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCgtcgtggt                                                                                                                                              >  7:495008/1‑149 (MQ=255)

                                                                                                                                                                                                                                       tAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTCGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTCATCGTCGCCCTGACCGCCGCAGCAACCGTCGTGGTCGTGATAGTGGTGGCAACCGTGCTCGCCCTGAACGTGCCCGTTGGCCGGTCCTTCTGCGTTCGCTCCGCTATGCCCCCCACCTTCACAGTTAGATTTCCGGTTCTGACCGGCGCCCATGTGGGTAc                >  4:68201/1‑274 (MQ=255)

                                                                                                                                                                                                                                       tAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATTACCGCCGCCGCAACCGTCGTGGTCGTGATCGTGGTGGTGATCGTGCGCGCCGTGAACGTGACCATGAGCCAGTTCTTCTTCATTCGCTTCCCGAATCGCCACAACTTCACCGTTTAATTTTGGTTTTTTACCGGCCCGCATGTTGTTAcc               >  4:541205/1‑275 (MQ=255)

                                                                                                                                                                                                                                       tAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGTGGTGATCGTGCGCGCCGTGAACGTGACCATGAGCCAGTTCTTCTTCAGTCGCTTCGCGAATCGCCACAACTTCAACGTTGAATTTCAGGTTCTGACCGGCCAGCATGTGGTTAcc               >  2:139302/1‑275 (MQ=255)

                                                                                                                                                                                                                                       tAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGTGGTGATCGTGCGCGCCGTGAACGTGACCATGAGCCAGTTCTTCTTCAGTCGCTTCGCGAATCGCCACAACTTCAACGTTGAATTTCAGGTTCTGACCGGCCAGCATGTGGTTAcc               >  4:420471/1‑275 (MQ=255)

                                                                                                                                                                                                                                        aaaaaaaGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGTGGTGATCGTGCGCGCCGTGAACGTGACCATGAGcc                                                                                               <  2:122256/194‑1 (MQ=255)

                                                                                                                                                                                                                                         aaaaaaGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCgtcgtggtcg                                                                                                                                            <  5:898203/148‑1 (MQ=255)

                                                                                                                                                                                                                                         aaaaaaGCGGGGAATCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCgtcgtggtcgt                                                                                                                                           >  7:468556/1‑149 (MQ=255)

                                                                                                                                                                                                                                          aaaaaGCGGGGAATCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCgtcgtggtcgt                                                                                                                                           >  5:419422/1‑148 (MQ=255)

                                                                                                                                                                                                                                             aaGCGGGGAATCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGAt                                                                                                                                        >  6:381363/1‑148 (MQ=255)

                                                                                                                                                                                                                                               gCGGGGATTCCCAGCTGTTTTGTCACTTTTTCGGTATTAGTGGCAACTGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGTGGTGATCGTGCGCGCCGTGAACGTGACCATGAGCCAGTTCTTCTTCAGTCGCTTCGCGAATCGCCACAACTTCAACGTTGAATTTCAGGTTCTGACCGGCCAGCATGTGGTTACCATCAACCa       <  2:596669/275‑1 (MQ=255)

                                                                                                                                                                                                                                                cGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCgtg                                                                                                                                    <  8:750485/149‑1 (MQ=255)

                                                                                                                                                                                                                                                   ggATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCgtg                                                                                                                                    >  7:725059/1‑146 (MQ=255)

                                                                                                                                                                                                                                                   ggATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCgtg                                                                                                                                    <  8:725059/146‑1 (MQ=255)

                                                                                                                                                                                                                                                    gATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCgtggtgg                                                                                                                                >  8:127215/1‑149 (MQ=255)

                                                                                                                                                                                                                                                    gATTCCCCGCTTTTTTGTCACTTTTTCGGTATGTGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGTGGTGATCGTGCGCGCCGTGAACGTGACCATGAGCCAGTTCTTCTTCAGTCGCTTCGCGAATCGCCACAACTTCAACGTTGAATTTCAGGTTCTGACCGGCCAGCATGTGGTTACCATCAACCacgacg  <  2:433250/275‑1 (MQ=255)

                                                                                                                                                                                                                                                     aTTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCgtggtggt                                                                                                                               >  7:551574/1‑149 (MQ=255)

GGTGAATCATGCAGGATTTATCATTGGAAGCACGCCTGGCAGAGCTGGAGAGCCGACTGGCTTTTCAGGAAATCACCATTGAAGAACTGAACGTCACGGTGACCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGTGGTGATCGTGCGCGCCGTGAACGTGACCATGAGCCAGTTCTTCTTCAGTCGCTTCGCGAATCGCCACAACTTCAACGTTGAATTTCAGGTTCTGACCGGCCAGCATGTGGTTACCATCAACCACGACG  >  NC_000913/3477632‑3478150

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