BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,680,329	A→G	13.0%	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,680,329	0	A	G	13.0%	424.5 / 44.6	182	intergenic (+33/+71)	PP_4137/chrR	outer membrane siderophore receptor/chromate reductase
Reads supporting (aligned to +/- strand): ref base A (39/102); new base G (5/16); total (48/134)
Fisher's exact test for biased strand distribution p-value = 7.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.37e-02

Only 100 of 185 total aligned reads displayed.

GCGCTATGCAGCGACCGTGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTA‑CCCG‑CGAATGCGTC‑ATTGAATTTGCCGACGCATTCGCGGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCATTG  >  NC_002947/4680182‑4680478

gcgcTATGCAGCGACCGTGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTA‑CCCG‑CGAATGCGTC‑Att                                                                                                                                                      <  2:4225177/150‑1 (MQ=255)

gcgcTATGCAGCGACCGTGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTA‑CCCG‑CGAATGCGTC‑Att                                                                                                                                                      <  1:6757816/150‑1 (MQ=255)

gcgcTATGCAGCGACCGTGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTA‑CCCG‑CGAATGCGTC‑Att                                                                                                                                                      <  2:1410411/150‑1 (MQ=255)

gcgcTATGCAGCGACCGTGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTA‑CCCG‑CGAATGCGTC‑Att                                                                                                                                                      <  1:1132366/150‑1 (MQ=255)

gcgcTATGCAGCGACCGTGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTA‑CCCG‑CGAATGCGTC‑Att                                                                                                                                                      <  1:1745767/150‑1 (MQ=255)

  gctcggCAGCCACCGTGCCGCCGCCCGACAACGACGCCGCCCAGGACACGGCCCGGGCGACCCCGGGCGACGGAATGGGCGTGTATGCCGGGGTGTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTA‑CCCG‑CGAATGCGTC‑ATTGa                                                                                                                                                    <  2:6592970/145‑1 (MQ=255)

  gcTATGCAGCGACCGTGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTA‑CCCG‑CGAATGCGTC‑ATTGa                                                                                                                                                    <  2:1331089/150‑1 (MQ=255)

   cTATGCAGCGACCGTGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTA‑CCCG‑CGAATGCGTC‑ATTGaa                                                                                                                                                   <  2:1146366/150‑1 (MQ=255)

    tATGCAGCGACCGTGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTA‑CCCG‑CGAATGCGTC‑ATTGAAt                                                                                                                                                  <  2:3961979/150‑1 (MQ=255)

     aTGCAGCGACCGTGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTA‑CCCG‑CGAATGCGTC‑ATTGAAtt                                                                                                                                                 <  2:3240686/150‑1 (MQ=255)

       gCAGCGACCGTGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGGTTT‑CCCG‑CGAATGCGGC‑ATTGAAtttt                                                                                                                                               >  2:5803546/1‑149 (MQ=255)

         aGCGACCGTGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGGTTTA‑ACC‑‑CGAAAGCGGCAATTGAAtttac                                                                                                                                              >  1:212075/1‑148 (MQ=255)

            gACCGCGACGCCGGCCCACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTGGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTA‑CCCG‑CGAATGCGTC‑ATTGAATTTGCCGAc                                                                                                                                          <  1:5627795/150‑1 (MQ=255)

                gTGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGGTTA‑ACCG‑CCGATTCGTG‑AATTAATTTGCCGCCGCAt                                                                                                                                      >  2:4987275/1‑150 (MQ=255)

                    cGCCGGCCTACACCGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTA‑CCCG‑CGAATGCGTC‑ATTGAATTTGCCGACGCATTcgc                                                                                                                                  <  2:4142955/150‑1 (MQ=255)

                    cGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTA‑CCCG‑CGAATGCGTC‑ATTGAATTTGCCGACGCATTcgc                                                                                                                                  <  2:3334653/150‑1 (MQ=255)

                     gCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTA‑CCCG‑CGAATGCGTC‑ATTGAATTTGCCGACGCATTcgcg                                                                                                                                 <  2:935962/150‑1 (MQ=255)

                      ccGGCCTACACCGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTA‑CCCG‑CGAATGCGTC‑GGCAAATTCAATGACGCATTCGCgg                                                                                                                                <  2:3934192/150‑1 (MQ=37)

                      ccGGCCTACAACGACGCCGGCCAGGACCATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTA‑CCCG‑CGAATGCGTC‑ATTGAATTTGCCGACGCATTCGCgg                                                                                                                                <  2:7384978/150‑1 (MQ=255)

                      ccGGCCTACAACGACGCCGGCCAGGACAATTCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTA‑CCCG‑CGAATGCGTC‑ATTGAATTTGCCGACGCATTCGCgg                                                                                                                                <  1:2955303/150‑1 (MQ=255)

                      ccGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTA‑CCCG‑CGAATGCGTC‑GGCAAATTCAATGACGCATTCGCgg                                                                                                                                <  2:4637637/150‑1 (MQ=37)

                      ccGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTA‑CCCG‑CGAATGCGTC‑GGCAAATTCAATGACGCATTCGCgg                                                                                                                                <  2:7335258/150‑1 (MQ=37)

                      ccGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTA‑CCCG‑CGAATGCGTC‑GGCAAATTCAATGACGCATTCGCgg                                                                                                                                <  2:7225851/150‑1 (MQ=37)

                      ccGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTA‑CCCG‑CGAATGCGTC‑GGCAAATTCAATGACGCATTCGCgg                                                                                                                                <  2:3934243/150‑1 (MQ=37)

                      ccGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTA‑CCCG‑CGAATGCGTC‑ATTGAATTTGCCGACGCATTCGCgg                                                                                                                                <  1:7009150/150‑1 (MQ=255)

                      ccGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTGCCGCTGACGGGTGGGCGC‑GGGTTTA‑CCCG‑CGAATGCGTC‑GGCAAATTCAATGACGCATTCGCgg                                                                                                                                <  2:3488774/150‑1 (MQ=37)

                       cGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTA‑CCCG‑CGAATGCGTC‑ATTGAATTTGCCGACGCATTCGCggg                                                                                                                               <  1:2188705/150‑1 (MQ=255)

                        ggCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTA‑CCCG‑CGAATGCGTC‑ATTGAATTTGCCGACGCATTCGCGGGt                                                                                                                              <  2:136245/150‑1 (MQ=255)

                         gCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTA‑CCCG‑CGAATGCGTC‑ATTGAATTTGCCGACGCATTCGCGGGTa                                                                                                                             <  2:5869333/150‑1 (MQ=255)

                           cTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTA‑CCCG‑CGAATGCGTC‑ATTGAATTTGCCGACGCATTCGCGGGTaaa                                                                                                                           <  1:1940611/150‑1 (MQ=255)

                           cTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGGTTT‑ACCG‑CGAAATCGTA‑AATGAATTTGGCGACCCAGTAGCGGGTaaa                                                                                                                           >  1:4888901/1‑150 (MQ=255)

                            tACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTA‑CCCG‑CGAATGCGTC‑ATTGAATTTGCCGACGCATTCGCGGGTAAA‑c                                                                                                                         <  1:5216451/150‑1 (MQ=255)

                            tACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTA‑CCCG‑CGAATGCGTC‑ATTGAATTTGCCGACGCATTCGCGGGTAAA‑c                                                                                                                         <  1:3639520/150‑1 (MQ=255)

                             aCAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTT‑CCCG‑CGCAGGCGGC‑ATTTAAATTTGCGACGCAATCGCGGGTaaaac                                                                                                                         >  2:4513022/1‑148 (MQ=255)

                             aCAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGGTTT‑CCCC‑CGCCTGCGTC‑CGTGGATTTTCCGCCGCCTTCGCGGGGTAA‑cc                                                                                                                        >  2:35886/1‑150 (MQ=255)

                                acgacgCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTA‑CCCG‑CGAATGCGTC‑GGCAAATTCAATGACGCATTCGCGGGTAAA‑CCCGc                                                                                                                     <  2:1925325/150‑1 (MQ=25)

                                  gacgCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTTTATGCCGGGGTTTCCTACAGCTGCCGCTGACGGGTGGGCGC‑GGTTTTA‑CCCG‑CGAATGCGTA‑ATTGAATTTGCCGACGATTTCGCGGTTAAA‑CCCGCTc                                                                                                                   <  2:4418511/150‑1 (MQ=255)

                                     gCCGGCCTGGACAATGCCGGTTCGACACCGGGCGATGGTTTTGGCGTGTTTGCGGGGGTTTCCTACAGCTTCGGCTGAGGGGTGGGCGG‑GGTTTTA‑CCCG‑CGAATGCGTC‑ATTGAATTTGCCGACGCTTTCGGGGTTAAA‑CCCGCTCCTa                                                                                                                <  2:5280214/150‑1 (MQ=255)

                                          ccAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTA‑ACCG‑CGAATGCGTC‑ATTTAATTTTCCGACGCATTCGCGGGTAAA‑CCCCCGCCAAAACAg                                                                                                           >  2:2432403/1‑150 (MQ=39)

                                              gACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTA‑ACCG‑CGCATGCGGC‑CTTGCATTTGGCGACGCATTCGCGGGTAAA‑CCCGCTCCTACACAggggg                                                                                                       >  2:7212382/1‑149 (MQ=38)

                                               aCAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACCGCTTCCGCTGACGGGTGGGCGC‑GGGTTTA‑CCCG‑CGAATGCGTC‑ATTTATTTTGCCGACGCATTCGCGGGTAAA‑CCCGCTCCTACACAGGGGCg                                                                                                      <  2:2566229/150‑1 (MQ=255)

                                                 aaggCCCGTTCGACACCGGGCGATGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTA‑CCCG‑CGAATGCGTC‑ATTGAATTTGCCGACGCATTCGCGGGTAAA‑CCCGCTCCTACACAGGGGCGtt                                                                                                    <  2:3609922/147‑1 (MQ=255)

                                                 aaTTCCCGTTCGACACGGGGGGAAGGATTTGGGGTGTATGCCGGGTTTTCCTACAGCTTCCGGTGACGGGTGGGGGC‑AGGTTGA‑CCCG‑CGAATGCGTC‑TTTGAATTTGCCGACGCATTCGCGGGTAAA‑CCCGCTCCTACACAGGGGCttt                                                                                                    <  2:952367/150‑1 (MQ=37)

                                                         ttCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTA‑CCCG‑CGAATGCGTC‑ATTGAATTTGCCGACGCATTCGCGGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCAcagca                                                                                            <  2:2961412/150‑1 (MQ=255)

                                                                ccGGGCTACGGTTTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGTTTTA‑CCCG‑CGATTGCGTC‑ATTGATTTTGCCGAGGCATTCGCGGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAg                                                                                     <  2:6212134/150‑1 (MQ=255)

                                                                 cGGGCGAAGGATTGGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTA‑CCCG‑CGAATGCGTC‑ATTGAATTTGCCGACGCATTCGCGGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAgg                                                                                    <  2:4477200/150‑1 (MQ=255)

                                                                  ggggGATGGTTTGGGCGTGTAGGCGGGGGTTCCCTCCCGCTGCCGCTGCGGGGGGGGCGG‑GGGTTCA‑CCCG‑CGACGGCGTC‑TTTGAATTGGCCGCCG‑ATTGGGGGGTAAACCCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAggg                                                                                   <  1:3537352/150‑1 (MQ=255)

                                                                  gggCGATGTATTGGGCGTGTAGGCGGGGGTTTCCTACCGCTTCCGCTGACGGGGGGGCGC‑GGGTTTA‑CCCG‑CGAATGCGTC‑ATTGAATTTGCCGACGCATTCGGGGTTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAggg                                                                                   <  1:7258350/150‑1 (MQ=37)

                                                                     ctacggttttGGGGTGTTTGCGGGGTTTTCCTCCAGCTTCCGCTGCGGGGGGGGCG‑‑GGGTTTACCCCG‑CTATTGCGTA‑TTTGATTTTGCCGACGATTTCGGGGTAAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGtt                                                                                <  1:5068731/143‑1 (MQ=255)

                                                                      gAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTT‑CCCG‑CGAATTCGTG‑ATTGAATTTGCCGACGCCGTCGCGGGTAAA‑CCCGCGCCGAACCAGGGGCGTTGCCCCGGCGGGTTAGGGGttg                                                                               >  2:447952/1‑149 (MQ=255)

                                                                      gAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTA‑ACCGCCAAAGGCGCC‑A‑TGCATTTGCCGACCCCATAGCCGGTAAA‑CCCCCTCCGCCACACGGGCGGTGCCCCGCACGGTTATGGGtta                                                                               >  2:5063974/1‑149 (MQ=255)

                                                                          gATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTA‑ACCG‑CGAATGCGGC‑ATTGAATTTGGCGACGCATTAGCGGGTTAA‑ACCGCAACTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTccggc                                                                           >  2:2367835/1‑146 (MQ=255)

                                                                           ttttGGCGTGTTTGCGGGGGTTTCCACCAGCTTCCGCTGAGGGGTGGGCGC‑GGGTTTA‑CCGG‑CTAATGCTTA‑ATTTAATTTGCCGACGATTTCGCGGTTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGAcc                                                                          <  1:5333217/149‑1 (MQ=255)

                                                                           aTTGGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTA‑CCCG‑CGAATGCGTC‑ATTGAATTTGCCGACGCATTCGCGGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGAcc                                                                          <  2:7185130/150‑1 (MQ=255)

                                                                             tgggCGTTTCTGCGGGGGTTTCCTACAGCTCCGGCTGCGGGGTGGGGGG‑GGGTTCA‑CCGG‑CTATGGCTTC‑TTTGATTTGGCCTCCGCATTGGGGGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGc                                                                        <  1:1472555/148‑1 (MQ=255)

                                                                              tGGCGTTTCTGCGGGGGTTTCCTCCCGCTCCCGCTGCGGGGGGGGCGC‑GGTTTTC‑CCCG‑CTTTGGCTTC‑TTTTATTTGGCCCACGAATTGGGGTTTAAC‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGcc                                                                       <  1:3994312/150‑1 (MQ=255)

                                                                              gggCGTGTATGCCGGGGTTACCTACAGCTTCCGCTGAGGGGTGGGGGC‑GGTTTTA‑CCCG‑CTAATGCTTA‑TTTTATTTTCCCGACGCTTTCGCGGTTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGcc                                                                       <  1:6421372/149‑1 (MQ=255)

                                                                                  gtgtATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTA‑CCCG‑CGATTGCTTC‑ATTTAATTTGCCGACGCATTCGCGGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGt                                                                   <  2:2845233/150‑1 (MQ=255)

                                                                                  gtgtATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTA‑CCCG‑CGAATGCGTC‑ATTGAATTTGCCGACGc                                                                                                                                        >  2:7106914/1‑82 (MQ=255)

                                                                                   tttATGCCGGGTTTACCAACAGATTCCGCTGACGGGTGGCCGC‑GGTATAA‑CCCG‑CGAATGCTTA‑ATTTAATTTCCCCACCAATTCGCGTTAAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGtt                                                                  <  1:2243376/148‑1 (MQ=255)

                                                                                    gttggcgggggTTTCCTACCGCTGCCGGTGCGGGGGGGGGGG‑GGGTTTC‑CCCG‑CGATGGCGTA‑TTTGAATTGGCCGCGGCATTGGGGGGTAAC‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTc                                                                 <  1:5782162/143‑1 (MQ=255)

                                                                                    gtATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTA‑CCCG‑CGAATGCGTC‑ATTGAATTTGCCGACGCATTCGCGGGTAAA‑CCCGCTCCTACACAGGGGTGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTc                                                                 <  1:7422755/150‑1 (MQ=255)

                                                                                     tATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTT‑CCCG‑CGAAATCGGG‑AGTTAATTTTGCGGCGCCTTAGCGGGGAAA‑ACCCCGGCGACACCGGGGGGGGGGGCCGCCGGGTTTGGGGGTGAGACCGCCCGGttga                                                                >  2:6296698/1‑148 (MQ=255)

                                                                                       tGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTA‑CCCG‑CGAATGCGTC‑ATTGAATTTGCCGACGCATTCGCGGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGtttcgg                                                              <  2:5566674/150‑5 (MQ=255)

                                                                                       tGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTA‑CCCG‑CGAATGCGTC‑ATTGAATTTGCCGACGCATTCGCGGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAAc                                                              <  2:7214037/150‑1 (MQ=255)

                                                                                        gCCGGGGTTTCCTACACCTACCGCTGACGGGTGGGCCC‑GCTTTTA‑CCCG‑CGAATCCGTC‑ATTTAATTTGCCGACGCATTCGCGGGAAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACt                                                             <  2:4692459/150‑1 (MQ=38)

                                                                                         ccGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTT‑ACCC‑CGAAAGCGGC‑CTTGAAATTTGCGAACCAGTAGCGGGTAAA‑ACCCCGCCAAAACACGGGCGGTTCACACCAGGCGTTGGGGGTAAGACCGCCCGGTTGaaaat                                                            >  2:3773788/1‑147 (MQ=255)

                                                                                          cGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTA‑CCCG‑CGAATGCGTC‑GGCAAATTCAATGACGCATTCGCgg                                                                                                                                >  1:1120291/1‑82 (MQ=14)

                                                                                          cGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTA‑CCCG‑CGAATGCGTC‑GGCAAATTCAATGACGCATTCGCgg                                                                                                                                <  2:1120291/82‑1 (MQ=14)

                                                                                             ggTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTA‑CCCG‑CGAATGCGTC‑ATTGAATTTGCCGACGCATTCGCGGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCAcc                                                        <  2:4149983/150‑1 (MQ=255)

                                                                                                    tACTGCTTCCGGTGTGGGGTGGGGGC‑GGTTTTA‑CCCG‑GTAATGCTTA‑ATTTAATTTGCCGACGCTTTGGCGGTTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACg                                                 <  1:6940421/150‑1 (MQ=255)

                                                                                                    tACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTA‑ACCG‑CGAATGCGTC‑ATTGAATTTGCCGACGCATTAGCGGGTAAA‑ACCCCTCCAACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCCCCCAGTTCAACATCACCCATGACg                                                 >  2:6737364/1‑150 (MQ=38)

                                                                                                        gCTTCCGCTGAGGGGTGGGCGC‑GGGTTTA‑CCCG‑CGAATGCGTC‑ATTGAATTTGCCGACGCATTCGCGGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCaaa                                             <  2:5591979/150‑1 (MQ=255)

                                                                                                         cTTCCGCTGACGGGTGGGCGC‑GGGTTTA‑CCCG‑CGAATGCGTC‑ATTGAATTTGCCGACGCATTCGCGGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAAc                                            <  2:1720166/150‑1 (MQ=255)

                                                                                                         cTTCCGCTGACGGGTGGGCGC‑GGGTTTA‑CCCG‑CGAATGCGTC‑ATTGAATTTGCCGACGCATTCGCGGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAAc                                            <  2:5492053/150‑1 (MQ=255)

                                                                                                           tCCGCTGACGGGTGGGCGC‑GGGTTTA‑CCCG‑CGAATGCGTC‑ATTGAATTTGCCGACGCATTCGCGGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCt                                          <  1:7361629/150‑1 (MQ=255)

                                                                                                              gCTGACGGGTGGGCGC‑GGGTTTA‑CCCG‑CGAATGCGTC‑ATTGAATTTGCCGACGCTTTCGCGGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTc                                       <  1:6072902/150‑1 (MQ=255)

                                                                                                                     ggTGGGCGC‑GGGTTTA‑CCCG‑CGACTGCGTC‑ATTGAATTTGCCGACGCTTTCGCGGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAAc                                <  1:2432403/150‑1 (MQ=255)

                                                                                                                        gggCGG‑GGGTTTC‑CCCG‑CTAATGGTTC‑TTTTATTTGGCCGACGATTTGGCGGTAAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGcc                             <  1:331813/150‑1 (MQ=255)

                                                                                                                         ggCGC‑GTGTTTA‑CCCG‑CGAATGCGTC‑ATTGAATTTGCCGACGCATTCGCGGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCt                            <  2:5304047/150‑1 (MQ=255)

                                                                                                                                  tttA‑CCCG‑CTATTGCGTC‑ATTTATTTTGCCGACGCTTTCGGGGTTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAAc                    <  1:737185/150‑1 (MQ=255)

                                                                                                                                  tttA‑CCCG‑CGAATGCGTC‑ATTGATTTTGCCGACGCATTCGCGGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAAc                    <  2:1376123/150‑1 (MQ=255)

                                                                                                                                        ccG‑CGAATGCGTC‑GGCAAATTCAATGACGCATTCGCGGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCAcc                                                        >  2:2517500/1‑109 (MQ=21)

                                                                                                                                        ccG‑CGAATGCGTC‑GGCAAATTCAATGACGCATTCGCGGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCAcc                                                        <  1:2517657/109‑1 (MQ=21)

                                                                                                                                        ccG‑CGAATGCGTC‑GGCAAATTCAATGACGCATTCGCGGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTcg               >  1:6725888/1‑150 (MQ=37)

                                                                                                                                        ccG‑CGAATGCGTC‑GGCAAATTCAATGACGCATTCGCGGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTcg               >  1:2981085/1‑150 (MQ=37)

                                                                                                                                        ccG‑CGAATGCGTC‑ATTGAATTTGCCGACGCATTCGCGGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTcg               >  1:5389858/1‑150 (MQ=255)

                                                                                                                                          g‑cgAATGCGTC‑ATTGAATTTGCCGACGCATTCGCGGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTcgcg             >  1:5514179/1‑150 (MQ=255)

                                                                                                                                                tGCTTC‑ATTTAATTTGCCGACGATTTCGCGGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTg        <  1:2787081/150‑1 (MQ=255)

                                                                                                                                                tGCGTC‑ATTGAATTTGCCGACGCATTCGCGGTTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTg        <  1:5848651/150‑1 (MQ=255)

                                                                                                                                                tGCGTC‑ATTGAATTTGCCGACGCATTCGCGGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTg        <  1:2367835/150‑1 (MQ=255)

                                                                                                                                                tGCGTC‑ATTGAATTTGCCGACGCATTCGCGGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTg        <  1:7280102/150‑1 (MQ=255)

                                                                                                                                                  cttC‑ATTTAATTCGCCGACGAATTCGCGGCTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTc      <  1:5942894/148‑1 (MQ=255)

                                                                                                                                                  cGTC‑ATTGATTTTGCCGACGCTTTGGGGGTAAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTc      <  1:4987275/150‑1 (MQ=255)

                                                                                                                                                  cGTC‑ATTGAATTTGCCGACGCATGCGGGGTTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTc      <  1:2192605/150‑1 (MQ=255)

                                                                                                                                                   gTC‑ATTGAATTTGCCGACGCATTCGCGGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCa     >  1:4221745/1‑150 (MQ=255)

                                                                                                                                                   gTC‑ATTGAATTTGCCGACGCATTCGCGGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCa     >  2:2330014/1‑150 (MQ=255)

                                                                                                                                                     c‑ATTGAATTTGCCGACGCATTCGCGGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCAtt   >  2:212282/1‑150 (MQ=255)

                                                                                                                                                     c‑ATTGAATTTGCCGACGCATTCGCGGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCAtt   >  2:209521/1‑150 (MQ=255)

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

GCGCTATGCAGCGACCGTGACGCCGGCCTACAACGACGCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGC‑GGGTTTA‑CCCG‑CGAATGCGTC‑ATTGAATTTGCCGACGCATTCGCGGGTAAA‑CCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCATTG  >  NC_002947/4680182‑4680478

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: `‑`