BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,820,823	C→A	5.2%	A97S (GCA→TCA) ‡	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,823	0	C	A	5.2%	1231.8 / 38.4	363	A97S (GCA→TCA) ‡	PP_3371	sensor histidine kinase
Reads supporting (aligned to +/- strand): ref base C (162/182); new base A (10/9); total (172/191)
Fisher's exact test for biased strand distribution p-value = 6.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 370 total aligned reads displayed.

CAGGTAAAAGCTGCACAACAGCGCGAACAGCGCGCTCATCACCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACGTCAAGCAACTGGCGGCCGGCCCGCTCGCCATGCCGTTCGAGCATTTGCAGCT  >  NC_002947/3820674‑3820972

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

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

   gTAAAAGCTGCACAACAGCGCGAACAGCGCGCTCATCACCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTg                                                                                                                                                    >  2:1964503/1‑150 (MQ=255)

   gTAAAAGCTGCACAACAGCGCGAACAGCGCGCTCATCACCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTg                                                                                                                                                    >  2:4236705/1‑150 (MQ=255)

         gCTGCACAACAGCGCGAACAGCGCGCTCATCACCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGaa                                                                                                                                              <  1:8114096/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                                                                                                                                             <  1:6006025/150‑3 (MQ=255)

          cTGCACAACAGCGCGAACAGCGCGCTCATCACCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGGCACTGCTCGacc                                                                                                                                             <  1:3395364/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:181116/1‑150 (MQ=255)

                aacagcgcgaacagcgcgCTCATCACCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGcc                                                                                                                                       >  1:3604146/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:1040199/1‑150 (MQ=255)

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

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

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

                                  cTCATCACCGCCCCAGTCACGAACGGCCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGgcca                                                                                                                     <  1:5683077/150‑1 (MQ=255)

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

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

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

                                                  tCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCtt                                                                                                     >  2:1464839/1‑150 (MQ=255)

                                                   cACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTg                                                                                                    <  2:699476/150‑1 (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)

                                                        aCGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGAtcgt                                                                                               <  1:5220899/150‑1 (MQ=255)

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

                                                            gCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTa                                                                                           >  1:8600807/1‑150 (MQ=255)

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

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

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

                                                                aGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAgg                                                                                       <  1:8178681/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)

                                                                        tGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTg                                                                               <  1:8307809/150‑1 (MQ=255)

                                                                        tGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTg                                                                               >  2:2787780/1‑150 (MQ=255)

                                                                        tGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTg                                                                               >  2:7262179/1‑150 (MQ=255)

                                                                          ggCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGca                                                                             >  2:5169865/1‑150 (MQ=255)

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

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

                                                                                   ccATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCg                                                                    >  1:8744011/1‑150 (MQ=255)

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

                                                                                     aTACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCgg                                                                                                         <  1:7678852/111‑1 (MQ=255)

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

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

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

                                                                                         cGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCtgtg                                                              >  1:7369828/1‑150 (MQ=255)

                                                                                         cGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCtgtg                                                              >  1:6565650/1‑150 (MQ=255)

                                                                                         cGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCtgtg                                                              >  2:323616/1‑150 (MQ=255)

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

                                                                                          gACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTgg                                                             >  2:241916/1‑150 (MQ=255)

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

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

                                                                                                ggTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGacac                                                       >  2:5630776/1‑150 (MQ=255)

                                                                                                  tGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACGt                                                     <  2:1759997/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                                                >  1:4525368/1‑150 (MQ=255)

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

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

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

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

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

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

                                                                                                                   agcaCCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACGTCAAGCAACTGGcggccg                                    <  1:8555558/150‑1 (MQ=255)

                                                                                                                     caCCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACGTCAAGCAACTGGcggccggc                                  >  1:2954712/1‑150 (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)

                                                                                                                              gcggcCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACGTCAAGCAACTGGCGGCCGGCCCGCTCGcc                         >  2:3728719/1‑150 (MQ=255)

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

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

                                                                                                                                  ccGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACGTCAAGCAACTGGCGGCCGGCCCGCTCGCCATGc                     >  2:2561695/1‑150 (MQ=255)

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

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

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

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

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

                                                                                                                                          ttCGAGCAGGCACTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGAGCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACGTCAAGCAACTGGCGGCCGGCCCGCTCGCCATGCCGTTCGAg             >  1:933931/2‑150 (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)

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

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

                                                                                                                                                     ccTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACGTCAAGCAACTGGCGGCCGGCCCGCTCGCCATGCCGTTCGAGCATTTGCAGCt  >  2:2098410/1‑150 (MQ=255)

CAGGTAAAAGCTGCACAACAGCGCGAACAGCGCGCTCATCACCGCCCCAGTCACGAACGGGCCAAGGTTGCGTGGGCGCGGTGCCATACCGACACGGGTGGCCACCGAAGCGCGCAGCACCAGTGGGCGGCCGTCACGGTCGAGCAGTGCCTGCTCGAAGCGCGCCACGTGCACCGTTTGGCCAGCCAGCAAGCGGTCTTGATCGTCGTAAAGGCCTACTTGCACATCCTGCGGCTGTGGCGACACGTCAAGCAACTGGCGGCCGGCCCGCTCGCCATGCCGTTCGAGCATTTGCAGCT  >  NC_002947/3820674‑3820972

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