BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,820,822	G→C	5.2%	A97G (GCA→GGA) ‡	PP_3371 ←	sensor histidine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,820,822	0	G	C	5.2%	1255.2 / 40.4	368	A97G (GCA→GGA) ‡	PP_3371	sensor histidine kinase
Reads supporting (aligned to +/- strand): ref base G (163/186); new base C (10/9); total (173/195)
Fisher's exact test for biased strand distribution p-value = 6.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

Only 100 of 369 total aligned reads displayed.

TCAGGTAAAAGCTGCACAACAGCGCGAACAGCGCGCTCATCACCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACGTCAAGCAACTGGCGGCCGGCCCGCTCGCCATGCCGTTCGAGCATTTGCA  >  NC_002947/3820673‑3820969

  aGGTAAAAGCTGCACAACAGCGCGAACAGCGCGCTCATCACCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGcc                                                                                                                                                   <  2:6082513/150‑1 (MQ=255)

     tAAAAGCTGCACAACAGCGCGAACAGCGCGCTCATCACCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGc                                                                                                                                                >  1:8685333/1‑150 (MQ=255)

          gCTGCACAACAGCGCGAACAGCGCGCTCATCACCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGaa                                                                                                                                           <  1:8113624/150‑1 (MQ=255)

          gCTGCACAACAGCGCGAACAGCGCGCTCATCACCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGaa                                                                                                                                           <  2:5189396/150‑1 (MQ=255)

          gCTGCACAACAGCGCGAACAGCGCGCTCATCACCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGAGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGaa                                                                                                                                           <  2:6378990/150‑1 (MQ=255)

           cTGCACACCAGCGCGAACAGCGCGCTCATCACCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAg                                                                                                                                          <  2:6647965/150‑1 (MQ=255)

           cTGCACAACAGCGCGAACAGCGCGCTCATCACCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAg                                                                                                                                          <  2:774822/150‑1 (MQ=255)

           cTGCACAACAGCGCGAACAGCGCGCTCATCACCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGGCACTGCTCGacc                                                                                                                                          <  2:1715928/150‑3 (MQ=255)

           cTGCACAACAGCGCGAACAGCGCGCTCATCACCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGGCACTGCTCGacc                                                                                                                                          <  1:3384809/150‑3 (MQ=255)

           cTGCACAACAGCGCGAACAGCGCGCTCATCACCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGGCACTGCTCGacc                                                                                                                                          <  1:6006025/150‑3 (MQ=255)

            tGCACAACAGCGCGAACAGCGCGCTCATCACCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAgc                                                                                                                                         >  2:5903804/1‑150 (MQ=255)

               acaacaGCGCGAACAGCGCGCTCATCACCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAgcgcg                                                                                                                                      <  2:8066638/150‑1 (MQ=255)

               acaacaGCGCGAACAGCGCGCTCATCACCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAgcgcg                                                                                                                                      <  2:3630860/150‑1 (MQ=255)

                 aacagcgcgaacagcgcgCTCATCACCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGcc                                                                                                                                    >  1:3604146/1‑150 (MQ=255)

                 aacagcgcgaacagcgcgCTCATCACCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGcc                                                                                                                                    >  1:181116/1‑150 (MQ=255)

                  acagcgcgaacagcgcgCTCATCACCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCa                                                                                                                                   >  1:2579593/1‑150 (MQ=255)

                        cgaacagcgcgCTCATCACCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCa                                                                                                                             <  2:5789076/150‑1 (MQ=255)

                         gaacagcgcgCTCATCACCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCAc                                                                                                                            <  1:2381248/150‑1 (MQ=255)

                             agcgcgCTCATCACCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGtt                                                                                                                        <  1:4750237/150‑1 (MQ=255)

                             agcgcgCTCATCACCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGtt                                                                                                                        >  1:1040199/1‑150 (MQ=255)

                              gcgcgcTCATCACCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGTGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGttt                                                                                                                       <  2:1580801/150‑1 (MQ=255)

                               cgcgcTCATCACCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTg                                                                                                                      <  2:330903/150‑1 (MQ=255)

                                  gcTCATCACCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGgcc                                                                                                                   <  2:3052440/150‑1 (MQ=255)

                                    tcatcaCCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGgccag                                                                                                                 >  1:7673019/1‑150 (MQ=255)

                                    tcatcaCCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGgccag                                                                                                                 >  2:3912848/1‑150 (MQ=255)

                                         aCCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCcagccagc                                                                                                            <  2:7522522/150‑1 (MQ=255)

                                         aCCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCcagccagc                                                                                                            <  1:6522191/150‑1 (MQ=255)

                                           cGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCaa                                                                                                          <  2:1826021/150‑1 (MQ=255)

                                                cAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGt                                                                                                     <  1:8037420/150‑1 (MQ=255)

                                                cAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGt                                                                                                     <  2:1623481/150‑1 (MQ=255)

                                                    cACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTg                                                                                                 <  2:699476/150‑1 (MQ=255)

                                                    cACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCACCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTg                                                                                                 >  1:8307338/1‑150 (MQ=255)

                                                     aCGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGa                                                                                                <  2:8424501/150‑1 (MQ=255)

                                                     aCGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGa                                                                                                >  1:3016789/1‑150 (MQ=255)

                                                      cGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGAt                                                                                               <  2:7347749/150‑1 (MQ=255)

                                                       gAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGTCAGCAAGCGGTCTTGAtc                                                                                              >  2:8260391/1‑150 (MQ=255)

                                                          cGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGAtcgtc                                                                                           <  2:792321/150‑1 (MQ=255)

                                                              ccAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTaa                                                                                       <  1:2311536/150‑1 (MQ=255)

                                                              ccAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTaa                                                                                       >  2:848452/1‑150 (MQ=255)

                                                              ccAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCATaa                                                                                       >  1:7251145/1‑150 (MQ=255)

                                                               cAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTaaa                                                                                      >  2:4380767/1‑150 (MQ=255)

                                                               cAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTaaa                                                                                      >  2:5645050/1‑150 (MQ=255)

                                                                 aGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAgg                                                                                    <  1:8178639/150‑1 (MQ=255)

                                                                 aGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAgg                                                                                    <  2:2404220/150‑1 (MQ=255)

                                                                     tGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTa                                                                                <  2:2288110/150‑1 (MQ=255)

                                                                     tGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTa                                                                                <  2:8615326/150‑1 (MQ=255)

                                                                       cGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACt                                                                              <  2:5375287/150‑1 (MQ=255)

                                                                              gcgGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACAt                                                                       >  1:8451707/1‑150 (MQ=255)

                                                                              gcgGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACAt                                                                       >  2:4685313/1‑150 (MQ=255)

                                                                                   gCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGc                                                                  <  1:8582868/150‑1 (MQ=255)

                                                                                     cATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCgg                                                                <  2:2080029/150‑1 (MQ=255)

                                                                                     cATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCgg                                                                <  1:4190029/150‑1 (MQ=255)

                                                                                      aTACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGc                                                               >  2:3413723/1‑150 (MQ=255)

                                                                                       tACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCt                                                              <  2:927392/150‑1 (MQ=255)

                                                                                         ccGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCtgt                                                            <  1:5633963/150‑1 (MQ=255)

                                                                                         ccGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCtgt                                                            <  1:918407/150‑1 (MQ=255)

                                                                                         ccGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCtgt                                                            >  2:5419098/1‑150 (MQ=255)

                                                                                           gACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTgg                                                          <  1:2935224/150‑1 (MQ=255)

                                                                                           gACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTgg                                                          <  1:7521634/150‑1 (MQ=255)

                                                                                              acGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGa                                                       >  1:3730429/1‑150 (MQ=255)

                                                                                              acGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGa                                                       >  2:5154800/1‑150 (MQ=255)

                                                                                                 ggTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGacac                                                    <  2:1084931/150‑1 (MQ=255)

                                                                                                 ggTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGacac                                                    <  2:3847365/150‑1 (MQ=255)

                                                                                                    ggCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACGTc                                                 <  2:4190418/150‑1 (MQ=255)

                                                                                                     gCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACGTCa                                                >  1:7802676/1‑150 (MQ=255)

                                                                                                      ccACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACGTCaa                                               <  2:6590810/150‑1 (MQ=255)

                                                                                                        aCCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGc                                                                                                                                                <  1:6709518/51‑1 (MQ=255)

                                                                                                        aCCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACGTCAAGc                                             >  2:1357830/1‑150 (MQ=255)

                                                                                                        aCCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACGTCAAGc                                             >  1:4525368/1‑150 (MQ=255)

                                                                                                        aCCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACGTCAAGc                                             >  2:21964/1‑150 (MQ=255)

                                                                                                           gAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACGTCAAGCAAc                                          >  1:4372731/1‑150 (MQ=255)

                                                                                                            aaGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACGTCAAGCAACt                                         >  2:2502366/1‑150 (MQ=255)

                                                                                                            aaGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACGTCAAGCAACt                                         >  1:3866005/1‑150 (MQ=255)

                                                                                                             aGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACGTCAAGCAACTg                                        <  1:50710/150‑1 (MQ=255)

                                                                                                              gcgcgcAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACGTCAAGCAACTgg                                       <  2:2878412/150‑1 (MQ=255)

                                                                                                               cgcgcAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACGTCAAGCAACTggc                                      <  2:4447831/150‑1 (MQ=255)

                                                                                                                gcgcAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACGTCAAGCAACTggcg                                     <  1:3998572/150‑1 (MQ=255)

                                                                                                                      caCCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACGTCAAGCAACTGGcggccggc                               >  1:1979373/1‑150 (MQ=255)

                                                                                                                       aCCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACGTCAAGCAACTGGCggccggcc                              >  1:3238304/1‑150 (MQ=255)

                                                                                                                       aCCAGTGGGCGGCCGTCACGGTCGAGCAGAGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGACAGTCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGAGGCGACACGTCAAGCAACTGGCggccgccc                              >  1:166258/1‑150 (MQ=255)

                                                                                                                         cAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACGTCAAGCAACTGGCGGCCGGCCCg                            >  2:4587388/1‑150 (MQ=255)

                                                                                                                                 ggcCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACGTCAAGCAACTGGCGGCCGGCCCGCTCGCCag                    >  1:6381362/1‑149 (MQ=255)

                                                                                                                                  gcCGTCACGGTCGAGCAGGGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACGTCAAGCAACTGGCGGCCGGCCCGCTCGCCATg                   <  1:5408643/150‑1 (MQ=255)

                                                                                                                                       gcgcttCGAGCAGGCACTGCTCGACGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACGTCAAGCAACTGGCGGCCGGCCCGCTCGCCATGCCGtt              >  1:1249173/6‑150 (MQ=255)

                                                                                                                                       cACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACGTCAAGCAACTGGCGGCCGGCCCGCTCGCCATGCCGtt              <  1:704447/150‑1 (MQ=255)

                                                                                                                                          ggTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACGTCAAGCAACTGGCGGCCGGCCCGCTCGCCATGCCGTTCGa           <  2:2644248/150‑1 (MQ=255)

                                                                                                                                          ggTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACGTCAAGCAACTGGCGGCCGGCCCGCTCGCCATGCCGTTCGa           <  1:4063953/150‑1 (MQ=255)

                                                                                                                                           ttCGAGCAGGCACTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACGTCAAGCAACTGGCGGCCGGCCCGCTCGCCATGCCGTTCGAg          >  2:2819256/2‑150 (MQ=255)

                                                                                                                                           ttCGAGCAGGCACTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACGTCAAGCAACTGGCGGCCGGCCCGCTCGCCATGCCGTTCGAg          >  1:2836304/2‑150 (MQ=255)

                                                                                                                                           gTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACGTCAAGCAACTGGCGGCCGGCCCGCTCGCCATGCCGTTCGAg          <  2:4317487/150‑1 (MQ=255)

                                                                                                                                            tCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACGTCAAGCAACTGGCGGCCGGCCCGCTCGCCATGCCGTTCGAGc         >  1:1515911/1‑150 (MQ=255)

                                                                                                                                            tCGAGCAGGCACTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACGTCAAGCAACTGGCGGCCGGCCCGCTCGCCATGCCGTTCGa           <  1:7468169/148‑1 (MQ=255)

                                                                                                                                               aGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACGTCAAGCAACTGGCGGCCGGCCCGCTCGCCATGCCGTTCGAGCAtt      <  2:6240990/150‑1 (MQ=255)

                                                                                                                                               aGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACGTCAAGCAACTGGCGGCCGGCCCGCTCGCCATGCCGTTCGAGCAtt      <  2:765522/150‑1 (MQ=255)

                                                                                                                                                gCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACGTCAAGCAACTGGCGGCCGGCCCGCTCGCCATGCCGTTCGAGCAttt     <  1:5552327/150‑1 (MQ=255)

                                                                                                                                                gCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACGTCAAGCAACTGGCGGCCGGCCCGCTCGCCATGCCGTTCGAGCAttt     >  2:4947874/1‑150 (MQ=255)

                                                                                                                                                 cAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACGTCAAGCAACTGGCGGCCGGCCCGCTCGCCATGCCGTTCGAGCATTTg    <  2:809366/150‑1 (MQ=255)

                                                                                                                                                  aGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACGTCAAGCAACTGGCGGCCGGCCCGCTCGCCATGCCGTTCGAGCATTTGc   <  2:7270675/150‑1 (MQ=255)

                                                                                                                                                   gTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCtt                                                                                                  <  2:6179267/54‑1 (MQ=255)

                                                                                                                                                   gTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCtt                                                                                                  >  1:6179267/1‑54 (MQ=255)

TCAGGTAAAAGCTGCACAACAGCGCGAACAGCGCGCTCATCACCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACGTCAAGCAACTGGCGGCCGGCCCGCTCGCCATGCCGTTCGAGCATTTGCA  >  NC_002947/3820673‑3820969

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