BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,581,045	T→A	5.9%	V150D (GTC→GAC)	ybhP‑II →	phosphohydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,581,045	0	T	A	5.9%	909.2 / 29.2	268	V150D (GTC→GAC)	ybhP‑II	phosphohydrolase
Reads supporting (aligned to +/- strand): ref base T (123/129); new base A (9/7); total (132/136)
Fisher's exact test for biased strand distribution p-value = 6.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 270 total aligned reads displayed.

CTACCCTGACGGCCACCACGGCAATGCGTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGAAAAGCCATGGCAACCGCACCCTTGGC  >  NC_002947/4580896‑4581193

  aCCCTGACGGCCACCACGGCAATGCGTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCt                                                                                                                                                    <  1:6039502/150‑1 (MQ=255)

   cccTGACGGCCACCACGGCAATGCGTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTg                                                                                                                                                   >  1:141453/1‑150 (MQ=255)

    ccTGACGGCCACCACGGCAATGCGTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGc                                                                                                                                                  <  1:2167786/150‑1 (MQ=255)

    ccTGACGGCCACCACGGCAATGCGTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGc                                                                                                                                                  <  2:3963954/150‑1 (MQ=255)

     cTGACGGCCACCACGGCAATGCGTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCg                                                                                                                                                 >  1:4473927/1‑150 (MQ=255)

         cGGCCACCACGGCAATGCGTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCa                                                                                                                                             >  2:2214653/1‑150 (MQ=255)

          ggCCACCACGGCAATGCGTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCAc                                                                                                                                            <  1:4326212/150‑1 (MQ=255)

          ggCCACCACGGCAATGCGTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGACAGGTGCACGCAGACCGCATGCAc                                                                                                                                            <  2:3116668/150‑1 (MQ=255)

            ccaccaTGGCAATGCGTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCt                                                                                                                                          <  2:2900753/150‑1 (MQ=255)

            ccaccaCGGCAATGCGTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCt                                                                                                                                          >  2:4738818/1‑150 (MQ=255)

            ccaccaCGGCAATGCGTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCt                                                                                                                                          <  1:1536489/150‑1 (MQ=255)

             caccacGGCAATGCGTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCACGCAGACCGCATGCACCTg                                                                                                                                         <  2:2014302/150‑1 (MQ=255)

              accacGGCAATGCGTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGt                                                                                                                                        >  1:40782/1‑150 (MQ=255)

              accacGGCAATGCGTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGt                                                                                                                                        >  1:5761370/1‑150 (MQ=255)

                   ggCAATGCGTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTg                                                                                                                                   <  1:3578063/150‑1 (MQ=255)

                    gCAATGCGTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGc                                                                                                                                  >  1:2965941/1‑150 (MQ=255)

                     cAATGCGTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCt                                                                                                                                 <  1:5771386/150‑1 (MQ=255)

                      aaTGCGTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTg                                                                                                                                <  2:5287748/150‑1 (MQ=255)

                      aaTGCGTTGCTGTCCAAATACCCGAGCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCGGGTGCACTGTGGGCTGGAGGGGCCGGGGCCAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTg                                                                                                                                >  2:6680757/1‑150 (MQ=255)

                          cGTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGaaa                                                                                                                            <  2:6350359/150‑1 (MQ=255)

                          cGTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGaaa                                                                                                                            <  1:2780116/150‑1 (MQ=255)

                          cGTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGaaa                                                                                                                            <  1:4010009/150‑1 (MQ=255)

                          cGTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGaaa                                                                                                                            <  1:3778482/150‑1 (MQ=255)

                          cGTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGaaa                                                                                                                            <  2:2606555/150‑1 (MQ=255)

                             tGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGcc                                                                                                                         <  2:4012805/150‑1 (MQ=255)

                                tGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCAcc                                                                                                                      >  2:2781520/1‑150 (MQ=255)

                                 gTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCg                                                                                                                     <  1:2068372/150‑1 (MQ=255)

                                  tCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCgg                                                                                                                    <  2:1902088/150‑1 (MQ=255)

                                         aCCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAagc                                                                                                             >  2:5617370/1‑150 (MQ=255)

                                          cccGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAagca                                                                                                            <  1:704095/150‑1 (MQ=255)

                                           ccGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAagcag                                                                                                           <  2:2554405/150‑1 (MQ=255)

                                           ccGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAagcag                                                                                                           >  2:6422584/1‑150 (MQ=255)

                                           ccGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAagcag                                                                                                           <  2:1999817/150‑1 (MQ=255)

                                           ccGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAagcag                                                                                                           >  1:2326855/1‑150 (MQ=255)

                                            cGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAagcagc                                                                                                          >  2:6346984/1‑150 (MQ=255)

                                            cGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAagcagc                                                                                                          <  1:2200961/150‑1 (MQ=255)

                                                  tcGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGt                                                                                                    >  2:3789224/1‑150 (MQ=255)

                                                   cGAGCACCGCAACCTTGACGTGTCGATCACCGGCCCCGAACTCCGCGGCCTGTTGCACTGTGTTCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTTTCGTTGCTTGAAAGCCACCGGAAAAAGCAGCTGAAGtt                                                                                                   <  2:1322359/150‑1 (MQ=255)

                                                   cGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGtt                                                                                                   >  2:3263939/1‑150 (MQ=255)

                                                   cGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGtt                                                                                                   <  2:4599197/150‑1 (MQ=255)

                                                       cACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTtgctg                                                                                               <  2:5037763/150‑1 (MQ=255)

                                                        aCCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTtgctgc                                                                                              <  1:4838090/150‑1 (MQ=255)

                                                        aCCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTtgctgc                                                                                              >  2:6039329/1‑150 (MQ=255)

                                                        aCCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTtgctgc                                                                                              <  1:6852061/150‑1 (MQ=255)

                                                          cGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGt                                                                                            >  2:4769198/1‑150 (MQ=255)

                                                          cGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGt                                                                                            >  2:3037706/1‑150 (MQ=255)

                                                           gCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTa                                                                                           <  2:913715/150‑1 (MQ=255)

                                                             aaCCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTaaa                                                                                         <  2:6176753/150‑1 (MQ=255)

                                                                       tgtCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACg                                                                               >  1:6021573/1‑150 (MQ=255)

                                                                        gtCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGc                                                                              >  2:1786189/1‑150 (MQ=255)

                                                                        gtCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGc                                                                              >  1:3715003/1‑150 (MQ=255)

                                                                           gATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCt                                                                           >  2:211691/1‑150 (MQ=255)

                                                                           gATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCt                                                                           >  1:4617570/1‑150 (MQ=255)

                                                                           gATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCt                                                                           >  1:178383/1‑150 (MQ=255)

                                                                              cACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCTGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGcc                                                                        <  1:645182/150‑1 (MQ=255)

                                                                                 cGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGc                                                                     >  2:5092726/1‑150 (MQ=255)

                                                                                 cGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGc                                                                     >  2:2390378/1‑150 (MQ=255)

                                                                                    ccccGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCaa                                                                  <  1:6394187/150‑1 (MQ=255)

                                                                                     cccGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAAc                                                                 >  1:3636221/1‑150 (MQ=255)

                                                                                            gccgcGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCgg                                                          >  2:6104892/1‑150 (MQ=255)

                                                                                             ccgcGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGt                                                         <  2:1700781/150‑1 (MQ=255)

                                                                                             ccgcGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGt                                                         >  1:6700388/1‑150 (MQ=255)

                                                                                               gcgGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGa                                                       >  1:2642887/1‑150 (MQ=255)

                                                                                                cgGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGat                                                      <  2:5957360/150‑1 (MQ=255)

                                                                                                 ggCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGCGatc                                                     <  1:25048/150‑1 (MQ=255)

                                                                                                    cTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGatcatc                                                  <  2:5741623/150‑1 (MQ=255)

                                                                                                    cTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGatcatc                                                  <  2:6336280/150‑1 (MQ=255)

                                                                                                     tGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCg                                                 >  1:3926381/1‑150 (MQ=255)

                                                                                                         gCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCgg                                             <  1:3141059/150‑1 (MQ=255)

                                                                                                         gCACCGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCgg                                             <  1:6544882/150‑1 (MQ=255)

                                                                                                           aCTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCg                                           >  2:421812/1‑150 (MQ=255)

                                                                                                            cTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGAGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCAGCAGTTGATGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGAGATCATCGCCGGCGa                                          >  2:4160429/1‑150 (MQ=255)

                                                                                                             tgtgtgCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGAc                                         >  1:155785/1‑150 (MQ=255)

                                                                                                                  gCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCAc                                                                                                                       <  1:1178369/67‑1 (MQ=255)

                                                                                                                  gCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCAc                                                                                                                       >  2:1178369/1‑67 (MQ=255)

                                                                                                                    tGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACg                                  >  1:2887912/1‑150 (MQ=255)

                                                                                                                        ggTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTg                              <  2:78802/150‑1 (MQ=255)

                                                                                                                           tttccagCAACGACAGGTGCACGCAGACCGCATGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGaa                           >  1:6678781/7‑150 (MQ=38)

                                                                                                                           gCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGaa                           <  2:357989/150‑1 (MQ=255)

                                                                                                                           gCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGaa                           >  2:6022084/1‑150 (MQ=255)

                                                                                                                            ccGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGaaa                          >  1:3590747/1‑150 (MQ=255)

                                                                                                                              ggggCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGCAAAg                        >  1:1699639/1‑150 (MQ=255)

                                                                                                                               cagCAACGACAGGTGCACGCAGACCGCATGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGAAAAGc                       >  2:728223/3‑150 (MQ=38)

                                                                                                                               cagCAACGACAGGTGCACGCAGACCGCATGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGAAAAGc                       >  2:717152/3‑150 (MQ=38)

                                                                                                                                  cAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGAAAAGCCAt                    <  2:1525922/150‑1 (MQ=255)

                                                                                                                                   aaGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGAACAGCCATg                   >  2:1124872/1‑150 (MQ=255)

                                                                                                                                   aaGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGACAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGAAAAGCCATg                   <  2:1836576/150‑1 (MQ=255)

                                                                                                                                     ggCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCAGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGAAAAGCCATGGc                 >  1:3441172/1‑150 (MQ=255)

                                                                                                                                        cAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGAAAAGCCATGGCAAc              <  2:6695383/150‑1 (MQ=255)

                                                                                                                                         aGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGAAAAGCCATGGCAAcc             <  1:5395063/150‑1 (MQ=255)

                                                                                                                                           gTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGAAAAGCCATGGCAACCGc           >  2:5905833/1‑150 (MQ=255)

                                                                                                                                           gTGCACGCAGACCGCATGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGAAAAGCCATGGCAACCGc           >  2:6559031/1‑150 (MQ=255)

                                                                                                                                           gTGCACGCAGACCGCATGCACCTGGCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGAAAAGCCATGGCAACCGc           >  1:746520/1‑150 (MQ=255)

                                                                                                                                             gCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGAAAAGCCATGGCAACCGCAc         >  1:1001665/1‑150 (MQ=255)

                                                                                                                                                tGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGAAAAGCCATGGCAACCGCACCCt      >  2:3399270/1‑150 (MQ=255)

                                                                                                                                                 gCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGAAAAGCCATGGCAACCGCACCCtt     <  1:1887516/150‑1 (MQ=255)

                                                                                                                                                 gCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGAAAAGCCATGGCAACCGCACCCtt     >  2:1265462/1‑150 (MQ=255)

                                                                                                                                                    gTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGAAAAGCCATGGCAACCGCACCCTTGGc  <  1:1310324/150‑1 (MQ=255)

                                                                                                                                                    gTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGAAAAGCCATGGCAACCGCACCCTTGGc  <  1:937203/150‑1 (MQ=255)

                                                                                                                                                    gTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGAAAAGCCATGGCAACCGCACCCTTGGc  <  2:141453/150‑1 (MQ=255)

CTACCCTGACGGCCACCACGGCAATGCGTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGAAAAGCCATGGCAACCGCACCCTTGGC  >  NC_002947/4580896‑4581193

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