BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,578,013	T→C	6.7%	H148H (CAT→CAC)	ybhP‑II →	phosphohydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,578,013	0	T	C	6.7%	407.8 / 12.7	133	H148H (CAT→CAC)	ybhP‑II	phosphohydrolase
Reads supporting (aligned to +/- strand): ref base T (61/63); new base C (3/6); total (64/69)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.82e-01

Only 100 of 133 total aligned reads displayed.

GTCTACCCTGACGGCCACCACGGCAATGCGTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGAAAAGCCATGGCAACCGCACCC  >  NC_002947/4577867‑4578161

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

      ccTGACGGCCACCACGGCAATGCGTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGc                                                                                                                                             >  1:1860106/1‑150 (MQ=255)

      ccTGACGGCCACCACGGCAATGCGTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGc                                                                                                                                             >  2:1267594/1‑150 (MQ=255)

        tGACGGCCACCACGGCAATGCGTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGt                                                                                                                                           <  1:1267594/150‑1 (MQ=255)

        tGACGGCCACCACGGCAATGCGTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGt                                                                                                                                           >  1:4135022/1‑150 (MQ=255)

        tGACGGCCACCACGGCAATGCGTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGt                                                                                                                                           <  1:2257393/150‑1 (MQ=255)

        tGACGGCCACCACGGCAATGCGTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGt                                                                                                                                           <  1:2107320/150‑1 (MQ=255)

        tGACGGCCACCACGGCAATGCGTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGACAGGTGCATGCGGTCTGCGt                                                                                                                                           >  2:1847815/1‑150 (MQ=255)

           cGGCCACCACGGCAATGCGTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCACGCAGACCGCATGCa                                                                                                                                        <  1:2639155/150‑1 (MQ=255)

            ggCCACCACGGCAATGCGTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGGCGGGGCAAGGCCAGGTGCACGCAGACCGCATGCAc                                                                                                                                       <  2:2751776/150‑1 (MQ=255)

            ggCCACCACGGCAATGCGTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCACGCAGACCGCATGCAc                                                                                                                                       <  2:2751529/150‑1 (MQ=255)

            ggCCACCACGGCAATGCGTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGACAGGTGCATGCGGTCTGCGTGCAc                                                                                                                                       <  1:1847815/150‑1 (MQ=255)

            ggCCACCACGGCAATGCGTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGAGCAACGACAGGTGCACGCAGACCGCATGCAc                                                                                                                                       <  1:1074244/150‑1 (MQ=255)

                 ccacGGCAATGCGTTGCTGTCCAAGTACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCACGCAGACCGCATGCACCTggc                                                                                                                                  <  2:1387278/150‑3 (MQ=255)

                 ccacGGCAATGCGTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTc                                                                                                                                  >  2:2097665/1‑150 (MQ=255)

                 ccacGGCAATGCGTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCATGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTc                                                                                                                                  >  1:4235061/1‑150 (MQ=255)

                  cacGGCAATGCGTTGCTGTCCAAATACCCGCTCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCg                                                                                                                                 <  2:1282101/150‑1 (MQ=255)

                  cacGGCAATGCGTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCg                                                                                                                                 <  2:1282000/150‑1 (MQ=255)

                  cacGGCAATGCGTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCg                                                                                                                                 <  2:1736882/150‑1 (MQ=255)

                    cGGCAATGCGTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGtt                                                                                                                               <  1:4265462/150‑1 (MQ=255)

                    cGGCAATGCGTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGtt                                                                                                                               <  2:21672/150‑1 (MQ=255)

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

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

                        aaTGCGTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTg                                                                                                                           >  2:3416481/1‑150 (MQ=255)

                           gCGTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGaa                                                                                                                        <  1:1995363/150‑1 (MQ=255)

                           gCGTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGaa                                                                                                                        >  2:2996462/1‑150 (MQ=255)

                           gCGTTGCTGTCCAAATACACGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGaa                                                                                                                        <  2:728932/150‑1 (MQ=255)

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

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

                                 cTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCAc                                                                                                                  <  2:116446/150‑1 (MQ=255)

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

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

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

                                       aaaTACCCAATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCaa                                                                                                            >  1:4654989/1‑150 (MQ=255)

                                           aCCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAagc                                                                                                        <  1:4894359/150‑1 (MQ=255)

                                            cccGATCATCGAGCACCGCAACCTCGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAagca                                                                                                       <  1:1690879/150‑1 (MQ=255)

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

                                               gATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCt                                                                                                    >  1:3690765/1‑150 (MQ=255)

                                                   atcGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGTGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAg                                                                                                >  2:2886341/1‑150 (MQ=255)

                                                     cGAGCACCGGAACCTGGTCGTGTCGATCACCGGCCCCGTACGCCGGGGCCTGGTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGtt                                                                                              <  2:3288000/150‑1 (MQ=255)

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

                                                                  cTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAActg                                                                                 >  2:2797514/1‑150 (MQ=255)

                                                                  cTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAActg                                                                                 >  2:3395898/1‑150 (MQ=255)

                                                                      aCGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTgg                                                                             <  1:2558977/150‑1 (MQ=255)

                                                                      aCGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTgg                                                                             <  2:3194241/150‑1 (MQ=255)

                                                                      aCGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTgg                                                                             >  1:1871331/1‑150 (MQ=255)

                                                                        gtgtCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGag                                                                           >  2:464808/1‑149 (MQ=255)

                                                                            cGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCCTTGCTGGAAAGCCACCGGCAAAACCAGCTGCAGTTGCTGCGTAAACTGCTGGACGccc                                                                       <  2:4135022/150‑1 (MQ=255)

                                                                                 aCCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCg                                                                  >  1:522893/1‑150 (MQ=255)

                                                                                 aCCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCg                                                                  <  1:4281417/150‑1 (MQ=255)

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

                                                                                     gCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCa                                                              >  1:1250241/1‑150 (MQ=255)

                                                                                           aaCGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGccc                                                        <  1:4703851/150‑1 (MQ=255)

                                                                                             cgccgcGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCg                                                      <  2:3502573/150‑1 (MQ=255)

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

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

                                                                                                      cTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGatcatc                                             >  2:246292/1‑150 (MQ=255)

                                                                                                      cTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGatcatc                                             >  2:306144/1‑150 (MQ=255)

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

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

                                                                                                         ttGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGcc                                          >  2:403801/1‑150 (MQ=255)

                                                                                                         ttGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGcc                                          >  1:3734343/1‑150 (MQ=255)

                                                                                                         ttGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCACCGGTGATCATCgtc                                          >  1:2513808/1‑148 (MQ=255)

                                                                                                         ttGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGATGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGcc                                          >  1:2058436/1‑150 (MQ=255)

                                                                                                            cACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGc                                       >  1:3023464/1‑150 (MQ=255)

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

                                                                                                               tgtgtgCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGAc                                    <  2:5049393/150‑1 (MQ=255)

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

                                                                                                               tgtgtgCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGAc                                    <  2:2876205/150‑1 (MQ=255)

                                                                                                                 tgtgCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACtt                                  >  1:5058177/1‑150 (MQ=255)

                                                                                                                 tgtgCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACtt                                  <  1:3416481/150‑1 (MQ=255)

                                                                                                                 tgtgCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACtt                                  >  1:5057745/1‑150 (MQ=255)

                                                                                                                  gtgCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCTCCGGCGACTTc                                 <  2:4235061/150‑1 (MQ=255)

                                                                                                                  gtgCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTc                                 <  2:4235377/150‑1 (MQ=255)

                                                                                                                   tgCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCa                                >  2:2417539/1‑150 (MQ=255)

                                                                                                                     cTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAAc                              <  1:3705234/150‑1 (MQ=255)

                                                                                                                     cTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAAc                              <  1:3919762/150‑1 (MQ=255)

                                                                                                                         tggTGCCGGGTCAAGGCCAGGTGAATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGTCAAAAGCAGCTGCAGTTGATGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACt                          <  1:2537041/149‑1 (MQ=255)

                                                                                                                               cGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGaaaa                    <  1:859584/150‑1 (MQ=255)

                                                                                                                               cGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGaaaa                    <  1:859608/150‑1 (MQ=255)

                                                                                                                               cGGGGCAAGGCCAGGTGAATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACTCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGaaaa                    <  2:4694656/150‑1 (MQ=255)

                                                                                                                               cGGGGCAAGGCCAGGGGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGaaaa                    <  2:2585965/150‑1 (MQ=255)

                                                                                                                                 gggCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGAAAAGc                  <  1:2097665/150‑1 (MQ=255)

                                                                                                                                   gCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGAAAAGCCa                >  2:1695464/1‑150 (MQ=255)

                                                                                                                                   gCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGAAAAGCCa                >  1:3913992/1‑150 (MQ=255)

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

                                                                                                                                     aaGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGAAAAGCCATg              <  2:1290476/150‑1 (MQ=255)

                                                                                                                                     aaGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGAAAAGCCATg              >  2:3231084/1‑150 (MQ=255)

                                                                                                                                     aaGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGAAAAGCCATg              >  2:447022/1‑150 (MQ=255)

                                                                                                                                          cAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGAAAAGCCATGGCAAc         >  2:813630/1‑150 (MQ=255)

                                                                                                                                          cAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGAAAAGCCATGGCAAc         >  1:4560050/1‑150 (MQ=255)

                                                                                                                                          cAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGAAAAGCCATGGCAAc         >  1:4560040/1‑150 (MQ=255)

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

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

                                                                                                                                           aGGTGCACGCAGACCGCATGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGAAAAGCCATGGCAAcc        >  1:4187966/1‑150 (MQ=255)

                                                                                                                                            ggTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGAAAAGCCATGGCAACCg       >  1:1724650/1‑150 (MQ=255)

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

                                                                                                                                             gTGCACGCAGACCGCATGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGa                                                  >  2:1974684/1‑106 (MQ=255)

                                                                                                                                             gTGCACGCAGACCGCATGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGa                                                  <  1:1974684/106‑1 (MQ=255)

                                                                                                                                                 aTGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGAAAAGCCATGGCAACCGCAccc  >  2:4794463/1‑150 (MQ=255)

GTCTACCCTGACGGCCACCACGGCAATGCGTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGAAAAGCCATGGCAACCGCACCC  >  NC_002947/4577867‑4578161

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 11 ≤` `ATCG/ATCG` `< 37 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`