BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,423,817	A→C	5.8%	V84G (GTA→GGA) ‡	PP_4758 ←	MFS transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,423,817	0	A	C	5.8%	1070.8 / 39.3	312	V84G (GTA→GGA) ‡	PP_4758	MFS transporter
Reads supporting (aligned to +/- strand): ref base A (139/155); new base C (13/5); total (152/160)
Fisher's exact test for biased strand distribution p-value = 5.14e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.01e-01

Only 100 of 315 total aligned reads displayed.

CTTCGGCCACCCCGAGCAGAAAACGCAGCACGTAGAACTGGGTTTCGTTCTGCACAAAGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGTGCGGGCGGACGAAGCCAATGTTGACGCGGTCGAGGTAGTTGACGATGAACATGATGACGAACAG  >  NC_002947/5423668‑5423965

cTTCGGCCACCCCGAGCAGAAAACGCAGCACGTAGAACTGGGTTTCGTTCTGCACAAAGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTa                                                                                                                                                      <  2:2862404/150‑1 (MQ=255)

cTTCGGCCACCCCGAGCAGAAAACGCAGCACGTAGAACTGGGTTTCGTTCTGCACAAAGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTa                                                                                                                                                      <  1:555457/150‑1 (MQ=255)

    ggCCACCCCGAGCAGAAAACGCAGCTCGTAGAACTGGGTTTCGTTCTGCACAAAGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTc                                                                                                                                                  <  2:7134076/150‑1 (MQ=255)

    ggCCACCCCGAGCAGAAAACGCAGCACGTAGAACTGGGTTTCGTTCTGCACAAAGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTc                                                                                                                                                  <  2:7241189/150‑1 (MQ=255)

    ggCCACCCCGAGCAGAAAACGCAGCACGTAGAACTGGGTTTCGTTCTGCACAAAGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTc                                                                                                                                                  >  2:1805057/1‑150 (MQ=255)

      ccACCCCGAGCAGAAAACGCAGCACGTAGAACTGGGTTTCGTTCTGCACAAAGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGTACCTTCGa                                                                                                                                                <  1:6378275/150‑1 (MQ=255)

      ccACCCCGAGCAGAAAACGCAGCACGTAGAACTGGGTTTCGTTCTGCACAAAGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGTACCTTCGa                                                                                                                                                <  1:3442386/150‑1 (MQ=255)

      ccACCCCGAGCAGAAAACGCAGCACGTAGAACTGGGTTTCGTTCTGCACAAAGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTGTTGCAGCAGCATGTTCGAAGGTACTTCGa                                                                                                                                                >  2:3336430/1‑150 (MQ=255)

        aCCCCGAGCAGAAAACGCAGCACGTAGAACTGGGTTTCGTTCTGCACAAAGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGAAc                                                                                                                                              <  2:2294450/150‑1 (MQ=255)

         ccccGAGCAGAAAACGCAGCACGTAGAACTGGGTTTCGTTCTGCACAAAGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGTACCTTCGAACa                                                                                                                                             <  2:6412306/150‑1 (MQ=255)

         ccccGAGCAGAAAACGCAGCACGTAGAACTGGGTTTCGTTCTGCACAAAGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGAACa                                                                                                                                             <  2:831612/150‑1 (MQ=255)

           ccGAGCAGAAAACGCAGCACGTAGAACTGGGTTTCGTTCTGCACAAAGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGAACAgc                                                                                                                                           <  1:2063700/150‑1 (MQ=255)

            cGAGCAGAAAACGCAGCACGTAGAACTGGGTTTCGTTCTGCACAAAGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGTACCTTCGAACatgc                                                                                                                                          <  1:4221605/150‑4 (MQ=255)

             gAGCAGAAAACGCAGCACGTAGAACTGGGTTTCGTTCTGCACAAAGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGAACAgcgc                                                                                                                                         <  2:4663150/150‑1 (MQ=255)

             gAGCAGAAAACGCAGCACGTAGAACTGGGTTTCGTTCTGCACAAAGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGAACAgcgc                                                                                                                                         <  1:3609604/150‑1 (MQ=255)

                  gAAAACGCAGCACGTAGAACTGGGTTTCGTTCTGCACAAAGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGAACAGCGCGTAAc                                                                                                                                    <  2:2109805/150‑1 (MQ=255)

                   aaaaCGCAGCACGTAGAACTGGGTTTCGTTCTGCACAAAGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGAACAGCGCGTAAcc                                                                                                                                   >  2:6648669/1‑150 (MQ=255)

                     aaCGCAGCACGTAGAACTGGGTTTCGTTCTGCACAAAGGCCATGGCCGTGGCCACCAG‑CCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGAACAGCGCGTAACCGAt                                                                                                                                <  2:4949602/150‑1 (MQ=255)

                      aCGCAGCACGTAGAACTGGGTTTCGTTCTGCACAAAGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGAACAGCGCGTAACCGAt                                                                                                                                <  1:662527/150‑1 (MQ=255)

                        gcagcaCGTAGAACTGGGTTTCGTTCTGCACAAAGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGAACAGCGCGTAACCGATga                                                                                                                              >  1:6709920/1‑150 (MQ=255)

                                  gAACTGGGTTTCGTTCTGCACAAAGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGAACAGCGCGTAACCGATGAAGAACAGccc                                                                                                                    <  2:5867570/150‑1 (MQ=255)

                                   aaCTGGGTTTCGTTCTGCACAAAGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCg                                                                                                                   <  1:3715013/150‑1 (MQ=255)

                                   aaCTGGGTTTCGTTCTGCACAAAGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCg                                                                                                                   <  1:148572/150‑1 (MQ=255)

                                   aaCTGGGTTTCGTTCTGCACAAAGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCg                                                                                                                   <  1:4922555/150‑1 (MQ=255)

                                    aCTGGGTTTCGTTCTGCACAAAGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCgg                                                                                                                  >  2:5933042/1‑150 (MQ=255)

                                      tGGGTTTCGTTCTGCACAAAGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCa                                                                                                                <  1:6042462/150‑1 (MQ=255)

                                          tttCGTTCTGCACAAAGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCAc                                                                                                               >  1:2183426/1‑147 (MQ=255)

                                           ttCGTTCTGCACAAAGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCTGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGaa                                                                                                           <  2:781137/150‑1 (MQ=255)

                                           ttCGTTCTGCACAAAGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGaa                                                                                                           <  2:3693306/150‑1 (MQ=255)

                                           ttCGTTCTGCACAAAGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGaa                                                                                                           <  1:526723/150‑1 (MQ=255)

                                           ttCGTTCTGCACAAAGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGaa                                                                                                           <  1:7234189/150‑1 (MQ=255)

                                           ttCGTTCTGCACAAAGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGaa                                                                                                           <  1:2389206/150‑1 (MQ=255)

                                             cGTTCTGCACAAAGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGc                                                                                                         <  2:5051483/150‑1 (MQ=255)

                                               ttCTGCACAAAGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGAACAGCGCGt                                                                                                                                       <  2:4811915/118‑1 (MQ=255)

                                                tCTGCACAAAGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCAt                                                                                                      <  2:1023262/150‑1 (MQ=255)

                                                tCTGCACAAAGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCAt                                                                                                      <  2:4850372/150‑1 (MQ=255)

                                                 cTGCACAAAGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATa                                                                                                     <  1:5055456/150‑1 (MQ=255)

                                                 cTGCACAAAGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATa                                                                                                     >  1:1223430/1‑150 (MQ=255)

                                                  tGCACAAAGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAg                                                                                                    >  1:4886803/1‑150 (MQ=255)

                                                  tGCACAAAGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAg                                                                                                    <  2:6533625/150‑1 (MQ=255)

                                                   gCACAAAGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGAACAGCGCGTAACCGATGAAGTACAGCCCGGCACCGAAGCCATAgg                                                                                                   >  2:4291923/1‑150 (MQ=255)

                                                    cacaAAGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAACCATAGGc                                                                                                  >  1:6992084/1‑150 (MQ=255)

                                                     acaAAGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGcc                                                                                                 >  1:2530238/1‑150 (MQ=255)

                                                          ggCCATGGCCGCGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGAACAGCGCGTAACCGa                                                                                                                                 <  2:3047428/113‑1 (MQ=255)

                                                           gCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAgcg                                                                                           >  1:108061/1‑150 (MQ=255)

                                                            ccATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAgcgc                                                                                          >  1:3071071/1‑150 (MQ=255)

                                                              aTGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTg                                                                                        >  2:2537792/1‑150 (MQ=255)

                                                                ggCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGAt                                                                                      <  2:6913427/150‑1 (MQ=255)

                                                                 gCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATg                                                                                     <  1:3404126/150‑1 (MQ=255)

                                                                        ccaccaGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGt                                                                              <  1:4623092/150‑1 (MQ=255)

                                                                         caccaGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTc                                                                             <  1:3210543/150‑1 (MQ=255)

                                                                           ccaGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGc                                                                           >  1:2970081/1‑150 (MQ=255)

                                                                               cccccAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTc                                                                       <  2:4527485/150‑1 (MQ=255)

                                                                                 cccAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCa                                                                     <  1:1240233/150‑1 (MQ=255)

                                                                                  ccAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAg                                                                    >  1:6920843/1‑150 (MQ=255)

                                                                                   cAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAgg                                                                   >  2:3485715/1‑150 (MQ=255)

                                                                                    aGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGt                                                                  >  2:6318221/1‑150 (MQ=255)

                                                                                     ggTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGTg                                                                 <  2:2226487/150‑1 (MQ=255)

                                                                                           cATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGTgcgggcg                                                           <  1:6108045/150‑1 (MQ=255)

                                                                                           cATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGTgcgggcg                                                           <  1:1587836/150‑1 (MQ=255)

                                                                                           cATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGTgcgggcg                                                           >  2:171531/1‑150 (MQ=255)

                                                                                           cATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGTgcgggcg                                                           >  2:4406861/1‑150 (MQ=255)

                                                                                              gATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGTGCGGGCGGAc                                                        <  1:1144336/150‑1 (MQ=255)

                                                                                               aTACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGTGCGGGCGGACg                                                       <  1:4522949/150‑1 (MQ=255)

                                                                                                 aCGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGTGCGGGCGGACGaa                                                     <  2:1668576/150‑1 (MQ=255)

                                                                                                 aCGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGTGCGGGCGGACGaa                                                     >  1:185767/1‑150 (MQ=255)

                                                                                                 aCGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGTGCGGGCGGACGaa                                                     <  1:2723303/150‑1 (MQ=255)

                                                                                                      tCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGTGCGGGCGGACGAAGCCaa                                                >  2:5118459/1‑150 (MQ=255)

                                                                                                          ccACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGTGCGGGCGGACGAAGCCAATGtt                                            >  2:3909745/1‑150 (MQ=255)

                                                                                                          ccACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGTGCGGGCGGACGAAGCCAATGtt                                            <  1:711198/150‑1 (MQ=255)

                                                                                                            acaGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGTGCGGGCGGACGAAGCCAATGTTGa                                          >  2:493667/1‑150 (MQ=255)

                                                                                                             caGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGTGCGGGCGGACGAAGCCAATGTTGAc                                         <  2:5268832/150‑1 (MQ=255)

                                                                                                                  gggCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGTGCGGGCGGACGAAGCCAATGTTGACGCGGt                                    >  2:3801995/1‑150 (MQ=255)

                                                                                                                   ggCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGTGCGGGCGGACGAAGCCAATGTTGACGCGGTc                                   >  1:4664341/1‑150 (MQ=255)

                                                                                                                     cACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGTGCGGGCGGACGAAGCCAATGTTGACGCGGTCGa                                 >  2:4100823/1‑150 (MQ=255)

                                                                                                                     cACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGTGCGGGCGGACGAAGCCAATGTTGACGCGGTCGa                                 <  2:1402061/150‑1 (MQ=255)

                                                                                                                     cACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGTGCGGGCGGACGAAGCCAATGTTGACGCGGTCGa                                 >  2:2949783/1‑150 (MQ=255)

                                                                                                                      aCCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGTGCGGGCGGACGAAGCCAATGTTGACGCGGTCGAg                                <  2:5357003/150‑1 (MQ=255)

                                                                                                                        cGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGTGCGGGCGGACGAAGCCAATGTTGACGCGGTCGAGGt                              >  2:4216818/1‑150 (MQ=255)

                                                                                                                        cGACCTATAGCAGCAGCATGATCGAAGGTACATCCAACATCGCGTAACCGATGAAGAACAGCCCGGCAACGAAGACATAGGCCGCAGCCCACATGCCCAGGACGCATTCCAGGTGCGGGCGGACGAAGACAAAGTTGACGATCTCGAGGt                              >  1:379264/1‑150 (MQ=255)

                                                                                                                          aCCTTTTGCAGCAGCATGTTCGAAGGTACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGTGCGGGCGGACGAAGCCAATGTTGACGCGGTCGAGGTAg                            <  2:7211456/150‑1 (MQ=255)

                                                                                                                          aCCTTTTGCAGCAGCATGTTCGAAGGTACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGTGCGGGCGGACGAAGCCAATGTTGACGCGGTCGAGGTAg                            <  2:6554276/150‑1 (MQ=255)

                                                                                                                          aCCTTTTGCAGCAGCATGTTCGAAGGTACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGTGCGGGCGGACGAAGCCAATGTTGACGCGGTCGAGGTAg                            <  1:5933042/150‑1 (MQ=255)

                                                                                                                          aCCTTTTGCAGCAGCATGTTCGAAGGTACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGTGCGGGCGGACGAAGCCAATGTTGACGCGGTCGAGGTAg                            <  1:1716051/150‑1 (MQ=255)

                                                                                                                           ccTTTTTCAGCAGCATGTTCGAAGGTACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGTGCGGGCGGACGAAGCCAATGTTGACGCGGTCGAGGTAGt                           <  2:1271244/150‑1 (MQ=255)

                                                                                                                             ttttGCAGCAGCATGTTCGAAGGTACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGTGCGGGCGGACGAAGCCAATGTTGACGCGGTCGAGGTAGTTg                         >  2:200824/1‑150 (MQ=255)

                                                                                                                             ttttGCAGCAGCATGTTCGAAGGTACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGTGCGGGCGGACGAAGCCAATGTTGACGCGGTCGAGGTAGTTg                         >  2:3209228/1‑150 (MQ=255)

                                                                                                                             ttttGCAGCAGCATGTTCGAAGGTACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGTGCGGGCGGACGAAGCCAATGTTGACGCGGTCGAGGTAGTTg                         >  2:6750741/1‑150 (MQ=255)

                                                                                                                              tttGCAGCAGCATGTTCGAAGGTACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGTGCGGGCGGACGAAGCCAATGTTGACGCGGTCGAGGTAGTTGa                        <  1:652680/150‑1 (MQ=255)

                                                                                                                                tGCAGCAGCATGTTCGAAGGTACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGTGCGGGCGGACGAAGCCAATGTTGACGCGGTCGAGGTAGTTGACg                      <  1:4279334/150‑1 (MQ=255)

                                                                                                                                 gcagcagcaTGTTCGAAGGTACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGTGCGGGCGGACGAAGCCAATGTTGACGCGGTCGAGGTAGTTGACGa                     >  1:5099730/1‑150 (MQ=255)

                                                                                                                                 gcagcagcaTGTTCGAAGGTACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGTGCGGGCGGACGAAGCCAATGTTGACGCGGTCGAGGTAGTTGACGa                     <  2:3569073/150‑1 (MQ=255)

                                                                                                                                        gctGTTCGAAGTACCTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGTGCGGGCGGACGAAGCCAATGTTGACGCGGTCGAGGTAGTTGACGATGAACat              >  1:3874296/3‑150 (MQ=255)

                                                                                                                                          tGTTCGAAGTACCTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGTGCGGGCGGACGAAGCCAATGTTGACGCGGTCGAGGTAGTTGACGATGAACatga            >  2:313839/1‑150 (MQ=255)

                                                                                                                                          tGTTCGAAGGTACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCCTAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGTGCGGGCGGACGAAGCCAATGTTGACGCGGTCGAGGTAGTTGACGATGAACatga            >  2:2722873/1‑150 (MQ=255)

                                                                                                                                          tGTTCGAAGGTACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGTGCGGGCGGACGAAGCCAATGTTGACGCGGTCGAGGTAGTTGACGATGAACatga            >  1:6600257/1‑150 (MQ=255)

                                                                                                                                           gTTCGAAGTACCTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGTGCGGGCGGACGAAGCCAATGTTGACGCGGTCGAGGTAGTTGACGATGAACatgat           >  1:479929/1‑150 (MQ=255)

                                                                                                                                             tCGAAGGTACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGTGCGGGCGGACGAAGCCAATGTTGACGCGGTCGAGGTAGTTGACGATGAACAtgatga         <  1:4431387/150‑1 (MQ=255)

                                                                                                                                              cGAAGGTACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGTGCGGGCGGACGAAGCCAATGTTGACGCGGTCGAGGTAGTTGACGATGAACATGATGAc        <  1:7399923/150‑1 (MQ=255)

                                                                                                                                                   gTACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGTGCGGGCGGACGAAGCCAATGTTGACGCGGTCGAGGTAGTTGACGATGAACATGATGACGAACa   >  1:1341385/1‑150 (MQ=255)

CTTCGGCCACCCCGAGCAGAAAACGCAGCACGTAGAACTGGGTTTCGTTCTGCACAAAGGCCATGGCCGTGGCCACCAGCCCCCAGGTGAACATGATACGGGTCAGCCACAGCCGGGCACCGACCTTTTGCAGCAGCATGTTCGAAGGTACTTCGAACAGCGCGTAACCGATGAAGAACAGCCCGGCACCGAAGCCATAGGCCGCAGCGCTGATGCCCAGGTCGCTTTCCAGGTGCGGGCGGACGAAGCCAATGTTGACGCGGTCGAGGTAGTTGACGATGAACATGATGACGAACAG  >  NC_002947/5423668‑5423965

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