BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,680,338	G→A	12.6%	intergenic (+42/+62)	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,338	0	G	A	12.6%	439.7 / 45.8	179	intergenic (+42/+62)	PP_4137/chrR	outer membrane siderophore receptor/chromate reductase
Reads supporting (aligned to +/- strand): ref base G (37/110); new base A (5/16); total (51/128)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.64e-01

Only 100 of 179 total aligned reads displayed.

CGACCGTGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTT‑A‑CCCG‑CGAATGCGTCATTGAATTTGCCGACGCATTCGC‑GGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCATTGAGCTTG  >  NC_002947/4680193‑4680484

cGACCGTGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTT‑A‑CCCG‑CGAATGCGTCATTGAATTTGCCGa                                                                                                                                                  <  2:414486/150‑1 (MQ=255)

 gACCGCGACGCCGGCCCACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTGGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTT‑A‑CCCG‑CGAATGCGTCATTGAATTTGCCGAc                                                                                                                                                 <  1:5627795/150‑1 (MQ=255)

  aCCGTGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGGTT‑T‑CCCC‑CGCATGCGGCCTTGGATTTGGCgccc                                                                                                                                                >  2:5068731/1‑147 (MQ=255)

     gTGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGGTT‑A‑ACCG‑CCGATTCGTGAATTAATTTGCCGCCGCAt                                                                                                                                             >  2:4987275/1‑150 (MQ=255)

       gACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTT‑A‑CCCG‑CGAATGCGTCATTGAATTTGCCGACGCATTc                                                                                                                                           <  1:6114672/150‑1 (MQ=255)

          gCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTT‑A‑CCCG‑CGAATGCGTCATTGAATTTGCCGACGCATTcgc‑g                                                                                                                                       <  2:935962/150‑1 (MQ=255)

           ccGGCCTACACCGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTT‑A‑CCCG‑CGAATGCGTCGGCAAATTCAATGACGCATTCGC‑gg                                                                                                                                      <  2:3934192/150‑1 (MQ=37)

           ccGGCCTACAACGACGCCGGCCAGGACCATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTT‑A‑CCCG‑CGAATGCGTCATTGAATTTGCCGACGCATTCGC‑gg                                                                                                                                      <  2:7384978/150‑1 (MQ=255)

           ccGGCCTACAACGACGCCGGCCAGGACAATTCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTT‑A‑CCCG‑CGAATGCGTCATTGAATTTGCCGACGCATTCGC‑gg                                                                                                                                      <  1:2955303/150‑1 (MQ=255)

           ccGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTT‑A‑CCCG‑CGAATGCGTCGGCAAATTCAATGACGCATTCGC‑gg                                                                                                                                      <  2:7225755/150‑1 (MQ=37)

           ccGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTT‑A‑CCCG‑CGAATGCGTCGGCAAATTCAATGACGCATTCGC‑gg                                                                                                                                      <  2:4637701/150‑1 (MQ=37)

           ccGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTT‑A‑CCCG‑CGAATGCGTCGGCAAATTCAATGACGCATTCGC‑gg                                                                                                                                      <  2:4637637/150‑1 (MQ=37)

           ccGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTT‑A‑CCCG‑CGAATGCGTCGGCAAATTCAATGACGCATTCGC‑gg                                                                                                                                      <  2:7335258/150‑1 (MQ=37)

           ccGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTT‑A‑CCCG‑CGAATGCGTCGGCAAATTCAATGACGCATTCGC‑gg                                                                                                                                      <  1:5627902/150‑1 (MQ=37)

           ccGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTT‑A‑CCCG‑CGAATGCGTCATTGAATTTGCCGACGCATTCGC‑gg                                                                                                                                      <  2:1328067/150‑1 (MQ=255)

           ccGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTT‑A‑CCCG‑CGAATGCGTCATTGAATTTGCCGACGCATTCGC‑gg                                                                                                                                      <  1:5333293/150‑1 (MQ=255)

           ccGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTT‑A‑CCCG‑CGAATGCGTCATTGAATTTGCCGACGCATTCGC‑gg                                                                                                                                      <  1:5937068/150‑1 (MQ=255)

           ccGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTT‑A‑CCCG‑CGAATGCGTCATTGAATTTGCCGACGCATTCGC‑gg                                                                                                                                      <  1:7009150/150‑1 (MQ=255)

           ccGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTGCCGCTGACGGGTGGGCGCGGGTTT‑A‑CCCG‑CGAATGCGTCGGCAAATTCAATGACGCATTCGC‑gg                                                                                                                                      <  2:3488774/150‑1 (MQ=37)

            cGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTT‑A‑CCCG‑CGAATGCGTCATTGAATTTGCCGACGCATTCGC‑ggg                                                                                                                                     <  1:2188705/150‑1 (MQ=255)

             ggCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTT‑A‑CCCG‑CGAATGCGTCATTGAATTTGCCGACGCATTCGC‑GGGt                                                                                                                                    <  2:136245/150‑1 (MQ=255)

              gCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTT‑A‑CCCG‑CGAATGCGTCATTGAATTTGCCGACGCATTCGC‑GGGTa                                                                                                                                   <  2:5869333/150‑1 (MQ=255)

               ccTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTT‑A‑CCCG‑CGAATGCGTCATTGAATTTGCCGACGCATTCGC‑GGGTaa                                                                                                                                  <  1:2028128/150‑1 (MQ=255)

                cTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTT‑A‑CCCG‑CGAATGCGTCATTGAATTTGCCGACGCATTCGC‑GGGTaaa                                                                                                                                 <  1:1940611/150‑1 (MQ=255)

                 tACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTT‑A‑CCCG‑CGAATGCGTCATTGAATTTGCCGACGCATTCGC‑GGGTAAA‑c                                                                                                                               <  1:5216451/150‑1 (MQ=255)

                 tACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTT‑A‑CCCG‑CGAATGCGTCATTGAATTTGCCGACGCATTCGC‑GGGTAAA‑c                                                                                                                               <  1:3639520/150‑1 (MQ=255)

                 tACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTT‑A‑CCCG‑CGAATGCGTCATTGAATTTGCCGACGCATTCGC‑GGGTAAA‑c                                                                                                                               <  2:32014/150‑1 (MQ=255)

                  aCAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTT‑T‑CCCG‑CGCAGGCGGCATTTAAATTTGCGACGCAATCGC‑GGGTaaaac                                                                                                                               >  2:4513022/1‑148 (MQ=255)

                  aCAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGGTT‑T‑CCCC‑CGCCTGCGTCCGTGGATTTTCCGCCGCCTTCGC‑GGGGTAA‑cc                                                                                                                              >  2:35886/1‑150 (MQ=255)

                    aaCGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTT‑A‑CCCG‑CGCATGCGGCATTGAATTTGGCGACGCAGTCGC‑GGGGCAA‑cccc                                                                                                                            >  2:3873266/1‑149 (MQ=255)

                     acgacgCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTT‑A‑CCCG‑CGAATGCGTCGGCAAATTCAATGACGCATTCGC‑GGGTAAA‑CCCGc                                                                                                                           <  2:1925325/150‑1 (MQ=25)

                     acgacgCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTT‑A‑CCCG‑CGAATGCGTCATTGAATTTGCCGACGCATTCGC‑GGGTaaaagcca                                                                                                                            >  2:1389598/1‑145 (MQ=255)

                           ccGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTAA‑ACCC‑CCAAT‑CGTCATTGTATTTTACGACCCAATAGCGGGGTAAA‑ccccaccgaa                                                                                                                      >  1:5442988/1‑143 (MQ=255)

                               ccAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTT‑A‑ACCG‑CGAATGCGTCATTTAATTTTCCGACGCATTCGC‑GGGTAAA‑CCCCCGCCAAAACAg                                                                                                                 >  2:2432403/1‑150 (MQ=39)

                                   gACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTT‑A‑ACCG‑CGCATGCGGCCTTGCATTTGGCGACGCATTCGC‑GGGTAAA‑CCCGCTCCTACACAggggg                                                                                                             >  2:7212382/1‑149 (MQ=38)

                                    aCAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACCGCTTCCGCTGACGGGTGGGCGCGGGTTT‑A‑CCCG‑CGAATGCGTCATTTATTTTGCCGACGCATTCGC‑GGGTAAA‑CCCGCTCCTACACAGGGGCg                                                                                                            <  2:2566229/150‑1 (MQ=255)

                                      acggCCGGTTCGACCCGGGGCGACGGATTGGGGGTGTTTGCGGGGGTTTCCTCCCGCTTCGGGGGCGGGGTGGGCGGGGTTTT‑A‑CCGG‑CGACTGCGTCATTGAATTTGCCGACGCATTCGG‑GGGTAAC‑CCCGCTCCTACACAGGGGCGtt                                                                                                          <  1:2327155/147‑1 (MQ=19)

                                      aaggCCCGTTCGACACCGGGCGATGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTT‑A‑CCCG‑CGAATGCGTCATTGAATTTGCCGACGCATTCGC‑GGGTAAA‑CCCGCTCCTACACAGGGGCGtt                                                                                                          <  2:3609922/147‑1 (MQ=255)

                                      aaggCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTT‑A‑CCCG‑CGAATGCGTCATTGAATTTGCCGACGCATTCGC‑GGGTAAA‑CCCGCTCCTACACAGGGGCGtt                                                                                                          <  2:6744737/147‑1 (MQ=255)

                                      aaTTCCCGTTCGACACGGGGGGAAGGATTTGGGGTGTATGCCGGGTTTTCCTACAGCTTCCGGTGACGGGTGGGGGCAGGTTG‑A‑CCCG‑CGAATGCGTCTTTGAATTTGCCGACGCATTCGC‑GGGTAAA‑CCCGCTCCTACACAGGGGCttt                                                                                                          <  2:952367/150‑1 (MQ=37)

                                              ttCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCCACAGCTTCCGCGGAGGGGTGGGCGGGGTTTT‑A‑CCCG‑CTAATGCGACATTGAATTTGCCGACGAATTCGC‑GGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCAcagca                                                                                                  <  1:2981155/150‑1 (MQ=38)

                                                       ggggGATGGTTTGGGCGTGTAGGCGGGGGTTCCCTCCCGCTGCCGCTGCGGGGGGGGCGGGGGTTC‑A‑CCCG‑CGACGGCGTCTTTGAATTGGCCGCCG‑ATTGGG‑GGGTAAACCCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAggg                                                                                         <  1:3537352/150‑1 (MQ=255)

                                                       gggCGATGTATTGGGCGTGTAGGCGGGGGTTTCCTACCGCTTCCGCTGACGGGGGGGCGCGGGTTT‑A‑CCCG‑CGAATGCGTCATTGAATTTGCCGACGCATTCGG‑GGTTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAggg                                                                                         <  1:7258350/150‑1 (MQ=37)

                                                       gggCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTT‑A‑CCCG‑CGAATGCGTCATTGAATTTGCCGACGCATTCGC‑GGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAggg                                                                                         <  2:3950683/150‑1 (MQ=255)

                                                          ctacggttttGGGGTGTTTGCGGGGTTTTCCTCCAGCTTCCGCTGCGGGGGGGGCG‑GGGTTT‑ACCCCG‑CTATTGCGTATTTGATTTTGCCGACGATTTCGG‑GGTAAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGtt                                                                                      <  1:5068731/143‑1 (MQ=255)

                                                           gAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTT‑T‑CCCG‑CGAATTCGTGATTGAATTTGCCGACGCCGTCGC‑GGGTAAA‑CCCGCGCCGAACCAGGGGCGTTGCCCCGGCGGGTTAGGGGttg                                                                                     >  2:447952/1‑149 (MQ=255)

                                                           gAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTT‑A‑ACCGCCAAAGGCGCCA‑TGCATTTGCCGACCCCATAGC‑CGGTAAA‑CCCCCTCCGCCACACGGGCGGTGCCCCGCACGGTTATGGGtta                                                                                     >  2:5063974/1‑149 (MQ=255)

                                                            aaGGATTTGGCGTGTAGGCCGGGGTTTCCTACAGCTCCCGCTGACGGGTGGGCGCGGGTTT‑A‑CCCG‑CGCATGCGTCATTGAATTTGCCGACGCATTCGC‑GGGTAAC‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCa                                                                                    <  2:4419906/150‑1 (MQ=255)

                                                               gATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTT‑A‑ACCG‑CGAATGCGGCATTGAATTTGGCGACGCATTAGC‑GGGTTAA‑ACCGCAACTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTccggc                                                                                 >  2:2367835/1‑146 (MQ=255)

                                                                aTTGGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTT‑A‑CCCG‑CGAATGCGTCATTGAATTTGCCGACGCATTCGC‑GGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGAcc                                                                                <  2:7185130/150‑1 (MQ=255)

                                                                   tGGCGTTTCTGCGGGGGTTTCCTCCCGCTCCCGCTGCGGGGGGGGCGCGGTTTT‑C‑CCCG‑CTTTGGCTTCTTTTATTTGGCCCACGAATTGGG‑GTTTAAC‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGcc                                                                             <  1:3994312/150‑1 (MQ=255)

                                                                   gggCGTGTATGCCGGGGTTACCTACAGCTTCCGCTGAGGGGTGGGGGCGGTTTT‑A‑CCCG‑CTAATGCTTATTTTATTTTCCCGACGCTTTCGC‑GGTTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGcc                                                                             <  1:6421372/149‑1 (MQ=255)

                                                                     gCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGCCGGGTGGGCGCGGGTTT‑A‑CCCG‑CGAATGCGTCATTGAATTTGCCGACGCATTCGC‑GGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCt                                                                           <  1:2406965/150‑1 (MQ=255)

                                                                       gtgtATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTT‑A‑CCCG‑CGATTGCTTCATTTAATTTGCCGACGCATTCGC‑GGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGt                                                                         <  2:2845233/150‑1 (MQ=255)

                                                                       gtgtATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTT‑A‑CCCG‑CGAATGCGTCATTGAATTTGCCGACGc                                                                                                                                               >  2:7106914/1‑82 (MQ=255)

                                                                         gtATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTT‑A‑CCCG‑CGAATGCGTCATTGAATTTGCCGACGCATTCGC‑GGGTAAA‑CCCGCTCCTACACAGGGGTGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTc                                                                       <  1:7422755/150‑1 (MQ=255)

                                                                            tGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTT‑A‑CCCG‑CGAATGCGTCATTGAATTTGCCGACGCATTCGC‑GGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGtttcgg                                                                    <  2:5566674/150‑5 (MQ=255)

                                                                               cGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTT‑A‑CCCG‑CGAATGCGTCGGCAAATTCAATGACGCATTCGC‑gg                                                                                                                                      >  1:1120291/1‑82 (MQ=14)

                                                                               cGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTT‑A‑CCCG‑CGAATGCGTCGGCAAATTCAATGACGCATTCGC‑gg                                                                                                                                      <  2:1120291/82‑1 (MQ=14)

                                                                               cGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTT‑A‑ACCC‑CGAAATCGTCATTTAATTTTCCGACGCAATAGC‑GGGTAAA‑ACCGCACCTACACAGGGGCGGTGCACACCAGGAATAGGGGTTCAGACCGCCCGGGTCAAcata                                                                 >  2:7317607/1‑147 (MQ=255)

                                                                                ggggTTTCCTCCTGCTTCCGGGGCCGGGTGGGCGGGGGTTT‑C‑CCCG‑GGATGGCGTCATTGAATTTGCCGACGCTTTGGG‑GGGTAAA‑CCCGCTCCGACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCa                                                                <  1:479015/150‑1 (MQ=255)

                                                                                     ttCCTACAGCTTCGGCTGTGGGGTGGGCGCGGGTTT‑A‑CCCG‑CGAATGCTTCTTTGATTTTGCCGACGATTTGGG‑GGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCAt                                                           <  1:5295272/150‑1 (MQ=255)

                                                                                        cTACAGCTTCCGCTGACGGGTGGGCGCGGGTTT‑A‑CCCG‑CGAATGCGTCATTGAATTTGCCGACGCATTCGC‑GGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGAc                                                        <  1:6737364/150‑1 (MQ=255)

                                                                                         tACTGCTTCCGGTGTGGGGTGGGGGCGGTTTT‑A‑CCCG‑GTAATGCTTAATTTAATTTGCCGACGCTTTGGC‑GGTTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACg                                                       <  1:6940421/150‑1 (MQ=255)

                                                                                             gCTTCCGCTGAGGGGTGGGCGCGGGTTT‑A‑CCCG‑CGAATGCGTCATTGAATTTGCCGACGCATTCGC‑GGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCaaa                                                   <  2:5591979/150‑1 (MQ=255)

                                                                                              cTTCCGCTGACGGGTGGGCGCGGGTTT‑A‑CCCG‑CGAATGCGTCATTGAATTTGCCGACGCATTCGC‑GGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAAc                                                  <  2:5492053/150‑1 (MQ=255)

                                                                                              cTTCCGCTGACGGGTGGGCGCGGGTTT‑A‑CCCG‑CGAATGCGTCATTGAATTTGCCGACGCATTCGC‑GGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAAc                                                  <  2:1720166/150‑1 (MQ=255)

                                                                                                tCCGCTGACGGGTGGGCGCGGGTTT‑A‑CCCG‑CGAATGCGTCATTGAATTTGCCGACGCATTCGC‑GGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCt                                                <  1:7361629/150‑1 (MQ=255)

                                                                                                 ccGCTGAGGGGTGGGCGCGGTTTT‑A‑CCGG‑CTATTGCTTCATTTAATTTGCCGACGCATTCGC‑GGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTg                                               <  1:5803546/150‑1 (MQ=255)

                                                                                                   gCTGACGGGTGGGCGCGGGTTT‑A‑CCCG‑CGAATGCGTCATTGAATTTGCCGACGCATTCGC‑GGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTc                                             <  2:6214972/150‑1 (MQ=255)

                                                                                                          ggTGGGCGCGGGTTT‑A‑CCCG‑CGACTGCGTCATTGAATTTGCCGACGCTTTCGC‑GGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAAc                                      <  1:2432403/150‑1 (MQ=255)

                                                                                                          ggTGGGCGCGGGTTT‑A‑CCCG‑CGAATGCGTCATTGAATTTGCCGACGCATTCGC‑GGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAAc                                      <  2:6705676/150‑1 (MQ=255)

                                                                                                             gggCGGGGGTTT‑C‑CCCG‑CTAATGGTTCTTTTATTTGGCCGACGATTTGGC‑GGTAAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGcc                                   <  1:331813/150‑1 (MQ=255)

                                                                                                              ggCGGGGGTTT‑C‑CCCG‑CGTATGCGTATTTTTTTTTGCCGACTAATCGGG‑GTGTAAC‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCt                                  <  1:2213157/150‑1 (MQ=255)

                                                                                                              ggCGCGTGTTT‑A‑CCCG‑CGAATGCGTCATTGAATTTGCCGACGCATTCGC‑GGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCt                                  <  2:5304047/150‑1 (MQ=255)

                                                                                                                 gcgGGTTT‑A‑CCGG‑CGAATGCGTCATTGAATTTGCCGACGCATTCGC‑GGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCa                               <  2:5621153/150‑1 (MQ=255)

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

                                                                                                                            cccG‑CGAATGCTTATTTTAATTTGCCGACGCTTTCGC‑GGTAAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTc                      <  1:1377852/150‑1 (MQ=255)

                                                                                                                             ccG‑CGAATGCGTCGGCAAATTCAATGACGCATTCGC‑GGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCAcc                                                              >  2:2517500/1‑109 (MQ=21)

                                                                                                                             ccG‑CGAATGCGTCGGCAAATTCAATGACGCATTCGC‑GGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCAcc                                                              <  1:2517657/109‑1 (MQ=21)

                                                                                                                             ccG‑CGAATGCGTCATTGAATTTGCCGACGCATTCGC‑GGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTcg                     >  2:2892444/1‑150 (MQ=255)

                                                                                                                             ccG‑CGAATGCGTCATTGAATTTGCCGACGCATTCGC‑GGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTcg                     >  2:3341215/1‑150 (MQ=255)

                                                                                                                                   aaTGCGTCATTGAATTTGCCGACGCATTCGC‑GGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCt                >  1:6029498/1‑150 (MQ=255)

                                                                                                                                     tGCTTCATTTAATTTGCCGACGATTTCGC‑GGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTg              <  1:2787081/150‑1 (MQ=255)

                                                                                                                                     tGCGTCATTGAATTTGCCGACGCATTCGC‑GGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTg              <  1:7280102/150‑1 (MQ=255)

                                                                                                                                       cGTCATTGATTTTGCCGACGCTTTGGG‑GGTAAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTc            <  1:4987275/150‑1 (MQ=255)

                                                                                                                                        gTCATTGAATTTGCCGACGCATTCGC‑GGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCa           >  1:4221745/1‑150 (MQ=255)

                                                                                                                                        gTCATTGAATTTGCCGACGCATTCGC‑GGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCa           >  2:2330014/1‑150 (MQ=255)

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

                                                                                                                                          cATTGAATTTGCCGACGCATTCGC‑GGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCAtt         >  2:209521/1‑150 (MQ=255)

                                                                                                                                          cATTGAATTTGCCGACGCATTCGC‑GGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCAtt         <  2:437106/150‑1 (MQ=255)

                                                                                                                                          cATTGAATTTGCCGACGCATTCGC‑GGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCAtt         >  2:212282/1‑150 (MQ=255)

                                                                                                                                           aTTGAATTTGCCGACGCATTCGC‑GGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCATTg        >  2:421366/1‑150 (MQ=255)

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

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

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

                                                                                                                                            ttGAATTTGCCGACGCATTCGC‑GGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCATTGa       >  2:3424817/1‑150 (MQ=255)

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

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

                                                                                                                                                aTTTGCCGACGCATTCGC‑GGGAAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCATTGAGCtt   <  2:374718/150‑1 (MQ=255)

CGACCGTGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTT‑A‑CCCG‑CGAATGCGTCATTGAATTTGCCGACGCATTCGC‑GGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCATTGAGCTTG  >  NC_002947/4680193‑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: `‑`