BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,814,860	A→G	5.1%	noncoding (59/78 nt)	micA →	sRNA antisense regulator of ompA, lamB, ompX, and phoP; Hfq‑dependent

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,814,860	0	A	G	5.1%	684.1 / 18.5	217	noncoding (59/78 nt)	micA	sRNA antisense regulator of ompA, lamB, ompX, and phoP; Hfq‑dependent
Reads supporting (aligned to +/- strand): ref base A (107/96); new base G (7/4); total (115/102)
Fisher's exact test for biased strand distribution p-value = 5.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.28e-01

Only 100 of 219 total aligned reads displayed.

CTGGCATCACTTCTTTGTTCGGCACGCAGAAGCGCAGATCGAACACGGTGATTGCGTCGCCATGCGGGGTGTTCATTGTTTTCGCCACCCGAACTGCAGGCGCTTCCATCCGGGTATGATCGACTGTGAAGCTATCTAACAACGGCATTTAGCCACCTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGC‑GTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCCTCTTCGCGCAGAATTTCCAGCGGCTCTTCACGCAGCAGATTACGGTAGGCTTCTGCAAATGCTTTGCCTGTTCCGCCAATACCAGTATCAATCATAGACCTTAAGATACGGGCAGAGAAAGTCAGATCGGGATTATCGAAGCA  >  NC_000913/2814586‑2815125

