BRESEQ :: Evidence

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

Only 100 of 102 total aligned reads displayed.

AGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCTGTAGGA‑GCGGGTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCAGGCCGC  >  NC_002947/955008‑955305

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

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

aGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTc                                                                                                                                                        <  2:2749605/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:2903532/150‑1 (MQ=38)

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

          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)

          cTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCTGTAGGA‑Gc                                                                                                                                             <  1:1158589/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                 >  1:4346146/1‑150 (MQ=255)

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

                                                                                                                                         ccGCTCCTACAAGGTTCCTGTAGGA‑GCGGGTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtct                 >  2:2525492/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        >  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)

                                                                                                                                                        tCCTTTGGGA‑GGGGTTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGCAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCAGgccgc  <  1:3998760/150‑1 (MQ=255)

                                                                                                                                                        tCCTGTAGGA‑GCGGGTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCAGgccgc  <  2:714139/150‑1 (MQ=255)

AGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCTGTAGGA‑GCGGGTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCAGGCCGC  >  NC_002947/955008‑955305

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