BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,630,667	A→G	10.4%	intergenic (‑65/+301)	yhiI ← / ← yhiJ	putative membrane fusion protein (MFP) of efflux pump/DUF4049 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,630,667	0	A	G	10.4%	465.6 / 34.1	155	intergenic (‑65/+301)	yhiI/yhiJ	putative membrane fusion protein (MFP) of efflux pump/DUF4049 family protein
Reads supporting (aligned to +/- strand): ref base A (70/68); new base G (8/8); total (79/76)
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.20e-01

Only 100 of 155 total aligned reads displayed.

CTTCTTTCACCAGAATGGTGTCGATACGCCCGGCAATTTTGCTGGCAATATCCACTTCCGTCGCTTCAATGCGCCCATTACTGACAGCAAAGCCTTCCGGCACACCTGCCGGGCGCAACAGCCACCACGCCACGATAGCCGCCACCGCCAGTAACCCGACAACCCACCACGCCAGATGGCGCTTACTCTTATCCATAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAA‑CTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATAC‑GACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGT  >  NC_000913/3630407‑3630926

                ggggctgttgcGCCCGGCAATTTTTCTGGCAATAACCACTTCCGTCGCTTCAATGCCCCCATTACTGACAGCAAAGCCTTCCGGCACACCTTCCGGGCGCACACGCCACCACCCCCCGCTAACCCCCACCCCCAGTCACCCGACCACCCACCACGCCACATGGCGCTTACCCCTCCTCATAAGCCACCCCCCATAATCCCCGTAAACCCACCCCACGTGTTACAACCCGCCCCCCAAC‑ACCTTCCGCGGATTCCTCCCCCCGGGATAAAACGGa                                                                                                                                                                                                                                         >  4:416905/11‑274 (MQ=255)

                       aTACGCCCGGCAATTTTGCTGGCAATATCCACTTCCGTCGCTTCAATGCGCCCATTACTGACAGCAAAGCCTTCCGGCACACCTGCCGGGCGCAACAGCCACCACGCCACGATAGCCGCCACCGCCAGTAACCCGACAACCCACCACGCCAGATGGCGCTTACTCTTATCCATAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAA‑CTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCa                                                                                                                                                                                                                                  >  2:280859/1‑274 (MQ=255)

                                gCAATTTTGCTGGCAATATCCACTTCCGTCGCTTCAATGCGCCCATTACTGACAGCAAAGCCTTCCGGCACACCTGCCGGGCGCAACAGCCACCACGCCACGATAGCCGCCACCGCCAGTAACCCGACAACCCACCACGCCAGATGGCGCTTACTCTTATCCATAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAA‑CTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATAC‑GACCa                                                                                                                                                                                                                       >  2:359167/1‑275 (MQ=255)

                                                       ttCCGTCGCTTCAATGCGCCCATTACTGACAGCAAAGCCTTCCGGCACACCTGCCGGGCGCAACAGCCACCACGCCACGATAGCCGCCACCGCCAGTAACCCGACAACCCACCACGCCAGATTGCGCTTACTCTTATCCATAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAA‑CTGCCTCCGAACGCCGTTCAGCGGATTCTTCCAACCGGGATAAGAAGGAAACGGCAAATAC‑GACCAGCGCGCATATCGTTGAGCGCCaa                                                                                                                                                                                                >  4:247158/1‑275 (MQ=255)

                                                            tCGCTTCAATGCGCCCATTACTGACAGCAAAGCCTTCCGGCACACCTGCCGGGCGCACCAGCCACCACGCCACGATAGCCGCCACCGCCAGTAACCCGACAACCCACCACGCCAGATGGCGCTTACTCTTATCCATAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAA‑CTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATAC‑GACCAGCGCGCAGAGCATGTGGCGCCAACgggg                                                                                                                                                                                           >  2:525852/1‑275 (MQ=255)

                                                                ttCAATGCGCCCATTACTGACAGCAAAGCCTTCCGGCACACCTGCCGGGCGCAACAGCCACCACGCCACGATAGCCGCCACCGCCAGTAACCCGACAACCCACCACGCCAGATGGCGCTTACTCTTATCCATAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAA‑CTGCCCCCGAAAGACGTTCAGCGGATTCTTCCCTCCGGGCTAAGCAGGAAACGGAAAATACACACAATCGC‑CAGAGCATTTAGAAACACCGGGGAAgg                                                                                                                                                                                       >  4:118355/1‑275 (MQ=255)

                                                                 tCAATGCGCCCATTACTGACAGCAAAGCCTTCCGGCACACCTGCCGGGCGCAACAGCCACCACGCCACGATAGCCGCCACCGCCAGTAACCCGACAACCCACCACGCCAGATGGCGCTTACTCTTATCCATAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAA‑CTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATAC‑GACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGc                                                                                                                                                                                      >  4:206500/1‑275 (MQ=255)

                                                                         cccATTACTGACAGCAAAGCCTTCCGGCACACCTGCCGGGCGCAACAGCCACCACGCCACGATAGCCGCCACCGCCAGTAACCCGACAACCCACCACGCCAGATGGCGCTTACTCTTATCCATAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAA‑CTGCCGCCGAACGACGTTCAGCGGATTCTTCCAACCGGGATAAGAAGGAAACGGCAAATAC‑GACCAGCGCTCAGAGCCTTGAGCGACAACGGGGCAGGCTTAg                                                                                                                                                                                  >  1:89427/1‑271 (MQ=255)

                                                                         cccATTACTGACAGCAAAGCCTTCCGGCACACCTGCCGGGCGCAACAGCCACCACGCCACGATAGCCGCCACCGCCAGTAACCCGACAACCCACCACGCCAGATGGCGCTTACTCTTATCCATAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAA‑CTGCCGCCGAACGACGTCCAGCGGATTCTTCCATCCGGGATAAGATGGAAACGGCAAATAC‑GCCCAGCCCGCTGAGCCTTACGCGCCAACGGGGAAGCCTTACTGaa                                                                                                                                                                              >  2:232131/1‑275 (MQ=255)

                                                                                      gCAAAGCCTTCCGGCACACCTGCCGGGCGCAACAGCCACCACGCCACGATAGCCGCCACCGCCAGTAACCCGACAACCCACCACGCCAGATGGCGCTTACTCTTATCCATAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAA‑CTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATAC‑GACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTaa                                                                                                                                                                  >  2:142629/1‑274 (MQ=255)

                                                                                       cAACGCCTCCCGCCACACCGCCCGGGCGAAAACGCCCCCACGCCACGATACCCGCCACCGCCAGTAACCCGACAACCCCCCACGCCAGATGGCGCTTACTCTTATCCATAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAA‑CTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATAC‑GACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAAtg                                                                                                                                                                <  2:476385/275‑1 (MQ=255)

                                                                                                                       cgCCCCCACGCCACCATGGCCCCCACCCCCAAACACCCGACAACCCACCCCCCCAGATGGCGCTTCCCCTTACCAAAAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAA‑CTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATAC‑GACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGa                                                                                                                                <  4:359166/274‑1 (MQ=255)

                                                                                                                       accccccaccccacgAGGCCCCCCACCGCCCGAAACCCCACAACCCACCCCGCCCGGTGTCCATTCCTCTTTTCCATAACCGCCCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAA‑CTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATAC‑GACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGa                                                                                                                                <  4:359775/270‑1 (MQ=255)

                                                                                                                         ccaccaCGCCACGATACCCGCCCCCCCCATTAACCCAACAACCCACCCCGCCAGATGGCGCTTACTCTTATCCATAACCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAA‑CTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATAC‑GACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGaaa                                                                                                                              <  2:24507/275‑1 (MQ=255)

                                                                                                                           accacGCCACGATAGCCGCCACCCCCAGTAACCCGACAACCCACCACGCCAGATGGCGCTTACTCTTATCCATAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAA‑CTGCCGCCGAACGACGTTCAGc                                                                                                                                                                                                                                                                    <  7:1114798/140‑1 (MQ=255)

                                                                                                                                gccacgATAGCCGCCACCGCAAGTACCCCGACAACCCACCAGGCCAAAGGGCGCTTACCCTAATCCATAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAA‑CTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATAC‑GACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCAt                                                                                                                       <  2:366074/275‑1 (MQ=255)

                                                                                                                                cccgcgATAGCCCCCACCGCCAGTAACCCGACAACCCACCACGCCAGATGGCGCTTACTCTTACCCAAAATCGCCCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAA‑CTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATAC‑GACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCAt                                                                                                                       <  2:143005/271‑1 (MQ=255)

                                                                                                                                     gacacccgccccccccaGTAACCCCACAACCCACCACGCCAGATGGCGCTTACTCTTATCCATAACCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAA‑CTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATAC‑GACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCa                                                                                                                  <  1:70421/270‑1 (MQ=255)

                                                                                                                                             cccccccccctAACCCCCCAACCCCCCACCCCGGTTGCCGTTAACCCTCATCCACAATCACCCCGCCATAATCCCTGTAAGCAAACGACACGTGTTACAA‑CTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATAC‑GACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGATAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGtgtg                                                                                                          <  4:222425/265‑1 (MQ=255)

                                                                                                                                             ccaccgcaaGTAACCCGACCACCCAACACGCCAGATTGCGCTTACTCTAACACATAACCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAA‑CTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATAC‑GACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGtgtg                                                                                                          <  2:258633/275‑1 (MQ=255)

                                                                                                                                                ccgccaGTAACCCGACAACCCACCACGCCAGATGGCGCTTACTCTTATCCATAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAA‑CTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATaaga                                                                                                                                                                                                                                             >  6:292904/1‑142 (MQ=255)

                                                                                                                                                ccgccaGTAACCCGACAACCCACCACGCCAGATGGCGCTTACTCTTATCCATAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAA‑CTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATaaga                                                                                                                                                                                                                                             <  5:292904/142‑1 (MQ=255)

                                                                                                                                                             gACAACCCACTACGCCAGATGGCGCTTACTCTTATCCATAATCGACCCGCCATAATCCGTGTAAGCAAACGACGCGTGTTACAA‑CTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATaaga                                                                                                                                                                                                                                             >  5:914964/1‑129 (MQ=255)

                                                                                                                                                             gACAACCCACCACGCCAGATGGCGCTTACTCTTATCCATAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAA‑CTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATaaga                                                                                                                                                                                                                                             <  6:914964/129‑1 (MQ=255)

                                                                                                                                                             gACAACCCACCACGCCAGATGGCGCTTACTCTTATCCATAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAA‑CTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATAC‑GACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCCTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACAAANNNACTCATAAGTTTAATCTATAAAACAggta                                                                                          >  3:155759/1‑273 (MQ=255)

                                                                                                                                                              aCAACCCACCACGCCAGATGGCGCTTACTCTTATCCATAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAA‑CTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATaaga                                                                                                                                                                                                                                             >  6:43995/1‑128 (MQ=255)

                                                                                                                                                              aCAACCCACCACGCCAAATGGCGCTTACTCTTATCCATAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAA‑CTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATaaga                                                                                                                                                                                                                                             <  5:43995/128‑1 (MQ=255)

                                                                                                                                                                aaCCCCCCAGGCCATATGGCGTTTACTTTTATCCATATTCCACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAA‑CTGCCGCCGAACGACGTTCAGCGGATTCTTACAGCCGGGATAAGAAGGAAACGGTAAATAC‑GACCAg                                                                                                                                                                                                                      <  5:158944/148‑1 (MQ=255)

                                                                                                                                                                       cacGCCAGATGGCGCTTACTCTTATCCATAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAA‑CTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATGC‑GACCAGCGCGCag                                                                                                                                                                                                               <  6:184268/148‑1 (MQ=255)

                                                                                                                                                                         cGCCAGATGGCGCTTACTCTTATCCATAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAA‑CTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATAC‑g                                                                                                                                                                                                                           >  6:730068/1‑134 (MQ=255)

                                                                                                                                                                             cgATTGCGTTTACTCTCACCAAGAACCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAA‑CTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATAC‑GACCAGCGCGCAGAGCATTGAGCGCCAACGGGGTAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTatcca                                                                          <  1:153135/274‑1 (MQ=255)

                                                                                                                                                                               aTGGCGCTTACTCTTATCCATAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAA‑CTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAAc                                                                                                                                                                                                                                      >  5:243218/1‑118 (MQ=255)

                                                                                                                                                                                     cTTACTCTTATCCATAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAA‑CTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATAC‑GACCAGCGCGCAGAGCATTGAGCGCCa                                                                                                                                                                                                 >  8:742378/1‑148 (MQ=255)

                                                                                                                                                                                           ctTATCCATAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAA‑CTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCg                                                                                                                                                                                                                                                     >  5:541147/1‑91 (MQ=255)

                                                                                                                                                                                           ctTACCCATAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAA‑CTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCg                                                                                                                                                                                                                                                     <  6:541147/91‑1 (MQ=255)

                                                                                                                                                                                             tATCCATAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAA‑CTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAgg                                                                                                                                                                                                                                          >  8:341904/1‑100 (MQ=255)

                                                                                                                                                                                             tATCCATAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAA‑CTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAgg                                                                                                                                                                                                                                          <  7:341904/100‑1 (MQ=255)

                                                                                                                                                                                             tATCCATAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAA‑CTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCGAATAC‑GACCAGCGCGCAGAGCATTGAGCGCCAACGGGGaa                                                                                                                                                                                         <  7:742378/148‑1 (MQ=255)

                                                                                                                                                                                             tATCCATAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAA‑CTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATAC‑GACCAg                                                                                                                                                                                                                      >  5:528762/1‑119 (MQ=255)

                                                                                                                                                                                             tATCCATAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAA‑CTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATAC‑GACCAGCGCGCAGAGCATTg                                                                                                                                                                                                        >  8:585850/1‑133 (MQ=255)

                                                                                                                                                                                             tATCCATAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAA‑CTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATAC‑GACCAGCGCGCAGAGCATTGAGCGCCAACGGGGaa                                                                                                                                                                                         >  8:1072944/1‑148 (MQ=255)

                                                                                                                                                                                                  aTAATCGACCCGCCATAATCCCTGTACGCAAACGACGCGTGTTACAA‑CTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCg                                                                                                                                                                                                                                                     >  6:771925/1‑84 (MQ=255)

                                                                                                                                                                                                  aTAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAA‑CTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGaa                                                                                                                                                                                                                                        >  8:13404/1‑97 (MQ=255)

                                                                                                                                                                                                  aTAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAA‑CTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGaa                                                                                                                                                                                                                                        <  7:13404/97‑1 (MQ=255)

                                                                                                                                                                                                  aTAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAA‑CTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATa                                                                                                                                                                                                                              <  6:1065970/107‑1 (MQ=255)

                                                                                                                                                                                                  aTAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAA‑CTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATa                                                                                                                                                                                                                              >  5:1065970/1‑107 (MQ=255)

                                                                                                                                                                                                  aTAACCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAA‑CTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCg                                                                                                                                                                                                                                                     <  5:771925/84‑1 (MQ=255)

                                                                                                                                                                                                   tAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAA‑CTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATAC‑GACCAGCGCGCAGAGCATTGAGCGCCAACggg                                                                                                                                                                                            >  5:597609/1‑139 (MQ=255)

                                                                                                                                                                                                      tCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAA‑CTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATaagaa                                                                                                                                                                                                                                            <  8:167985/89‑1 (MQ=255)

                                                                                                                                                                                                      tCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAA‑CTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATAC‑GACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGgtagt                                                                                                                                                                    >  2:554261/1‑160 (MQ=255)

                                                                                                                                                                                                         aCCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAA‑CTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATAC‑GACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGtt                                              <  2:222195/275‑1 (MQ=255)

                                                                                                                                                                                                          cccGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAA‑CTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATAC‑GACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAAtg                                                                                                                                                                <  4:554260/160‑1 (MQ=255)

                                                                                                                                                                                                           ccGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAA‑CTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATAC‑GACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGaaaa                                                                                                                                                                            >  6:363895/1‑147 (MQ=255)

                                                                                                                                                                                                           ccGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAA‑CTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATAC‑GACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCa                                            >  2:278005/1‑275 (MQ=255)

                                                                                                                                                                                                           ccGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAA‑CTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATAC‑GACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCa                                            <  2:206501/275‑1 (MQ=255)

                                                                                                                                                                                                            cGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAA‑CTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATAC‑GACCAGCGCGCAGAGCATTGAGCGCCAACGGGGa                                                                                                                                                                                          <  6:469882/132‑1 (MQ=255)

                                                                                                                                                                                                            cGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAA‑CTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATAC‑GACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCtt                                                                                                                                                                                    <  7:615963/138‑1 (MQ=255)

                                                                                                                                                                                                            cGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAA‑CTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATAC‑GACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCtt                                                                                                                                                                                    >  8:615963/1‑138 (MQ=255)

                                                                                                                                                                                                             gCCATAATCCCTGTAAGCAAACGACGCGTGTTACAA‑CTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATaag                                                                                                                                                                                                                                              >  7:368093/1‑80 (MQ=255)

                                                                                                                                                                                                             gCCATAATCCCTGTAAGCAAACGACGCGTGTTACAA‑CTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATaag                                                                                                                                                                                                                                              <  8:368093/80‑1 (MQ=255)

                                                                                                                                                                                                              ccATAATCCCTGTAAGCAAACGACGCGTGTTACAA‑CTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCggg                                                                                                                                                                                                                                                   <  6:1116363/74‑1 (MQ=255)

                                                                                                                                                                                                              ccATAATCCCTGTAAGCAAACGACGCGTGTTACAA‑CTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCggg                                                                                                                                                                                                                                                   >  5:1116363/1‑74 (MQ=255)

                                                                                                                                                                                                                aTAATCCCCTTAAGCAAACGACGCGTGTTACAA‑CTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATAC‑GACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGc                                       <  2:104816/275‑1 (MQ=255)

                                                                                                                                                                                                                    tCCCTGTAAGCAAACGACGCGTGTTACAA‑CTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATAC‑GACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTaa                                                                                                                                                                  <  7:1072944/148‑1 (MQ=255)

                                                                                                                                                                                                                      ccTGTAAGCAAACGACGCGTGTTACAA‑CTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATAC‑GACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACg                                                                                                                                                                          >  8:1085101/1‑138 (MQ=255)

                                                                                                                                                                                                                      ccTGTAAGCAAACGACGCGTGTTACAA‑CTGCCGCCGAACGACGATCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATAC‑GACCAGCGCGCAGAGc                                                                                                                                                                                                            >  8:702424/1‑104 (MQ=255)

                                                                                                                                                                                                                        tGTAAGCAAACGACGCGTGTTACAA‑CTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATAC‑GACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAg                                                                                                              >  4:5683/1‑196 (MQ=255)

                                                                                                                                                                                                                         gTAAGCAAACGACGCGTGTTACAA‑CTGCCGCCGAACGACGTTCAGCGGAttct                                                                                                                                                                                                                                                             <  5:743130/53‑1 (MQ=255)

                                                                                                                                                                                                                         gTAAGCAAACGACGCGTGTTACAA‑CTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATAC‑GACCa                                                                                                                                                                                                                       <  8:751556/90‑1 (MQ=255)

                                                                                                                                                                                                                         gTAAGCAAACGACGCGTGTTACAA‑CTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATAC‑GACCAGCGCGCAGAGCATTGAGCGCCAAc                                                                                                                                                                                               >  4:26965/1‑114 (MQ=255)

                                                                                                                                                                                                                         gTAAGCAAACGACGCGTGTTACAA‑CTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATAC‑GACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACa                                >  4:384898/1‑273 (MQ=255)

                                                                                                                                                                                                                             gcaaaacgcgcgGGTGAAAA‑CGGCCGCCGAACGCCGGTTAGCGGATTCCTCCATCCGGGATAAAAAGGAAACGGCAAATAC‑GAGCAGCGCGCAGGGCATTGAGCGCCAACGGGGAAGGCGTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCAt                          <  4:398029/267‑1 (MQ=255)

                                                                                                                                                                                                                             gCAAACGACGCGTGTTACAA‑CTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATAC‑GACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGtcat                                                                                                                                                         >  5:431386/1‑148 (MQ=255)

                                                                                                                                                                                                                                aaCGACGCGTGTTACAA‑CTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATAC‑GACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGtcatcat                                                                                                                                                      >  6:672510/1‑148 (MQ=255)

                                                                                                                                                                                                                                       cgTGTTACAA‑CTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATAC‑GACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCAtt                                                                                                                                               <  5:342301/148‑1 (MQ=255)

                                                                                                                                                                                                                                            tACAA‑CTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATAC‑GACCAgcgcgc                                                                                                                                                                                                                 >  7:20140/1‑77 (MQ=255)

                                                                                                                                                                                                                                            tACAA‑CTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATAC‑GACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAAtgttgt                                                                                                                                                            >  6:133160/1‑130 (MQ=255)

                                                                                                                                                                                                                                            tACAA‑CTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATAC‑GACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAAtgttgt                                                                                                                                                            <  5:133160/130‑1 (MQ=255)

                                                                                                                                                                                                                                             aCAA‑CTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATAC‑GACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGc                                                                                                                                         >  7:490062/1‑148 (MQ=255)

                                                                                                                                                                                                                                              cAA‑CTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATAC‑GACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGcc                                                                                                                                        <  8:490062/148‑1 (MQ=255)

                                                                                                                                                                                                                                                a‑CTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATAC‑GACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGtcatca                                                                                                                                                       >  6:244798/1‑131 (MQ=255)

                                                                                                                                                                                                                                                a‑CTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATAC‑GACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGt                                                                                                                                      <  5:672510/148‑1 (MQ=255)

                                                                                                                                                                                                                                                  cTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGAAAATAC‑GACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCaa      <  2:416906/275‑1 (MQ=255)

                                                                                                                                                                                                                                                    tccgctgAACGTCGTTCGGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATAC‑GACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGtcat                                                                                                                                                         >  5:265863/2‑126 (MQ=255)

                                                                                                                                                                                                                                                     ccgctgAACGTCGTTCGGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATAC‑GACCAGCGCGCag                                                                                                                                                                                                               >  7:778419/1‑71 (MQ=38)

                                                                                                                                                                                                                                                     ccgctgAACGTCGTTCGGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATAC‑GACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTa                                                                                                                                                                                   <  6:551073/99‑1 (MQ=255)

                                                                                                                                                                                                                                                     ccgctgAACGTCGTTCGGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATAC‑GACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAg                                                                                                                                                                                  >  5:129050/1‑100 (MQ=255)

                                                                                                                                                                                                                                                     ccgctgAACGTCGTTCGGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATAC‑GACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAg                                                                                                                                                                                  <  6:129050/100‑1 (MQ=255)

                                                                                                                                                                                                                                                     ccgctgAACGTCGTTCGGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATAC‑GACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGgt                                                                                                                                                                       <  5:164718/111‑1 (MQ=255)

                                                                                                                                                                                                                                                     ccgctgAACGTCGTTCGGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATAC‑GACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGgt                                                                                                                                                                       >  6:164718/1‑111 (MQ=255)

                                                                                                                                                                                                                                                      cgctgAACGTCGTTCGGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATAC‑GACCAGCGCGCAGAGCATTGAgcg                                                                                                                                                                                                    <  6:311755/81‑1 (MQ=255)

                                                                                                                                                                                                                                                       gctgAACGTCGTTCGGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATAC‑GACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGgtagt                                                                                                                                                                    <  7:608569/109‑1 (MQ=255)

                                                                                                                                                                                                                                                       gccgAACGACGTTCAGCGGATTCTTCCATCCGGGAGAAAAAGGAAACGCGAAAGAC‑GACCAGCGCGCAGAGCATTGAGCACCAACGGGGAAGGCTTAGTGGAAACGAGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGt  <  4:232130/274‑1 (MQ=255)

                                                                                                                                                                                                                                                        ccgAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATAC‑GACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAAtgttgt                                                                                                                                                            <  7:723362/119‑1 (MQ=255)

                                                                                                                                                                                                                                                         cgAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATAC‑GACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCtt                                                                                                                                                                                    <  7:817387/94‑1 (MQ=255)

                                                                                                                                                                                                                                                         cgAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATAC‑GACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCtt                                                                                                                                                                                    >  8:817387/1‑94 (MQ=255)

                                                                                                                                                                                                                                                               tcgTTCGGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATAC‑GACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGCCTTAg                                                                                                                                                                                  >  6:362915/2‑90 (MQ=255)

                                                                                                                                                                                                                                                               tcgTTCGGCCGATTCTTCCATCCGGGAGAAGAAGGAAACGGCAAATAC‑GACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAg                                                                                                                                                                                  <  5:362915/89‑1 (MQ=255)

                                                                                                                                                                                                                                                                cgTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATAC‑GACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGgtagta                                                                                                                                                                   <  8:222694/104‑1 (MQ=255)

                                                                                                                                                                                                                                                                cgTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATAC‑GACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGgtagta                                                                                                                                                                   >  7:222694/1‑104 (MQ=255)

CTTCTTTCACCAGAATGGTGTCGATACGCCCGGCAATTTTGCTGGCAATATCCACTTCCGTCGCTTCAATGCGCCCATTACTGACAGCAAAGCCTTCCGGCACACCTGCCGGGCGCAACAGCCACCACGCCACGATAGCCGCCACCGCCAGTAACCCGACAACCCACCACGCCAGATGGCGCTTACTCTTATCCATAATCGACCCGCCATAATCCCTGTAAGCAAACGACGCGTGTTACAA‑CTGCCGCCGAACGACGTTCAGCGGATTCTTCCATCCGGGATAAGAAGGAAACGGCAAATAC‑GACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGT  >  NC_000913/3630407‑3630926

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