BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,426,592	T→A	6.1%	intergenic (+35/+36)	rplI → / ← yjfZ	50S ribosomal subunit protein L9/uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,426,592	0	T	A	6.1%	812.8 / 27.6	232	intergenic (+35/+36)	rplI/yjfZ	50S ribosomal subunit protein L9/uncharacterized protein
Reads supporting (aligned to +/- strand): ref base T (104/112); new base A (5/9); total (111/121)
Fisher's exact test for biased strand distribution p-value = 4.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.82e-01

Only 100 of 241 total aligned reads displayed.

TGTTACCATCGCGTCTAAAGCTGGCGACGAAGGTAAACTGTTCGGTTCCATCGGTACTCGCGACATCGCTGACGCTGTAACTGCAGCTGGCGTTGAAGTGGCTAAGAGCGAAGTTCGTCTGCCGAACGGCGTTCTGCGTACCACTGGCGAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCGCTCCTGCGGCACTCGCTCGAAACCCAAACCCGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAA  >  NC_000913/4426338‑4426865

     ccATCGCGTCTAAAGCTGGCGACGAAGGTAAACTGTTCGGTTCCATCGGTACTCGCGACATCGCTGACGCTGTAACTGCAGCTGGCGTTGAAGTGGCTAAGAGCGAAGTTCGTCTGCCGAACGGCGTTCTGCGTACCACTGGCGAACACGAAGTGCGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGGTATTCAACGAGACGTAAAAACGCCGACCATTGGACGGAGTTTTGGTTTTTAtttt                                                                                                                                                                                                                                                          >  4:360141/1‑275 (MQ=255)

         cgcgTCTAAAGCTGGCGACGAAGGTAAACTGTTCGGTTCCATCGGTACTCGCGACATCGCTGACGCTGTAACTGCAGCTGGCGTTGAAGTGGCTAAGAGCGAAGTTCGTCTGCCGAACGGCGTTCTGCGTACCACTGGNGAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTACACGTAGTAGCTGAATAATTNGTTATTCAACGCGACGTAGAAACGACGACCNNTGTGCGGGGTTTGTCTTTCTAgttttttt                                                                                                                                                                                                                                                      >  3:79332/1‑272 (MQ=255)

            gTCTAAAGCTGGCGACGAAGGTAAACTGTTCGGTTCCATCGGTACTCGCGACATCGCTGACGCTGTAACTGCAGCTGGCGTTGAAGTGGCTCAGCGCGAAGTTCGTCTGCAGAACGGCGTTCTGCGTACCACTTGCGANCACGAAGTGAGCCTCCAGGTTCACAGCGCAGTATTCTCGCAAGTACTGGTAAACGTAGTTGCTGAATCATTAGTTATTCANTGAGCAGGAAAAACGACTAACATACNTCGTGGATTTTTTTTCTttttttgtgga                                                                                                                                                                                                                                                    >  3:140577/1‑269 (MQ=255)

               tAAAGCTGGCGACGAAGGTAAACTGTTCGGTTCCATCGGTACTCGCGACATCGCTGACGCTGTAACTGCAGCTGGCGTTGAAGTGGCTAAGAGCGAAGTTCGTCTGCCGAACGGCGTTCTGCGTACCACTGGCGAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTa                                                                                                                                                                                                                                                <  2:640367/275‑1 (MQ=255)

                                       gTTCGGTTCCAGCGGTACTCTCGCAATCGCTGACGCTGTAACTGCAGCTTGGGTTTAAGTGGCTAAGAGCGAAGTTCGGCTGCCGAACGGCGTTCTGCGTACCACTGGCGAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCAtt                                                                                                                                                                                                                         <  4:285256/274‑1 (MQ=255)

                                                cATCGGTACTCGCGACATCGCTGACGCTGTAACTGCAGCTGGCGTTGAAGTGGCTAAGAGCGAAGTTCGTCTGCCGAACGGCGTTCTGCGTACCACTGGCGAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATtct                                                                                                                                                                                                                >  2:174016/1‑274 (MQ=255)

                                                  tCGGTACTCGCGACATCGCTGACGCTGTAACTGCAGCTGGCGTTGAAGTGGCTAAGAGCGAAGTTCGGCTGCCGAACGGCGTTCTGCGTACCACTGGCGAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCtt                                                                                                                                                                                                             <  2:379947/275‑1 (MQ=255)

                                                                aTCGCTGACGCTGTAACTGCAGCTGGCGTTGAAGTGGCTAAGAGCGAAGTTCGTCTGCCGAACGGCGTTCTGCGTACCACTGGCGAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGGCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTATATTCGCAAGTAGATTCATTCACCATTCTCTTGGTTCCAGAAGTCa                                                                                                                                                                                               >  4:628228/1‑275 (MQ=255)

                                                                                cgggaGCGGGCGTTGACGTGGGTAAGGGGGAAGTTCAGCTACCGAACGGCGGTCTGCGTACCACTGGCGAACACGAAGTGAGCTTCCAGGTTCACAGAGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTAGTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCg                                                                                                                                                                               <  4:464450/271‑1 (MQ=255)

                                                                                   cAGCTGGCGTTGAAGTGGCTAAGAGCGAAGTTCGTCTGCCGAACGGCGTTCTGCGTACCACTGGCGAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGtt                                                                                                                                                                             <  2:410793/274‑1 (MQ=255)

                                                                                               aaGTGGTGAAGAGGGAAGTTCGTCGGCCGAACGGCGTTTGGCGTACCACTGGCGAAAACGAAGTGAGCTTCCAGGTTAACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAg                                                                                                                                                                <  4:8809/275‑1 (MQ=255)

                                                                                                          agCGAAGTTCGTCTGCCGAACGGCGTTCTGCGTACCACTGGCGAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCAt                                                                                                                                                                                                                                                                                   >  7:295475/1‑149 (MQ=255)

                                                                                                           gCGAAGTTCGTCTGCCGAACGGCGTTCTGCGTACCACTGGCGAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCAtt                                                                                                                                                                                                                                                                                  >  8:523206/1‑149 (MQ=255)

                                                                                                              aaGTTCGTCTGCCGAACGGCGTTCTGCGTACCACTGGCGAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGtct                                                                                                                                                  <  2:600009/274‑1 (MQ=255)

                                                                                                                gTTCGTCTGCCGAACGGCGTTCTGCGTACCACTGGCGAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTc                                                                                                                                                                                                                                                                              <  6:769549/148‑1 (MQ=255)

                                                                                                                 ttCGTCTGCCGAACGGCGTTCTGCGTACCACTGGCGAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTTGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCgg                                                                                                                                                                                                                                                                            >  5:110753/1‑149 (MQ=255)

                                                                                                                   cGTCTGCCGAACGGCGTTCTGCGTACCACTGGCGAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCAATGGTCGGCg                                                                                                                                                                                                                                                                          <  7:917961/149‑1 (MQ=255)

                                                                                                                        gCCGAACGGCGTTCTGCGTACCACTGGCGAACACGAAGTGAGCTTACAGGTTCACTGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCAATGGTCGGCGttttt                                                                                                                                                                                                                                                                     <  7:542928/149‑2 (MQ=255)

                                                                                                                            aaCGGCGTTCTGCGTACAACTGGCGAAAACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCGCTCCTGCGGCAct                                                                                                                                   <  4:209087/275‑1 (MQ=255)

                                                                                                                             aCGGCGTTCTGCGTACCACTGGCGAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCGCTCCTGCGGCAct                                                                                                                                   <  4:174016/274‑1 (MQ=255)

                                                                                                                                   ttCTGCGTACCACTGGCGAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCa                                                                                                                                                                                                                                     >  4:244883/1‑170 (MQ=255)

                                                                                                                                     cTGCGTACCACTGGCGAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTc                                                                                                                                                                                                                                                        <  6:581291/149‑1 (MQ=255)

                                                                                                                                       gCGTACCACTGGCGAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTa                                                                                                                                                                                                                                 <  2:244883/170‑1 (MQ=255)

                                                                                                                                         gTACCACTGGCGAACACGAAGTGAGCTTTCAGGTTCACAGCGAAGGATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGTTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCGCTCCTGCGGCACTCGCTCGaaaccca                                                                                                                      <  2:427579/275‑1 (MQ=255)

                                                                                                                                         gTACCACTGGCGAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCGCTCCTGCGGCACTCGCTCGaaaccca                                                                                                                      <  2:634537/275‑1 (MQ=255)

                                                                                                                                          tACCACTGGCGAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCGCTCCTGCGGCACTCGCTCGaaaccca                                                                                                                      <  2:356838/274‑1 (MQ=255)

                                                                                                                                             cACTGGCGAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTa                                                                                                                                                                                                                                                >  6:363342/1‑149 (MQ=255)

                                                                                                                                                  gCGAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTTTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATAAATACGAGCGAAACTGATGTCTCGCTCCTGCGGCACTCCCTCGAAACCCAAACCCGAGt                                                                                                             >  2:560171/1‑275 (MQ=255)

                                                                                                                                                   cGAACACGAAGTGAGCTTCCAGGTTCACAGCCAAGTATTCGCGAAAGTGATCGTACACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGTGTTTTGCTTTCTATTTTTCGTCAGGTAGTAGtt                                                                                                                                                                                                                                          <  5:315764/149‑1 (MQ=255)

                                                                                                                                                    gAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTTGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGttt                                                                                                                                                                                                                                         >  8:205889/1‑149 (MQ=255)

                                                                                                                                                    gAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGttt                                                                                                                                                                                                                                         >  5:566856/1‑149 (MQ=255)

                                                                                                                                                    gAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGttt                                                                                                                                                                                                                                         <  7:323872/149‑1 (MQ=255)

                                                                                                                                                       cacGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGt                                                                                                                                                                                                                                  >  1:79012/1‑153 (MQ=255)

                                                                                                                                                       cacGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCGCTCCTGCGGCACTCGCTCGaaacccaa                                                                                                                     >  4:56548/1‑262 (MQ=255)

                                                                                                                                                        acGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGACGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCa                                                                                                                                                                                                                                     >  5:852683/1‑149 (MQ=255)

                                                                                                                                                         cGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCaa                                                                                                                                                                                                                                    <  8:153531/149‑1 (MQ=255)

                                                                                                                                                          gatgTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCACTGGTCGGCGTTTTGCTTTCTAttt                                                                                                                                                                                                                                                           >  8:355233/4‑125 (MQ=255)

                                                                                                                                                              tGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGa                                                                                                                                                                                                                               >  5:761294/1‑149 (MQ=255)

                                                                                                                                                                aGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGtt                                                                                                                                                                             <  4:384506/197‑1 (MQ=255)

                                                                                                                                                                       aGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCAc                                                                                                                                                                                                                      >  6:750199/1‑149 (MQ=255)

                                                                                                                                                                          ttCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCAt                                                                                                                                                                                                                   >  6:676939/1‑149 (MQ=255)

                                                                                                                                                                               aGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTctct                                                                                                                                                                                                              <  7:515226/149‑1 (MQ=255)

                                                                                                                                                                               aGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTctct                                                                                                                                                                                                              >  7:752763/1‑149 (MQ=255)

                                                                                                                                                                                      tATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTctct                                                                                                                                                                                                              <  7:175301/142‑1 (MQ=255)

                                                                                                                                                                                      tATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTcc                                                                                                                                                                                                       >  5:672123/1‑149 (MQ=255)

                                                                                                                                                                                      tATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTcc                                                                                                                                                                                                       <  7:832057/149‑1 (MQ=255)

                                                                                                                                                                                        ttCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCAt                                                                                                                                                                                                     >  5:772472/1‑149 (MQ=255)

                                                                                                                                                                                          cgcgAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCTACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATaa                                                                                                                                                                                                   <  8:298596/149‑1 (MQ=255)

                                                                                                                                                                                          cgcgAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATaa                                                                                                                                                                                                   <  8:959439/149‑1 (MQ=255)

                                                                                                                                                                                           gcgAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCGCTCCTGCGGCActcg                                                                                                                                 >  2:498986/1‑214 (MQ=255)

                                                                                                                                                                                               aaGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAg                                                                                                                                                                                              <  8:408693/149‑1 (MQ=255)

                                                                                                                                                                                               aaGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCGCTCCTGCGGCActcgctcg                                                                                                                             <  4:498986/214‑1 (MQ=255)

                                                                                                                                                                                                aGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCa                                                                                                                                                                                                                           >  1:106351/1‑119 (MQ=255)

                                                                                                                                                                                                 gTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGa                                                                                                                                                                                            <  7:140450/149‑1 (MQ=255)

                                                                                                                                                                                                  ttaTCGAAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAAAGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCGCTCCTGCGGCACTCGCTCGAAACCCAAACCCGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCa                                                              <  1:167612/272‑1 (MQ=255)

                                                                                                                                                                                                  tGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGAt                                                                                                                                                                                           <  6:649376/149‑1 (MQ=255)

                                                                                                                                                                                                  tGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCGCTCCTGCGGCACTCGCTCGAAACCCAAACCCGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCa                                                              <  1:99065/274‑1 (MQ=255)

                                                                                                                                                                                                     tCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGt                                                                                                                                                                                        <  6:163623/149‑1 (MQ=255)

                                                                                                                                                                                                      cGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGtt                                                                                                                                                                                       >  5:831294/1‑149 (MQ=255)

                                                                                                                                                                                                       gTAAACGTTTTTGCTGAATAATTCGTTATTCAACGAGTCGTAAAAACTCCATGCATTTGTCGGCGTTTTGCTTTCTATTTTTTGTCAAGTATTAGTTTCGCATGTAAATCCATTCACCATTCTCTTGTGTCCATAAATCCGGATGGTTTGCTGTCGTTAGTTACATACGAGCGCAACTGATGTCTCGCTCCTACGGCACTCGCTCGAAACCCAAACCCGAGTAGATTGCGGCTGATCTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGt                                                        <  4:428192/275‑1 (MQ=255)

                                                                                                                                                                                                       gTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGa                                                                                                                                                                                            >  5:446341/1‑143 (MQ=255)

                                                                                                                                                                                                       gTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGa                                                                                                                                                                                            <  6:446341/143‑1 (MQ=255)

                                                                                                                                                                                                        tAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTg                                                                                                                                                                                     >  7:680655/1‑149 (MQ=255)

                                                                                                                                                                                                          aaCGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCt                                                                                                                                                                                   >  8:23147/1‑149 (MQ=255)

                                                                                                                                                                                                           aCGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCCTTCTCtt                                                                                                                                                                                                             <  6:964109/122‑1 (MQ=255)

                                                                                                                                                                                                           aCGTAGGAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCATGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCGCTCCTGCGGCACTCGCTCGAAACCCAAACCCGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAAt                                                    <  2:301319/275‑1 (MQ=255)

                                                                                                                                                                                                             gtagtaGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCAt                                                                                                                                                                                                     <  7:21732/128‑1 (MQ=255)

                                                                                                                                                                                                                 tagCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATACGTCAGGATGGTTTGCt                                                                                                                                                                                   >  2:671838/1‑142 (MQ=255)

                                                                                                                                                                                                                   gCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGa                                                                                                                                                                          <  5:658190/149‑1 (MQ=255)

                                                                                                                                                                                                                    cTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCg                                                                                                                                                                                                                                                       <  3:32055/71‑1 (MQ=255)

                                                                                                                                                                                                                    cTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCAc                                                                                                                                                                                                                      >  2:12658/1‑104 (MQ=255)

                                                                                                                                                                                                                     tGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCg                                                                                                                                                                               <  4:671838/142‑1 (MQ=255)

                                                                                                                                                                                                                      gAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTGTGCTTTCTATTTTTCATCAGGTACTAGTTTCGCAAGTAGATCCATTCACCACTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGatac                                                                                                                                                                       <  5:294210/149‑1 (MQ=255)

                                                                                                                                                                                                                       aataatTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGataca                                                                                                                                                                      >  6:180322/1‑149 (MQ=255)

                                                                                                                                                                                                                       aataatTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCGCTCCTGCGGCActcg                                                                                                                                 >  4:735304/1‑186 (MQ=255)

                                                                                                                                                                                                                        ataatTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCAtt                                                                                                                                                                                                                  <  4:12658/104‑1 (MQ=255)

                                                                                                                                                                                                                         taatTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTgg                                                                                                                                                                                 >  3:176202/1‑136 (MQ=255)

                                                                                                                                                                                                                           atTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCtt                                                                                                                                                                                                             <  4:32530/106‑1 (MQ=255)

                                                                                                                                                                                                                           atTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGc                                                                                                                                                                                    <  2:507333/131‑1 (MQ=255)

                                                                                                                                                                                                                           atTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCGCTCCTGCGGCActcgctcg                                                                                                                             <  2:735304/186‑1 (MQ=255)

                                                                                                                                                                                                                            ttCGTTATTCAACGAGACGTAAAAACGCCGACCATCGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGa                                                                                                                                                                 <  8:704585/149‑1 (MQ=255)

                                                                                                                                                                                                                             tCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGtt                                                                                                                                                                             <  1:176202/136‑1 (MQ=255)

                                                                                                                                                                                                                                ttATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGtttt                                                                                                                                                                                                                                                                      >  8:189133/1‑44 (MQ=255)

                                                                                                                                                                                                                                 tATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCa                                                                                                                                                            <  7:352720/149‑1 (MQ=255)

                                                                                                                                                                                                                                 tATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCGCTCCTGCGGCACTCGCTCGAAACCCAAACCCGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTACAACATCCTGGGGTAAATCAGAATACGCGTTTAGCAATg                               >  4:417197/1‑274 (MQ=255)

                                                                                                                                                                                                                                  aTTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCaa                                                                                                                                                           >  7:378989/1‑149 (MQ=255)

                                                                                                                                                                                                                                       aCGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTTAGGTATTAGTTTCGCAAGTAGATCCATTCACCGTTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGAt                                                                                                                                                      <  5:1069568/149‑1 (MQ=255)

                                                                                                                                                                                                                                        cGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATg                                                                                                                                                     >  7:207537/1‑149 (MQ=255)

                                                                                                                                                                                                                                         gagaCGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGt                                                                                                                                                    >  6:708387/1‑149 (MQ=255)

                                                                                                                                                                                                                                           gaCGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGtct                                                                                                                                                  >  8:580969/1‑149 (MQ=255)

                                                                                                                                                                                                                                             cGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCg                                                                                                                                                <  6:272212/149‑1 (MQ=255)

                                                                                                                                                                                                                                              gTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCGc                                                                                                                                               >  5:360097/1‑149 (MQ=255)

                                                                                                                                                                                                                                              gTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCGc                                                                                                                                               <  7:580969/149‑1 (MQ=255)

                                                                                                                                                                                                                                              gTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCGc                                                                                                                                               <  7:713136/149‑1 (MQ=255)

                                                                                                                                                                                                                                                aaaaaCGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCGCTc                                                                                                                                             >  6:494603/1‑149 (MQ=255)

                                                                                                                                                                                                                                                  aaaCGCCGACCAATGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCCGGCGTTAGATACATACGAGCGCAACTGATGTCTCGCTCCt                                                                                                                                           >  5:944575/1‑149 (MQ=255)

                                                                                                                                                                                                                                                   aaCGCCGACCAATGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTc                                                                                                                                                                                                                                        >  8:842440/1‑55 (MQ=255)

                                                                                                                                                                                                                                                      gCCGACCAATGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCGCTCCTGCgg                                                                                                                                       <  5:512063/149‑1 (MQ=255)

                                                                                                                                                                                                                                                              ttGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGa                                                                                                                                                                                            >  6:314249/1‑88 (MQ=255)

                                                                                                                                                                                                                                                              ttGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGa                                                                                                                                                                                            >  6:825844/1‑88 (MQ=255)

TGTTACCATCGCGTCTAAAGCTGGCGACGAAGGTAAACTGTTCGGTTCCATCGGTACTCGCGACATCGCTGACGCTGTAACTGCAGCTGGCGTTGAAGTGGCTAAGAGCGAAGTTCGTCTGCCGAACGGCGTTCTGCGTACCACTGGCGAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCGCTCCTGCGGCACTCGCTCGAAACCCAAACCCGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAA  >  NC_000913/4426338‑4426865

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