BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,680,337	T→C	12.8%	intergenic (+41/+63)	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,680,337	0	T	C	12.8%	450.9 / 45.6	178	intergenic (+41/+63)	PP_4137/chrR	outer membrane siderophore receptor/chromate reductase
Reads supporting (aligned to +/- strand): ref base T (45/101); new base C (5/17); total (50/128)
Fisher's exact test for biased strand distribution p-value = 6.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.13e-01

Only 100 of 182 total aligned reads displayed.

GCAGCGACCGTGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTT‑A‑CCCGCGAATGCGTC‑A‑TTGAATTTGCCGACGCATTCGC‑GGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCATTGAGCTTG  >  NC_002947/4680189‑4680484

  aGCGACCGTGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGGTTT‑A‑ACC‑CGAAAGCGGCAA‑TTGAAtttac                                                                                                                                                     >  1:212075/1‑148 (MQ=255)

  aGCGACCGTGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTT‑T‑CCCGCGCAAGCGTC‑AGTTGAttttcc                                                                                                                                                     >  2:658200/1‑148 (MQ=255)

      aCCGTGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGGTT‑T‑CCCCCGCATGCGGC‑C‑TTGGATTTGGCgccc                                                                                                                                                >  2:5068731/1‑147 (MQ=255)

         gTGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTT‑A‑CCCGCGAATGCGTC‑A‑TTGAATTTGCCGACGCAt                                                                                                                                             <  1:516811/150‑1 (MQ=255)

           gACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTT‑A‑CCCGCGAATGCGTC‑A‑TTGAATTTGCCGACGCATTc                                                                                                                                           <  1:6114672/150‑1 (MQ=255)

             cGCCGGCCTACACCGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTT‑A‑CCCGCGAATGCGTC‑A‑TTGAATTTGCCGACGCATTcgc                                                                                                                                         <  2:4142955/150‑1 (MQ=255)

             cGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTT‑A‑CCCGCGAATGCGTC‑A‑TTGAATTTGCCGACGCATTcgc                                                                                                                                         <  2:3334653/150‑1 (MQ=255)

              gCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTT‑A‑CCCGCGAATGCGTC‑A‑TTGAATTTGCCGACGCATTcgc‑g                                                                                                                                       <  2:935962/150‑1 (MQ=255)

               cctgCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTT‑A‑CCCGCGAATGCGTC‑A‑TTGAATTTGCCGACGCATTCGC‑gg                                                                                                                                      <  1:3644874/147‑1 (MQ=255)

               ccGGCCTACACCGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTT‑A‑CCCGCGAATGCGTC‑G‑GCAAATTCAATGACGCATTCGC‑gg                                                                                                                                      <  2:3934192/150‑1 (MQ=37)

               ccGGCCTACAACGACGCCGGCCAGGACCATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTT‑A‑CCCGCGAATGCGTC‑A‑TTGAATTTGCCGACGCATTCGC‑gg                                                                                                                                      <  2:7384978/150‑1 (MQ=255)

               ccGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTT‑A‑CCCGCGAATGCGTC‑G‑GCAAATTCAATGACGCATTCGC‑gg                                                                                                                                      <  1:5627902/150‑1 (MQ=37)

               ccGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTT‑A‑CCCGCGAATGCGTC‑G‑GCAAATTCAATGACGCATTCGC‑gg                                                                                                                                      <  2:7225851/150‑1 (MQ=37)

               ccGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTT‑A‑CCCGCGAATGCGTC‑G‑GCAAATTCAATGACGCATTCGC‑gg                                                                                                                                      <  2:7225755/150‑1 (MQ=37)

               ccGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTT‑A‑CCCGCGAATGCGTC‑G‑GCAAATTCAATGACGCATTCGC‑gg                                                                                                                                      <  2:7335258/150‑1 (MQ=37)

               ccGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTT‑A‑CCCGCGAATGCGTC‑G‑GCAAATTCAATGACGCATTCGC‑gg                                                                                                                                      <  2:3934243/150‑1 (MQ=37)

               ccGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTT‑A‑CCCGCGAATGCGTC‑G‑GCAAATTCAATGACGCATTCGC‑gg                                                                                                                                      <  2:3971678/150‑1 (MQ=37)

               ccGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTT‑A‑CCCGCGAATGCGTC‑G‑GCAAATTCAATGACGCATTCGC‑gg                                                                                                                                      <  2:4637701/150‑1 (MQ=37)

               ccGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTT‑A‑CCCGCGAATGCGTC‑G‑GCAAATTCAATGACGCATTCGC‑gg                                                                                                                                      <  1:2171907/150‑1 (MQ=37)

               ccGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTT‑A‑CCCGCGAATGCGTC‑A‑TTGAATTTGCCGACGCATTCGC‑gg                                                                                                                                      <  2:1328067/150‑1 (MQ=255)

               ccGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTT‑A‑CCCGCGAATGCGTC‑A‑TTGAATTTGCCGACGCATTCGC‑gg                                                                                                                                      <  1:5937068/150‑1 (MQ=255)

               ccGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTT‑A‑CCCGCGAATGCGTC‑A‑TTGAATTTGCCGACGCATTCGC‑gg                                                                                                                                      <  1:5333293/150‑1 (MQ=255)

               ccGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTT‑A‑CCCGCGAATGCGTC‑A‑TTGAATTTGCCGACGCATTCGC‑gg                                                                                                                                      <  1:7009150/150‑1 (MQ=255)

               ccGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTGCCGCTGACGGGTGGGCGC‑GGGTTT‑A‑CCCGCGAATGCGTC‑G‑GCAAATTCAATGACGCATTCGC‑gg                                                                                                                                      <  2:3488774/150‑1 (MQ=37)

                 ggCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTT‑A‑CCCGCGAATGCGTC‑A‑TTGAATTTGCCGACGCATTCGC‑GGGt                                                                                                                                    <  2:136245/150‑1 (MQ=255)

                 ggCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTT‑A‑CCCGCGAATGCGTC‑A‑TTGAATTTGCCGACGCATTCGC‑GGGt                                                                                                                                    <  2:3574636/150‑1 (MQ=255)

                  gCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTT‑A‑CCCGCGAATGCGTC‑A‑TTGAATTTGCCGACGCATTCGC‑GGGTa                                                                                                                                   <  2:5869269/150‑1 (MQ=255)

                    cTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTT‑A‑CCCGCGAATGCGTC‑A‑TTGAATTTGCCGACGCATTCGC‑GGGTaaa                                                                                                                                 <  1:1940611/150‑1 (MQ=255)

                         acgacgCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTT‑A‑CCCGCGAATGCGTC‑A‑TTGAATTTGCCGACGCATTCGC‑GGGTaaaagcca                                                                                                                            >  2:1389598/1‑145 (MQ=255)

                           gacgCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGAGGGGTGGGCGC‑GGGTTT‑A‑CCCGCGAATGCGTC‑A‑TTTAATTTGCCGACGCATTCGC‑GGGTAAA‑CCCGCTc                                                                                                                         <  2:1031901/150‑1 (MQ=255)

                             cgCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTT‑A‑CCCGCGAATGCGTC‑A‑TTGAATTTGCCGACGCATTCGC‑GGGTAAA‑CCCGCTCCt                                                                                                                       <  2:1760007/150‑1 (MQ=255)

                               ccGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTAA‑ACCCCCAAT‑CGTC‑A‑TTGTATTTTACGACCCAATAGCGGGGTAAA‑ccccaccgaa                                                                                                                      >  1:5442988/1‑143 (MQ=255)

                                       gACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTT‑A‑ACCGCGCATGCGGC‑C‑TTGCATTTGGCGACGCATTCGC‑GGGTAAA‑CCCGCTCCTACACAggggg                                                                                                             >  2:7212382/1‑149 (MQ=38)

                                        aCAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACCGCTTCCGCTGACGGGTGGGCGC‑GGGTTT‑A‑CCCGCGAATGCGTC‑A‑TTTATTTTGCCGACGCATTCGC‑GGGTAAA‑CCCGCTCCTACACAGGGGCg                                                                                                            <  2:2566229/150‑1 (MQ=255)

                                          aaTTCCCGTTCGACACGGGGGGAAGGATTTGGGGTGTATGCCGGGTTTTCCTACAGCTTCCGGTGACGGGTGGGGGC‑AGGTTG‑A‑CCCGCGAATGCGTC‑T‑TTGAATTTGCCGACGCATTCGC‑GGGTAAA‑CCCGCTCCTACACAGGGGCttt                                                                                                          <  2:952367/150‑1 (MQ=37)

                                                cGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTT‑A‑CCCGCGAATGCGTC‑A‑TTGAATTTGCCGACGCATTcg                                                                                                                                          <  2:557622/114‑1 (MQ=255)

                                                cGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTT‑A‑CCCGCGAATGCGTC‑A‑TTGAATTTGCCGACGCATTcg                                                                                                                                          >  1:557622/1‑114 (MQ=255)

                                                  ttCGACACCGGGCGAAGGATTTGGCGTTTTTGCCGGGGTTCCCTACCGCTTCCGCTGACGGGTGGGCGC‑GGGTTT‑A‑CCCGCGAATGCGTC‑A‑TTTAATTTGCCGCCGCATTCGC‑GGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCAcagca                                                                                                  <  2:1999740/150‑1 (MQ=38)

                                                  ttCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTT‑A‑CCCGCGAATGCGTC‑A‑TTGAATTTGCCGACGCATTCGC‑GGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCAcagca                                                                                                  <  2:2961412/150‑1 (MQ=255)

                                                         ccGGGCTACGGTTTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGTTTT‑A‑CCCGCGATTGCGTC‑A‑TTGATTTTGCCGAGGCATTCGC‑GGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAg                                                                                           <  2:6212134/150‑1 (MQ=255)

                                                          cGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGGTT‑A‑CCCCCCAAAGCGGC‑A‑ATTAATTTGCCGACGCATTAGC‑GGGGTAA‑ACCGCACCAACACACGGGCGGTGCACCGCAGGATTAgg                                                                                          >  2:4511288/1‑150 (MQ=255)

                                                          cGGGCGAAGGATTGGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTT‑A‑CCCGCGAATGCGTC‑A‑TTGAATTTGCCGACGCATTCGC‑GGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAgg                                                                                          <  2:4477200/150‑1 (MQ=255)

                                                           ggggGATGGTTTGGGCGTGTAGGCGGGGGTTCCCTCCCGCTGCCGCTGCGGGGGGGGCGG‑GGGTTC‑A‑CCCGCGACGGCGTC‑T‑TTGAATTGGCCGCCG‑ATTGGG‑GGGTAAACCCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAggg                                                                                         <  1:3537352/150‑1 (MQ=255)

                                                           gggCGATGTATTGGGCGTGTAGGCGGGGGTTTCCTACCGCTTCCGCTGACGGGGGGGCGC‑GGGTTT‑A‑CCCGCGAATGCGTC‑A‑TTGAATTTGCCGACGCATTCGG‑GGTTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAggg                                                                                         <  1:7258350/150‑1 (MQ=37)

                                                              ctacggttttGGGGTGTTTGCGGGGTTTTCCTCCAGCTTCCGCTGCGGGGGGGGCG‑‑GGGTTT‑ACCCCGCTATTGCGTA‑T‑TTGATTTTGCCGACGATTTCGG‑GGTAAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGtt                                                                                      <  1:5068731/143‑1 (MQ=255)

                                                               gAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTT‑A‑ACCGCGAATTCGTC‑A‑TTGAATTTGGCGACGCATTCGC‑GGGTAAA‑CCCGCGCCGACACAGGGGCGTTGCACAGCCGGCTTAGGGGTTc                                                                                     >  2:2367909/1‑150 (MQ=38)

                                                                aaGGATTTGGCGTGTAGGCCGGGGTTTCCTACAGCTCCCGCTGACGGGTGGGCGC‑GGGTTT‑A‑CCCGCGCATGCGTC‑A‑TTGAATTTGCCGACGCATTCGC‑GGGTAAC‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCa                                                                                    <  2:4419906/150‑1 (MQ=255)

                                                                 ttcgtaaggCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGGTTT‑A‑CCC‑CGAAATCGTC‑A‑ATTAAATTTCCGGCGCAATAGC‑GGGTAAA‑ACCCCACCAAAACACGGGCGTTTCACAGCAGGATTAGGGGTTCAg                                                                                   >  2:3598080/8‑150 (MQ=255)

                                                                    aTTGGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTT‑A‑CCCGCGAATGCGTC‑A‑TTGAATTTGCCGACGCATTCGC‑GGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGAcc                                                                                <  2:7185130/150‑1 (MQ=255)

                                                                      ttGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTT‑A‑CCCGCGAATGCGTC‑A‑TTGAATTTGCCGACGCATTCGC‑GGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGc                                                                              <  2:2443603/150‑1 (MQ=255)

                                                                       gggCGTGTATGCCGGGGTTACCTACAGCTTCCGCTGAGGGGTGGGGGC‑GGTTTT‑A‑CCCGCTAATGCTTA‑T‑TTTATTTTCCCGACGCTTTCGC‑GGTTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGcc                                                                             <  1:6421372/149‑1 (MQ=255)

                                                                           gtgtATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTT‑A‑CCCGCGATTGCTTC‑A‑TTTAATTTGCCGACGCATTCGC‑GGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGt                                                                         <  2:2845233/150‑1 (MQ=255)

                                                                           gtgtATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTT‑A‑CCCGCGAATGCGTC‑A‑TTGAATTTGCCGACGc                                                                                                                                               >  2:7106914/1‑82 (MQ=255)

                                                                            tttATGCCGGGTTTACCAACAGATTCCGCTGACGGGTGGCCGC‑GGTATA‑A‑CCCGCGAATGCTTA‑A‑TTTAATTTCCCCACCAATTCGC‑GTTAAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGtt                                                                        <  1:2243376/148‑1 (MQ=255)

                                                                             gtATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTT‑A‑CCCGCGAATGCGTC‑A‑TTGAATTTGCCGACGCATTCGC‑GGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTc                                                                       <  2:1126246/150‑1 (MQ=255)

                                                                              tATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTT‑T‑CCCGCGAAATCGGG‑A‑GTTAATTTTGCGGCGCCTTAGC‑GGGGAAA‑ACCCCGGCGACACCGGGGGGGGGGGCCGCCGGGTTTGGGGGTGAGACCGCCCGGttga                                                                      >  2:6296698/1‑148 (MQ=255)

                                                                                tGCCGGGGTTTCCTACCGCTTCCGCTGACGGGTGGGCGC‑GGGTTT‑A‑CCCGCGAATGCGTC‑A‑TTGAATTTGCCGACGCATTCGC‑GGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAAc                                                                    <  2:4519889/150‑1 (MQ=255)

                                                                                tGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTT‑A‑CCCGCGAATGCGTC‑A‑TTGAATTTGCCGACGCATTCGC‑GGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGtttcgg                                                                    <  2:5566674/150‑5 (MQ=255)

                                                                                  ccGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTT‑T‑ACCCCGAAAGCGGC‑C‑TTGAAATTTGCGAACCAGTAGC‑GGGTAAA‑ACCCCGCCAAAACACGGGCGGTTCACACCAGGCGTTGGGGGTAAGACCGCCCGGTTGaaaat                                                                  >  2:3773788/1‑147 (MQ=255)

                                                                                   cGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTT‑A‑CCCGCGAATGCGTC‑G‑GCAAATTCAATGACGCATTCGC‑gg                                                                                                                                      >  1:1120291/1‑82 (MQ=14)

                                                                                    ggggTTTCCTCCTGCTTCCGGGGCCGGGTGGGCGG‑GGGTTT‑C‑CCCGGGATGGCGTC‑A‑TTGAATTTGCCGACGCTTTGGG‑GGGTAAA‑CCCGCTCCGACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCa                                                                <  1:479015/150‑1 (MQ=255)

                                                                                      ggTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTT‑A‑CCCGCGAATGCGTC‑A‑TTGACTTTGCCGACGCATTCGC‑GGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCAcc                                                              <  2:4930318/150‑1 (MQ=255)

                                                                                      ggTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTT‑A‑CCCGCGAATGCGTC‑A‑TTGAATTTGCCGACGCATTCGC‑GGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCAcc                                                              <  2:4149983/150‑1 (MQ=255)

                                                                                         ttCCTACAGCTTCGGCTGTGGGGTGGGCGC‑GGGTTT‑A‑CCCGCGAATGCTTC‑T‑TTGATTTTGCCGACGATTTGGG‑GGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCAt                                                           <  1:5295272/150‑1 (MQ=255)

                                                                                            cTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTT‑A‑CCCGCGAATGCGTC‑A‑TTGAATTTGCCGACGCATTCGC‑GGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGAc                                                        <  1:6737364/150‑1 (MQ=255)

                                                                                             tACTGCTTCCGGTGTGGGGTGGGGGC‑GGTTTT‑A‑CCCGGTAATGCTTA‑A‑TTTAATTTGCCGACGCTTTGGC‑GGTTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACg                                                       <  1:6940421/150‑1 (MQ=255)

                                                                                             tACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTT‑A‑ACCGCGAATGCGTC‑A‑TTGAATTTGCCGACGCATTAGC‑GGGTAAA‑ACCCCTCCAACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCCCCCAGTTCAACATCACCCATGACg                                                       >  2:6737364/1‑150 (MQ=38)

                                                                                                  cTTCGGCTGAGGGGTGGGCGG‑GGGTTT‑C‑CCGGCTATGGCTTC‑T‑TTTATTTGGCCGACGCTTTGGG‑GGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAAc                                                  <  1:2726974/150‑1 (MQ=255)

                                                                                                  cTTCCGCTGACGGGTGGGCGC‑GGGTTT‑A‑CCCGCGAATGCGTC‑A‑TTGAATTTGCCGACGCATTCGC‑GGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAAc                                                  <  2:5492053/150‑1 (MQ=255)

                                                                                                    tCCGCTGACGGGTGGGCGC‑GGGTTT‑A‑CCCGCGAATGCGTC‑A‑TTGAATTTGCCGACGCATTCGC‑GGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCt                                                <  1:7361629/150‑1 (MQ=255)

                                                                                                     ccGCTGAGGGGTGGGCGC‑GGTTTT‑A‑CCGGCTATTGCTTC‑A‑TTTAATTTGCCGACGCATTCGC‑GGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTg                                               <  1:5803546/150‑1 (MQ=255)

                                                                                                              ggggggCGG‑GGTTTC‑A‑CCCGCTAATGGGTC‑T‑TTTAATTGGCCGCCGTTTTGGG‑GGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAAc                                      <  1:35886/147‑1 (MQ=255)

                                                                                                              ggTGGGCGC‑GGGTTT‑A‑CCCGCGACTGCGTC‑A‑TTGAATTTGCCGACGCTTTCGC‑GGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAAc                                      <  1:2432403/150‑1 (MQ=255)

                                                                                                                 gggCGG‑GGGTTT‑C‑CCCGCTAATGGTTC‑T‑TTTATTTGGCCGACGATTTGGC‑GGTAAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGcc                                   <  1:331813/150‑1 (MQ=255)

                                                                                                                  ggCGC‑GTGTTT‑A‑CCCGCGAATGCGTC‑A‑TTGAATTTGCCGACGCATTCGC‑GGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCt                                  <  2:5304047/150‑1 (MQ=255)

                                                                                                                     gc‑gGGTTT‑A‑CCGGCGAATGCGTC‑A‑TTGAATTTGCCGACGCATTCGC‑GGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCa                               <  2:5621153/150‑1 (MQ=255)

                                                                                                                          gTTT‑A‑CCCGCGAATGCGTC‑A‑TTGAATTTGCCGACGCATTCGC‑GGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGaa                           >  1:7377177/1‑150 (MQ=255)

                                                                                                                                 cccGCGAATGCTTA‑T‑TTTAATTTGCCGACGCTTTCGC‑GGTAAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTc                      <  1:1377852/150‑1 (MQ=255)

                                                                                                                                  ccGCGAATGCGTC‑G‑GCAAATTCAATGACGCATTCGC‑GGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCAcc                                                              >  2:2517657/1‑109 (MQ=21)

                                                                                                                                  ccGCGAATGCGTC‑G‑GCAAATTCAATGACGCATTCGC‑GGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCAcc                                                              <  1:2517657/109‑1 (MQ=21)

                                                                                                                                  ccGCGAATGCGTC‑G‑GCAAATTCAATGACGCATTCGC‑GGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTcg                     >  1:2981085/1‑150 (MQ=37)

                                                                                                                                  ccGCGAATGCGTC‑A‑TTGAATTTGCCGACGCATTCGC‑GGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTcg                     >  2:500970/1‑150 (MQ=255)

                                                                                                                                       aaTGCGTC‑A‑TTGAATTTGCCGACGCATTCGC‑GGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCt                >  1:6029498/1‑150 (MQ=255)

                                                                                                                                        aTGCGTC‑A‑TTGAATTTGCCGACGCATTCGC‑GGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCtt               <  2:3639385/150‑1 (MQ=255)

                                                                                                                                         tGCGTC‑A‑TTGAATTTGCCGACGCATTCGC‑GGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTg              <  1:7280102/150‑1 (MQ=255)

                                                                                                                                         tGCGTC‑A‑TTGAATTTGCCGACGCATTCGC‑GGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTg              <  1:2367835/150‑1 (MQ=255)

                                                                                                                                           cGTC‑A‑TTGATTTTGCCGACGCTTTGGG‑GGTAAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTc            <  1:4987275/150‑1 (MQ=255)

                                                                                                                                           cGTC‑A‑TTGAATTTGCCGACGCATGCGG‑GGTTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTc            <  1:2192605/150‑1 (MQ=255)

                                                                                                                                            gTC‑A‑TTGAATTTGCCGACGCATTCGC‑GGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCa           >  1:2900942/1‑150 (MQ=255)

                                                                                                                                            gTC‑A‑TTGAATTTGCCGACGCATTCGC‑GGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCa           >  2:2330014/1‑150 (MQ=255)

                                                                                                                                              taa‑TTGAATTTGCCGACGCATTCGC‑GGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCAt          <  1:1659559/148‑1 (MQ=255)

                                                                                                                                              c‑A‑TTGAATTTGCCGACGCATTCGC‑GGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCAtt         >  2:210178/1‑150 (MQ=255)

                                                                                                                                              c‑A‑TTGAATTTGCCGACGCATTCGC‑GGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCAtt         <  2:437106/150‑1 (MQ=255)

                                                                                                                                              c‑A‑TTGAATTTGCCGACGCATTCGC‑GGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCAtt         >  2:212282/1‑150 (MQ=255)

                                                                                                                                                ctttGAATTTGCCGACGCTTTCGC‑GGTTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCAtt         <  1:4207697/148‑1 (MQ=255)

                                                                                                                                                a‑TTGAATTTGCCGACGCATTCGC‑GGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTAGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCATTg        >  1:5516308/1‑150 (MQ=255)

                                                                                                                                                  ttGAATTTGCCGACGCATTCGC‑GGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCATTGa       <  2:4888901/150‑1 (MQ=255)

                                                                                                                                                   ttaaTTGGCCGACGCATTCGC‑GGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCATTGAg      <  1:7135055/148‑1 (MQ=255)

                                                                                                                                                      aTTTGCCGACGCATTCGC‑GGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCATTGAGCtt   >  1:1157106/1‑150 (MQ=255)

                                                                                                                                                       tttGCCGACGCATTCGC‑GGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCATTGAGCTTg  >  2:3342972/1‑150 (MQ=255)

GCAGCGACCGTGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTT‑A‑CCCGCGAATGCGTC‑A‑TTGAATTTGCCGACGCATTCGC‑GGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCATTGAGCTTG  >  NC_002947/4680189‑4680484

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