BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,378,143	T→C	8.0%	P464P (CCA→CCG)	dnaK ←	chaperone protein DnaK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,378,143	0	T	C	8.0%	647.7 / 30.1	213	P464P (CCA→CCG)	dnaK	chaperone protein DnaK
Reads supporting (aligned to +/- strand): ref base T (93/103); new base C (5/12); total (98/115)
Fisher's exact test for biased strand distribution p-value = 2.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.89e-01

Only 100 of 217 total aligned reads displayed.

CGCTCGATCTCTTCGTCCGACAGGCCGGAGTTGGCCTTGATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTT  >  NC_002947/5377994‑5378292

 gCTCGATCTCTTCGTCCGACAGGCCGGAGTTGGCCTTGATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCAAGCCGCGCGGTGCTg                                                                                                                                                      >  2:4321262/1‑150 (MQ=255)

  cTCGATCTCTTCGTCCGACAGGCCGGAGTTGGCCTTGATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTgg                                                                                                                                                     <  2:5091592/150‑1 (MQ=255)

  cTCGATCTCTTCGTCCGACAGGCCGGAGTTGGCCTTGATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTgg                                                                                                                                                     >  2:578811/1‑150 (MQ=255)

   tCGATCTCTTCGTCCGACAGGCCGGAGTTGGCCTTGATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGc                                                                                                                                                    >  1:1037691/1‑150 (MQ=255)

   tCGATCTCTTCGTCCGACAGGCCGGAGTTGGCCTTGATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGc                                                                                                                                                    <  1:4363801/150‑1 (MQ=255)

    cGATCTCTTCGTCCGACAGGCCGGAGTTGGCCTTGATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCg                                                                                                                                                   <  2:1081509/150‑1 (MQ=255)

           ttCGTCCGACAGGCCGGAGTTGGCCTTGATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTc                                                                                                                                            >  1:4493529/1‑150 (MQ=255)

           ttCGTCAGACAGGCCGGAGTTGGCCTTGATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTc                                                                                                                                            <  1:3981082/150‑1 (MQ=255)

            tCGTCCGACAGGCCGGAGTTGGCCTTGATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCa                                                                                                                                           <  1:1571527/150‑1 (MQ=255)

            tCGTCCGACAGGCCGGAGTTGGCCTTGATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCa                                                                                                                                           <  2:4318133/150‑1 (MQ=255)

               tCCGACAGGCCGGAGTTGGCCTTGATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGcc                                                                                                                                        >  1:4134574/1‑150 (MQ=255)

                   aCAGGCCGGAGTTGGCCTTGATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGt                                                                                                                                    >  2:4612775/1‑150 (MQ=255)

                       gCCGGAGTTGGCCTTGATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGaa                                                                                                                                <  1:1734416/150‑1 (MQ=255)

                       gCCGGAGTTGGCCTTGATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGaa                                                                                                                                >  2:476317/1‑150 (MQ=255)

                       gCCGGAGTTGGCCTTGATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGaa                                                                                                                                >  1:3175843/1‑150 (MQ=255)

                       gCCGGAGTTGGCCTTGATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGaa                                                                                                                                >  1:324049/1‑150 (MQ=255)

                        ccGGAGTTGGCCTTGATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAAc                                                                                                                               >  1:591486/1‑150 (MQ=255)

                         cGGAGTTGGCCTTGATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACt                                                                                                                              >  1:2556631/1‑150 (MQ=255)

                          ggAGTTGGCCTTGATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACtt                                                                                                                             <  1:1462890/150‑1 (MQ=255)

                           gAGTTGGCCTTGATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTg                                                                                                                            <  2:593756/150‑1 (MQ=255)

                            aGTTGGCCTTGATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGATGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGc                                                                                                                           <  1:268872/150‑1 (MQ=255)

                                 gCCTTGATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGc                                                                                                                      <  1:3209244/150‑1 (MQ=255)

                                 gCCTTGATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGc                                                                                                                      >  2:2967868/1‑150 (MQ=255)

                                   cTTGATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGa                                                                                                                    >  1:2844680/1‑150 (MQ=255)

                                    ttGATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGAc                                                                                                                   >  2:2811455/1‑150 (MQ=255)

                                      gATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACtt                                                                                                                 <  2:2087899/150‑1 (MQ=255)

                                         cACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGtt                                                                                                              <  2:2880804/150‑1 (MQ=255)

                                         cACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGtt                                                                                                              >  2:2643237/1‑150 (MQ=255)

                                         cACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGtt                                                                                                              <  1:4281932/150‑1 (MQ=255)

                                         cACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGtt                                                                                                              >  2:729601/1‑150 (MQ=255)

                                           cGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCt                                                                                                            <  2:1256471/150‑1 (MQ=255)

                                             aTCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTgc                                                                                                          <  2:1660977/150‑1 (MQ=255)

                                                  cTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGc                                                                                                     <  2:4571235/150‑1 (MQ=255)

                                                  cTGAGTCTTGCCGGTGGCCTTGCCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGc                                                                                                     >  2:3821940/1‑150 (MQ=255)

                                                   tGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGATGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGcc                                                                                                    <  1:505282/150‑1 (MQ=255)

                                                   tGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGcc                                                                                                    >  1:1260680/1‑150 (MQ=255)

                                                   tGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGcc                                                                                                    >  1:4645724/1‑150 (MQ=255)

                                                        cTTGCCTGTGCCCTTGTCCTTTGCGGACACGTTCACGATACCGTTGGCGTCGATGTCGAAGGTCACTTCGATCTGTGGAACGCCACGCGGTGCCGGCGGAATGCCAACGAGGTCGAACTGGCCGAGGATCTTGTTGTCGGCAGccatttc                                                                                               <  2:3224403/150‑6 (MQ=255)

                                                        cTTGCCTGTGCCCTTGTCCTTTGCGGACACGTTCACGATACCGTTGGCGTCGATGTCGAAGGTCACTTCGATCTGTGGAACGCCACGCGGTGCCGGCGGAATGCCAACGAGGTCGAACTGGCCGAGGATCTTGTTGTCGGCAGccatttc                                                                                               >  1:4242162/1‑145 (MQ=255)

                                                         ttGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTa                                                                                              <  1:983126/150‑1 (MQ=255)

                                                          tGCCTGTGCCCTTGTCCTTTGCGGACACGTTCACGATACCGTTGGCGTCGATGTCGAAGGTCACTTCGATCTGTGGAACGCCACGCGGTGCCGGCGGAATGCCAACGAGGTCGAACTGGCCGAGGATCTTGTTGTCGGCAGccatttcac                                                                                             <  1:2042803/150‑8 (MQ=255)

                                                             cGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCt                                                                                          >  2:1471894/1‑150 (MQ=255)

                                                                    cTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTg                                                                                   >  1:799997/1‑150 (MQ=255)

                                                                    cTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTg                                                                                   <  2:4063415/150‑1 (MQ=255)

                                                                    cTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTg                                                                                   >  2:1772708/1‑150 (MQ=255)

                                                                     ttGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTgc                                                                                  >  1:4528725/1‑150 (MQ=255)

                                                                     ttGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTgc                                                                                  >  2:3502261/1‑150 (MQ=255)

                                                                     ttGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTgc                                                                                  >  1:410824/1‑150 (MQ=255)

                                                                     ttGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGGTCTGCGCAGCCTGCTTACGCTCACCTTgc                                                                                  >  2:2791644/1‑150 (MQ=255)

                                                                      tGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTgca                                                                                 >  1:1226198/1‑150 (MQ=255)

                                                                          tttCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCAc                                                                             >  1:4581226/1‑150 (MQ=255)

                                                                            tCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGt                                                                           <  2:112320/150‑1 (MQ=255)

                                                                            tCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGt                                                                           <  1:134568/150‑1 (MQ=255)

                                                                              gCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGa                                                                         >  2:530346/1‑150 (MQ=255)

                                                                              gCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGa                                                                         <  2:883537/150‑1 (MQ=255)

                                                                                 ccTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATg                                                                      <  2:1333852/150‑1 (MQ=255)

                                                                                    aCGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTAGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTc                                                                   >  2:5019528/1‑150 (MQ=255)

                                                                                        tcacgATACCGTTGGCGTCGATGTCGAAGGTCACTTCGATCTGTGGAACGCCACGCGGTGCCGGCGGAATGCCAACGAGGTCGAACTGGCCGAGGATCTTGTTGTCGGCAGCCATTTCACGCTCGCCCTGGAACACGCGGATCGTCACgg                                                               <  2:248699/146‑1 (MQ=255)

                                                                                          aGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGgcg                                                             <  1:3568253/150‑1 (MQ=255)

                                                                                              tGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGctcg                                                         >  2:3369739/1‑149 (MQ=255)

                                                                                               gCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTg                                                        <  1:5000267/150‑1 (MQ=255)

                                                                                                  gTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAACGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGtt                                                     <  2:1246397/150‑1 (MQ=255)

                                                                                                     ggCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTgtc                                                  >  1:2403828/1‑150 (MQ=255)

                                                                                                        gTCGATGTCGAAGGTCACTTCGATCTGTGGAACGCCACGCGGTGCCGGCGGAATGCCAACGAGGTCGAACTGGCCGAGGATCTTGTTGTCGGCAGCCATTTCACGCTCGCCCTGGAACACGCGGATCGTCACGGCCGACTGATTgtcttc                                               <  2:4370864/150‑1 (MQ=255)

                                                                                                         tCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGtcgtcg                                              <  2:5097330/150‑1 (MQ=255)

                                                                                                         tCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGtcgtcg                                              >  2:4666834/1‑150 (MQ=255)

                                                                                                         tCGATGTCGAAGGTCACTTCGATCTGTGGAACGCCACGCGGTGCCGGCGGAATGCCAACGAGGTCGAACTGGCCGAGGATCTTGTTGTCGGCAGCCATTTCACGCTCGCCCTGGAACACGCGGATCGTCACGGCCGACTGATTgtcttca                                              <  2:4439874/150‑2 (MQ=255)

                                                                                                           gATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGc                                            >  1:469152/1‑150 (MQ=255)

                                                                                                                 gAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCga                                      >  1:4344091/1‑150 (MQ=255)

                                                                                                                 gAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCga                                      >  2:237251/1‑150 (MQ=255)

                                                                                                                 gAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCga                                      <  1:1772708/150‑1 (MQ=255)

                                                                                                                 gAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCga                                      <  1:1247385/150‑1 (MQ=255)

                                                                                                                   aGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCgaga                                    <  1:3607536/150‑1 (MQ=255)

                                                                                                                       cACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAacac                                <  2:4134574/150‑1 (MQ=255)

                                                                                                                        aCTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACAcc                               >  2:1713056/1‑150 (MQ=255)

                                                                                                                         cTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCt                              >  2:4296584/1‑150 (MQ=255)

                                                                                                                          ttCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTg                             >  2:607731/1‑150 (MQ=255)

                                                                                                                          ttCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTg                             >  1:1822627/1‑150 (MQ=255)

                                                                                                                          ttCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTg                             >  2:2929225/1‑150 (MQ=255)

                                                                                                                          ttCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTg                             >  1:48050/1‑150 (MQ=255)

                                                                                                                           tCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGc                            <  2:1404324/150‑1 (MQ=255)

                                                                                                                           tCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGc                            <  1:2107363/150‑1 (MQ=255)

                                                                                                                           tCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGc                            <  2:4805448/150‑1 (MQ=255)

                                                                                                                              ccgcgcggcgTGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGAc                         >  2:232066/5‑150 (MQ=255)

                                                                                                                               tttGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGAct                        >  1:347614/1‑150 (MQ=255)

                                                                                                                                 tGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGActtc                      <  2:4850466/150‑1 (MQ=255)

                                                                                                                                 tGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGActtc                      <  2:866904/150‑1 (MQ=255)

                                                                                                                                 tGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGAGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGActtc                      <  2:2280047/150‑1 (MQ=255)

                                                                                                                                  gCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGActtct                     <  2:4503400/150‑1 (MQ=255)

                                                                                                                                     gCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTgg                  >  1:3025280/1‑150 (MQ=255)

                                                                                                                                       aCGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTc                <  1:4227781/150‑1 (MQ=255)

                                                                                                                                         gccgcGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCgg              <  1:2318081/150‑1 (MQ=255)

                                                                                                                                             cgcgGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGAtg          <  1:4466815/150‑1 (MQ=255)

                                                                                                                                             cgcgGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGAtg          <  2:1856094/150‑1 (MQ=255)

                                                                                                                                              gcgGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGAtgg         <  1:4612775/150‑1 (MQ=255)

                                                                                                                                              gcgGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGAtgg         >  1:4424048/1‑150 (MQ=255)

                                                                                                                                               cgGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGAtggt        <  2:1368673/150‑1 (MQ=255)

                                                                                                                                                ggTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGAtggtg       <  2:5059510/150‑1 (MQ=255)

                                                                                                                                                ggTGCCGGCGGAATGCCAACGAGGTCGAACTGGCCGAGGATCTTGTTGTCGGCAGCCATTTCACGCTCGCCCTGGAACACGCGGATCGTCACGGCCGACTGATTGTCTTCAGCGGTCGAGAAGGTCTGCGACTTCTTGGTCGGGAtcgtc       <  1:5059538/150‑2 (MQ=255)

                                                                                                                                                     tGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGtt  >  1:4380889/1‑150 (MQ=255)

CGCTCGATCTCTTCGTCCGACAGGCCGGAGTTGGCCTTGATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTT  >  NC_002947/5377994‑5378292

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