BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,479,779	T→G	5.9%	Q60P (CAG→CCG) ‡	PP_3968 ←	sensor histidine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,479,779	0	T	G	5.9%	700.3 / 25.1	202	Q60P (CAG→CCG) ‡	PP_3968	sensor histidine kinase
Reads supporting (aligned to +/- strand): ref base T (92/97); new base G (5/7); total (98/104)
Fisher's exact test for biased strand distribution p-value = 7.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

Only 100 of 216 total aligned reads displayed.

AGAGCGATAGAACACTCGAGCGCAGCGCTCGGTACTCAGAGACCAATTGATCCATCGAAAAGCCCGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAGTCCGTCTTGACCGAACGTGCGAAATCCTCCCAAGCCTGGAGAATGCTTCCCATA  >  NC_002947/4479630‑4479928

agagCGATAGAACACTCGAGCGCAGCGCTCGGTACTCAGAGACCAATTGATCCATCGAAAAGCCCGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGActgct                                                                                                                                                       >  2:3871121/1‑150 (MQ=255)

   gCGATAGAACACTCGAGCGCAGCGCTCGGTACTCAGAGACCAATTGATCCATCGAAAAGCCCGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGAt                                                                                                                                                    <  1:85278/150‑1 (MQ=255)

     gATAGAACACTCGAGCGCAGCGCTCGGTACTCAGAGACCAATTGATCCATCGAAAAGCCCGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATgc                                                                                                                                                  <  1:2632276/150‑1 (MQ=255)

          aaCACTCGAGCGCAGCGCTCGGTACTCAGAGACCAATTGATCCATCGAAAAGCCCGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGt                                                                                                                                             >  2:3112610/1‑150 (MQ=255)

          aaCACTCGAGCGCAGCGCTCGGTACTCAGAGACCAATTGATCCATCGAAAAGCCCGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGt                                                                                                                                             >  1:449511/1‑150 (MQ=255)

           acacTCGAGCGCAGCGCTCGGTACTCAGAGACCAATTGATCCATCGAAAAGCCCGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTc                                                                                                                                            >  2:3256629/1‑150 (MQ=255)

            cacTCGAGCGCAGCGCTCGGTACTCAGAGACCAATTGATCCATCGAAAAGCCCGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCt                                                                                                                                           <  1:688362/150‑1 (MQ=255)

             acTCGAGCGCAGCGCTCGGTACTCAGAGACCAATTGATCCATCGAAAAGCCCGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTg                                                                                                                                          >  1:3669067/1‑150 (MQ=255)

             acTCGAGCGCAGCGCTCGGTACTCAGAGACCAATTGATCCATCGAAAAGCCCGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTg                                                                                                                                          >  2:2298415/1‑150 (MQ=255)

             acTCGAGCGCAGCGCTCGGTACTCAGAGACCAATTGATCCATCGAAAAGCCCGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTg                                                                                                                                          >  2:1051384/1‑150 (MQ=255)

              cTCGAGCGCAGCGCTCGGTACTCAGAGACCAATTGATCCATCGAAAAGCCCGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGa                                                                                                                                         <  1:3212192/150‑1 (MQ=255)

              cTCGAGCGCAGCGCTCGGTACTCAGAGACCAATTGATCCATCGAAAAGCCCGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCGCATCAGCAGTCTGa                                                                                                                                         <  1:1680126/150‑1 (MQ=255)

               tCGAGCGCAGCGCTCGGTACTCAGAGACCAATTGATCCATCGAAAAGCCCGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAg                                                                                                                                        <  2:2056372/150‑1 (MQ=255)

               tCGAGCGCAGCGCTCGGTACTCAGAGACCAATTGATCCATCGAAAAGCCCGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAg                                                                                                                                        <  1:3866588/150‑1 (MQ=255)

               tCGAGCGCAGCGCTCGGTACTCAGAGACCAATTGATCCATCGAAAAGCCCGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAg                                                                                                                                        >  1:4272181/1‑150 (MQ=255)

               tCGAGCGCAGCGCTCGGTACTCAGAGACCAATTGATCCATCGAAAAGCCCGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAg                                                                                                                                        >  2:1393677/1‑150 (MQ=255)

               tCGAGCGCAGCGCTCGGTACTCAGAGACCAATTGATCCATCGAAAAGCCCGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAg                                                                                                                                        <  2:223671/150‑1 (MQ=255)

               tCGAGCGCAGCGCTCGGTACTCAGAGACCAATTGATCCATCGAAAAGCCCGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAg                                                                                                                                        >  1:1132039/1‑150 (MQ=255)

                 gAGCGCAGCGCTCGGTACTCAGAGACCAATTGATCCATCGAAAAGCCCGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCGCATCAGCAGTCTGAgaa                                                                                                                                      <  1:4453172/150‑3 (MQ=255)

                    cgcagcgcTCGGTACTCAGAGACCAATTGATCCATCGAAAAGCCCGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTGCTCAGACTGCGCATCAGCAGTCTGAGAAggc                                                                                                                                   <  1:1481641/150‑3 (MQ=255)

                       tgcgcTCGGTACTCAGAGACCAATTGATCCATCGAAAAGCCCGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGc                                                                                                                                >  2:1460180/2‑150 (MQ=255)

                        gcgcTCGGTACTCAGAGACCAATTGATCCATCGAAAAGCCCGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCa                                                                                                                               <  2:2442478/150‑1 (MQ=255)

                               gTACTCAGAGACCAATTGATCCATCGAAAAGCCCGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGc                                                                                                                        <  2:1240951/150‑1 (MQ=255)

                                    cAGAGACCAATTGATCCATCGAAAAGCCCGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCaa                                                                                                                   <  2:4431068/150‑1 (MQ=255)

                                         aCCAATTGATCCATCGAAAAGCCCGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTg                                                                                                              >  2:2922641/1‑150 (MQ=255)

                                          ccAATTGATCCATCGAAAAGCCCGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATGAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCGTATCAGCGGCAACGGTgc                                                                                                             <  2:328003/150‑1 (MQ=255)

                                           cAATTGATCCATCGAAAAGCCCGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTgcg                                                                                                            >  1:1087112/1‑150 (MQ=255)

                                           cAATTGATCCATCGAAAAGCCCGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTgcg                                                                                                            >  1:88521/1‑150 (MQ=255)

                                           cAATTGATCCATCGAAAAGCCCGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTgcg                                                                                                            >  2:3091029/1‑150 (MQ=255)

                                            aaTTGATCCATCGAAAAGCCCGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTgcgc                                                                                                           >  2:2844841/1‑150 (MQ=255)

                                             aTTGATCCATCGAAAAGCCCGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCa                                                                                                          <  2:449511/150‑1 (MQ=255)

                                                gATCCATCGAAAAGCCCGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGaa                                                                                                       <  2:428606/150‑1 (MQ=255)

                                                 aTCCATCGAAAAGCCCGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAAt                                                                                                      >  1:720643/1‑150 (MQ=255)

                                                    cATCGAAAAGCCCGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAAtctt                                                                                                   >  1:3219327/1‑150 (MQ=255)

                                                     aTCGAAAAGCCCGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAAtcttc                                                                                                  >  1:3160193/1‑150 (MQ=255)

                                                     aTCGAAAAGCCCGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAAtcttc                                                                                                  <  1:2127325/150‑1 (MQ=255)

                                                        gAAAAGCCCGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCa                                                                                               <  2:4030332/150‑1 (MQ=255)

                                                               ccGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCACAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGAt                                                                                        <  2:1986635/150‑1 (MQ=255)

                                                                cGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGAtt                                                                                       >  1:158355/1‑150 (MQ=255)

                                                                cGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGAtt                                                                                       >  2:2518566/1‑150 (MQ=255)

                                                                cGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGAtt                                                                                       >  1:4514769/1‑150 (MQ=255)

                                                                cGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGAtt                                                                                       <  1:291662/150‑1 (MQ=255)

                                                                 gCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTc                                                                                      <  1:3292259/150‑1 (MQ=255)

                                                                 gCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTc                                                                                      >  1:2808994/1‑150 (MQ=255)

                                                                    aTCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTg                                                                                   <  1:773528/150‑1 (MQ=255)

                                                                     tCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGa                                                                                  <  1:2257226/150‑1 (MQ=255)

                                                                     tCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGa                                                                                  <  1:201646/150‑1 (MQ=255)

                                                                     tCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGAATGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGa                                                                                  <  2:506908/150‑1 (MQ=255)

                                                                      cAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAg                                                                                 <  1:958795/150‑1 (MQ=255)

                                                                      cAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAg                                                                                 <  2:2509606/150‑1 (MQ=255)

                                                                       aaaCGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAgg                                                                                >  2:4216167/1‑150 (MQ=255)

                                                                       aaaCGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATGAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCAGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCACCGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAgg                                                                                >  2:2992883/1‑150 (MQ=255)

                                                                         aCGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGcc                                                                              <  1:1051384/150‑1 (MQ=255)

                                                                           gAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTGGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCgg                                                                            <  2:3962797/150‑1 (MQ=255)

                                                                           gAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCgg                                                                            >  1:3725740/1‑150 (MQ=255)

                                                                             gTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCt                                                                          <  2:3213332/150‑1 (MQ=255)

                                                                              tAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTg                                                                         >  2:2888113/1‑150 (MQ=255)

                                                                               aaTGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGa                                                                        <  1:3028057/150‑1 (MQ=255)

                                                                                 tGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAAt                                                                      >  1:2219311/1‑150 (MQ=255)

                                                                                    cATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCt                                                                   >  2:1744268/1‑150 (MQ=255)

                                                                                    cATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCt                                                                   >  2:3806623/1‑150 (MQ=255)

                                                                                     aTGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTg                                                                  >  2:3683488/1‑150 (MQ=255)

                                                                                     aTGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTg                                                                  >  2:3413684/1‑150 (MQ=255)

                                                                                      tGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGa                                                                 <  1:1181560/150‑1 (MQ=255)

                                                                                      tGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCGGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGa                                                                 <  2:310466/150‑1 (MQ=255)

                                                                                      tGGGTCTGAGCCGCCGTATCGCTGCACATGAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGACTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCAGAGGCCGGCTGAATCCTGa                                                                 <  2:295063/150‑1 (MQ=255)

                                                                                       gggTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAg                                                                >  1:4431822/1‑150 (MQ=255)

                                                                                            tGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGt                                                           >  2:4373872/1‑150 (MQ=255)

                                                                                               gccgccGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAg                                                        <  2:3972897/150‑1 (MQ=255)

                                                                                                ccgccgTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAGt                                                       <  2:1143890/150‑1 (MQ=255)

                                                                                                  gccgTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAGTcc                                                     <  1:4352945/150‑1 (MQ=255)

                                                                                                   ccgTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAGTCCg                                                    <  1:225279/150‑1 (MQ=255)

                                                                                                      tATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAGTCCGTCt                                                 <  1:2420432/150‑1 (MQ=255)

                                                                                                         cGCTGCTCCTTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAGTCCGTCTTGa                                              <  2:3746475/150‑1 (MQ=255)

                                                                                                            tgctCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAGTCCGTCTTGACCg                                           >  2:3575079/1‑150 (MQ=255)

                                                                                                                cATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAGTCCGTCTTGACCGAACg                                       >  2:3203405/1‑150 (MQ=255)

                                                                                                                  ttAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAGTCCGTCTTGACCGAACGTg                                     >  1:4191209/1‑150 (MQ=255)

                                                                                                                  ttAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAGTCCGTCTTGACCGAACGTg                                     >  1:1237321/1‑150 (MQ=255)

                                                                                                                         gCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAGTCCGTCTTGACCGAACGTGCGAAAtc                              <  1:4046576/150‑1 (MQ=255)

                                                                                                                         gCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAGTCCGTCTTGACCGAACGTGCGAAAtc                              <  1:1914678/150‑1 (MQ=255)

                                                                                                                           ctctGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAGTCCGTCTTGACCGAACGTGCGAAAtcct                            >  2:1758660/1‑150 (MQ=255)

                                                                                                                            tctGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAGTCCGTCTTGACCGAACGTGCGAAAtcctc                           >  1:470329/1‑150 (MQ=255)

                                                                                                                             ctGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAGTCCGTCTTGACCGAACGTGCGAAAtcctcc                          <  1:1624710/150‑1 (MQ=255)

                                                                                                                                  ttGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAGTCCGTCTTGACCGAACGTGCGAAATCCTCCCAAGc                     >  2:3787124/1‑150 (MQ=255)

                                                                                                                                  ttGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAGTCCGTCTTGACCGAACGTGCGAAATCCTCCCAAGc                     >  1:1090088/1‑150 (MQ=255)

                                                                                                                                       ccTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAGTCCGTCTTGACCGAACGTGCGAAATCCTCCCAAGCCTGga                <  1:1460180/150‑1 (MQ=255)

                                                                                                                                        cTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAGTCCGTCTTGACCGAACGTGCGAAATCCTCCCAAGCCTGgag               >  1:3888475/1‑150 (MQ=255)

                                                                                                                                            tcAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAGTCCGTCTTGACCGAACGTGCGAAATCCTCCCAAGCCTGGAGAATg           >  1:2979075/1‑150 (MQ=255)

                                                                                                                                            tcAGACTGCGGATCAGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAGTCCGTCTTGACCGAACGTGCGAAATCCTCCCAAGCCTGGAGAATg           >  2:1577816/1‑150 (MQ=255)

                                                                                                                                             cAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAGTCCGTCTTGACCGAACGTGCGAAATCCTCCCAAGCCTGGAGAATGc          >  1:1365186/1‑150 (MQ=255)

                                                                                                                                               gACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAGTCCGTCTTGACCGAACGTGCGAAATCCTCCCAAGCCTGGAGAATGCtt        >  2:2119558/1‑150 (MQ=255)

                                                                                                                                                 ctgctgATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAGTCCGTCTTGACCGAACGTGCGAAATCCTCCCAAGCCTGGAGAATGCTTcc      <  2:1159354/150‑1 (MQ=255)

                                                                                                                                                 ctgatgATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACTGTGCGCAGAATATTCTCTGAGTGATTCCTGATGCCGTCTGTATCCTGAGTGTGTAAGTCCTTCTTGACCGAACGTGCGAAATACTCCAAAGCCTGGATAATTCttgc      <  2:881035/150‑3 (MQ=255)

                                                                                                                                                  tgctgATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAGTCCGTCTTGACCGAACGTGCGAAATCCTCCCAAGCCTGGAGAATGCTTccc     >  1:3458860/1‑150 (MQ=255)

                                                                                                                                                    ctgATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAGTCCGTCTTGACCGAACGTGCGAAATCCTCCCAAGCCTGGAGAATGCTTCCCat   >  2:412229/1‑150 (MQ=255)

                                                                                                                                                     tgATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTggg                                                            >  1:169009/1‑92 (MQ=255)

                                                                                                                                                     tgATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTggg                                                            <  2:169009/92‑1 (MQ=255)

                                                                                                                                                     tgATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAGTCCGTCTTGACCGAACGTGCGAAATCCTCCCAAGCCTGGAGAATGCTTCCCcta  >  2:3947774/1‑150 (MQ=255)

                                                                                                                                                     tgATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAGTCCGTCTTGACCGAACGTGCGAAATCCTCCCAAGCCTGGAGAATGCTTCCCata  >  2:1837544/1‑150 (MQ=255)

                                                                                                                                                     tgATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAGTCCGTCTTGACCGAACGTGCGAAATCCTCCCAAGCCTGGAGAATGCTTCCCata  >  1:3244922/1‑150 (MQ=255)

AGAGCGATAGAACACTCGAGCGCAGCGCTCGGTACTCAGAGACCAATTGATCCATCGAAAAGCCCGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAGTCCGTCTTGACCGAACGTGCGAAATCCTCCCAAGCCTGGAGAATGCTTCCCATA  >  NC_002947/4479630‑4479928

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