BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,888,030	A→T	5.0%	intergenic (‑39/+31)	rnd ← / ← fadD	ribonuclease D/acyl‑CoA synthetase (long‑chain‑fatty‑acid‑‑CoA ligase)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,888,030	0	A	T	5.0%	855.1 / 25.2	242	intergenic (‑39/+31)	rnd/fadD	ribonuclease D/acyl‑CoA synthetase (long‑chain‑fatty‑acid‑‑CoA ligase)
Reads supporting (aligned to +/- strand): ref base A (115/112); new base T (7/5); total (124/118)
Fisher's exact test for biased strand distribution p-value = 7.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.82e-01

Only 100 of 249 total aligned reads displayed.

CGGATCGCGCAGGATCGCTTTCAGCGGTGACCAGTCGGTGATCCCGAGTGGATCGATTAGCGCCAGATGCTCGCCATCGAAAAGTTGAATCAACCCCAGCTGCGGGTAATAAGTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGA‑CTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTACTTTGTATCCCGTGAGCTGACGGCGGCAAAAAGTCACCAGTGACTCTTCGGTAAGCGATGGATCTTTTTTCACTACGAAGATTTTCACCGCTTCACCACTGGAGCCGGAAGGTACGCCAACAGCCGCGAC  >  NC_000913/1887773‑1888300

                gCTTTCAGCGGTGACCAGTCGGTGATCCCGAGTGGATCGATTAGCGCCAGATGCTCGCCATCGAAAAGTTGAATCAACCCCAGCTGCGGGTAATAAGTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGA‑CTGACTTAACGCTc                                                                                                                                                                                                                                                <  2:354794/274‑1 (MQ=255)

                                            cGAGTGGATAGATTAGCGCCAGATGCTCGCCATCGAAAAGTTGAATCAACCCCAGCTGCGGGTAATAAGTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGA‑CTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCt                                                                                                                                                                                                                   <  2:102370/275‑1 (MQ=255)

                                                   aTCGATTAGCGCCAGATGCTCGCCATCGAAAATTTGAATCAACTCCGGCTGCGGGTAATAAGTGCGCGTAAGAACTAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGTACGGCTTCACAAAAAGAAGCCAACGCATCGTCCGTGGTAATCATTTTGTAATTCCATTTTTTTTCTGTGTAGTAGGCTTC‑AAAAAAAACGCCGGATTAACCGGCGTATGACGA‑CTGGATTTACGGCTAGGATTTTGTTTCCATTTTGCCCCACGCTTCGTCaa                                                                                                                                                                                                            >  4:650010/1‑274 (MQ=255)

                                                   aTCGATTAGCGCCAGATGCTCGCCATCGAAAAGTTGAATCAACCCCAGCTGCGGGTAATAAGTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGA‑CTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCAc                                                                                                                                                                                                            >  2:473446/1‑275 (MQ=255)

                                                        ttAGCGCCAGATGCTCGCCATCGAAAAGTTGAATCAACCCCAGCTGCGGGTAATAAGTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCCCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGAACA‑CTGACTTAACGCCCAGGCTTTATTGTCCACTTTGCCGCCCGCTTCGCCACGTAAt                                                                                                                                                                                                       >  4:615936/1‑275 (MQ=255)

                                                                             cGAAAAGTTGAATCAACCCCAGCTGCGGGTAATAAGTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGA‑CTGACTTAACGCTc                                                                                                                                                                                                                                                >  1:174792/1‑213 (MQ=255)

                                                                                        aTCAACCCCAGCTGCGGGTAATAAGTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGA‑CTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGt                                                                                                                                                                   >  1:291100/1‑279 (MQ=255)

                                                                                              cccAGCGGCGGGTAATAAGTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGA‑CTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGt                                                                                                                                                                   <  2:151007/273‑1 (MQ=255)

                                                                                                  gCTGCGGGTAATAAGTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGA‑CTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGAt                                                                                                                                                                         <  4:118586/263‑1 (MQ=255)

                                                                                                             taaGTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGa                                                                                                                                                                                                                                                                                 <  8:1557956/149‑1 (MQ=255)

                                                                                                             taaGTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGa                                                                                                                                                                                                                                                                                 <  8:352505/149‑1 (MQ=255)

                                                                                                             taaGTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGa                                                                                                                                                                                                                                                                                 >  7:550395/1‑149 (MQ=255)

                                                                                                                 tGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTaa                                                                                                                                                                                                                                                                             <  7:467403/149‑1 (MQ=255)

                                                                                                                     cgTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTCATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTga                                                                                                                                                                                                                                                                  >  4:313290/1‑156 (MQ=255)

                                                                                                                     cgTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCgg                                                                                                                                                                                                                                                                         <  6:793511/149‑1 (MQ=255)

                                                                                                                       tACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGA‑CTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCtt                                                                                                                                         <  2:457990/274‑1 (MQ=255)

                                                                                                                        aCGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGTTAATCCGGCGt                                                                                                                                                                                                                                                                      <  6:781391/149‑1 (MQ=255)

                                                                                                                        aCGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGt                                                                                                                                                                                                                                                                      <  5:1320995/149‑1 (MQ=255)

                                                                                                                         cGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGA‑CTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCg                                                                                                                                       <  2:615940/274‑1 (MQ=255)

                                                                                                                           aaCAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGTTAATCCGGCGttt                                                                                                                                                                                                                                                                    <  5:662541/148‑3 (MQ=255)

                                                                                                                           aaCAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTg                                                                                                                                                                                                                                                                   >  7:1311611/1‑149 (MQ=255)

                                                                                                                            aCAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGGAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTga                                                                                                                                                                                                                                                                  <  7:349139/149‑1 (MQ=255)

                                                                                                                              aaaTTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTgacg                                                                                                                                                                                                                                                                >  6:119685/1‑149 (MQ=255)

                                                                                                                                aTTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTgacga‑c                                                                                                                                                                                                                                                             <  5:1106419/149‑1 (MQ=255)

                                                                                                                                 ttCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACga‑ct                                                                                                                                                                                                                                                            >  7:833515/1‑149 (MQ=255)

                                                                                                                                   cAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACga‑ctga                                                                                                                                                                                                                                                          >  7:268778/1‑149 (MQ=255)

                                                                                                                                        tCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGA‑CTGACTTaa                                                                                                                                                                                                                                                     <  6:1295640/149‑1 (MQ=255)

                                                                                                                                           aGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGA‑CTGACTTAACGc                                                                                                                                                                                                                                                  >  7:1393435/1‑149 (MQ=255)

                                                                                                                                                tATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGGC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGA‑CTGACTTACCGCGCATGCTTTATTCTCCCATTTGCCCCCCGCTGCGTCACGTAATTCTCTCCCCCACAGTTTTCCGACGTTCGATGTCGTTCAACCATCAAGACACTCCAACACCTTCGGTACTTTCAATACCGTGAGCTGAc                                                                                                               >  4:101094/1‑275 (MQ=255)

                                                                                                                                                  tCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGA‑CTGACTTAACGCt                                                                                                                                                                                                                                                 >  3:18867/1‑143 (MQ=255)

                                                                                                                                                    gCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTACGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGA‑CTGACTTAACGCTCAGGCttt                                                                                                                                                                                                                                         <  5:502834/149‑1 (MQ=255)

                                                                                                                                                     ccGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGA‑CTGACTTAACGCTCAGGCTTTa                                                                                                                                                                                                                                        <  8:183470/149‑1 (MQ=255)

                                                                                                                                                               acggacggCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGA‑CTGACTTAACGCTCAGGCTTTATTGTCCACtt                                                                                                                                                                                                                              >  8:722907/1‑149 (MQ=255)

                                                                                                                                                                   acggCTTCACAAAAAGAAGCCAGCGCATCGTCCGTAGTAATCATTTGGTAATTCAAATTGTTTTCGCTTTAATGGCCGAC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACAA‑CTGACGTAACGCTCAGGCTTTATTGTCCACTTTGcc                                                                                                                                                                                                                          <  8:3233/149‑1 (MQ=255)

                                                                                                                                                                       cTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAAc                                                                                                                                                                                                                                                                            >  5:971068/1‑96 (MQ=255)

                                                                                                                                                                        ttCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGa                                                                                                                                                                                                                                                                                 <  6:1602273/90‑1 (MQ=255)

                                                                                                                                                                        ttCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGA‑CTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCcgcgcg                                                                                                                                                                                                                     >  5:1002753/1‑149 (MQ=255)

                                                                                                                                                                          cacacaAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGA‑CTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCt                                                                                                                                                                                                                   <  7:1520234/149‑1 (MQ=255)

                                                                                                                                                                            cacaAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGA‑CTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTc                                                                                                                                                                                                                 <  7:1225403/149‑1 (MQ=255)

                                                                                                                                                                             acaAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGA‑CTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCg                                                                                                                                                                                                                >  6:1687341/1‑149 (MQ=255)

                                                                                                                                                                             acaAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGA‑CTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCg                                                                                                                                                                                                                >  6:1631417/1‑149 (MQ=255)

                                                                                                                                                                              caAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGA‑CTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGt                                                                                                                                                                                                               <  8:1139655/149‑1 (MQ=255)

                                                                                                                                                                               aaaGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGA‑CTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTc                                                                                                                                                                                                              >  6:522336/1‑149 (MQ=255)

                                                                                                                                                                                 agaagCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGA‑CTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCAc                                                                                                                                                                                                            >  8:418549/1‑149 (MQ=255)

                                                                                                                                                                                  gaagCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGA‑CTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACg                                                                                                                                                                                                           >  6:1375171/1‑149 (MQ=255)

                                                                                                                                                                                   aagCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGA‑CTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGt                                                                                                                                                                                                          >  7:1259330/1‑149 (MQ=255)

                                                                                                                                                                                     gCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGA‑CTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTaa                                                                                                                                                                                                        >  8:873593/1‑149 (MQ=255)

                                                                                                                                                                                      ccAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGA‑CTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAAt                                                                                                                                                                                                       >  7:990431/1‑149 (MQ=255)

                                                                                                                                                                                       cAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGA‑CTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAAtt                                                                                                                                                                                                      <  5:1281003/149‑1 (MQ=255)

                                                                                                                                                                                            cATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGA‑CTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCtcgt                                                                                                                                                                                                 <  8:1613301/149‑1 (MQ=255)

                                                                                                                                                                                             aTCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGt                                                                                                                                                                                                                                                                      >  5:73661/1‑80 (MQ=255)

                                                                                                                                                                                                   cGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGA‑CTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTACTTTGTATCCCGTGAGCTGACGGCGGCAAAAAGTCACCAGTGACTCTTCGGTAAGCGATg                                                                        <  2:22829/263‑1 (MQ=255)

                                                                                                                                                                                                    gTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGTTAATCCGGCGtt                                                                                                                                                                                                                                                                     >  5:1573054/1‑73 (MQ=39)

                                                                                                                                                                                                       gTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACga‑ct                                                                                                                                                                                                                                                            >  8:1625483/1‑79 (MQ=255)

                                                                                                                                                                                                        tAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGA‑CTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAAttttt                                                                                                                                                                                     <  8:119831/149‑1 (MQ=255)

                                                                                                                                                                                                          aTCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGA‑CTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAAATCCACCAGCTTCGGTACTTTGTATCCCGTGAGCTGACGGCGGCAAAAAGTCACCAGTGACTCTGCGGTCCGCTATGGATCGTTTTTTTATACGACGAt                                                  >  1:16381/1‑278 (MQ=255)

                                                                                                                                                                                                              tttGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGA‑CTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACg                                                                                                                                                                               >  8:689058/1‑149 (MQ=255)

                                                                                                                                                                                                                   tAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGA‑CTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGa                                                                                                                                                                          >  8:165101/1‑149 (MQ=255)

                                                                                                                                                                                                                     aTTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGA‑CTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGAtt                                                                                                                                                                        <  8:1311611/149‑1 (MQ=255)

                                                                                                                                                                                                                      ttCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGA‑CTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAACTTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTACTTTGTATCCCGTGAGCTGACGGCGGCAAAAAGTCACCGGTGACTCTTCGGTCAGCGAGGGATCTTTTTTACATACGAAACTTTTCACCg                                          >  4:422649/1‑274 (MQ=255)

                                                                                                                                                                                                                      ttCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGA‑CTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGAttt                                                                                                                                                                       >  7:936001/1‑149 (MQ=255)

                                                                                                                                                                                                                       tCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAACAACGCCGGATTAACCGGCGTCTGACGA‑CTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTc                                                                                                                                                                      <  8:526291/149‑1 (MQ=255)

                                                                                                                                                                                                                          aaTTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGA‑CTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGt                                                                                                                                                                   >  7:1665934/1‑149 (MQ=255)

                                                                                                                                                                                                                           aTTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGc                                                                                                                                                                                                                                                                        <  6:671788/48‑1 (MQ=255)

                                                                                                                                                                                                                           aTTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGA‑CTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTaa                                                                                                                                                                                                        <  7:586177/111‑1 (MQ=255)

                                                                                                                                                                                                                           aTTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGA‑CTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTa                                                                                                                                                                  <  8:936001/149‑1 (MQ=255)

                                                                                                                                                                                                                              gTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGA‑CTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACt                                                                                                                                                               <  7:1606983/149‑1 (MQ=255)

                                                                                                                                                                                                                               ttttCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGA‑CTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACtc                                                                                                                                                              >  5:416335/1‑149 (MQ=255)

                                                                                                                                                                                                                                 ttCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGA‑CTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACtcat                                                                                                                                                            >  5:812204/1‑149 (MQ=255)

                                                                                                                                                                                                                                  tctcTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGA‑CTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACtcatc                                                                                                                                                           >  5:293685/1‑149 (MQ=255)

                                                                                                                                                                                                                                      tttAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGA‑CTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCAc                                                                                                                                                                                                            >  8:1292004/1‑96 (MQ=255)

                                                                                                                                                                                                                                      tttAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGA‑CTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCg                                                                                                                                                                     >  5:1146169/1‑135 (MQ=255)

                                                                                                                                                                                                                                      tttAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGA‑CTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGa                                                                                                                                                       >  6:230095/1‑149 (MQ=255)

                                                                                                                                                                                                                                        tAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGA‑CTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGaaa                                                                                                                                                     <  6:994824/149‑1 (MQ=255)

                                                                                                                                                                                                                                        tAGTGAGCGTCAAAAAAAAACGCCGGATTAACCGGCGTCTGACGA‑CTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGaa                                                                                                                                                      <  6:718734/149‑1 (MQ=255)

                                                                                                                                                                                                                                         aGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGA‑CTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTACTTTGTATCCCGTGAGCTGACGGCGGCAAAAAGTCACCAGTGACTCTTCGGTAAGCGATGGATCTTTTTTCACTACGAAGATTTTCACCGCTTCACCACTGGAGCCGGAAGGt                   >  1:207537/1‑278 (MQ=255)

                                                                                                                                                                                                                                           tGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGA‑CTGACTTa                                                                                                                                                                                                                                                      >  2:757230/1‑49 (MQ=255)

                                                                                                                                                                                                                                             ggCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGA‑CTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCNCGAAACTCCACCAGCTTCGGTCATTTGTATCCCGTGAGCTGACGGCCGCAAAAAGTCACACGTGACTCTCCGGTAAGCGATgnt                                                                      >  3:243347/1‑221 (MQ=255)

                                                                                                                                                                                                                                             ggCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGA‑CTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTACTTTGTATCCCGTGAGCTGACGGCGGCAAAAAGTCACCAGTGACTCTTCGGTAAGCGATGGa                                                                      <  1:243347/223‑1 (MQ=255)

                                                                                                                                                                                                                                              gCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGA‑CTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTccacc                                                                                                                                              >  1:138643/1‑150 (MQ=255)

                                                                                                                                                                                                                                               cGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGA‑CTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTccacc                                                                                                                                              <  7:722907/149‑1 (MQ=255)

                                                                                                                                                                                                                                                gTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGA‑CTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTccacc                                                                                                                                              <  1:30949/148‑1 (MQ=255)

                                                                                                                                                                                                                                                gTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGA‑CTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTc                                                                                                                                                  >  4:525747/1‑144 (MQ=255)

                                                                                                                                                                                                                                                 tC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGA‑CTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAg                                                                                                                                            >  8:1057415/1‑149 (MQ=255)

                                                                                                                                                                                                                                                 tC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGA‑CTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTACTTTGTATCCCGTGAGCTGACGGCGGCAAAAAGTCACCAGTGACTCTTCGGTAAGCGATGGATCTTTTTTCACTACGAAGATTTTCACCGCTTCACCACTGGAGCCGGAAGGTACGcc              <  2:422653/275‑1 (MQ=255)

                                                                                                                                                                                                                                                  tcgaaaaaaaCGCCGGATTAACCGGCGTCTGACGA‑CTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCATAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAg                                                                                                                                            >  8:555243/4‑149 (MQ=255)

                                                                                                                                                                                                                                                     aaaaaaaCGCCGGATTAACCGGCGTCTGACGA‑CTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTccacc                                                                                                                                              <  2:525751/144‑1 (MQ=255)

                                                                                                                                                                                                                                                     aaaaaaaCGCCGGATTAACCGGCGTCTGACGA‑CTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTccac                                                                                                                                               <  6:925692/143‑1 (MQ=255)

                                                                                                                                                                                                                                                      aaaaaaCGCCGGATTAACCGGCGTCTGACGA‑CTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACt                                                                                                                                                               >  5:1134923/1‑126 (MQ=255)

                                                                                                                                                                                                                                                      aaaaaaCGCCGGATTAACCGGCGTCTGACGA‑CTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACt                                                                                                                                                               <  6:1134923/126‑1 (MQ=255)

                                                                                                                                                                                                                                                      aaaaaaCGCCGGATTAACCGGCGTCTGACGA‑CTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTc                                                                                                                                        <  6:410940/149‑1 (MQ=255)

                                                                                                                                                                                                                                                       tatacctccGGAGTAACGGGGCTCTGACGA‑CTGACTTAACACTCAGGTTTTATTGTCCACTTTTCCGCGCGCTTCGTCACGAAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGGAACTCATCACGAAACTCCACCAGCTTCGGTAATTTGTATCCCGTGAGCTGACGGCGGCAAAAAGTCACCAGTGAATCTTCGGTAAGCGATGGATCTTTTTTCACTACGAAGATTTTCACCGCTTCACCACTGGAGCCGGAAGGTACGCCAACAg         <  4:555389/268‑1 (MQ=255)

                                                                                                                                                                                                                                                       aaaaaCGCCGGATTAACCGGCGTCTGACGA‑CTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCAc                                                                                                                                                                                                            >  5:1467181/1‑80 (MQ=255)

                                                                                                                                                                                                                                                        aaaaCCCCGTATTTACCGGCGGC‑GANNACCTTNTTTAACCCTCCGGCTTTATAGTCCACTTTTCCGCGCGCTTCTTAACGTAAGTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTACTTTGTATCCCGTGAGCTGACGGCGGCAAAAAGTCACCAGTGACTCTTCGGTAAGCGATGGATCTTTTTTCACTACGAAGATTTTCACCGCTTCACCACTGGAGCCGGAAGGTACGCCAACAGCcgc     <  3:15193/278‑1 (MQ=255)

                                                                                                                                                                                                                                                         agaCGCCGGTTAATCCGGCGTCTGACGA‑CTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGt                                                                                                                                     >  7:921732/3‑149 (MQ=255)

                                                                                                                                                                                                                                                          gaCGCCGGTTAATCCGGCGTCTGACGA‑CTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTa                                                                                                                                    >  6:607385/2‑149 (MQ=255)

                                                                                                                                                                                                                                                            cGCCGGATTAACCGGCGTCTGACGA‑CTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTACTTTGTATCCCGTGAGCTGACGGCGGCAAAAAGTCACCAGTGACTCTTCGGTAAGCGATGGATCTTTTTTCACTACGAAGATTTTCACCGCTTCACCACTGGAGCCGGAAGGTACGCCAACAGCcgc     <  2:21378/274‑1 (MQ=255)

                                                                                                                                                                                                                                                              caggATTAACCGGCGTCTGACGA‑CTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGGCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAGCTCCACCAGCTTCGGTACttt                                                                                                                                <  8:1447263/147‑1 (MQ=255)

                                                                                                                                                                                                                                                               cGGATTAACCGGCGTCTGACGA‑CTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAAttttt                                                                                                                                                                                     >  8:494478/1‑95 (MQ=255)

                                                                                                                                                                                                                                                                  aTTAACCGGCGTCTGACGA‑CTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTACTTTGTAt                                                                                                                            <  6:1016004/149‑1 (MQ=255)

CGGATCGCGCAGGATCGCTTTCAGCGGTGACCAGTCGGTGATCCCGAGTGGATCGATTAGCGCCAGATGCTCGCCATCGAAAAGTTGAATCAACCCCAGCTGCGGGTAATAAGTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTC‑AAAAAAAACGCCGGATTAACCGGCGTCTGACGA‑CTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTACTTTGTATCCCGTGAGCTGACGGCGGCAAAAAGTCACCAGTGACTCTTCGGTAAGCGATGGATCTTTTTTCACTACGAAGATTTTCACCGCTTCACCACTGGAGCCGGAAGGTACGCCAACAGCCGCGAC  >  NC_000913/1887773‑1888300

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: `‑`