cTGGCATCACTTCTTTGTTCGGCACGCAGAAGCGCAGATCGAACACGGTGATTGCGTCGCCATGCGGGGTGTTCATTGTTTTCGCCACCCGAACTGCAGGCGCTTCCATCCGGGTATGATCGACTGTGAAGCTATCTAACAACGGCATTTAGCCACCTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCa                                                                                                                                                                                                                                                                            <  2:290372/275‑1 (MQ=255)

                                                                   ggTGTTCATTGTTTTCGCCACCCGAACTGCAGCCGCTTCCCTCCGGGTATGATCGACTGTGAAGCTATCTAACAACGGCATTTAGCCACCTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCACAGGCGAATCTGAGTATCTGAAAGACGCGCATTTGTTATCATCCTCCCTGAATTCAGCAATGAAATTTTGGCCACTAACGAGTGGCCTTTTTGTTTTATGTCAGGCGGTGTTTTCCCA‑CCACGCCTCACACGGTTCGTGATCacc                                                                                                                                                                                                        >  4:210263/1‑273 (MQ=255)

                                                                                                     gcTTCCATCCGGGTATGATCGACTGTGAAGCTATCTAACAACGGCATTTAGCCACCTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGC‑GTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGcc                                                                                                                                                                      <  2:666386/275‑1 (MQ=255)

                                                                                                          cATCCGGGTATGATCGACTGTGAAGCTATCTAACAACGGCATTTAGCCACCTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGC‑GTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTcgcg                                                                                                                                                                 <  4:102395/275‑1 (MQ=255)

                                                                                                             ccGGGTATGATCGACTGTGAAGCTATCTAACAACGGCATTTAGCCACCTCCGGTAATTTTTTTAAAACTTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTCGTTACCATCATCCCTGAATCCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCCGTC‑GTGTTTTTCCCGCAACATCGCCACTCGTTCGGTATCGGCTTCCTCCATTTCTTGCTGCCGCCGCTCatgcccccccagcc                                                                                                                                                               >  4:7926/1‑261 (MQ=255)

                                                                                                                               gAAGCTATCTAACAACGGCATTTAGCCACCTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCAc                                                                                                                                                                                                                                                                           <  6:1084662/149‑1 (MQ=255)

                                                                                                                                 aGCTATCTAACAACGGCATTTAGCCACCTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGa                                                                                                                                                                                                                                                                         <  6:1524122/149‑1 (MQ=255)

                                                                                                                                  gCTATCTAACAACGGCATTTAGCCACCTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAg                                                                                                                                                                                                                                                                        >  7:698061/1‑149 (MQ=255)

                                                                                                                                     aTCTAACAACGGCATTTAGCCACCTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTgg                                                                                                                                                                                                                                                                     >  6:81731/1‑149 (MQ=255)

                                                                                                                                     aTCTAACAACGGCATTTAGCCACCTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGC‑GTGTTTTTGCATCACAACCGCCAAGGGTTCGGTATCCGCGGCTTCCCTTTCCTGCTGACGGCCTTCAAACGCCCGGCGCCCCGCTACAAAATCCCCTTCACGCa                                                                                                                                       >  4:129894/1‑274 (MQ=255)

                                                                                                                                      tCTAACAACGGCATTTAGCCACCTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGC‑gtg                                                                                                                                                                                                                                            <  3:56394/172‑1 (MQ=255)

                                                                                                                                       cTAACAACGGCATTTAGCCACCTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCGTGAGTGGcc                                                                                                                                                                                                                                                                   <  5:408585/149‑1 (MQ=255)

                                                                                                                                        tAACAACGGCATTTAGCCACCTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCt                                                                                                                                                                                                                                                                  >  5:30522/1‑149 (MQ=255)

                                                                                                                                         aacaacGGCATTTAGCCACCTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCGTGAGTGGccaa                                                                                                                                                                                                                                                                 <  5:1201959/149‑3 (MQ=255)

                                                                                                                                            aacGGCATTTAGCCACCTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCttttt                                                                                                                                                                                                                                                              >  6:1300725/1‑149 (MQ=255)

                                                                                                                                            aacGGCATTTAGCCACCTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCttttt                                                                                                                                                                                                                                                              >  8:1084078/1‑149 (MQ=255)

                                                                                                                                                  aTTTAGCCACCTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTttttcttttc                                                                                                                                                                                                                                                        <  8:221116/149‑1 (MQ=255)

                                                                                                                                                  aTTTAGCCACCTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTttttcttttc                                                                                                                                                                                                                                                        <  7:587023/149‑1 (MQ=255)

                                                                                                                                                  aTTTAGCCACCTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTttttcttttc                                                                                                                                                                                                                                                        >  5:1356908/1‑149 (MQ=255)

                                                                                                                                                  aTTTAGCCACCTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTttttcttttc                                                                                                                                                                                                                                                        >  8:1540301/1‑149 (MQ=255)

                                                                                                                                                  aTTTAGCCACCTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTttttcttttc                                                                                                                                                                                                                                                        <  7:1137924/149‑1 (MQ=255)

                                                                                                                                                    ttAGCCACCTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTg                                                                                                                                                                                                                                                      <  8:142921/149‑1 (MQ=255)

                                                                                                                                                     tAGCCACCTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGt                                                                                                                                                                                                                                                     <  8:1689976/149‑1 (MQ=255)

                                                                                                                                                      aGCCACCTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCa                                                                                                                                                                                                                                                                            >  4:282068/1‑125 (MQ=255)

                                                                                                                                                        ccACCTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTttttctt                                                                                                                                                                                                                                                           >  6:1131684/1‑140 (MQ=255)

                                                                                                                                                          aCCTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAg                                                                                                                                                                                                                                                                        <  2:282071/125‑1 (MQ=255)

                                                                                                                                                          aCCTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGc                                                                                                                                                                                                                                                <  5:1300725/149‑1 (MQ=255)

                                                                                                                                                            cTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGC‑gt                                                                                                                                                                                                                                             >  8:473566/1‑149 (MQ=255)

                                                                                                                                                            cTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGC‑gt                                                                                                                                                                                                                                             >  8:1674406/1‑149 (MQ=255)

                                                                                                                                                             tCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGC‑gtg                                                                                                                                                                                                                                            >  5:514979/1‑149 (MQ=255)

                                                                                                                                                             tCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGC‑GTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCCTCTTCGCGCAGAATTTCCAGCGGCTCTTCACgcag                                                                                                              <  2:622284/275‑1 (MQ=255)

                                                                                                                                                             tCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGC‑GTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCCTCTTCGCGCAGAATTTCCAGCGGCTCTTCACgcag                                                                                                              <  2:614660/275‑1 (MQ=255)

                                                                                                                                                             tCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTAACGAGTGGCCTTTTTCTTTTCTGTCAGGC‑GTGTTTTGCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCCTCTTCGCGCAGAATTTCCAGCGGCTCTTCACgcag                                                                                                              <  2:545234/275‑1 (MQ=255)

                                                                                                                                                               cGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGC‑GTGtt                                                                                                                                                                                                                                          >  8:729999/1‑149 (MQ=255)

                                                                                                                                                                     tttttttAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGC‑GTGTTTTTCCa                                                                                                                                                                                                                                    <  5:1502246/149‑1 (MQ=255)

                                                                                                                                                                     tttttttAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGC‑GTGTTTTTCCa                                                                                                                                                                                                                                    <  5:117004/149‑1 (MQ=255)

                                                                                                                                                                      ttttttAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGC‑GTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCCTCTTCGCGCAGAATTTCCAGCGGCTCTTCACGCAGCAGATTACg                                                                                                     >  2:131333/1‑275 (MQ=255)

                                                                                                                                                                       tttttAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGC‑GTGTTTTTCCAGc                                                                                                                                                                                                                                  <  8:1460349/149‑1 (MQ=255)

                                                                                                                                                                            aaaaaTTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGC‑GTGTTTTTCCAGCCacac                                                                                                                                                                                                                             >  7:241240/1‑149 (MQ=255)

                                                                                                                                                                                aTTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGC‑GTGTTTTTCCAGCCACACCGCa                                                                                                                                                                                                                         >  8:409375/1‑149 (MQ=255)

                                                                                                                                                                                 ttttCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGC‑GTGTTTTTCCAGCCACACCGc                                                                                                                                                                                                                          <  5:600465/147‑1 (MQ=255)

                                                                                                                                                                                      tGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGC‑GTGTTTTTCCAGCCACACCGCAAACGGt                                                                                                                                                                                                                   <  6:862180/149‑1 (MQ=255)

                                                                                                                                                                                       gAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGC‑GTGTTTTTCCAGCCACACCGCAAACGGtt                                                                                                                                                                                                                  <  8:598432/149‑1 (MQ=255)

                                                                                                                                                                                             tttcttCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGC‑GTGTTTTTCCAGCCACACCGCAAACGGTTCGGTAt                                                                                                                                                                                                            >  6:281831/1‑149 (MQ=255)

                                                                                                                                                                                             tttcttCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGC‑GTGTTTTTCCAGCCACACCGCAAACGGTTCGGTAt                                                                                                                                                                                                            >  7:1646565/1‑149 (MQ=255)

                                                                                                                                                                                              ttcttcCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGC‑GTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATc                                                                                                                                                                                                           >  5:748739/1‑149 (MQ=255)

                                                                                                                                                                                              tnttgcccTGGCGAGTCGGATTATATGAACGTCNNTATTTTGTTATCCTCCTCCCTTAATTCAGANATGCATTTTTGGCCACTCACGTGGGGCCTTTTTCTTTTCGGTCAGGC‑GTGTTTTTTCAGCCACACCGCAAACGGNTCGGTATCAGCTGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCCTCTTCGCGCAGAATTTCCAGCGGCTCTTCACGCAGCAGATTACGGTAGGCTTCTGCAAATGCTTTGCCTGtt                                                                         <  3:125494/274‑1 (MQ=255)

                                                                                                                                                                                               tcttcCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGC‑GTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCa                                                                                                                                                                                                          <  8:698061/149‑1 (MQ=255)

                                                                                                                                                                                               tcttcCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGC‑GTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCa                                                                                                                                                                                                          <  5:1591671/149‑1 (MQ=255)

                                                                                                                                                                                                    ccAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCtt                                                                                                                                                                                                                                                                 <  6:133260/90‑1 (MQ=255)

                                                                                                                                                                                                      aGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGC‑GTGTTTTTCCAGCCACACCGCAAACGGTTc                                                                                                                                                                                                                 >  6:537121/1‑135 (MQ=255)

                                                                                                                                                                                                       ggCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGCGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGC‑GTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCTTCGCGCTCGACTACCAAATCCTCTTCGCGCAGCATTTCCAGCGGCTCTTCACGCACCAGATTACGGTATGCTTCTGCTACTGCTTTGCCTGTTCCGcc                                                                    >  2:344045/1‑275 (MQ=255)

                                                                                                                                                                                                             tCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCCCGAGTGGCCTTTTTCTTTTCTGTCAGGC‑GTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCCTCTTCGCGCAGAATTTCCAGCGGCTCTTCACGCAGCAGATTACGGTAGGCTTCTGCAAATGCTTTGCCTGTTc                                                                        >  2:123637/1‑265 (MQ=255)

                                                                                                                                                                                                             tCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGC‑GTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTc                                                                                                                                                                                            >  5:834713/1‑149 (MQ=255)

                                                                                                                                                                                                             tCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGC‑GTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTc                                                                                                                                                                                            <  6:992234/149‑1 (MQ=255)

                                                                                                                                                                                                              cTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGC‑GTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTcc                                                                                                                                                                                           <  5:1315285/149‑1 (MQ=255)

                                                                                                                                                                                                               tGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGAAGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGC‑GTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCct                                                                                                                                                                                          <  6:1415409/149‑1 (MQ=255)

                                                                                                                                                                                                                gAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGC‑GTG‑TTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCctgc                                                                                                                                                                                        <  7:729999/149‑1 (MQ=255)

                                                                                                                                                                                                                  gTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGC‑GTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCctgct                                                                                                                                                                                       >  6:570753/1‑149 (MQ=255)

                                                                                                                                                                                                                   tataTGAAAGACGCGCGTTTGTTATCATCATCCCTGAATTCAGAGATGAAATGTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGC‑GTGTTTTTCCAGCCACACCGCaaa                                                                                                                                                                                                                       >  7:903848/1‑116 (MQ=255)

                                                                                                                                                                                                                   tataTGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGC‑GTGTTTTTCCAGCCACACCGCaaa                                                                                                                                                                                                                       <  8:903848/116‑1 (MQ=255)

                                                                                                                                                                                                                   tataTGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGC‑GTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCctgctg                                                                                                                                                                                      >  5:784956/1‑149 (MQ=255)

                                                                                                                                                                                                                    atatGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGC‑GTGTTTTTCCAGCCACACCGCAAACGGTTc                                                                                                                                                                                                                 <  6:1437261/121‑1 (MQ=255)

                                                                                                                                                                                                                     tatGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGC‑GTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGAc                                                                                                                                                                                    <  5:281831/149‑1 (MQ=255)

                                                                                                                                                                                                                      atGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGC‑GTGTTTTTCCAGCCACACCGCAAACGGtt                                                                                                                                                                                                                  <  7:767905/118‑1 (MQ=255)

                                                                                                                                                                                                                      atGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGC‑GTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCCCGCTCGGCTACAAAATCCTCTTCGCGCAGAATTTCCAGCGACTCTTCTCGCCGCAGATCACGGTAGCCTTCTGCAAATGCTTTGCCTGTTCCGCCATTCCCCGTatccat                                                     >  4:65024/1‑275 (MQ=255)

                                                                                                                                                                                                                        gAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGC‑GTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGc                                                                                                                                                                                 >  6:759025/1‑149 (MQ=255)

                                                                                                                                                                                                                          agagacaccATTTGTTTGCTTCAGCCCTGAATTCAGAGATGAAATTTTGGCCACCCACGAGTGGCCTTTTTCTTTTCTGTCAGGC‑GTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCCTCTTCGCGCAGAATTTCCAGCGGCTCTTCACGCAGCAGATTACGGTAGGCTTCTGCAAATGCTTTGCCTGGTCCGCCAATACCAGTATCAATCATa                                                 <  4:617182/267‑1 (MQ=255)

                                                                                                                                                                                                                          aaGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGC‑GTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGt                                                                                                                                                                               <  6:784956/149‑1 (MQ=255)

                                                                                                                                                                                                                           aGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGC‑GTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCCTCTTCGCGCAGAATTTCCAGCGGCTCTTCACGCAGCAGATTACGGTAGGCTTCTGCAAATGCTTTGCCTGTTCCGCCAATACCAGTATCAATCATa                                                 <  2:391277/274‑1 (MQ=255)

                                                                                                                                                                                                                             aCGCGCATTTGTTATCATCATCCCTGAATTCAGGGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGC‑GTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCCTCTTCGCGCAGAATTTCCAGCGGCTCTTCACGCAGCAGATTACGGTAGGCTTCTGCAAATGCTTTGCCTGTTCCGCCAATACCAGTat                                                         <  4:357592/264‑1 (MQ=255)

                                                                                                                                                                                                                              ctcgcATTTGTTATCATCATCCCTGAATTCAGGGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGC‑GTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCGGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCCTCTTCGCGCAGAATTTCCAGCGGCTCTTCACGCAGCAGATTACGGTAGGCTTCTGCAAATGCTTTGCCTGTTCCGCCAATACCAGTat                                                         <  1:73464/261‑1 (MQ=255)

                                                                                                                                                                                                                               gcgcATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGC‑GTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGTCGGCGTTCAGa                                                                                                                                                                          >  5:1421637/1‑149 (MQ=255)

                                                                                                                                                                                                                                 gcATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGC‑GTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCgg                                                                                                                                                                                                      <  2:736376/119‑1 (MQ=255)

                                                                                                                                                                                                                                   aTTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGC‑GTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGAc                                                                                                                                                                                    <  2:140215/135‑1 (MQ=255)

                                                                                                                                                                                                                                        ttATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGC‑GTGTTTTTCCAGCCACACCGCAAACGGtt                                                                                                                                                                                                                  <  5:655769/100‑1 (MQ=255)

                                                                                                                                                                                                                                        ttATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGC‑GTGTTTTTCCAGCCACACCGCAAACGGtt                                                                                                                                                                                                                  >  6:655769/1‑100 (MQ=255)

                                                                                                                                                                                                                                        ttATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGC‑GTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTcgcg                                                                                                                                                                 >  5:1081568/1‑149 (MQ=255)

                                                                                                                                                                                                                                        ttATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGC‑GTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTcgcg                                                                                                                                                                 >  6:954110/1‑149 (MQ=255)

                                                                                                                                                                                                                                              tcatcCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGC‑GTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGc                                                                                                                                                           >  7:1477748/1‑149 (MQ=255)

                                                                                                                                                                                                                                              tcatcCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGC‑GTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTAAGACGCCTCGCGCTCGGc                                                                                                                                                           >  8:11898/1‑149 (MQ=255)

                                                                                                                                                                                                                                                  cccTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGC‑GTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACa                                                                                                                                                       >  6:1638097/1‑149 (MQ=255)

                                                                                                                                                                                                                                                   ccTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGC‑GTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCt                                                                                                                                                               <  4:545753/140‑1 (MQ=255)

                                                                                                                                                                                                                                                    cTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGC‑GTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACaaa                                                                                                                                                     >  8:1409785/1‑149 (MQ=255)

                                                                                                                                                                                                                                                            agagATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGC‑GTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCCTCtt                                                                                                                                             <  5:1481916/149‑1 (MQ=255)

                                                                                                                                                                                                                                                            agagATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGC‑GTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCCTCTTCGCGCAGAATTTCCAGCGGCTCTTCACGCAGCAGATTACGGTAGGCTTCTGCAAATGCTTTGCCTGTTCCGCCAATCCCAGTATCAATCATAGACCTTAAGATACGGGCAGATAAAGTCAGAACgg               >  2:458980/1‑275 (MQ=255)

                                                                                                                                                                                                                                                                 tGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGC‑GTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCt                                                                                                                                                                                                    <  7:1279154/89‑1 (MQ=255)

                                                                                                                                                                                                                                                                 tGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGC‑GTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCCTCTTCgcgc                                                                                                                                        >  7:20885/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                  gAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGC‑GTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCCTCTTCGCGCa                                                                                                                                       <  5:658497/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                  gAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGC‑GTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCCTCTTCGCGCa                                                                                                                                       >  7:1534037/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                  gAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGC‑GTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCCTCTTCGCGCAGAATTTCCAGCGGCTCTTCACGCAGCAGATTACGGTAGGCTTCTGCAAATGCTTTGCCTGTTCCGCCAATACCAGTATCAATCATAGACCTTAAGATACGGGCAGAGAAAGTCAGATCGGGATTATCGa      <  1:226834/278‑1 (MQ=255)

                                                                                                                                                                                                                                                                     aTTTTGGCCACTCACGAGTTGCCTTTTTCTTTTCTGTCAGGC‑GTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCCTCTTCGCGCAGAATTTCCAGCGGCTCTTCACGCAGCAGATTACGGTAGGCTTCTGCAAATGCTTTGCCTGTTCCGCCAATACCAGTATCAATCATAGACCTTAAGATACGGGCAGAGAAAGTCAGATCGGGATTATCGAAGCa  <  1:15807/279‑1 (MQ=255)

                                                                                                                                                                                                                                                                      ttttGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGC‑GTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCCTCTTCGCGCAGAAt                                                                                                                                   >  8:1310168/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                      ttttGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGC‑GTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCCTCTTCGCGCAGAAt                                                                                                                                   >  5:46765/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                       tttGGCCACTCACGAGTGGCCTTTTTTTTTTCTGTCAGGC‑GTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCCTCTTCGCGCAGAATTTCCAGCGGCTCTTCACGCAGCAGATTACGGTAGGCTTCTGCAAATGCTTTGCCTGTTCCGCCAATACCAGTATCAATCATAGACCTTAAGATACGGGCAGAGAAAGTCAGATCGGGATTATCGAAg    <  4:132807/275‑1 (MQ=255)

                                                                                                                                                                                                                                                                        aaggCCACTCGTGAGTGGCCTTTTTCTTTTCTGTCAGGC‑GTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCCTCTTCGCGCAGAAttt                                                                                                                                 >  5:1032874/3‑149 (MQ=255)

                                                                                                                                                                                                                                                                         agtccACTCGTGAGTGGCCTTTTTCTTATGTGTCAGGC‑GTGTTTTTCCAGCCACACCGCAAACGGTTCGGAATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCCTCTTCGCGCAGAATTTc                                                                                                                                >  6:304485/4‑149 (MQ=255)

                                                                                                                                                                                                                                                                         aggCCACTCGTGAGTGGCCTTTTTCTTTTCTGTCAGGC‑GTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTc                                                                                                                                                              <  7:443675/118‑1 (MQ=255)

                                                                                                                                                                                                                                                                         aggCCACTCGTGAGTGGCCTTTTTCTTTTCTGTCAGGC‑GTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCCTCTTCGCGCAGAATTTc                                                                                                                                >  8:438646/2‑149 (MQ=255)

                                                                                                                                                                                                                                                                            ccACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGC‑GTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCCTCTTCGCGCAGAATTTCCAg                                                                                                                             >  5:13893/1‑149 (MQ=255)

CTGGCATCACTTCTTTGTTCGGCACGCAGAAGCGCAGATCGAACACGGTGATTGCGTCGCCATGCGGGGTGTTCATTGTTTTCGCCACCCGAACTGCAGGCGCTTCCATCCGGGTATGATCGACTGTGAAGCTATCTAACAACGGCATTTAGCCACCTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGC‑GTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCCTCTTCGCGCAGAATTTCCAGCGGCTCTTCACGCAGCAGATTACGGTAGGCTTCTGCAAATGCTTTGCCTGTTCCGCCAATACCAGTATCAATCATAGACCTTAAGATACGGGCAGAGAAAGTCAGATCGGGATTATCGAAGCA  >  NC_000913/2814586‑2815125

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