BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,407,309	Δ1 bp	5.3%	coding (50/330 nt)	PP_2103 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,407,309	0	T	.	5.3%	1715.7 / 81.2	486	coding (50/330 nt)	PP_2103	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (234/225); new base . (14/12); total (248/238)
Fisher's exact test for biased strand distribution p-value = 8.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 494 total aligned reads displayed.

CGATGCCCCACAGTGGCGCATGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATTCGTCGAACACCGCCGCCTTCAGGCGCACGCCGCAGTACA  >  NC_002947/2407160‑2407458

 gATGCCCCACAGTGGCGCATGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCtt                                                                                                                                                      >  2:794378/1‑150 (MQ=255)

 gATGCCCCACAGTGGCGCATGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCtt                                                                                                                                                      <  1:4113288/150‑1 (MQ=255)

 gATGCCCCACAGTGGCGCATGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCtt                                                                                                                                                      <  2:3792978/150‑1 (MQ=255)

   tGCCCCACAGTGGCGCATGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTgc                                                                                                                                                    <  2:3884229/150‑1 (MQ=255)

     ccccACAGTGGCGCATGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTgcgc                                                                                                                                                  <  1:2269595/150‑1 (MQ=255)

     ccccACAGTGGCGCATGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTgcgc                                                                                                                                                  <  1:5543035/150‑1 (MQ=255)

      cccACAGTTGCGCATGCTTGATACCCAGGCCTTCATGTAGCGAGCGGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCTGCTGCGCC‑TGCGCat                                                                                                                                                <  1:2672899/150‑2 (MQ=255)

      cccACAGTGGCGCATGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCa                                                                                                                                                 >  2:6731603/1‑150 (MQ=255)

      cccACAGTGGCGCATGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCa                                                                                                                                                 <  1:3301104/150‑1 (MQ=255)

       ccACAGTGGCGCATGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAg                                                                                                                                                >  2:3027233/1‑150 (MQ=255)

       ccACAGTGGCGCATGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAg                                                                                                                                                >  2:1803874/1‑150 (MQ=255)

        cacaGTGGCGCATGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGt                                                                                                                                               <  2:1158771/150‑1 (MQ=255)

             tGGCGCATGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTcgcg                                                                                                                                          <  1:4438662/150‑1 (MQ=255)

               gcgcATGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTcgcgcg                                                                                                                                        >  1:1976295/1‑150 (MQ=255)

                cgcATGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCgg                                                                                                                                       <  1:2732639/150‑1 (MQ=255)

                    tGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTa                                                                                                                                   <  1:2168552/150‑1 (MQ=255)

                       ttGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTAcgc                                                                                                                                >  2:6203938/1‑150 (MQ=255)

                       ttGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTAcgc                                                                                                                                >  2:369957/1‑150 (MQ=255)

                       ttGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTAcgc                                                                                                                                >  1:402991/1‑150 (MQ=255)

                        tGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCGGTTCGCGCGGGTTAcgcg                                                                                                                               >  2:7289685/1‑150 (MQ=255)

                             cccAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGAGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGttt                                                                                                                          >  1:3874944/1‑150 (MQ=255)

                              ccAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTg                                                                                                                         >  1:3147459/1‑150 (MQ=255)

                              ccAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTg                                                                                                                         >  1:3765105/1‑150 (MQ=255)

                              ccAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTg                                                                                                                         >  2:3225279/1‑150 (MQ=255)

                               cAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTgg                                                                                                                        >  1:5705626/1‑150 (MQ=255)

                               cAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTgg                                                                                                                        >  2:3383887/1‑150 (MQ=255)

                                 ggCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCa                                                                                                                      <  2:6625955/150‑1 (MQ=255)

                                     ttCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTc                                                                                                                  >  1:6585817/1‑150 (MQ=255)

                                     ttCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTc                                                                                                                  >  2:147209/1‑150 (MQ=255)

                                        aTGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTgg                                                                                                               >  1:3507863/1‑150 (MQ=255)

                                            aGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAgg                                                                                                           <  2:6782377/150‑1 (MQ=255)

                                               gAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTc                                                                                                        <  1:3668160/150‑1 (MQ=255)

                                                 gCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCAt                                                                                                      <  1:2540873/150‑1 (MQ=255)

                                                  cTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGGGTTTGGCAGTTCTGGCAGGCTCATg                                                                                                     >  2:1244234/1‑150 (MQ=255)

                                                     cccATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCa                                                                                                  >  1:3039906/1‑150 (MQ=255)

                                                       cATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAAc                                                                                                >  1:3427195/1‑150 (MQ=255)

                                                          aCCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAgg                                                                                             <  1:6995423/150‑1 (MQ=255)

                                                            cgcgCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTg                                                                                                                         >  2:2620331/1‑120 (MQ=255)

                                                            cgcgCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTg                                                                                                                         <  1:2620331/120‑1 (MQ=255)

                                                            cgcgCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGgcg                                                                                           >  2:551037/1‑150 (MQ=255)

                                                            cgcgCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGgcg                                                                                           >  1:6615408/1‑150 (MQ=255)

                                                            cgcgCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGgcg                                                                                           >  1:4486157/1‑150 (MQ=255)

                                                             gcgcACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGgcgc                                                                                          <  1:6185879/150‑1 (MQ=255)

                                                              cgcACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCt                                                                                         >  1:4222605/1‑150 (MQ=255)

                                                                cACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTcc                                                                                       <  2:2561665/150‑1 (MQ=255)

                                                                        gcAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCtt                                                                               >  1:5888164/1‑150 (MQ=255)

                                                                        gcAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCtt                                                                               <  1:5071405/150‑1 (MQ=255)

                                                                        gcAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCtt                                                                               >  2:2889296/1‑150 (MQ=255)

                                                                        gcAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCtt                                                                               >  2:5149000/1‑150 (MQ=255)

                                                                         cAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTg                                                                              <  2:1710518/150‑1 (MQ=255)

                                                                         cAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTg                                                                              <  1:1000632/150‑1 (MQ=255)

                                                                            cGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCa                                                                           <  1:3935796/150‑1 (MQ=255)

                                                                             ggTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCaa                                                                          <  2:1602527/150‑1 (MQ=255)

                                                                               ttCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGt                                                                        >  1:799485/1‑150 (MQ=255)

                                                                                   aGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTtcgt                                                                    >  2:4674912/1‑150 (MQ=255)

                                                                                    gCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTtcgtc                                                                   <  1:5438267/150‑1 (MQ=255)

                                                                                     cAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTtcgtcg                                                                  >  1:6455783/1‑150 (MQ=255)

                                                                                       aCTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAg                                                                <  2:7171140/150‑1 (MQ=255)

                                                                                        cTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAgg                                                               <  1:1145256/150‑1 (MQ=255)

                                                                                         tGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAg                                                                >  2:2257793/1‑148 (MQ=255)

                                                                                           aCCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTg                                                            >  1:2825368/1‑150 (MQ=255)

                                                                                           aCCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTg                                                            <  2:5888164/150‑1 (MQ=255)

                                                                                           aCCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTg                                                            <  1:5949038/150‑1 (MQ=255)

                                                                                            ccGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTgg                                                           <  2:273817/150‑1 (MQ=255)

                                                                                             cGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGgc                                                          <  2:1939894/150‑1 (MQ=255)

                                                                                             cGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGgc                                                          >  2:3323653/1‑150 (MQ=255)

                                                                                               gACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTc                                                                                  <  1:4755333/124‑1 (MQ=255)

                                                                                               gACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGgcgc                                                        <  2:2474514/150‑1 (MQ=255)

                                                                                                     acaCGGATTTACTCGCGGAGCATTTCAATGCGC‑AGGCGCAGCAGCGCATAGCGCAGTTCGCGCGGGATACGCGAGTTTCGAAGTTCTGGCAGGCTCATTGCAACAGGCGCACCATCTACTTGGAAAGTACGTCGAGGATCGAGCCGAAgg                                                  >  2:5626532/3‑150 (MQ=255)

                                                                                                      gACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAggg                                                 >  1:7426269/1‑150 (MQ=255)

                                                                                                       aCGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAg                                                                <  1:674578/134‑1 (MQ=255)

                                                                                                         ggATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCga                                              <  1:1634839/150‑1 (MQ=255)

                                                                                                         ggATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCga                                              <  2:4788467/150‑1 (MQ=255)

                                                                                                            tttCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCgatga                                           <  1:4992893/150‑1 (MQ=255)

                                                                                                              tCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGAtt                                         >  1:2407583/1‑150 (MQ=255)

                                                                                                               ccTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATtc                                        <  2:4214372/150‑1 (MQ=255)

                                                                                                                 tGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCga                                              <  2:4892846/142‑1 (MQ=255)

                                                                                                                  ggCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATtcgtc                                     <  1:7262090/150‑1 (MQ=255)

                                                                                                                     ggTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATTCGTCGaa                                  <  1:1942831/150‑1 (MQ=255)

                                                                                                                         cATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTc                                                                                        >  2:594425/1‑92 (MQ=255)

                                                                                                                           gttCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTTGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATTCGTCGAACAccgc                            <  1:5191152/149‑1 (MQ=255)

                                                                                                                             tCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATTCGTCGAACAccgccg                          >  2:3538100/1‑150 (MQ=255)

                                                                                                                             tCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATTCGTCGAACAccgccg                          >  1:6911651/1‑150 (MQ=255)

                                                                                                                             gaactGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATTCGTCGAACAccgccg                          >  2:5257187/5‑150 (MQ=255)

                                                                                                                               cATGCGC‑AGGAGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGTCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATTCGTCGAACACCgccgcc                        <  1:5626532/150‑1 (MQ=255)

                                                                                                                                ctGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAg                                                                >  1:1334675/2‑109 (MQ=255)

                                                                                                                                ctGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAg                                                                <  2:1334678/108‑1 (MQ=255)

                                                                                                                                aTGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCgat                                             >  2:6480527/1‑128 (MQ=255)

                                                                                                                                 tGCGCAAGGCGCTGCTGCGCC‑TGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATTCGTCGAACACCGCCGCCtg                      >  2:6348596/1‑149 (MQ=255)

                                                                                                                                 tGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATTCGTCGAACACCGCCGCCtt                      <  1:6509858/150‑1 (MQ=255)

                                                                                                                                   cgc‑aggcgcagCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATTCGTCGAACACCGCCGCCTTCa                    <  2:5419485/150‑1 (MQ=255)

                                                                                                                                   cgc‑aggcgcagCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATTCGTCGAACACCGCCGCCTTCa                    <  2:4461682/150‑1 (MQ=255)

                                                                                                                                    gcaaggcgctgCTGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATTCGTCGAACACCGCCGCCTTCa                    >  1:7316657/4‑150 (MQ=255)

                                                                                                                                    gc‑aggcgcagCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATTCGTCGAACACCGCCGCCTTCAg                   >  1:1797246/1‑150 (MQ=255)

                                                                                                                                     caaggcgctgCTGCGCC‑TGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATTCGTCGAACACCGCCGCCTTCAgg                  <  1:5369080/148‑1 (MQ=255)

                                                                                                                                       aggcgcagCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATTCGTCGAACACCGCCGCCTTCAGgc                 <  2:346380/150‑1 (MQ=255)

                                                                                                                                           gcagcaGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCCAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATTCGTCGAACACCACCGCCTTCAGGCGCAc             >  2:707413/1‑150 (MQ=255)

                                                                                                                                              gcagcGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATTCGTCGAACACCGCCGCCTTCAGGCGCAcgcc          <  2:5328873/150‑1 (MQ=255)

                                                                                                                                                 gcgcCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATTCGTCGAACACCGCCGCCTTCAGGCGCACGCCGCa       <  1:4634540/150‑1 (MQ=255)

                                                                                                                                                    ccTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATTCGTCGAACACCGCCGCCTTCAGGCGCACGCCGCAGTa    >  1:2962759/1‑150 (MQ=255)

CGATGCCCCACAGTGGCGCATGCTTGATACCCAGGCCTTCATGTAGCGAGCTGCCCATACCGCGCACGCGGCGCAGCGGTTCAAGCAACTGACCGGACTCATGACGGATTTCCTGGCGGTGCATTTCCATGCGC‑AGGCGCAGCAGCGCCTTGCGCAGTTCGCGCGGGTTACGCGTGTTTGGCAGTTCTGGCAGGCTCATGGCAACAGGCGCTCCCTGTCCTTGGCAAGTTCGTCGAGGGTGGCGCCGAAGGGCGATGATTCGTCGAACACCGCCGCCTTCAGGCGCACGCCGCAGTACA  >  NC_002947/2407160‑2407458

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