BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,479,784	C→A	6.3%	A58A (GCG→GCT)	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,784	0	C	A	6.3%	719.8 / 22.3	212	A58A (GCG→GCT)	PP_3968	sensor histidine kinase
Reads supporting (aligned to +/- strand): ref base C (99/98); new base A (7/8); total (106/106)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 220 total aligned reads displayed.

GATAGAACACTCGAGCGCAGCGCTCGGTACTCAGAGACCAATTGATCCATCGAAAAGCCCGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAGTCCGTCTTGACCGAACGTGCGAAATCCTCCCAAGCCTGGAGAATGCTTCCCATATCACT  >  NC_002947/4479635‑4479933

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

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

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

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

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

        acTCGAGCGCAGCGCTCGGTACTCAGAGACCAATTGATCCATCGAAAAGCCCGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCGCATCAGCAGTCTg                                                                                                                                               <  1:1379429/150‑1 (MQ=255)

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

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

         cTCGAGCGCAGCGCTCGGTACTCAGAGACCAATTGATCCATCGAAAAGCCCGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTGCTCAGACTGCGCATCAGCAGTCTGa                                                                                                                                              <  2:3430363/150‑1 (MQ=255)

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

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

          tCGAGCGCAGCGCTCGGTACTCAGAGACCAATTGATCCATCGAAAAGCCCGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAg                                                                                                                                             <  2:223706/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)

                         ggTACTCAGAGACCAATTGATCCATCGAAAAGCCCGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAg                                                                                                                              >  2:1679445/1‑150 (MQ=255)

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

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

                             cTCAGAGACCAATTGATCCATCGAAAAGCCCGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGc                                                                                                                          >  1:881035/1‑150 (MQ=255)

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

                                agagACCACTTGATCCATCGAAAAGCCCGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAAc                                                                                                                       <  1:1935033/150‑1 (MQ=255)

                                   gaCCAATTGATCCATCGAAAAGCCCGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGt                                                                                                                    <  2:3846524/150‑1 (MQ=255)

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

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

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

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

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

                                             tCCATCGAAAAGCCCGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAAtc                                                                                                          <  1:3498058/150‑1 (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:2507443/1‑150 (MQ=255)

                                                   gAAAAGCCCGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCa                                                                                                    >  2:2066528/1‑150 (MQ=255)

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

                                                   gAAAAGCCCGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCa                                                                                                    >  2:3795192/1‑150 (MQ=255)

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

                                                       aGCCCGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGt                                                                                                <  2:393366/150‑1 (MQ=255)

                                                         cccGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGa                                                                                              >  2:2920746/1‑150 (MQ=255)

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

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

                                                           cGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGAtt                                                                                            <  1:291662/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)

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

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

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

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

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

                                                                      gAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCgg                                                                                 <  2:3962750/150‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

                                                                                  gggTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAg                                                                     <  2:3780638/150‑1 (MQ=255)

                                                                                  gggTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAg                                                                     <  2:1724088/150‑1 (MQ=255)

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

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

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

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

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

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

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

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

                                                                                                     gctgctCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAGTCCGTCTTGAc                                                  >  2:2606324/1‑150 (MQ=255)

                                                                                                      ctgctCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAGTCCGTCTTGAcc                                                 >  2:4279916/1‑150 (MQ=255)

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

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

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

                                                                                                               agagGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAGTCCGTCTTGACCGAACGTGCg                                        <  2:2398908/150‑1 (MQ=255)

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

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

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

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

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

                                                                                                                                     tctcAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAGTCCGTCTTGACCGAACGTGCGAAATCCTCCCAAGCCTGGAGaa                  <  1:1683381/150‑1 (MQ=255)

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

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

                                                                                                                                       tcAGACTGCGCATCAGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAGTCCGTCTTGACCGAACGTGCGAAATCCTCCCAAGCCTGGAGAATg                >  2:2986975/1‑150 (MQ=255)

                                                                                                                                         aGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCACCGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAGTCCGTCTTGACCGAACGTGCGAAATCCTCCCAAGCCTGGAGAATGCt              >  1:1724394/1‑150 (MQ=255)

                                                                                                                                         aGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAGTCCGTCTTGACCGAACGTGCGAAATCCTCCCAAGCCTGGAGAATGCt              >  2:3098704/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          >  2:2931958/1‑150 (MQ=255)

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

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

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

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

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

                                                                                                                                                   tGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAGTCCGTCTTGACCGAACGTGCGAAATCCTCCCAAGCCTGGAGAATGCTTCCCATATCa    <  1:3256823/150‑1 (MQ=255)

                                                                                                                                                   tGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAGTCCGTCTTGACCGAACGTGCGAAATCCTCCCAAGCCTGGAGAATGCTTCCCATATCa    >  2:1931665/1‑150 (MQ=255)

                                                                                                                                                    gcgcAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAGTCCGTCTTGACCGAACGTGCGAAATCCTCCCAAGCCTGGAGAATGCTTCCCATATCAc   <  1:3649846/150‑1 (MQ=255)

GATAGAACACTCGAGCGCAGCGCTCGGTACTCAGAGACCAATTGATCCATCGAAAAGCCCGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAGTCCGTCTTGACCGAACGTGCGAAATCCTCCCAAGCCTGGAGAATGCTTCCCATATCACT  >  NC_002947/4479635‑4479933

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