BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,102,219	G→C	7.2%	intergenic (+213/+79)	PP_mr66 → / ← acsA‑I	P6/acetyl‑CoA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,102,219	0	G	C	7.2%	527.8 / 25.4	166	intergenic (+213/+79)	PP_mr66/acsA‑I	P6/acetyl‑CoA synthetase
Reads supporting (aligned to +/- strand): ref base G (79/74); new base C (2/10); total (81/85)
Fisher's exact test for biased strand distribution p-value = 3.22e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.72e-01

Only 100 of 166 total aligned reads displayed.

TGCCTTGATGCACAAAATGCTGACTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCA  >  NC_002947/5102071‑5102365

tGCCTTGATGCACAAAATGCTGACTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAGTGCAAAAGCGGGCTTGCCCGCGAATGCGTGa                                                                                                                                                     <  1:4726008/150‑1 (MQ=255)

tGCCTTGATGCACAAAATGCTGACTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGa                                                                                                                                                     <  2:4738157/150‑1 (MQ=255)

tGCCTTGATGCACAAAATGCTGACTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGa                                                                                                                                                     <  2:6449569/150‑1 (MQ=255)

tGCCTTGATGCACAAAATGCTGACTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGa                                                                                                                                                     <  1:6279044/150‑1 (MQ=255)

tGCCTTGATGCACAAAATGCTGACTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGa                                                                                                                                                     <  2:3662882/150‑1 (MQ=255)

tGCCTTGATGCACAAAATGCTGACTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCAGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGa                                                                                                                                                     <  2:4341693/150‑1 (MQ=255)

tGCCTTGATGCACAAAATGCTGACTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCAGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGa                                                                                                                                                     <  2:3127969/150‑1 (MQ=255)

 gCCTTGATGCACAAAATGCTGACTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAg                                                                                                                                                    <  1:1135269/150‑1 (MQ=255)

 gCCTTGATGCACAAAATGCTGACTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAg                                                                                                                                                    <  2:2534571/150‑1 (MQ=255)

   cTTGATGCACAAAATGCTGACTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCCCGACAGCGAGAggg                                                                                                                                                  >  1:5482529/1‑149 (MQ=255)

    ttGATGCACAAAATGCTGACTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTc                                                                                                                                                 <  2:5065000/150‑1 (MQ=255)

     tGATGCACAAAATGCTGACTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCa                                                                                                                                                <  1:6139263/150‑1 (MQ=255)

     tGATGCACAAAATGCTGACTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCa                                                                                                                                                <  1:6868160/150‑1 (MQ=255)

       aTGCACAAAATGCTGACTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAgg                                                                                                                                              <  2:1223055/150‑1 (MQ=255)

       aTGCACAAAATGCTGACTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAgg                                                                                                                                              <  1:3932825/150‑1 (MQ=255)

       aTGCACAAAATGCTGACTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAgg                                                                                                                                              >  2:1550896/1‑150 (MQ=255)

       aTGCACAAAATGCTGACTCCGATGCACCGAAACGGGCCGGAATTGTTCCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAgg                                                                                                                                              <  2:4093835/150‑1 (MQ=255)

        tGCACAAAATGCTGACTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTTACAGAACAGTAACCGATGATTGTCTTGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGATGTCAggt                                                                                                                                             <  2:2377209/150‑2 (MQ=255)

        tGCACAAAATGCTGACTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAggg                                                                                                                                             <  2:4805538/150‑1 (MQ=255)

          cacaAAATGCTGACTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAggggc                                                                                                                                           >  2:3117134/1‑149 (MQ=255)

          cacaAAATGCTGACTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGGGAGGTCAGGGGa                                                                                                                                           >  2:5278987/1‑150 (MQ=255)

             aaaaTGCTGACTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAg                                                                                                                                        <  1:1734015/150‑1 (MQ=255)

             aaaaTGCTGACTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAg                                                                                                                                        <  1:1734860/150‑1 (MQ=255)

                 tGCTGACTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑Ta                                                                                                                                   <  2:4947621/150‑1 (MQ=255)

                    tGACTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGcc                                                                                                                                <  1:6114791/150‑1 (MQ=255)

                      aCTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTg                                                                                                                              <  1:5936346/150‑1 (MQ=255)

                        tCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACCGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTAGCCCGCGCATGCGGGAGGTCAGGGGCCAGCG‑TAGCCTGaa                                                                                                                            >  2:6430789/1‑150 (MQ=255)

                        tCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATTCGTGAGGTGAGGGGACAGCT‑TAGCCTGaa                                                                                                                            >  2:3163627/1‑150 (MQ=255)

                        tCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGaa                                                                                                                            <  2:5148412/150‑1 (MQ=255)

                          cGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGGCCGCGCAAGCGTGAGGTCAGGGGACAGCT‑TGTCCGGAACt                                                                                                                          >  1:2328048/1‑150 (MQ=255)

                             tGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGAc                                                                                                                       <  1:1470480/150‑1 (MQ=255)

                             tGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGAc                                                                                                                       <  2:244196/150‑1 (MQ=255)

                             tGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGAc                                                                                                                       <  1:5790312/150‑1 (MQ=255)

                             tGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGAc                                                                                                                       <  2:4340414/150‑1 (MQ=255)

                                aCCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCa                                                                                                                    <  2:3896620/150‑1 (MQ=255)

                                aCCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCGTGCCCGCGCATGCGTGAGGTCAGGGGCCAGCT‑TAGCCTGAACTGACGCa                                                                                                                    >  2:6778825/1‑150 (MQ=255)

                                      aCGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTCAGTTCAGGCTA‑AGCTGTCCCCTGACCTCACGCATTcgcg                                                                                                              <  1:3005047/150‑1 (MQ=21)

                                      aCGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTCAGTTCAGGCTA‑AGCTGTCCCCTGACCTCACGCATTcgcg                                                                                                              <  1:6001794/150‑1 (MQ=21)

                                       cGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCgg                                                                                                             <  1:6425583/150‑1 (MQ=255)

                                       cGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCgg                                                                                                             <  2:6639158/150‑1 (MQ=255)

                                       cGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCgg                                                                                                             >  2:3578442/1‑150 (MQ=255)

                                       cGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCgg                                                                                                             <  2:4732994/150‑1 (MQ=255)

                                       cGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTCAGTTCAGGCTA‑AGCTGTCCCCTGACCTCACGCATTCGCgg                                                                                                             <  1:4371051/150‑1 (MQ=21)

                                       cGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTCAGTTCAGGCTA‑AGCTGTCCCCTGACCTCACGCATTCGCgg                                                                                                             <  2:2122742/150‑1 (MQ=21)

                                       cGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTCAGTTCAGGCTA‑AGCTGTCCCCTGACCTCACGCATTCGCgg                                                                                                             <  1:2184440/150‑1 (MQ=21)

                                        gggCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCggg                                                                                                            >  2:5997961/1‑150 (MQ=255)

                                           ccGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCaa                                                                                                         <  2:1996520/150‑1 (MQ=255)

                                           ccGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCaa                                                                                                         <  2:2287833/150‑1 (MQ=255)

                                           ccGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTCAGTTCAGGCTA‑AGCTGTCCCCTGACCTCACGCATTCGCGGGCaa                                                                                                         <  1:342006/150‑1 (MQ=21)

                                            cGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCG‑TAGCCTGAACTGAAGCATTAGCGGGCAA‑g                                                                                                       >  2:1444565/1‑150 (MQ=255)

                                              gAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAA‑Gcc                                                                                                     <  2:1870415/150‑1 (MQ=255)

                                                    tttCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAAATGAAGCATTCGCGGGCAA‑GCCCGCGcc                                                                                               >  2:699231/1‑150 (MQ=255)

                                                    tttCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCCAATTCGGGAGGTGAGGGGACAGCT‑TATCCGGAACGTACGCAATCGCGGGCAA‑GCCCGCGcc                                                                                               >  1:330212/1‑150 (MQ=255)

                                                       tcaccggcCCCGCACAGGTTCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCGTGCCCCCGCATGCGTGTGGTCAGGGGGCAGCT‑TTGCCTGAACAGAAGCATTCGCGGGCAA‑GCCCCCTCCGGc                                                                                            >  1:3707524/7‑150 (MQ=255)

                                                       ccTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGc                                                                                            >  2:5610052/1‑150 (MQ=255)

                                                         tGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCCCGAATGCGTGAGGTCAGGGGACAGCG‑TAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAg                                                                                          >  2:738482/1‑150 (MQ=255)

                                                           gcgcCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCGTGCCCCCGAAAGCGTGACGTGAGGGGACAGCT‑TAGCCTGAACTGACGCAATAGCGGGGAA‑GCCCGCGCCAGCAggg                                                                                        >  1:4142424/1‑150 (MQ=255)

                                                                ccGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTTCCCGCGCAAGGGGGAGGTCCGGGGCCAGCG‑TAGGCTGAACAGAACCATTAGCGGGCAA‑GCCCGCGCCAGCCggggggcg                                                                                   >  1:621663/1‑147 (MQ=255)

                                                                     ccTGTGCGGGGCGGGTGACAGAACATTAACCGATGATTGTCGGTTAAAAATGCAAACGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCg                                                                              <  2:6229551/150‑1 (MQ=39)

                                                                         tgCGGGGCCGGTGACAGAACAGTACCCGATGTTTGGCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGGGTGAGGTCAGGGGACAGCT‑TAGCCTGAACGGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGaa                                                                          <  2:3952833/150‑1 (MQ=255)

                                                                             gggCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCGTTGCCGCGAATGCGTGAGGGGCGGGGGCAGCG‑TAGGCGGGACTGACGCAGTCGCGGGGAA‑GCCCCCGCCGGCCGGGGGGGTTGCCGGGGAgggg                                                                      >  1:3626055/1‑149 (MQ=255)

                                                                             ggcccGGTTAAAGAACCTTAACCTATTATTGTCGGTAAACAATGAAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTATGGGGACAGAT‑TACCCTGAACTGACGCATTCGCGGGAAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGt                                                                      <  2:1612100/147‑1 (MQ=255)

                                                                                cgggTGACATAACATTAACCGTTGTTTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGTGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCTTTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTa                                                                   <  2:2740430/148‑1 (MQ=255)

                                                                                ccGGTTACACAACAGTACCCGATGATTGTCGGTTAAAAATGCAAACGGGGGCTGGCCCGCTAATGCTTGAGGTCAGGGGACAGCT‑TCGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTa                                                                   <  2:2740326/150‑1 (MQ=255)

                                                                                        aGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTTCCCGCCAATTCGTGAGGTCAGGGGGCACCT‑TATCCGGAACGGACGCAATAGCGGGCCA‑CCCCCCGCCTGCAGGGGGTATTGCCGGGGAGGGTTTACGCCGAAg                                                           >  1:111352/1‑150 (MQ=255)

                                                                                        aGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGCAAGCGTGAGGTCAGGGGACAGCT‑TAGCCGGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGGATTGCCGGGACGGGTTTACGCCGAAg                                                           >  1:5057228/1‑150 (MQ=255)

                                                                                        aGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAg                                                           >  2:1560024/1‑150 (MQ=255)

                                                                                        aGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCG‑TAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTCTTGCCGGGCAGGGTTTAAGCCGTAg                                                           >  2:6500934/1‑150 (MQ=39)

                                                                                            ctgTACCCTATGTTTGTCGGGTAACAATGCAAACGCGGGCGTGCCCGCGATTGCGTGCGTTCAGGTTACAGCT‑TCGCCTGACCTGCCGCTTTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAg                                                       <  2:1057310/148‑1 (MQ=255)

                                                                                               tAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGtt                                                    >  2:3837473/1‑150 (MQ=255)

                                                                                               tAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGtt                                                    >  2:5042795/1‑150 (MQ=255)

                                                                                               tAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCAATCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGtt                                                    >  2:518363/1‑150 (MQ=255)

                                                                                                 aCCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCa                                                  >  2:6090232/1‑150 (MQ=255)

                                                                                                  ccGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAAATCGTGAGGTAAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCCA‑GCCCCCTCCTGCAGGGGGTATTGCCGGGAAGGGGTTAAGCCGCAGCCAGGTTCaa                                                 >  2:1834126/1‑149 (MQ=255)

                                                                                                   cGATGATTGTCGGGTAACAATGCAAAAGAGGGCTTGCCCGCGAATGCGTGTGGTCAGGGGACAGCA‑TAGCCAGAACTGACGCATTCGCGGGCAA‑GCCCGCACCAGCAGGGGGTCTGGCCGGGAAGGGTTTACGACGAAACCAGGGTCATc                                                >  1:5679547/1‑150 (MQ=255)

                                                                                                    gatgatTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCAtag                                               >  2:933672/1‑148 (MQ=255)

                                                                                                    gatgatTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCG‑TAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCCGGTTCATCg                                               >  2:1484756/1‑150 (MQ=255)

                                                                                                    gatgatTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGGGAGGTCAGGGGACAGCT‑TAGCCTGAACTGAAGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTAAGCCGAAGCCAGGTTAAtag                                               >  2:229483/1‑148 (MQ=255)

                                                                                                        atTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCCGGGGGTATTGCCGGGGAGGGTTTACGCCGAAGCCCGGTTCATCGCCtt                                           >  2:4092368/1‑150 (MQ=255)

                                                                                                           gtgggttAACAATGCAAACGGGGGCTTGCCCGCGAATGCGTGCGGTCTGGGGACAGCT‑TAGCCTCAACGGACGAATCCGGGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTtgtg                                        <  2:4347900/144‑1 (MQ=255)

                                                                                                           gTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGTCAGCT‑TAGCCTGAACGGTCGCGTTCGCGGGCAA‑GCCCGCTGCGGCAGGGGGTGTTGCCGGGCAGGGTTTACGCCGAAGCCAGGGTAATCGCCTtgtg                                        >  2:5357049/1‑150 (MQ=255)

                                                                                                            tCGGTTAACATTGCAAACGCGGGCTTGCCCGCTAATGCGTTAGGTCTGGGTACAGAT‑TAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTgg                                       <  2:2119238/150‑1 (MQ=255)

                                                                                                               ggTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAAATGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTa                                    >  2:2202833/1‑150 (MQ=255)

                                                                                                                          attcgCGGGCAAGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTg                         >  2:6105863/5‑150 (MQ=255)

                                                                                                                          aaaaGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTg                         >  2:3707831/1‑150 (MQ=255)

                                                                                                                                    ttGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACcaca               >  1:966987/1‑150 (MQ=255)

                                                                                                                                       cccGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACAccc            >  1:5444030/1‑150 (MQ=255)

                                                                                                                                       cccGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACAccc            >  2:5131755/1‑150 (MQ=255)

                                                                                                                                       cccGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACAccc            >  1:4213730/1‑150 (MQ=255)

                                                                                                                                       cccGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACAccc            >  2:4581113/1‑150 (MQ=255)

                                                                                                                                       cccGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACAccc            >  1:2707981/1‑150 (MQ=255)

                                                                                                                                       cccGCGAATGCGTCAGTTCAGGCTA‑AGCTGTCCCCTGACCTCACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACcaacaag            >  1:4882449/1‑145 (MQ=18)

                                                                                                                                        ccGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCg           >  1:5949125/1‑150 (MQ=255)

                                                                                                                                        ccGCGAAGCCGTGAGGTCAGGGGACCGCT‑CACCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCg           <  2:5318345/150‑1 (MQ=255)

                                                                                                                                         cgcgAATGCGTGAGGTCTGGGGACAGCT‑TCGCCTGAACTGCCGCTTTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCgg          <  2:6588364/150‑1 (MQ=255)

                                                                                                                                          gcgAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCggg         >  1:1639150/1‑150 (MQ=255)

                                                                                                                                           gttttGCTTGTGGTCGGGGGCCTGCT‑TCGCCCCAAGTGCCGC‑TTCCGGGGCAACCCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGt        <  1:1061555/146‑1 (MQ=255)

                                                                                                                                            gAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTc       >  1:2014613/1‑150 (MQ=255)

                                                                                                                                                gCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAgg                                                                        <  2:550024/81‑1 (MQ=255)

                                                                                                                                                 cGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCa  >  2:2893058/1‑150 (MQ=255)

TGCCTTGATGCACAAAATGCTGACTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCA  >  NC_002947/5102071‑5102365

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