BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,423,231	T→A	6.3%	intergenic (+37/+192)	yhdZ → / ← rrfF	putative amino acid ABC transporter ATPase/5S ribosomal RNA of rrnD operon

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,423,231	0	T	A	6.3%	583.4 / 20.8	171	intergenic (+37/+192)	yhdZ/rrfF	putative amino acid ABC transporter ATPase/5S ribosomal RNA of rrnD operon
Reads supporting (aligned to +/- strand): ref base T (85/75); new base A (6/5); total (91/80)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.51e-01

Only 100 of 174 total aligned reads displayed.

AAAGAGGTGCTGGATACGATGATTGGGCTGGCGCAGTCGGGTATGACAATGTTGTGTGTAACACATGAGATGGGGTTTGCACGAACCGTCGCTGACCGGGTAATTTTTATGGATCGTGGGGAAATAGTGGAGCAAGCTGCACCTGA‑TGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGC  >  NC_000913/3422985‑3423502

                                                 tgttgtGTGTAACACATGAGATGGGGTTTGCACGAACCGTCGCTGACCGGGTAATTTTTATGGATCGTGGGGAAATAGTGGAGCAAGCTGCACCTGA‑TGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTTCCCGCCACCCCCCGGGCCGTTTCtt                                                                                                                                                                                                    >  4:520817/1‑275 (MQ=255)

                                                  gttgtGTGTAACACATGAGATGGGGTTTGCACGAACCGTCGCTGACCGGGTAATTTTTATGGATCGTGGGGAAATAGTGGAGCAAGCTGCACCTGA‑TGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCAACCTCCGGGCCGTTTCtt                                                                                                                                                                                                    >  2:337228/1‑274 (MQ=255)

                                                      tgtgcAGCACATGAGATGGGGTTGGCACGAACCGTCGCTGACCGGGTAATTTTGATTGATCGTGTGGAAATAGTGGAGCAAGCTGCACCTGA‑TGAATTTTGTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCa                                                                                                                                                                                                <  1:1432/274‑1 (MQ=255)

                                                             cacaTGAGATGGGGTTTGCACGAACCGTCGCTGACCGGGTAATTTTTATGGATCGTGGGGAAATAGTGGAGCAAGCTGCACCTGA‑TGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAAAGTTCaa                                                                                                                                                                                        >  2:558981/1‑275 (MQ=255)

                                                                       ggggTTTGCACGAACCGTCGCTGACCGGGTAATTTTTATGGATCGTGGGGAAATAGTGGAGCAAGCTGCACCTGA‑TGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTccc                                                                                                                                                                              >  1:87900/1‑275 (MQ=255)

                                                                            ttGCACGAACCGTCGCTGACCGGGTAATTTTTATGGATCGTGGGGAAATAGTGGAGCAAGCTGCACCTGA‑TGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCCCCCCcggcgg                                                                                                                                                                         >  2:594224/1‑275 (MQ=255)

                                                                             tGCACGAACCGTCGCTGACTGGGTAATTTTTAGGGATCGCGGGGAAATAGTGGAGCAAGCTGCACCTGA‑TGAATTTATTGCGCATCCTAAATCAGCGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATAGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTCGATGTCTGGCAGTCAATGGCGGGGGTCCTGCCCGCCACCCCACGGGCCGGTGCTTGGCTACGGTGAAACCCGCCACcggacgg                                                                                                                                                                        >  4:593825/1‑271 (MQ=255)

                                                                                       gtccgtGACAGGTTAATTTTTAGGGTTGAGGGGGAAAATGTTGAGCAAACGGACCCTGA‑TGAATTTTTTTCGCATCCTAAATCATAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGCTTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTa                                                                                                                                                               <  4:294022/269‑1 (MQ=25)

                                                                                         tgCTGACCGGGGAATTTTTATGGATCGTGGGGAAATAGTGGAGCAAGCTGCACCTGA‑TGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACt                                                                                                                                                             <  2:418854/273‑1 (MQ=255)

                                                                                             gcacGGGTAATTTTTATGGTTCGTGGGGAAATAGTGGGGCAACCTGCAACTGA‑TGAATTTTTTGCGCATCCTAAATCAGAGGGTACGGGGGCATTTTTAGCGCAGGTAATCCATTAATTAAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCggg                                                                                                                                                         <  4:558981/271‑1 (MQ=255)

                                                                                              aCCGGGTAATTTTTATGGATCGTGGGGAAATAGTGGAGCAAGCTGCACCTGA‑TGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCataaactgccag                                                                                                                                                                                                           >  2:486372/1‑212 (MQ=255)

                                                                                                      aTTTTTATGGATCGTGGGGAAATAGTGGAGCAAGCTGCACCTGA‑TGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGa                                                                                                                                                                                                                                                                              <  6:536943/149‑1 (MQ=255)

                                                                                                          ttATGGATCGTGGGGAAATAGTGGAGCAAGCTGCACCTGA‑TGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGc                                                                                                                                                                                                                                                                          <  8:835006/149‑1 (MQ=255)

                                                                                                             tGGATCGTGGGGAAATAGTGGAGCAAGCTGCACCTGA‑TGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCTTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCCGTTTATGGCGGGCGTCCTTCCCTCCACCCTCCGGGCCGTTGCTTTGCAACGTTCAAAGCCGCTCCCCGCGGGTTTGTCCTCCTCGCTCAAGTGTTCACCGaaca                                                                                                                                        >  4:593970/1‑272 (MQ=255)

                                                                                                             tGGATCGTGGGGAAATAGTGGAGCAAGCTGCACCTGA‑TGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCGAAAGGATGGCC‑tt                                                                                                                                                                                                                                                                      <  5:305452/149‑1 (MQ=255)

                                                                                                               gATCGTGGGGAAATAGTGGAGCAAGCTGCACCTGA‑TGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTc                                                                                                                                                                                                                                                                    <  8:418133/149‑1 (MQ=255)

                                                                                                                 tCGTGGGGAAATAGTGGAGCAAGCTGCACCTGA‑TGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGc                                                                                                                                                                                                                                                                  <  6:763247/149‑1 (MQ=255)

                                                                                                                 tCGTGGGGAAATAGTGGAGCAAGCTGCACCTGA‑TGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGc                                                                                                                                                                                                                                                                  >  8:126180/1‑149 (MQ=255)

                                                                                                                  cGTGGGGAAATAGTGGAGCAAGCTGCACCTGA‑TGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATAGCC‑TTTCGCt                                                                                                                                                                                                                                                                 <  8:311236/149‑1 (MQ=255)

                                                                                                                   agagacggaTAGTGGAGCAAGCTGAACCTTAGTGATTTTTTTG‑GAACCCTAAATAAGAGGGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAGAAGGCCATCATTTCGGATGGCC‑TTTCGCTTGATTTGAGGTCTGGCAGTTTATGGCGGGCGTCCTGTCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTAACCGACAaacaac                                                                                                                                   <  4:283140/266‑1 (MQ=21)

                                                                                                                    tGGGGAAATAGTGGAGCAAGCTGCACCTGA‑TGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCttg                                                                                                                                                                                                                                                               >  7:591964/1‑149 (MQ=255)

                                                                                                                          aaTAGTGGAGCAAGCTGCACCTGA‑TGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCttgatttga                                                                                                                                                                                                                                                         <  5:572356/149‑1 (MQ=255)

                                                                                                                           aTAGTGGAGCAAGCTGCACCTGA‑TGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCttgatttgat                                                                                                                                                                                                                                                        <  7:633928/149‑1 (MQ=255)

                                                                                                                           aTAGTGGAGCAAGCTGCACCTGA‑TGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCttgatttgat                                                                                                                                                                                                                                                        <  8:362829/149‑1 (MQ=255)

                                                                                                                              gTGGAGCAAGCTGCACCTGA‑TGAATTTTTTGCGCATCCTAAATCAGAGCGTACGGGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCGAAAGGATGGCC‑ttt                                                                                                                                                                                                                                                                     <  7:755456/133‑1 (MQ=255)

                                                                                                                              gTGGAGCAAGCTGCACCTGA‑TGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCGAAAGGATGGCC‑ttt                                                                                                                                                                                                                                                                     >  8:755456/1‑133 (MQ=255)

                                                                                                                                ggAGCAAGCTGCACCTGA‑TGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTg                                                                                                                                                                                                                                                   <  5:282269/149‑1 (MQ=255)

                                                                                                                                  aGCAAGCTGCACCTGA‑TGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGc                                                                                                                                                                                                                                                 >  5:1003534/1‑149 (MQ=255)

                                                                                                                                   gCAAGCTGCACCTGA‑TGAATTTTTTGCGCATCCTAAATCAGGGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCa                                                                                                                                                                                                <  4:305626/197‑1 (MQ=255)

                                                                                                                                   gCAAGCTGCACCTGA‑TGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGc                                                                                                                                                                                                                                                                          >  7:106151/1‑124 (MQ=255)

                                                                                                                                    cAAGCTGCACCTGA‑TGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTACTCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTTAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACCGATAACACAAAAGAcca                                                                                                                 >  2:145361/1‑274 (MQ=255)

                                                                                                                                          gCACCTGA‑TGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATg                                                                                                                                                                                                                                         <  7:908095/149‑1 (MQ=255)

                                                                                                                                              cTGA‑TGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATggcgg                                                                                                                                                                                                                                     <  6:1003534/149‑1 (MQ=255)

                                                                                                                                              cTGA‑TGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATggcgg                                                                                                                                                                                                                                     >  6:418727/1‑149 (MQ=255)

                                                                                                                                               tga‑tgaATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGa                                                                                                       <  2:520817/274‑1 (MQ=255)

                                                                                                                                                ga‑tgaATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATggcgggc                                                                                                                                                                                                                                   >  7:999221/1‑149 (MQ=255)

                                                                                                                                                ga‑tgaATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATggcgggc                                                                                                                                                                                                                                   >  8:591272/1‑149 (MQ=255)

                                                                                                                                                    gaATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCttgatttga                                                                                                                                                                                                                                                         >  6:981818/1‑124 (MQ=255)

                                                                                                                                                    gaATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCttgatttga                                                                                                                                                                                                                                                         <  5:981818/124‑1 (MQ=255)

                                                                                                                                                     aaTTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGc                                                                                                                                                                                                                                                                          <  5:132601/107‑1 (MQ=255)

                                                                                                                                                                 aTCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCAccc                                                                                                                                                                                                                   >  6:716192/1‑149 (MQ=255)

                                                                                                                                                                    cTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTGGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTcc                                                                                                                                                                                                                >  5:1062505/1‑149 (MQ=255)

                                                                                                                                                                    cTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGt                                                                                                                                                                                                                                                      >  7:203006/1‑111 (MQ=255)

                                                                                                                                                                    cTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGt                                                                                                                                                                                                                                                      <  8:203006/111‑1 (MQ=255)

                                                                                                                                                                     tAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCg                                                                                                                                                                                                               >  5:832907/1‑149 (MQ=255)

                                                                                                                                                                        aTCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCC‑CCGGGcc                                                                                                                                                                                                           <  5:418727/149‑1 (MQ=255)

                                                                                                                                                                            gagCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGtt                                                                                                                                                                                                        <  7:162410/149‑1 (MQ=255)

                                                                                                                                                                            gagCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGtt                                                                                                                                                                                                        <  8:440588/149‑1 (MQ=255)

                                                                                                                                                                             agCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATaaa                                                                                                                            >  4:194432/1‑224 (MQ=255)

                                                                                                                                                                             agCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTgg                                                                          >  2:598184/1‑274 (MQ=255)

                                                                                                                                                                             agCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCAACTGAGCCTTTTGTTTTATTTGATTTCTgg                                                                          >  2:273027/1‑274 (MQ=255)

                                                                                                                                                                              gCGTACGGGGGTTTTTTTTTCGCAGGTAATCCATTAATTTAATTGTAGTTCGAAAAGCAAAAAGGTCATCCGTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCa                                                                        <  4:145361/275‑1 (MQ=255)

                                                                                                                                                                              gCGTACGAGGGCATTTTTATCGCAGGTAATGCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGc                                                                                                                                                                                                      <  5:716192/149‑1 (MQ=255)

                                                                                                                                                                              gCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGACAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGc                                                                                                                                                                                                      <  8:562208/149‑1 (MQ=255)

                                                                                                                                                                              gCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGc                                                                                                                                                                                                      <  6:135587/149‑1 (MQ=255)

                                                                                                                                                                                 tACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑tt                                                                                                                                                                                                                                                                      >  1:5901/1‑82 (MQ=255)

                                                                                                                                                                                 tACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGgaga                                                                                                                                                      <  1:126340/194‑1 (MQ=255)

                                                                                                                                                                                 tACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACaa                                                                                                                        <  2:194432/224‑1 (MQ=255)

                                                                                                                                                                                    ggtggCATTTTTATCGCAGGTAATCCATTAATTGAATGTTCGTTCGAAAAGCAAAAAGGCCCTCCTGTCGTGTGGCC‑TTTCCCTTGATTTGATGTCTGTCAGTTTATGGCGGGCGTCCTGCCCTCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTccc                                                                  <  4:598184/272‑1 (MQ=255)

                                                                                                                                                                                     aGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGc                                                                                                                                                                                                                                                                  <  3:5901/82‑1 (MQ=255)

                                                                                                                                                                                      gggCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAAc                                                                                                                                                                                              <  8:839344/149‑1 (MQ=255)

                                                                                                                                                                                        gCATTTTTATTGCAGGGAATACATTAATTGAATGTTAGGTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTAct                                                              <  1:22400/275‑1 (MQ=255)

                                                                                                                                                                                         cATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGtt                                                                                                                                                                                           <  5:173899/149‑1 (MQ=255)

                                                                                                                                                                                          aTTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTc                                                                                                                                                                                          >  7:772629/1‑149 (MQ=255)

                                                                                                                                                                                          aTTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTc                                                                                                                                                                                          >  7:377376/1‑149 (MQ=255)

                                                                                                                                                                                           tttttATCGCAGGTAATCCATTAATTGCATGTTAGTTCGAAAAGCAAAAAGCCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGCAGTCCTGCCCGCCACCCTCCTGGCCGGCTCTTCCCAACGTTTAAATTCCCTCCCGTCCGTGTCTTCCTTCCCGGGGGAGTGTTCACCGACCAACAACCGCTAACACAAAAGCCCCCGTCTTCCGCTTTCACCTTCTGTTTTATTCTATTTCTGGCCGTTCCctcctcc                                                            >  4:300299/1‑273 (MQ=255)

                                                                                                                                                                                                    cAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTc                                                                                                                                                                                >  5:5755/1‑149 (MQ=255)

                                                                                                                                                                                                      ggTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGgag                                                 >  2:432835/1‑274 (MQ=255)

                                                                                                                                                                                                       gTAATACATTAATTAAATGTTAGTGCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGgaga                                                <  1:165294/274‑1 (MQ=255)

                                                                                                                                                                                                           tCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcggcgg                                                                                                                                                                         >  8:272142/1‑149 (MQ=255)

                                                                                                                                                                                                              aTTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAtt                                                                                                                                                                      >  7:699050/1‑149 (MQ=255)

                                                                                                                                                                                                              aTTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAtt                                                                                                                                                                      >  6:627153/1‑149 (MQ=255)

                                                                                                                                                                                                              aTTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAtt                                                                                                                                                                      >  8:662272/1‑149 (MQ=255)

                                                                                                                                                                                                               ttAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAttt                                                                                                                                                                     <  5:1102537/149‑1 (MQ=255)

                                                                                                                                                                                                               ttAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAacaaca                                                                                                                                  >  2:288687/1‑184 (MQ=255)

                                                                                                                                                                                                               ttAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATggg                                                    >  2:467755/1‑262 (MQ=255)

                                                                                                                                                                                                               ttAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCacacc                                       >  1:69112/1‑274 (MQ=255)

                                                                                                                                                                                                                   ttGAATGTTAGTTCGAAAAGCAAAAATGCCATCCTTACGGATGGCC‑TTTCGCTTGATTTGATATCTGGCAGTTTATGGCGGACGGCCTTCCCGCCACCTTTCGGGCCGTTGCTTCTCAACGTTCAAATCCGCTCCCCCCCGCTTTGTcc                                                                                                                                                                 >  8:929686/1‑149 (MQ=18)

                                                                                                                                                                                                                   ttGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTcc                                                                                                                                                                 >  6:933647/1‑149 (MQ=255)

                                                                                                                                                                                                                   ttGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGgaga                                                <  4:467755/262‑1 (MQ=255)

                                                                                                                                                                                                                   ttGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGACAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTcc                                                                                                                                                                 >  7:215526/1‑149 (MQ=255)

                                                                                                                                                                                                                          ttAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCgg                                                                                                                                                          >  6:824060/1‑149 (MQ=255)

                                                                                                                                                                                                                            aGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGga                                                                                                                                                        <  6:5755/149‑1 (MQ=255)

                                                                                                                                                                                                                               tagAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGagag                                                                                                                                                     <  6:651389/147‑1 (MQ=255)

                                                                                                                                                                                                                                   aaaGCAAAAAGGCCATCCTGTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGATATTGtt                                                                                                                                                 <  8:105140/149‑1 (MQ=255)

                                                                                                                                                                                                                                    aaGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTc                                                                                                                                                >  7:832040/1‑149 (MQ=255)

                                                                                                                                                                                                                                    aaGCAAAAAGGCCATCCTTTCGGATGGCC‑CTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGttt                  <  4:420096/275‑1 (MQ=255)

                                                                                                                                                                                                                                     aagcgaaaGGCCATCCGAAAGGATGGCCTTTTTGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAgtg                                                                                                                                                   <  6:488484/141‑1 (MQ=11)

                                                                                                                                                                                                                                     aagcgaaaGGCCATCCGAAAGGATGGCCTTTTTGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAgtg                                                                                                                                                   >  5:488484/6‑146 (MQ=11)

                                                                                                                                                                                                                                      gCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAg                                                                                                                     <  2:671337/174‑1 (MQ=255)

                                                                                                                                                                                                                                         aaaaGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGa                                                                                                                                           >  7:992165/1‑149 (MQ=255)

                                                                                                                                                                                                                                            aGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGacaa                                                                                                                                        >  8:515901/1‑149 (MQ=255)

                                                                                                                                                                                                                                              gCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGacaaac                                                                                                                                      <  8:992165/149‑1 (MQ=255)

                                                                                                                                                                                                                                               ccATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGacaaa                                                                                                                                       <  8:543978/147‑1 (MQ=255)

                                                                                                                                                                                                                                               ccATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGacaaa                                                                                                                                       >  7:543978/1‑147 (MQ=255)

                                                                                                                                                                                                                                                cATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGacaaacaa                                                                                                                                    >  7:179256/1‑149 (MQ=255)

                                                                                                                                                                                                                                                 aTCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAaacaac                                                                                                                                   >  6:443767/1‑149 (MQ=255)

                                                                                                                                                                                                                                                 aTCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTTGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCCCTACGGCGTTTACCTTCTGAGTTc     >  2:371001/1‑275 (MQ=255)

                                                                                                                                                                                                                                                    cTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGa                                                                                                                                <  6:285677/149‑1 (MQ=255)

                                                                                                                                                                                                                                                     tttGGGATGGCC‑TTTCGCTTGATTTGGGGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGc  <  4:377037/274‑1 (MQ=255)

AAAGAGGTGCTGGATACGATGATTGGGCTGGCGCAGTCGGGTATGACAATGTTGTGTGTAACACATGAGATGGGGTTTGCACGAACCGTCGCTGACCGGGTAATTTTTATGGATCGTGGGGAAATAGTGGAGCAAGCTGCACCTGA‑TGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCC‑TTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGC  >  NC_000913/3422985‑3423502

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