BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,477,878	T→A	5.9%	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.9%	622.0 / 20.7	185	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 (92/81); new base A (4/7); total (96/89)
Fisher's exact test for biased strand distribution p-value = 3.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.86e-01

Only 100 of 189 total aligned reads displayed.

AATCCGCGGCGACCATTCGACAAAAGAGGTGAATCATGCAGGATTTATCATTGGAAGCACGCCTGGCAGAGCTGGAGAGCCGACTGGCTTTTCAGGAAATCACCATTGAAGAAC‑TGAACGTCACGGTGACCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGC‑CTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGTGGTGATCGTGCGCGCCGTGAACGTGACCATGAGCCAGTTCTTCTTCAGTCGCTTCGCGAATCGCCACAACTTCAACGTTGAATTTCAGGTTCTGACCGGCCAGCATGTGGTTACCATCAACCACGACG  >  NC_000913/3477605‑3478150

                                                   tGGAAGCACGCCTGGCAGAGCTGGAGAGCCGACTGGCTTTTCAGGAAATCACCATTGAAGAAC‑TGAACGTCACGGTGACCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCCCTTTTTCGGTATTAGTGGCAACCGCAACCGCCGtt                                                                                                                                                                                                                               >  2:324501/1‑275 (MQ=255)

                                                     gAAGCCGGCCTGGCAGAGCTGGAGAGCCGCCTGGCTTTTCAGGAAATCACCATTGAAGAAC‑TGAACGTCACGGTGACCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCGGACCGAGAAGTTAAAAGCCAGCCCGCCGTCGACCATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCCGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCttt                                                                                                                                                                                                                         <  1:213778/279‑1 (MQ=255)

                                                          aCGCCTGGCAGAGCTGGAGAGCCGACTGGCTTTTCAGGAAATCACCATTGAAGAAC‑TGAACGTCACGGTGACCGCTCATGAAATGGAGATGGCGACACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACAACGCGTCGCAGGCCGAAGAAACGCCACACCCACATTCTTGAGCCTTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCCTTTTTTCGGATTTAGGTGCAACCGCAACCGCCTTTGCCcttg                                                                                                                                                                                                                        >  4:311514/1‑274 (MQ=255)

                                                                        tggaggggcggatGGCTTTGAAGGAAATAACCATTGAAGAAA‑TGAACGTCACGGTGACCGCGCATGAAATGGAGAGGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTGAAAAGCCAGCCAGCAGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACagcagc                                                                                                                                                                                                           <  2:538255/262‑1 (MQ=255)

                                                                                  cttgctttttAAGGAAATAACC‑TTGAAAAACTTAAACGCCACGGTGGCCGTCCTGGAAATGGGGATGGCGAACCGTGGCGATACTCGGCGTCTGCTGACCGGGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAACCGCCCCCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGC‑CTTCGCCAcc                                                                                                                                                                                                <  4:120556/268‑1 (MQ=255)

                                                                                   cgggCTTTTAAGGAAATAACCATTGAAGAAC‑TGAACGTCACGGTGACCGCTCATGAAATGGGGATGGCGAAACGGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAGGCCAGCCAGCCGGCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGC‑CTTCGCCAcc                                                                                                                                                                                                <  2:445340/271‑1 (MQ=255)

                                                                                    tGGCTTTTCAGGAAATCACCATTGAAGAAC‑TGAACGTCACGGTGACCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCGGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGC‑CTTCGCCAcc                                                                                                                                                                                                <  2:425745/272‑1 (MQ=255)

                                                                                      gCTTTTCAGGAAAGAACCATTGAAGAAC‑TGAACGTCACGGTGACCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGC‑CTTCGCCACCGTGTTCAt                                                                                                                                                                                        <  1:273506/278‑1 (MQ=255)

                                                                                                                         tCACGGTGACCGCTCATGAAAGGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGC‑CTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCgcagca                                                                                                                                                         <  2:151478/275‑1 (MQ=255)

                                                                                                                          cccgttgacttcTCATGAAATGGACATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACAGCGACGGTAAAAGCCAGCCAGCCGTCGACCATCGCGTCGTAGGGTGAAGGAACGCCACCGCCACATTATTGACGCCGAACCAGAGCGGGGATCGCCCCCGTTCTTTGCACTGTTTCGGTATCACTTGGAACAGCCACAGCCTGTGCCCTTACCGACCCAACTGCACTTC‑CCACCCTTCTCATGACCGTTTTCGTCGCCATGACCCCCccaccac                                                                                                                                                        >  4:241578/12‑274 (MQ=255)

                                                                                                                               tGACCGCTAATGAAAGGGAGATGGCGAAACTGCGCGATCATCTGCGTATGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATGGCGTGGCAGGCTGAAGAAACGCCACCGCCACATTATGGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGC‑CTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACcg                                                                                                                                                     <  2:672270/273‑1 (MQ=255)

                                                                                                                                   cGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAATGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGAATcccc                                                                                                                                                                                                                                                                              >  7:1441774/1‑149 (MQ=255)

                                                                                                                                   cGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTcccc                                                                                                                                                                                                                                                                              <  5:213700/149‑1 (MQ=255)

                                                                                                                                   cGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTcccc                                                                                                                                                                                                                                                                              >  5:1447179/1‑149 (MQ=255)

                                                                                                                                     cTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGC‑CTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGCCGTGGTCGTGa                                                                                                                                         >  1:14840/1‑279 (MQ=255)

                                                                                                                                      tCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCt                                                                                                                                                                                                                                                                           <  8:636034/149‑1 (MQ=255)

                                                                                                                                      tCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCt                                                                                                                                                                                                                                                                           >  7:1602617/1‑149 (MQ=255)

                                                                                                                                       cATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCtt                                                                                                                                                                                                                                                                          <  5:1072176/149‑1 (MQ=255)

                                                                                                                                       cATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCtt                                                                                                                                                                                                                                                                          <  8:1687314/149‑1 (MQ=255)

                                                                                                                                       cATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGAATCCCCGCtt                                                                                                                                                                                                                                                                          <  5:842729/149‑1 (MQ=255)

                                                                                                                                       cATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAATTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTCAAGAAGCGCCACCGCCACCTTATTGAGGCGTAAAAAAAGCGGGGAATCCCCGCtt                                                                                                                                                                                                                                                                          <  7:1002935/149‑1 (MQ=255)

                                                                                                                                        aTGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCCTAAAAAAAGCGGGGAATCCCCGCttt                                                                                                                                                                                                                                                                         <  8:1441774/149‑1 (MQ=255)

                                                                                                                                          gAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCttttt                                                                                                                                                                                                                                                                       >  5:1375130/1‑149 (MQ=255)

                                                                                                                                          gAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCAAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCATAAAAAAAGCGGGGATTCCCCGCttttt                                                                                                                                                                                                                                                                       <  7:1510793/149‑1 (MQ=255)

                                                                                                                                             aTGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGt                                                                                                                                                                                                                                                                    >  7:769224/1‑149 (MQ=255)

                                                                                                                                              tGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTc                                                                                                                                                                                                                                                                   <  8:554188/149‑1 (MQ=255)

                                                                                                                                                    tGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACtttt                                                                                                                                                                                                                                                             >  6:1555047/1‑149 (MQ=255)

                                                                                                                                                     ggCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACttttt                                                                                                                                                                                                                                                            <  7:1269214/149‑1 (MQ=255)

                                                                                                                                                     ggCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGC‑CTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCgttgtggtta                                                                                                                             >  2:473035/1‑273 (MQ=255)

                                                                                                                                                      gCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTc                                                                                                                                                                                                                                                           <  8:1602617/149‑1 (MQ=255)

                                                                                                                                                      gCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGC‑CTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGTGGTGATCGTg                                                                                                                        >  1:22874/1‑279 (MQ=255)

                                                                                                                                                       cGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGACCATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTCTCg                                                                                                                                                                                                                                                          <  6:1208335/149‑1 (MQ=255)

                                                                                                                                                        gAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCgg                                                                                                                                                                                                                                                         <  6:1375130/149‑1 (MQ=255)

                                                                                                                                                        gAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCgg                                                                                                                                                                                                                                                         >  6:1089469/1‑149 (MQ=255)

                                                                                                                                                           aCTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTAt                                                                                                                                                                                                                                                      >  7:469830/1‑149 (MQ=255)

                                                                                                                                                             tGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTa                                                                                                                                                                                                                                                    <  5:61856/149‑1 (MQ=255)

                                                                                                                                                               cgcgATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGt                                                                                                                                                                                                                                                  <  6:1447179/149‑1 (MQ=255)

                                                                                                                                                               cgcgATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGC‑CTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGTGGTGATCGTgcgcgc                                                                                                                   >  4:162799/1‑275 (MQ=255)

                                                                                                                                                               cgcgATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGC‑CTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGTGGTGATCGTgcgcgc                                                                                                                   >  2:306107/1‑275 (MQ=255)

                                                                                                                                                                 cgATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTgg                                                                                                                                                                                                                                                <  6:1140475/149‑1 (MQ=255)

                                                                                                                                                                  gATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGgc                                                                                                                                                                                                                                               >  5:573593/1‑149 (MQ=255)

                                                                                                                                                                   atcatcTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCtt                                                                                                                                                                                                                                                                          >  8:874864/1‑121 (MQ=255)

                                                                                                                                                                   atcatcTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCtt                                                                                                                                                                                                                                                                          <  7:874864/121‑1 (MQ=255)

                                                                                                                                                                    tcatcTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGgcaa                                                                                                                                                                                                                                             <  6:897378/149‑1 (MQ=255)

                                                                                                                                                                       tcTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGgcaaccg                                                                                                                                                                                                                                          >  7:1588095/1‑149 (MQ=255)

                                                                                                                                                                       tcTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGgcaaccg                                                                                                                                                                                                                                          >  7:508639/1‑149 (MQ=255)

                                                                                                                                                                       tcTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGgcaaccg                                                                                                                                                                                                                                          >  5:1222325/1‑149 (MQ=255)

                                                                                                                                                                              ctgctgACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGC‑CTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTTGTGGTGATCGTGCGCGCCGTGAACGTTACCAt                                                                                                    >  2:573827/1‑275 (MQ=255)

                                                                                                                                                                                gctgACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGC‑CTTCGc                                                                                                                                                                                                    >  2:24999/1‑177 (MQ=255)

                                                                                                                                                                                  tgACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGtt                                                                                                                                                                                                                               >  8:20723/1‑149 (MQ=255)

                                                                                                                                                                                   gACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGAATCCCCGCttttt                                                                                                                                                                                                                                                                       >  7:1357493/1‑108 (MQ=255)

                                                                                                                                                                                    aCCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGc                                                                                                                                                                                                                             <  5:1555047/149‑1 (MQ=255)

                                                                                                                                                                                     ccGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGcc                                                                                                                                                                                                                            <  5:1504564/149‑1 (MQ=255)

                                                                                                                                                                                      cGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGgcaaccg                                                                                                                                                                                                                                          >  6:758133/1‑134 (MQ=255)

                                                                                                                                                                                         gaAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTa                                                                                                                                                                                                                        >  5:68039/1‑149 (MQ=255)

                                                                                                                                                                                             ttAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGc                                                                                                                                                                                                                    >  7:773370/1‑149 (MQ=255)

                                                                                                                                                                                             ttAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGc                                                                                                                                                                                                                    >  6:573105/1‑149 (MQ=255)

                                                                                                                                                                                                 aaGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCTCCGTTGCCTTTACCGCcaca                                                                                                                                                                                                                >  7:981813/1‑149 (MQ=255)

                                                                                                                                                                                                     cagccagccGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCCCCAcagcag                                                                                                                                                                                                            >  8:1245831/1‑149 (MQ=255)

                                                                                                                                                                                                       gccagccGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGc‑c                                                                                                                                                                                                         <  7:736080/149‑1 (MQ=255)

                                                                                                                                                                                                         cagccGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCtttttt                                                                                                                                                                                                                                                                      >  6:357870/1‑87 (MQ=255)

                                                                                                                                                                                                         cagccGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCtttttt                                                                                                                                                                                                                                                                      <  5:357870/87‑1 (MQ=255)

                                                                                                                                                                                                                 gAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGC‑CTTCGCCACCg                                                                                                                                                                                               >  5:1203419/1‑149 (MQ=255)

                                                                                                                                                                                                                  aaaacggcgTCGAGGGCTGAAGAACCGCCACCCCCACTTTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACCGCAGC‑CTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGTGGTGATCGTGCGCGCCGTGAACGTGACCATGAGCCAGTTCTTCTTCAGTCGCTTCGCGAATCGc                                                                  <  2:174732/267‑1 (MQ=255)

                                                                                                                                                                                                                  aaCATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGC‑CTTCGc                                                                                                                                                                                                    >  4:138732/1‑143 (MQ=255)

                                                                                                                                                                                                                      tCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGC‑CTTCGCCAcc                                                                                                                                                                                                <  2:138734/143‑1 (MQ=255)

                                                                                                                                                                                                                      tCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGC‑CTTCGCCACCGTGTTc                                                                                                                                                                                          <  6:68039/149‑1 (MQ=255)

                                                                                                                                                                                                                         cgTCGCAGGCTGAAGAAACGCCACCGCCACATTATT‑AGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGC‑CTTCGCCACCGTGTTCATGa                                                                                                                                                                                      <  8:940024/149‑1 (MQ=255)

                                                                                                                                                                                                                          gTCGGAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGC‑CTTCGCCACCGTGTTCATGa                                                                                                                                                                                      <  8:773370/149‑1 (MQ=255)

                                                                                                                                                                                                                                 gCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGC‑CTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGTGGTGATCGTGCGCGTCGGGAACGTGACCATGAGCCAGTGTTTCTTCAGTCGCTGCGCGAGTTGCCACAAGTTCAACGGTGACtt                                              >  1:100084/1‑278 (MQ=255)

                                                                                                                                                                                                                                  cTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGC‑CTTCGCCCCCGTGTTCATGACCGTGATc                                                                                                                                                                              <  6:1514007/149‑1 (MQ=255)

                                                                                                                                                                                                                                   tGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGgc                                                                                                                                                                                                                                               <  8:231572/84‑1 (MQ=255)

                                                                                                                                                                                                                                   tGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGgc                                                                                                                                                                                                                                               >  7:231572/1‑84 (MQ=255)

                                                                                                                                                                                                                                       gaaACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTGTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGC‑TTTCGCCACCGTGTTCATGACCGTGATCGTGGc                                                                                                                                                                         >  7:1419904/1‑149 (MQ=255)

                                                                                                                                                                                                                                                       attattGAGGCGTAAAAAAAGCTGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCCCAGCAGC‑CTTCGCCACCGTGTTCATGACCGTGACCGTGGCCATGACCGCCgcagca                                                                                                                                                         >  8:149582/1‑149 (MQ=255)

                                                                                                                                                                                                                                                        ttattGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGACAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGC‑CTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCaa                                                                                                                                                        >  6:737548/1‑149 (MQ=255)

                                                                                                                                                                                                                                                              aGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGC‑CTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACcgtcg                                                                                                                                                  >  6:251757/1‑149 (MQ=255)

                                                                                                                                                                                                                                                              aGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGC‑CTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGTGGTGATCg                                                                                                                          >  2:52154/1‑173 (MQ=255)

                                                                                                                                                                                                                                                                gCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGC‑CTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCgtcgtggtcgt                                                                                                                                           >  2:362993/1‑154 (MQ=255)

                                                                                                                                                                                                                                                                gCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGC‑CTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGTGGTGATCGTGCGCGCCGTGAACGTGACCATGAGCCAGTTCTTCTTCAGTCGCTTTGCGAATCGCCACAACTTCGACGTTTAATTTCAGGTTCTGACCGGCTAGCATGTGGt                   >  2:221772/1‑274 (MQ=255)

                                                                                                                                                                                                                                                                  gTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACagcagc                                                                                                                                                                                                           <  7:34128/89‑1 (MQ=255)

                                                                                                                                                                                                                                                                  gTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACagcagc                                                                                                                                                                                                           >  8:34128/1‑89 (MQ=255)

                                                                                                                                                                                                                                                                  gTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGC‑CTTCGCCACCgt                                                                                                                                                                                              <  4:734727/101‑1 (MQ=255)

                                                                                                                                                                                                                                                                  gTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGC‑CTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCgtcgtggt                                                                                                                                              >  7:437405/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                  gTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGC‑CTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCgtcgtggt                                                                                                                                              >  7:644868/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                   tAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGC‑CTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGTGGTGATCGTgcgcg                                                                                                                    >  3:77442/1‑174 (MQ=255)

                                                                                                                                                                                                                                                                   tAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGC‑CTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGTGGTGATCGTGCGCGCCGTGAACGTGACCATGAGCCAGTTCTTCTTCAGTCGCTTCGCGAATCGCCACAACTTCAACGTTGAATTTCAGGTTCTGACCGGCCAGCATGTGGTTAcc               >  2:620860/1‑275 (MQ=255)

                                                                                                                                                                                                                                                                   tAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGC‑CTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGTGGTGATCGTGCGCGCCGTGAACGTGACCATGAGCCAGTTCTTCTTCAGTCGCTTCGCGAATCGCCACAACTTCAACGTTGAATTTCAGGTTCTGACCGGCCAGCATGTGGTTAcc               >  2:550617/1‑275 (MQ=255)

                                                                                                                                                                                                                                                                    aaaaaaaGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGC‑CTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATc                                                                                                                                       <  4:362989/154‑1 (MQ=255)

                                                                                                                                                                                                                                                                     aaaaaaGCGGGGAATCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGC‑CTTCGCCACCGTGTTCATGACCGTGATc                                                                                                                                                                              >  5:1264276/1‑114 (MQ=255)

                                                                                                                                                                                                                                                                       aaaaGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGC‑CTTCGCCACCGTGTTCATGACCGTGATc                                                                                                                                                                              <  2:754564/112‑1 (MQ=255)

                                                                                                                                                                                                                                                                          aGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGC‑CTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCg                                                                                                                                      >  6:102215/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                           gCGGGGTTCCCCCGCTTTTTTGTCACTTTTTCGTTATTATTGGCAACCGCAACCGCCTGTGCCTTTACCGCCACAGCAGC‑CTTCGCCACCGTTTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGTGGTGATCGTGCGCGCCGTGAACGTGACCATGAGCCAGTTCTTCTTCAGTCGCTTCGCGAATCGCCACAACTTCAACGTTGAATTTCAGGTTCTGACCGGCCAGCATGTGGTTACCATCAACCa       <  4:550613/275‑1 (MQ=255)

                                                                                                                                                                                                                                                                           gCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGC‑CTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCgt                                                                                                                                     >  6:1498254/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                             ggggATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGC‑CTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAAcc                                                                                                                                                      >  4:281626/1‑130 (MQ=255)

                                                                                                                                                                                                                                                                                gtttCCCCGTTGTTTTGTAAGGTTTTCGGTATAGGGGGAAAACCAACCCGCGTTTGCGTTTCCCGCCACTGCAGC‑CTTCGCCACCGTGTGCATGACCGTGATGTTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGGTCGGGGTGGTGATCGGGCGCGCCGTGAACGTGACCATGAGCCAGTTCTTCTTCAGTCGCTTCGCGAATCGCCACAACTTCAACGTTGAATTTCAGGTTCTGACCGGCCAGCATGTGGTTACCATCAACCacgacg  <  4:652449/273‑1 (MQ=255)

                                                                                                                                                                                                                                                                                 aTTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGC‑CTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACcgtcg                                                                                                                                                  <  2:281629/130‑1 (MQ=255)

                                                                                                                                                                                                                                                                                 aTTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGC‑CTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCgtggtggt                                                                                                                               <  8:482351/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                                 aTTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGC‑CTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCgtggtggt                                                                                                                               <  8:744047/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                                  ttCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGACAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGC‑CTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCgtggtggtg                                                                                                                              <  5:737548/149‑1 (MQ=255)

AATCCGCGGCGACCATTCGACAAAAGAGGTGAATCATGCAGGATTTATCATTGGAAGCACGCCTGGCAGAGCTGGAGAGCCGACTGGCTTTTCAGGAAATCACCATTGAAGAAC‑TGAACGTCACGGTGACCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGC‑CTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGTGGTGATCGTGCGCGCCGTGAACGTGACCATGAGCCAGTTCTTCTTCAGTCGCTTCGCGAATCGCCACAACTTCAACGTTGAATTTCAGGTTCTGACCGGCCAGCATGTGGTTACCATCAACCACGACG  >  NC_000913/3477605‑3478150

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