BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,677,302	A→G	10.1%	intergenic (+33/+71)	PP_4137 → / ← chrR	outer membrane siderophore receptor/chromate reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,677,302	0	A	G	10.1%	237.3 / 13.4	102	intergenic (+33/+71)	PP_4137/chrR	outer membrane siderophore receptor/chromate reductase
Reads supporting (aligned to +/- strand): ref base A (21/60); new base G (4/5); total (29/73)
Fisher's exact test for biased strand distribution p-value = 2.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.37e-01

Only 100 of 103 total aligned reads displayed.

CAGCGCTATGCAGCGACCGTGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTC  >  NC_002947/4677153‑4677442

caGCGCTATGCAGCGACCGTGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑GTTTCCCCCGCAAGCGTCa                                                                                                                                               >  2:3237869/1‑150 (MQ=255)

        tGCAGCGACCGTGACGCCGGCCTACAACGACGCCGGTCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGCATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAAttt                                                                                                                                      <  2:2827482/150‑1 (MQ=255)

        tGCAGCGACCGTGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAAttt                                                                                                                                      <  1:221941/150‑1 (MQ=255)

        tGCAGCGACCGTGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAAttt                                                                                                                                      <  1:1116646/150‑1 (MQ=255)

        tGCAGCGACCGTGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAAttt                                                                                                                                      <  1:265930/150‑1 (MQ=255)

        tGCAGCGACCGTGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAAttt                                                                                                                                      <  2:3213715/150‑1 (MQ=255)

          cAGCGACCGTGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGc                                                                                                                                    <  1:546905/150‑1 (MQ=255)

              gACCGTGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGAc                                                                                                                                <  1:4779669/150‑1 (MQ=255)

               aCCGTGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGCGGGG‑TTTACCC‑CGAATTCGTCAGT‑GAATTTTCCGAcc                                                                                                                               >  1:437016/1‑149 (MQ=255)

               aCCGTGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCCCGCATTCCGCATT‑TAAATTTCAGACg                                                                                                                               >  1:436750/1‑150 (MQ=255)

                ccGTGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGACGc                                                                                                                              <  2:1309216/150‑1 (MQ=255)

                  gTGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGACGCAt                                                                                                                            <  2:3665596/150‑1 (MQ=255)

                      cGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGACGCATTcgc                                                                                                                        <  1:3763219/150‑1 (MQ=255)

                      cGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGACGCATTcgc                                                                                                                        <  1:3763054/150‑1 (MQ=255)

                      cGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGACGCATTcgc                                                                                                                        <  1:4707577/150‑1 (MQ=255)

                      cGCCGGCCGACCCCGACGCCGGCCAGGACCCGGCCCGTTCGACCCCGGGCGAAGGATGTGGCGTGTAGGCCGGGGTTT‑CCTACCGCTTCCGCTGACGGGTGGGCGC‑GGG‑TGTACCCGCGCATGCGTCATT‑GAATGTGCCGACGCATTcgc                                                                                                                        <  2:1556516/150‑1 (MQ=255)

                       gCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCGGC‑AAATTCAATGACGCATTcgcg                                                                                                                       <  2:1367917/150‑1 (MQ=37)

                       gCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑Gaa                                                                                                                                         <  1:4680640/132‑1 (MQ=255)

                       gCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑Gaa                                                                                                                                         >  2:4680640/1‑132 (MQ=255)

                        ccGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCGGC‑AAATTCAATGACGCATTCGCgg                                                                                                                      <  2:2634065/150‑1 (MQ=37)

                        ccGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCGGC‑AAATTCAATGACGCATTCGCgg                                                                                                                      <  2:443606/150‑1 (MQ=37)

                        ccGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGACGCATTCGCgg                                                                                                                      <  2:3446300/150‑1 (MQ=255)

                        ccGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGACGCATTCGCgg                                                                                                                      <  2:2052829/150‑1 (MQ=255)

                        ccGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGACGCATTCGCgg                                                                                                                      <  2:4565245/150‑1 (MQ=255)

                        ccGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGACGCATTCGCgg                                                                                                                      <  2:5300967/150‑1 (MQ=255)

                        ccGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGAGC‑GGG‑TTTACCCGCGAATGCGTCGGC‑AAATTCAATGACGCATTCGCgg                                                                                                                      <  2:4215963/150‑1 (MQ=37)

                          ggCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGACGCATTCGCGGGt                                                                                                                    >  2:5341912/1‑150 (MQ=255)

                             cTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTAACCCCCAAAGCGTCAAT‑GAATTTTACGACGCATTAGCCGGGaaa                                                                                                                 >  1:4484483/1‑150 (MQ=255)

                              tACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGACGCATTCGCGGGTAAAc                                                                                                                <  1:1996649/150‑1 (MQ=255)

                               aCAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTTCCCGCGCATGCGTCATT‑GAATTTGCCGGCGCCTTAGCGGGTaaaac                                                                                                               >  2:1663163/1‑148 (MQ=255)

                                cAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGACGCATTCGCGGGTAAAccc                                                                                                              <  1:4063022/150‑1 (MQ=255)

                                 aaCGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGAc                                                                                                                                <  1:4728946/131‑1 (MQ=255)

                                 aaCGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGAc                                                                                                                                >  2:1989606/1‑131 (MQ=255)

                                 aaCGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGAc                                                                                                                                <  1:1989606/131‑1 (MQ=255)

                                 aaCGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGAc                                                                                                                                >  2:4728946/1‑131 (MQ=255)

                                   cgacgCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGACGCATTCGCGGGTAAACCCGCt                                                                                                           <  2:4166837/150‑1 (MQ=255)

                                    gccgCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACCGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCCACGCATTCGCGGGTAAACCCGCTc                                                                                                          <  2:4597795/148‑1 (MQ=255)

                                              aGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAggg                                                                                                <  2:1748011/150‑1 (MQ=255)

                                                   aaTGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGACGAATTCGCGGGTAAACCCGCTCCTACACAGGGGCGtt                                                                                           <  2:1642017/150‑1 (MQ=255)

                                                   aaTGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTGCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGtt                                                                                           <  2:3268618/150‑1 (MQ=255)

                                                    aTGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTg                                                                                          <  2:1770353/150‑1 (MQ=255)

                                                    aTGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTg                                                                                          <  2:1770174/150‑1 (MQ=255)

                                                     tGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTCCCGCTGACGGGTGGGCGC‑GGT‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGc                                                                                         <  1:4932882/150‑1 (MQ=255)

                                                      gCCCGTTCGACACCGGGCGAAGTTTTGGGGGTGTATGCGGGGGTTC‑CCTACAGCTTCCGGTGACGGGTGGGCGC‑GGT‑TTTACCCGCGAATGCGTATTT‑GAATTTGCCGACGAATTCGCGGTTAAACCCGCTCCTACACAGGGGCGTTGca                                                                                        <  2:1430412/150‑1 (MQ=255)

                                                      gCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGca                                                                                        <  2:4329158/150‑1 (MQ=255)

                                                        ccGTTCGACACCGGGCGATGGTTTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGcaca                                                                                      <  1:4953213/150‑1 (MQ=255)

                                                         cGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑GGTCCCCGCGCATGCGGCCTG‑GCATTTGCCCCCGCCGGCCCGGGGAAACCCGCGCCGACACCGGGGCGGTGCAcag                                                                                     >  2:3573696/1‑150 (MQ=255)

                                                          gTTCGACACCGGGCGATGGATTGGGCGTGTATGCCGGGGTTT‑CCTCCAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCAcagc                                                                                    <  2:3735245/150‑1 (MQ=255)

                                                           ttCGACACCGGGCTAAGTATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTCCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCTAATGCTTCATT‑TAATTGGCCCACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCAcagca                                                                                   <  2:622517/150‑1 (MQ=38)

                                                                cacCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTTACCCCCGAAT‑CGTCATT‑GAATTTTCCCACGCAGTCGCGGGGAAACCCCCGCCTACACAGGGGGGGTGGACAGCAGGGtt                                                                              >  1:3597797/1‑150 (MQ=255)

                                                                 acCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTAACCCCGAATTCGTCATT‑GTATTTTGCGACCCATTCGCCGGTAAAACCCCTCCAAAAAAGGGGCGGTTCCCAGCAGGCTTa                                                                             >  2:3405312/1‑150 (MQ=255)

                                                                 acCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑GTTACCCGCGACAGCGGCCAT‑GCAGTTTGCGCCGCAGTCGCGGGGGAACCCGCGCCGACACCGGGGCGGTGCCCAGCCggggtt                                                                             >  2:331027/1‑146 (MQ=255)

                                                                  ccGGGCGATGGATTTGGCGTTTATGCCGGGTTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGG‑GGT‑TTTACCCGCGAATGCTTCATT‑GAATTTGCCGACGCTTTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAg                                                                            <  1:801518/150‑1 (MQ=255)

                                                                    gggCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGCGGGG‑TTAACCCCCCAAT‑CGTCATT‑GAAATTTCCGACCCAATAGCGGGTAAAACCCCAACAACACACGGGCGGTGCACAGCCGGGATAggg                                                                          >  2:4205153/1‑150 (MQ=255)

                                                                    gggCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTAACCCCCAAAGCGGGAAT‑GTATTTTGCCACCCAATAGCGGGTAAAACCCCACCAAAACACGGGCGGTGCACAGCACGATTAggg                                                                          >  2:4205146/1‑150 (MQ=255)

                                                                        gacggATTGGGCGTGTCTGCCGGGGTGT‑CCTCCCGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTc                                                                      <  2:583681/147‑1 (MQ=38)

                                                                         ttgtttttGGCGTTTTTGCGGGGTTTT‑CCTACTGCTTCCGGTGCGGGGGGGGGGG‑GGT‑TTTACCCGCTAATGCTTCTTT‑TATTTTGCCGACGAATGGGGGGTTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCa                                                                     <  1:1119245/145‑1 (MQ=255)

                                                                              tttGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTTACCGCGAATTCGTCATT‑GAATTTTCCGACGCATTCGCGGGTAAACCCGCGCCTAAAAAGGGGCGTTGCACAGCAGGGTTAGGGGTTGAGACCg                                                                >  2:329190/1‑150 (MQ=255)

                                                                                tGGCGCGTAGGCCGGGGTTT‑CCAACCGCGCCCGCTGACGGGTGGGCGC‑GGG‑TGGCCCCGCGAATGCGACATG‑GAATTTGCCGACGCATGCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGcc                                                              <  2:5097005/150‑1 (MQ=255)

                                                                                 ggCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑GTTTCCCCCGAAAGCGGCCTTGGAATTTGCC‑ACGCAATAGCGGGTAAAACCCCGCCGACACAAGGGCGGTGCACCCCCGGATTAGGGGTTGAGAcccccc                                                             >  2:1187038/1‑150 (MQ=255)

                                                                                  gCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTAACCCCGAAAGCGTCATT‑GAATTTTGCGACCCATTAGCGGGTAAAACCGCACCTAAAAAAGGGCGTTGCACAGCAGGATTAGGGGTTCAGAcccccca                                                            >  2:3502255/1‑149 (MQ=255)

                                                                                  gCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑GTTAACCCCGAATTGGGGATT‑TAATTTTGCGACGCAATGGCGGGTTAACCCCCACCAAAACACGGGCGGTTCACAGCAGGATTATGGGTTAAGAACccccc                                                            >  2:3502333/1‑149 (MQ=255)

                                                                                    gtgtATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑GTTACCCGCGAATGGGGCCTT‑TGGTTTGGCGGCGCATTCGCGGGGTAACCCCCTCCGACACAGGGGCGGTGCACCGCCGGGTTAGGGGTTGAGACCCCCCGGt                                                          >  2:4798119/1‑150 (MQ=255)

                                                                                      gtATGCGGGGGTTT‑CCTACCGCTTCCGCTGTGGGGTGGGCGC‑GGG‑TTTACCCGCGATTGCTTATTT‑TATTTTGCCGACGCTTTGGGGGTTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTc                                                        <  1:4205153/150‑1 (MQ=255)

                                                                                      gtATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTc                                                        <  2:4972664/150‑1 (MQ=255)

                                                                                      gtATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTc                                                        <  2:435105/150‑1 (MQ=255)

                                                                                       tttggcGGGGTTTGCCTACCGCTGCGGCTGCGGGGGGGGCGC‑GGG‑TTTCCCCGGTATTGCTTATTT‑GATTTGGCCGACGATTTCGGGGTTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTc                                                        <  1:4205146/145‑1 (MQ=255)

                                                                                        ttGCGGGGTTTT‑CCTCCTGCTTCCGCTTCGGGGTGGGGGG‑GGT‑TTTCCCGGCTAGTGCTTCATT‑TATTTTGCCGCCGCTTTGGGGGTAAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCaa                                                      <  1:4674234/149‑1 (MQ=255)

                                                                                        tggcgggggTTT‑CCTACCGCTTCCGCTGAGGGGTGGGCGG‑GGT‑TTTCCCCGCTAATGCGTATTT‑TAATTGGCCGACGCATTCGGGGGAAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCaa                                                      <  1:3237869/145‑1 (MQ=255)

                                                                                        aTGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCaa                                                      <  1:4955390/150‑1 (MQ=255)

                                                                                         tGCCGGGGTTT‑CCTACAGCTGCCGCTGACGGGTGGGCGC‑GGG‑TTTACCGGCGAATGCGTCATT‑GAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAAc                                                     <  2:4656945/150‑1 (MQ=255)

                                                                                               ggTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑TAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCAcc                                               <  2:4707938/150‑1 (MQ=255)

                                                                                               ggTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCAcc                                               <  2:4354594/150‑1 (MQ=255)

                                                                                                  tt‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCAt                                            <  2:1590183/150‑1 (MQ=255)

                                                                                                     ccTACTGCTTCCGCTGACGGGTGGGCGG‑GGG‑TTTACCCGCTATTGCGTAATT‑TAATTTGCCGACGCTTTCGCGGTTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGa                                          <  1:3502255/150‑1 (MQ=255)

                                                                                                     ccTACAGCTCCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCTAATGCGTCATT‑TATTTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGa                                          <  1:3502333/150‑1 (MQ=255)

                                                                                                      atACTGCTTCGGGTGCGGGGTGGGGGG‑GGT‑TTTACCCGCTATTGCTTCTTT‑TATTTTGCCGCCGATTTCGGGGTTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGAc                                         <  1:189902/149‑1 (MQ=255)

                                                                                                            cTTCCGCTGAGGGGTGGGCGC‑GGG‑TTTACCCGCGATTGCGTCATT‑GAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAAc                                   <  2:2017390/150‑1 (MQ=255)

                                                                                                              tCCGCTGCGGGGGGGGCGC‑GGG‑TTGACCCGCGAAGGCGTCATT‑GAATTTGCCGACGCTTTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCt                                 <  1:1187038/150‑1 (MQ=255)

                                                                                                                cggggccGGGTGGGCGC‑GGG‑TTTACCCGCGATTGCGTCTTT‑GAATTTGCCGCCGCTTTGGGGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGt                               <  2:2274546/144‑1 (MQ=255)

                                                                                                                   ggACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGACGAATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGa                            <  2:1587192/149‑1 (MQ=255)

                                                                                                                            ggCGC‑GGG‑TTTACCCGCGAATGCGTAATT‑TATTTTGCCGACGCTTTCGCGGTTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCt                   <  1:3281242/150‑1 (MQ=255)

                                                                                                                             gcgc‑GGG‑TTTACCCGCTAATGCGTCATT‑GAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTg                  <  2:436750/150‑1 (MQ=255)

                                                                                                                             gcgc‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTg                  <  2:436768/150‑1 (MQ=255)

                                                                                                                             gcgc‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTg                  <  2:437016/150‑1 (MQ=255)

                                                                                                                               ggggggttttACCCGCGAATGCGTCATT‑GAATTTGCCGACGCATTCGCGGTTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTg                  <  2:2341525/143‑1 (MQ=255)

                                                                                                                               gc‑gGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCa                >  1:2443995/1‑150 (MQ=255)

                                                                                                                                    g‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGt                                                          >  1:1633197/1‑104 (MQ=255)

                                                                                                                                    g‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGt                                                          <  2:1633197/104‑1 (MQ=255)

                                                                                                                                    g‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGaa            <  2:1689060/150‑1 (MQ=255)

                                                                                                                                          cccGCGAATGCGTCATT‑GAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTc       <  2:2443995/150‑1 (MQ=255)

                                                                                                                                           ccGCGAATGCGTCGGC‑AAATTCAATGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTcg      >  1:2396191/1‑150 (MQ=37)

                                                                                                                                           ccGCGAATGCGTCGGC‑AAATTCAATGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTcg      >  2:237113/1‑150 (MQ=37)

                                                                                                                                           ccGCGAATGCGTCGGC‑AAATTCAATGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATCACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTcg      >  2:5112792/1‑150 (MQ=37)

                                                                                                                                           ccGCGAATGCGTCATT‑GAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTcg      >  2:2874174/1‑150 (MQ=255)

                                                                                                                                            cgcgAATGCGTCTTT‑GAATTGGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTcgc     <  1:4296739/150‑1 (MQ=255)

                                                                                                                                            cgcgAATGCGTCATT‑GAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTcgc     <  2:821310/150‑1 (MQ=255)

                                                                                                                                            cgcgAATGCGTCATT‑GAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTcgc     >  1:5266800/1‑150 (MQ=255)

                                                                                                                                             gcgACTGCGTCATT‑GAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACgggcatt    <  2:2145972/150‑6 (MQ=255)

                                                                                                                                               gAATGCGTCATT‑GAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTc  >  2:1947230/1‑150 (MQ=255)

CAGCGCTATGCAGCGACCGTGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTT‑CCTACAGCTTCCGCTGACGGGTGGGCGC‑GGG‑TTTACCCGCGAATGCGTCATT‑GAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTC  >  NC_002947/4677153‑4677442

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