BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	955,155	T→A	9.1%	intergenic (+106/‑207)	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,155	0	T	A	9.1%	270.2 / 17.9	99	intergenic (+106/‑207)	cyoE/alaC	protoheme IX farnesyltransferase/aminotransferase
Reads supporting (aligned to +/- strand): ref base T (42/45); new base A (3/6); total (48/51)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.85e-01

Only 100 of 101 total aligned reads displayed.

ACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCTGTAGGA‑GCGGGTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCAGGC  >  NC_002947/955006‑955302

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

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

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

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

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

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

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

     cAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCTg                                                                                                                                                  <  1:2561845/150‑1 (MQ=255)

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

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

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

      aGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCTGt                                                                                                                                                 <  1:719996/150‑1 (MQ=255)

        aCCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAAACCGCTCCGACAAGGTTGCTGTAg                                                                                                                                               >  2:4261689/1‑150 (MQ=255)

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

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

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

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

            cTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAGGAACCTTGTAGGA‑Gc                                                                                                                                          <  2:246154/150‑1 (MQ=38)

            cTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAGGAACCTTGTAGGA‑Gc                                                                                                                                          <  2:2903532/150‑1 (MQ=38)

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

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

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

            cTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCTGTAGGA‑Gc                                                                                                                                          <  1:3093212/150‑1 (MQ=255)

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

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

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

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

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

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

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

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

                 gATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCTGTAGGA‑GCGGGtt                                                                                                                                     <  1:3050452/150‑1 (MQ=255)

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

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

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

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

                              ccTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAATCCCGCTCCTACAAGGTTCCTGGAGGT‑GCGGGTTTTCCCGCGAagag                                                                                                                        >  2:3363419/1‑150 (MQ=255)

                              ccTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGGGAACCCGCGCCAACAAGGGTGCGGTAGGA‑GCGGGGTTTCCCCCGCagag                                                                                                                        >  1:2253098/1‑150 (MQ=255)

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

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

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

                                          gAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAAACCGCGCCTACAAGGTTCCTGGAGGA‑GCGGGTTTACCCGCGCAGAGGCCAGCGCC‑Gaa                                                                                                           >  2:1746550/1‑150 (MQ=38)

                                             aCTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCTGTAGGA‑GCGGGTTTACCCGCGAAGAGGCCAGCACA‑Gaaaaa                                                                                                        <  2:2560407/150‑1 (MQ=255)

                                             aCTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTGT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCTGTAGGA‑GCGGGTTTACCCGCGAAGAGGCCAGCACA‑Gaaaaa                                                                                                        <  2:2635085/150‑1 (MQ=255)

                                             aCTGCTGCCTGTCGCACGAAACCCCCGGCCAGTGT‑GCGGGGT‑TTTTTTTTGGGGGGTATTTTT‑CGGGGGCCGGCCTCTGGGC‑GGGTAAACCCGCTCCTCCCCGGTTCCTGTAGGA‑GCGGGTTTACCCGCGAAGAGGCCAGCACA‑Gaaaaa                                                                                                        <  2:3810503/150‑1 (MQ=25)

                                                  tgcCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGGTAACCCCCTCCGACAAGGTTACTGTAGGA‑GCGGGTTTTCCCGCGGAGAGGCCAGCCCC‑Gacaaataaaa                                                                                                   >  2:245371/1‑150 (MQ=255)

                                                    ccTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTAGC‑GGGTTAACACCCGCCGACAAGGTTGCAGGTGGAGGGGGGTTTACCC‑CGAAGGGGGCCGCACA‑Cgcaaataaaaat                                                                                                 >  1:1793927/1‑150 (MQ=255)

                                                    ccTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTATTCGC‑GGGTAAACCCGCTCCTACCAGGTTCCGGTAGGA‑GCGGGTTTACCCGCGTAGTGGCCAGCAAA‑Gaaaaataaaaat                                                                                                 >  2:364766/1‑150 (MQ=38)

                                                      tGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTAGCGGGGTATACCC‑CGCCCACCACGTTGGGGGTGGG‑GGGGGGGTTCCCGCGGAGAGGCCCGCCCA‑Gcaaaaaaaaaaaaa                                                                                               >  1:3129213/1‑147 (MQ=255)

                                                        tttCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAAACCGCACCTACAAGGTTGCTGTATGA‑GCGGGTTTTCCCGCGAAGAGGCCAGCAAA‑GAAAAATAAAAATACCt                                                                                             >  2:4223411/1‑150 (MQ=38)

                                                        ccaaaaaaaaaCCCCGG‑CACAGTGGCCGGGGTTTTTTTTTTGGTGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCGCCTTCAAGGTTCCGGTAGGG‑GCGGGGTTACCCGCGAAGAGGCCAGCACA‑GGAAAATAAAAATAcc                                                                                              >  2:2233506/8‑150 (MQ=21)

                                                          tCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCACCTACAAAGTTCCGGTAGGA‑GCGGGTGTACCCCCCCACAGGCCCGCAAC‑GAAAAATAAAAATAccagg                                                                                           >  2:789220/1‑147 (MQ=37)

                                                            aCGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAAACCCCTCCGACCAGGTTGCGGGTGGG‑GCGGGGTTAACCGCGCAGCGGCCAGCACC‑GAAAAATAAAAAAACCTGCtt                                                                                         >  2:4159108/1‑150 (MQ=255)

                                                                 aaCCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCGTCGC‑GGGGAAACCCGCACCCACAAGGGTCCGGGAGGC‑GCGGGGGTACCCGCGCAGAGGCCCCCCCC‑GCAAAAAAAAAAAACCGGATTTTCaa                                                                                    >  2:1956709/1‑150 (MQ=255)

                                                                   ccccGGCCCCTGT‑GCGGGTG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACACGGTTCCTTTTGGA‑GCGGGTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCaaaa                                                                                  <  1:734189/150‑1 (MQ=37)

                                                                           aCTGT‑GCCGGGG‑ATTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAAACCCCTCCTACAAGGGTGCTGTAGGG‑GCGGGTTTTCCCGCGACGAGGCCAGCACA‑GAAAAAAAAAAATACCTGATTTTGAAAAAATAGAta                                                                          >  2:2791069/1‑149 (MQ=37)

                                                                            cTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCTGTAGGA‑GCGGGTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCt                                                                         <  2:3381908/150‑1 (MQ=255)

                                                                                       gttttttttttGGGGTTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGTTAAACCCGCTCCTACACGGTTCCTGTAGGA‑GCGGGTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTc                                                                <  2:429023/148‑1 (MQ=255)

                                                                                        ttttttttttGGGGTTATTTTT‑CCGGTGCCGGCCTCTTCGC‑GGTTAAACCCGCTCCTACAAGGTTCCTGTCGGA‑GCGGGTTAACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCt                                                               <  2:1877266/149‑1 (MQ=38)

                                                                                            ttttttGGGGGTTTTTTTCCCGGTGCCGCCCTCTTCGC‑GGGAACCCCCGCTCCTACACGGTACGTGTAGGA‑GCGGGTTCACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCtaat                                                            <  1:2702163/150‑1 (MQ=37)

                                                                                                 tGGGGGTATTTTC‑CCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCTGTAGGA‑GCGGTTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACa                                                      <  2:1028763/150‑1 (MQ=255)

                                                                                                     ggTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTTAACCCGCTCCTACAAGGTTGCTGGAGGA‑GCGGGTTTAACCGCGGAGGGGCCAGCACC‑GGAAAATAAAAATACCTGTTTTTAAAAAAATAGGTATTAAAAGTTGTAATAAAACAGGaa                                                  >  2:3510506/1‑150 (MQ=25)

                                                                                                     ggTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAAGCCGCTCCTACAATGTTGCTGTAGGT‑GCGGGTTTAGCCGCGTAGAGGCCAGCACA‑GAAAAATTAAAATACCTGATTTTGTAAAAATAGATCTGAAAAGTTGTAATAAAAGAGGaa                                                  >  2:1441133/1‑150 (MQ=25)

                                                                                                        aTTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCTGTAGGA‑GCGGTTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAAtt                                               <  2:3075426/150‑1 (MQ=255)

                                                                                                        aTTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCTGTAGGA‑GCGGGTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAAtt                                               <  2:3075433/150‑1 (MQ=255)

                                                                                                               cgtgtgCCGCCCTCTGGGG‑GGGTAACCCCGCTCCTCCACGGTTCCTGTTGGC‑GCGGTTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCtt                                         <  1:3216641/148‑1 (MQ=255)

                                                                                                                             ttCGC‑GGGTAAACCCGCGCCTACAAGGTTCCTGTAGGA‑GCGGGTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGTTCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGt                           >  2:1183987/1‑150 (MQ=255)

                                                                                                                               cgc‑gGGTAAACCCGCTCCTACAAGGTTCCTGTAGGA‑GCGGGTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGttt                         >  2:3702639/1‑150 (MQ=255)

                                                                                                                                gc‑gGGTAAACCCGCTCCTACACGGTTCCTGTAGGA‑GCGGGTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTc                        <  2:77860/150‑1 (MQ=255)

                                                                                                                                gc‑gGGTAAACCCGCTCCTACAAGGTTCCTGTAGGA‑GCGGGTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAGTACAGGAAATTATCCTTTACAGGCATCCTGTTTc                        >  2:3148245/1‑150 (MQ=255)

                                                                                                                                gc‑gGGAAAACCCGCTCCTACAAGGTTCCTGTTGGA‑GCGGGTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTc                        <  1:713829/150‑1 (MQ=255)

                                                                                                                                 ggggttcaaCCCGATCCTACACGGTTCCTGTTGGT‑GGGTGTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTc                        <  1:2790275/143‑1 (MQ=255)

                                                                                                                                     gtaatgccGCTCCTACAAGGTTCGTGTAGGA‑GCGGGTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAATTAGATCTGAATAGTTCTAATAATACAGGTAATTATCGTTTACAGGCATGCTGTTTCGGct                    >  2:2341787/7‑149 (MQ=38)

                                                                                                                                       aaaCCCGCTCCTACAAGGTTCCTGTAGGA‑GCGGGTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACt                  >  2:3899452/1‑150 (MQ=255)

                                                                                                                                        aaCCCGCTCCTACAAGGTTCCTGTAGGA‑GCGGGTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTa                 >  1:1241004/1‑150 (MQ=255)

                                                                                                                                          cccGCTCCTACAGGAACCTTGTAGGA‑GCGGGTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtc               >  1:1961007/1‑150 (MQ=38)

                                                                                                                                           ccGCTCCTACAGGAACCTTGTAGGA‑GCGGGTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtct              >  1:3496614/1‑150 (MQ=38)

                                                                                                                                           ccGCTCCTACAGGAACCTTGTAGGA‑GCGGGTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtct              >  2:1238897/1‑150 (MQ=38)

                                                                                                                                           ccGCTCCTACAAGGTTCCTGTAGGA‑GCGGGTTTACCCGCGAAGAGGCCAGCACAGGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtc               >  1:1363848/1‑150 (MQ=255)

                                                                                                                                           ccGCTCCTACAAGGTTCCTGTAGGA‑GCGGGTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtct              >  1:1156226/1‑150 (MQ=255)

                                                                                                                                           ccGCTCCTACAAGGTTCCTGTAGGA‑GCGGGTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtct              >  2:378955/1‑150 (MQ=255)

                                                                                                                                           ccGCTCCTACAAGGTTCCTGTAGGA‑GCGGGTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtct              >  2:2525492/1‑150 (MQ=255)

                                                                                                                                           ccGCTCCTACAAGGTTCCTGTAGGA‑GCGGGTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtct              >  1:4346146/1‑150 (MQ=255)

                                                                                                                                            cGCTCCTACAAGGTTCCTGTAGGA‑GCGGGTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtctt             >  2:2354765/1‑150 (MQ=255)

                                                                                                                                             gcacacacacGGTCCCGGTCGGC‑GGGGGTTTCCCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtcttc            <  1:3832093/144‑1 (MQ=255)

                                                                                                                                             gCTCCTACACGGTTCCTGTAGGC‑GCGGGTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtcttc            <  2:1331475/150‑1 (MQ=255)

                                                                                                                                             gCTCCTACACGGTACGTGTAGTC‑GGGGGTTTCCCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtcttc            <  1:3937322/150‑1 (MQ=255)

                                                                                                                                             gCTCCTACAAGGTTCCTGTAGGC‑GCGGGTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtcttc            <  2:2814860/150‑1 (MQ=255)

                                                                                                                                             gCTCCTACAAGGTTCCTGTAGGA‑GCGGGTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtcttc            >  1:1556563/1‑150 (MQ=255)

                                                                                                                                                 cTACAAGGTTCCTGTAGGA‑GCGGGTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGAt        >  2:643587/1‑150 (MQ=255)

                                                                                                                                                  tACAAGGTTCCTGTAGGA‑GCGGGTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGAtt       >  1:1209208/1‑150 (MQ=255)

                                                                                                                                                  tACAAGGTTCCTGTAGGA‑GCGGGTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGAtt       >  1:266595/1‑150 (MQ=255)

                                                                                                                                                  tACAAGGTTCCTGTAGGA‑GCGGGTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGAta       >  2:2952856/1‑149 (MQ=255)

                                                                                                                                                   aCAAGGTTCCTGTAGGA‑GCGGGTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTGCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTATGATTc      >  1:3136085/1‑150 (MQ=255)

                                                                                                                                                    cAAGGTTCCTGTAGGA‑GCGGGTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCa     >  2:2285089/1‑150 (MQ=255)

                                                                                                                                                    cAAGGTTCCTGTAGGA‑GCGGGTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCa     >  1:1119233/1‑150 (MQ=255)

                                                                                                                                                     aaGGTTCCTGTAGGA‑GCGGGTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCAg    >  2:3241215/1‑150 (MQ=255)

                                                                                                                                                       ggTTCCTGTAGGA‑GCGGGTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCAGgc  >  2:2462501/1‑150 (MQ=255)

                                                                                                                                                       ggTTCCTGTAGGA‑GCGGGTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCAGgc  >  2:3757116/1‑150 (MQ=255)

                                                                                                                                                       ggTTCCTGTAGGA‑GCGGGTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATCATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCAGgc  >  2:2858249/1‑150 (MQ=255)

ACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCTGTAGGA‑GCGGGTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCAGGC  >  NC_002947/955006‑955302

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