BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	955,085:1	+G	6.6%	intergenic (+36/‑277)	cyoE → / → alaC	protoheme IX farnesyltransferase/aminotransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	955,085	1	.	G	6.6%	585.7 / 34.6	137	intergenic (+36/‑277)	cyoE/alaC	protoheme IX farnesyltransferase/aminotransferase
Reads supporting (aligned to +/- strand): ref base . (36/92); new base G (2/7); total (38/99)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 137 total aligned reads displayed.

TGGGCCCGCAAGGTGTTCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCTGTAGGA‑GCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCT  >  NC_002947/954939‑955231

tGGGCCCGCAAGGTGTTCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑Gcc                                                                                                                                                    <  2:2070693/150‑1 (MQ=255)

  ggCCCGCAAGGTGTTCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCgg                                                                                                                                                  <  2:1435479/150‑1 (MQ=255)

   gCCCGCAAGGTGTTCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCggg                                                                                                                                                 <  1:786708/150‑1 (MQ=255)

   gCCCGCAAGGTGTTCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCggg                                                                                                                                                 <  1:3669958/150‑1 (MQ=255)

     ccGCAAGGTGTTCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑t                                                                                                                                              >  1:272561/1‑150 (MQ=255)

     ccGCAAGGTGTTCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑t                                                                                                                                              <  1:3342949/150‑1 (MQ=255)

      cGCAAGGTGTTCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGG‑CACAGTGGCCGGGG‑tt                                                                                                                                             <  2:2462807/150‑1 (MQ=255)

        cAAGGTGTTCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑tttt                                                                                                                                           <  1:1301559/150‑1 (MQ=255)

        cAAGGTGTTCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGG‑CACAGTGGCCGGGG‑tttt                                                                                                                                           <  2:3019714/150‑1 (MQ=255)

         aaGGTGTTCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑ttttt                                                                                                                                          <  1:2375462/150‑1 (MQ=255)

         aaGGTGTTCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑ttttt                                                                                                                                          <  1:4139813/150‑1 (MQ=255)

         aaGGTGTTCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGG‑CACAGTGGCCGGGG‑ttttc                                                                                                                                          <  2:643539/150‑2 (MQ=255)

         aaGGTGTTCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGG‑CACAGTGGCCGGGG‑ttttc                                                                                                                                          <  1:3438949/150‑2 (MQ=255)

           ggTGTTCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑ttttttt                                                                                                                                        >  1:2560407/1‑150 (MQ=255)

           ggTGTTCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGG‑CACAGTGGCCGGGG‑ttttcgt                                                                                                                                        <  1:1739923/150‑4 (MQ=255)

             tgtTCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑ttttttttt                                                                                                                                      >  2:4345819/1‑150 (MQ=255)

             tgtTCGGCTTCTCCATCCTCACCGTCACGGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCGGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑ttttttttt                                                                                                                                      <  2:286482/150‑1 (MQ=255)

              gtTCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTg                                                                                                                                     <  1:235877/150‑1 (MQ=255)

              gtTCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTg                                                                                                                                     <  2:3409219/150‑1 (MQ=255)

              gtTCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTg                                                                                                                                     <  1:3283938/150‑1 (MQ=255)

              gtTCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTg                                                                                                                                     <  2:3409125/150‑1 (MQ=255)

              gtTCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTg                                                                                                                                     <  2:2383051/150‑1 (MQ=255)

               ttCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTgg                                                                                                                                    <  1:2518310/150‑1 (MQ=255)

               ttCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTgg                                                                                                                                    <  1:2812879/150‑1 (MQ=255)

               ttCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTgg                                                                                                                                    <  2:2973894/150‑1 (MQ=255)

               ttCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGG‑CACAGTGGCCGGGG‑tt                                                                                                                                             <  2:1316148/141‑1 (MQ=255)

               gtCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCGGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTgg                                                                                                                                    <  2:3818914/149‑1 (MQ=255)

                tCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTggg                                                                                                                                   <  1:4146347/150‑1 (MQ=255)

                tCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTggg                                                                                                                                   <  1:3089424/150‑1 (MQ=255)

                     ttCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGACACTGCTGCCTGTTTCACGAAAACCCCGGCTACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTAt                                                                                                                              >  1:3075433/1‑150 (MQ=255)

                     ttCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCTACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTAt                                                                                                                              >  1:3075426/1‑150 (MQ=255)

                     ttCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTAt                                                                                                                              >  1:2635085/1‑150 (MQ=255)

                             ccTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCtg                                                                                                                      <  2:1511969/150‑1 (MQ=255)

                               tCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCtgtg                                                                                                                    >  1:1877266/1‑150 (MQ=255)

                                  ccGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCg                                                                                                                 <  2:1748906/150‑1 (MQ=255)

                                   cGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTgggg                                                                                                                                  <  1:1076085/132‑1 (MQ=255)

                                       aCTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTACTGTGCTGGCctc                                                                                                            >  1:2415202/1‑150 (MQ=255)

                                        cTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCctct                                                                                                           >  2:2702163/1‑150 (MQ=255)

                                            ccTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTcgc                                                                                                       <  1:2422145/150‑1 (MQ=255)

                                             cTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTcgc‑g                                                                                                     <  1:3686209/150‑1 (MQ=255)

                                              tGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑gg                                                                                                    <  1:4408450/150‑1 (MQ=255)

                                               gAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑ggg                                                                                                   <  1:3411251/150‑1 (MQ=255)

                                               gAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑ggg                                                                                                   <  2:131095/150‑1 (MQ=255)

                                               gAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑ggg                                                                                                   <  2:131173/150‑1 (MQ=255)

                                                      atgatgGGGGTGGACAGCCAGACCGCTGCGGATGTGCGGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAAccc                                                                                            <  2:1829358/150‑1 (MQ=255)

                                                        gatgGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGc                                                                                          <  2:24304/150‑1 (MQ=255)

                                                         atgGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTGGC‑GGGTTAACCCGcg                                                                                         >  1:3589542/1‑149 (MQ=255)

                                                          tgGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTc                                                                                        <  2:407320/150‑1 (MQ=255)

                                                            ggTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑‑TTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTccga                                                                                     >  2:3998760/1‑148 (MQ=255)

                                                             gtgtGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCtc                                                                                     <  1:3352901/150‑2 (MQ=255)

                                                                tgGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACaa                                                                                  <  2:2543490/150‑1 (MQ=255)

                                                                 ggACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAg                                                                                 <  2:2268638/150‑1 (MQ=255)

                                                                   aCAGCCAGACCGCTGCGGTTGTGCTGTTGAACTTCGCACGCTGAAACTGCTGCCTTTTTCACGAAAACCCGGGCAACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAAACCGCTCCTACAAGGt                                                                               <  1:155555/150‑1 (MQ=38)

                                                                     aGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAATCCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTc                                                                             <  2:3263366/150‑1 (MQ=255)

                                                                     aGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTc                                                                             <  2:2749605/150‑1 (MQ=255)

                                                                     aGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTc                                                                             <  2:3263324/150‑1 (MQ=255)

                                                                     aGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCGCTTCGC‑GGGGGAACCCGCGCCGACCAgggtg                                                                             >  1:1735055/1‑147 (MQ=255)

                                                                      gCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAAACCGCTCCTACAAGGttac                                                                            >  2:2790275/1‑148 (MQ=255)

                                                                       ccAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCt                                                                           <  1:2918796/150‑1 (MQ=255)

                                                                        cAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCTg                                                                          <  1:2852808/150‑1 (MQ=255)

                                                                        cAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAAACCGCTCCAACAAGGTTccgg                                                                          >  2:3937322/1‑148 (MQ=255)

                                                                         aGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAAGTGCTGCCTGTTTCACGAAAACCCCGGCCTCTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCGGGGCTCTTCGC‑GGGTAAACCCGCTCCTACCCGGTTCgtgt                                                                         <  2:494172/150‑1 (MQ=255)

                                                                         aGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCTGt                                                                         <  2:2839473/150‑1 (MQ=255)

                                                                            ccGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCTGTAgg                                                                      <  2:1485955/150‑1 (MQ=255)

                                                                             cGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCTGTAGGa                                                                     <  2:3010037/150‑1 (MQ=255)

                                                                             cGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCTGTAGGa                                                                     <  2:4312527/150‑1 (MQ=255)

                                                                              gCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CTGTGCTGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAGGAACCTTGTAGGA‑Gc                                                                  <  1:1886367/150‑1 (MQ=25)

                                                                               cTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCTGTAGGA‑Gc                                                                  <  1:1158589/150‑1 (MQ=255)

                                                                               cTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCTGTAGGA‑Gc                                                                  <  1:475127/150‑1 (MQ=255)

                                                                               cTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCTGTAGGA‑Gc                                                                  <  1:475235/150‑1 (MQ=255)

                                                                                tGCGTATGAGCTGATGACCTACGCACGATGAAACTGCTGCCTGTTTCACTAAAACCCCGCCCACTGT‑GCCGGGG‑TATTTTATTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCTGTAGGA‑GCg                                                                 <  2:2569575/150‑1 (MQ=38)

                                                                                tGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CTGTGCTGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAGGAACCTTGTAGGA‑GCgg                                                                <  1:3794698/150‑1 (MQ=25)

                                                                                 gCGGTTGTGCTGATGCCCTACGCACGCTGAAACTGCGGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGT‑TTTTTTTTTGGGGGTATTTTTTCTGTGCTGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAGGAACCTTGTAGGA‑GCgg                                                                <  1:4326311/150‑1 (MQ=25)

                                                                                 gCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCTGTAGGA‑GCgg                                                                <  2:3590650/150‑1 (MQ=255)

                                                                                 gCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCTGTAGGA‑GCgg                                                                <  1:3296377/150‑1 (MQ=255)

                                                                                  cGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCTGTAGGA‑GCggg                                                               <  1:4274752/150‑1 (MQ=255)

                                                                                  cGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCTGTAGGA‑GCggg                                                               <  1:2237250/150‑1 (MQ=255)

                                                                                   ggATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCTGTAGGA‑GCGGGt                                                              <  2:613526/150‑1 (MQ=255)

                                                                                         gCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGACAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCTGTAGGA‑GCGGGTTTAccc                                                        <  2:2231241/150‑1 (MQ=255)

                                                                                           tgatgaCCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAGGAACCTTGTAGGA‑GCGGGTTTACCcgc                                                      <  2:272561/150‑1 (MQ=37)

                                                                                           tgatgaCCGACGCACGCTGAAACTGCTGCCTGTATCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TATTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCTGTAGGA‑GCGGGTTTACCcgc                                                      <  2:3358492/150‑1 (MQ=255)

                                                                                            gatgaCCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAAACCGCGCCGACAAGGTTCCGGTAGGA‑GCGGGGTTTACcgcg                                                     >  2:3832093/1‑150 (MQ=255)

                                                                                                 ccTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTGTGGGGGTATTTTTCCAGTGCCGGCCTCTTCGC‑GGGGAAAACCCCAACAACGAAGTTACGGGAGGG‑GCGGGTGTTCCcgcgtaaag                                                >  1:2252846/1‑145 (MQ=255)

                                                                                                  cTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTGAACCCCCGCCGACAAGGTTCCGGTAGGT‑GCGGGGTTTCCCGCGGAGAgg                                               >  1:714139/1‑150 (MQ=255)

                                                                                                  cTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCGGTAGGG‑GCGGGTTTACCCGCGAAGAgg                                               >  2:713829/1‑150 (MQ=255)

                                                                                                             gAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAAACCGCGCCTACAAGGTTCCTGGAGGA‑GCGGGTTTACCCGCGCAGAGGCCAGCGCCGaa                                    >  2:1746550/1‑150 (MQ=38)

                                                                                                                aCTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTGTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCTGTAGGA‑GCGGGTTTACCCGCGAAGAGGCCAGCACAGaaaaa                                 <  2:2635085/150‑1 (MQ=255)

                                                                                                                aCTGCTGCCTGTCGCACGAAACCCCCGGCCAGTGT‑GCGGGGT‑TTTTTTTTGGGGGGTATTTTTCGGGGGCCGGCCTCTGGGC‑GGGTAAACCCGCTCCTCCCCGGTTCCTGTAGGA‑GCGGGTTTACCCGCGAAGAGGCCAGCACAGaaaaa                                 <  2:3810503/150‑1 (MQ=25)

                                                                                                                     tgcCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGGTAACCCCCTCCGACAAGGTTACTGTAGGA‑GCGGGTTTTCCCGCGGAGAGGCCAGCCCCGacaaataaaa                            >  2:245371/1‑150 (MQ=255)

                                                                                                                       ccTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTAGC‑GGGTTAACACCCGCCGACAAGGTTGCAGGTGGAGGGGGGTTTACCC‑CGAAGGGGGCCGCACACgcaaataaaaat                          >  1:1793927/1‑150 (MQ=255)

                                                                                                                       ccTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTATTCGC‑GGGTAAACCCGCTCCTACCAGGTTCCGGTAGGA‑GCGGGTTTACCCGCGTAGTGGCCAGCAAAGaaaaataaaaat                          >  2:364766/1‑150 (MQ=38)

                                                                                                                         tGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTAGCGGGGTATACCC‑CGCCCACCACGTTGGGGGTGGG‑GGGGGGGTTCCCGCGGAGAGGCCCGCCCAGcaaaaaaaaaaaaa                        >  1:3129213/1‑147 (MQ=255)

                                                                                                                         tGTTTCACGAAAACCCCGGCCACGGT‑GCCGGGG‑TTTTTTTTGGGGGGTATTTTTCCGGTGCCGGCCTCTCCGC‑GGGTAAACCCGCCCCACCCCGGTTCCTGTAGGA‑GCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATAc                        <  2:856846/150‑1 (MQ=37)

                                                                                                                           tttCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAAACCGCACCTACAAGGTTGCTGTATGA‑GCGGGTTTTCCCGCGAAGAGGCCAGCAAAGAAAAATAAAAATACCt                      >  2:4223411/1‑150 (MQ=38)

                                                                                                                           ccaaaaaaaaaCCCCGG‑CACAGTGGCCGGGGTTTTTTTTTTGGTGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCGCCTTCAAGGTTCCGGTAGGG‑GCGGGGTTACCCGCGAAGAGGCCAGCACAGGAAAATAAAAATAcc                       >  2:2233506/8‑150 (MQ=21)

                                                                                                                             tCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCACCTACAAAGTTCCGGTAGGA‑GCGGGTGTACCCCCCCACAGGCCCGCAACGAAAAATAAAAATAccagg                    >  2:789220/1‑147 (MQ=37)

                                                                                                                               aCGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAAACCCCTCCGACCAGGTTGCGGGTGGG‑GCGGGGTTAACCGCGCAGCGGCCAGCACCGAAAAATAAAAAAACCTGCtt                  >  2:4159108/1‑150 (MQ=255)

                                                                                                                                    aaCCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCGTCGC‑GGGGAAACCCGCACCCACAAGGGTCCGGGAGGC‑GCGGGGGTACCCGCGCAGAGGCCCCCCCCGCAAAAAAAAAAAACCGGATTTTCaa             >  2:1956709/1‑150 (MQ=255)

                                                                                                                                      ccccGGCCCCTGT‑GCGGGTG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACACGGTTCCTTTTGGA‑GCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCaaaa           <  1:734189/150‑1 (MQ=37)

                                                                                                                                               cTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCTGTAGGA‑GCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCt  <  2:3381908/150‑1 (MQ=255)

TGGGCCCGCAAGGTGTTCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCTGTAGGA‑GCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCT  >  NC_002947/954939‑955231

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