BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,299,504	T→C	5.5%	E861E (GAA→GAG)	sbcC ←	exonuclease subunit SbcC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,299,504	0	T	C	5.5%	576.7 / 13.8	181	E861E (GAA→GAG)	sbcC	exonuclease subunit SbcC
Reads supporting (aligned to +/- strand): ref base T (78/92); new base C (3/7); total (81/100)
Fisher's exact test for biased strand distribution p-value = 5.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.46e-01

Only 100 of 189 total aligned reads displayed.

TGGCCTTGAGCTCGCTGGCCAGTTGCACGCCTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTCCAGTTCGTC  >  NC_002947/2299356‑2299653

  gCCTTGAGCTCGCTGGCCAGTTGCACGCCTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCg                                                                                                                                                    <  1:506058/150‑1 (MQ=255)

  gCCTTGAGCTCGCTGGCCAGTTGCACGCCTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCg                                                                                                                                                    >  2:4141364/1‑150 (MQ=255)

   ccTTGAGCTCGCTGGCCAGTTGCACGCCTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAACTGCTCgg                                                                                                                                                   <  1:5205940/150‑1 (MQ=255)

    cTTGAGCTCGCTGGCCAGTTGCACGCCTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGc                                                                                                                                                  >  1:4411757/1‑150 (MQ=255)

    cTTGAGCTCGCTGGCCAGTTGCACGCCTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGc                                                                                                                                                  <  2:2607970/150‑1 (MQ=255)

    cTTGAGCTCGCTGGCCAGTTGCACGCCTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGc                                                                                                                                                  >  1:988215/1‑150 (MQ=255)

    cTTGAGCTCGCTGGCCAGTTGCACGCCTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGc                                                                                                                                                  <  1:1001075/150‑1 (MQ=255)

    cTTGAGCTCGCTGGCCAGTTGCACGCCTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGc                                                                                                                                                  <  2:1224329/150‑1 (MQ=255)

       gAGCTCGCTGGCCAGTTGCACGCCTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAACTGCTCGGCGaa                                                                                                                                               <  2:3212145/150‑2 (MQ=255)

        aGCTCGCTGGCCAGTTGCACGCCTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGt                                                                                                                                              <  2:5370134/150‑1 (MQ=255)

        aGCTCGCTGGCCAGTTGCACGCCTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAACTGCTCGGCGaac                                                                                                                                              <  1:4084054/150‑3 (MQ=255)

         gCTCGCTGGCCAGTTGCACGCCTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGGGTCAAGGGTGCGGGCGTGCTCCAGCGGGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCAGGTGCGCCGAGCGGTTCGGCGCGTg                                                                                                                                             >  2:5022344/1‑150 (MQ=255)

           tCGCTGGCCAGTTGCACGCCTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAACTGCTCGGCGaacatg                                                                                                                                           <  1:1322508/150‑6 (MQ=255)

             gCTGGCCAGTTGCACGCCTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGc                                                                                                                                         >  2:4446154/1‑150 (MQ=255)

               tGGCCAGTTGCACGCCTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCa                                                                                                                                       >  2:482997/1‑150 (MQ=255)

                ggCCAGTTGCACGCCTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAg                                                                                                                                      <  1:1660314/150‑1 (MQ=255)

                ggCCAGTTGCACGCCTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAg                                                                                                                                      <  1:1460080/150‑1 (MQ=255)

                 gCCAGTTGCACGCCTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAgg                                                                                                                                     >  2:2191313/1‑150 (MQ=255)

                  ccAGTTGCACGCCTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAggg                                                                                                                                    <  2:216498/150‑1 (MQ=255)

                   cAGTTGCACGCGTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAgggg                                                                                                                                   >  2:155140/1‑149 (MQ=255)

                   cAGTTGCACGCCTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGc                                                                                                                                   >  1:4465892/1‑150 (MQ=255)

                    aGTTGCACGCCTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCt                                                                                                                                  <  2:4515129/150‑1 (MQ=255)

                       tGCACGCCTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGc                                                                                                                               <  1:12579/150‑1 (MQ=255)

                        gCACGCCTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCt                                                                                                                              >  2:4191823/1‑150 (MQ=255)

                            gCCTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCTCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCGTGCTGGCg                                                                                                                          >  2:2944295/1‑150 (MQ=255)

                               ttGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTg                                                                                                                       >  2:3244974/1‑150 (MQ=255)

                                 gggTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCt                                                                                                                     >  1:916767/1‑150 (MQ=255)

                                   gTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCtcg                                                                                                                   >  2:1648706/1‑150 (MQ=255)

                                         aGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAgg                                                                                                             <  1:1026282/150‑1 (MQ=255)

                                         aGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAgg                                                                                                             <  2:1953591/150‑1 (MQ=255)

                                             cTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGggcgg                                                                                                         >  2:5444731/1‑150 (MQ=255)

                                                  aGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGt                                                                                                    <  2:3407477/150‑1 (MQ=255)

                                                    gcgcTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGttt                                                                                                  <  2:2832923/150‑1 (MQ=255)

                                                    gcgcTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGttt                                                                                                  <  2:1056882/150‑1 (MQ=255)

                                                    gcgcTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGttt                                                                                                  >  2:839477/1‑150 (MQ=255)

                                                     cgcTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTc                                                                                                 <  2:2303212/150‑1 (MQ=255)

                                                      gcTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTct                                                                                                >  1:1709575/1‑150 (MQ=255)

                                                         gggCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTctgct                                                                                             <  1:5253519/150‑1 (MQ=255)

                                                            cGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTa                                                                                          >  2:4950174/1‑150 (MQ=255)

                                                             ggTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTaa                                                                                         <  2:1120101/150‑1 (MQ=255)

                                                                   ggCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCa                                                                                                      >  1:3165777/1‑131 (MQ=255)

                                                                       tCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTctg                                                                                               >  2:4695392/1‑134 (MQ=255)

                                                                       tCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTgg                                                                               <  2:5181936/150‑1 (MQ=255)

                                                                            ggTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTg                                                                          >  1:4115320/1‑150 (MQ=255)

                                                                                 gggCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCAcgc                                                                     >  2:761182/1‑150 (MQ=255)

                                                                                  ggCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAggg                                                                                                            >  2:2405601/1‑110 (MQ=255)

                                                                                   gCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCtt                                                                           <  1:5353415/142‑1 (MQ=255)

                                                                                   gCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCtt                                                                           >  2:5353415/1‑142 (MQ=255)

                                                                                   gCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCAcgcgc                                                                   <  2:361595/150‑1 (MQ=255)

                                                                                   gCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCAcgcgc                                                                   >  1:2390717/1‑150 (MQ=255)

                                                                                     gTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGt                                                                                                                                              >  1:14610/1‑73 (MQ=255)

                                                                                     gTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGt                                                                                                                                              <  2:14610/73‑1 (MQ=255)

                                                                                     gTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACgcgcgc                                                                 <  1:931198/150‑1 (MQ=255)

                                                                                     gTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACgcgcgc                                                                 <  1:5022344/150‑1 (MQ=255)

                                                                                         tCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCctgc                                                             >  2:3828169/1‑150 (MQ=255)

                                                                                         tCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCctgc                                                             >  2:1939212/1‑150 (MQ=255)

                                                                                           cAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCctgctg                                                           >  2:2068214/1‑150 (MQ=255)

                                                                                            aGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGt                                                          >  1:22960/1‑150 (MQ=255)

                                                                                            aGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGt                                                          <  1:5071649/150‑1 (MQ=255)

                                                                                             gcgcGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGTTGGCCTTGCACGCGCGCctgcaaga                                                         <  2:1325286/150‑5 (MQ=255)

                                                                                             gcgcGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGtt                                                         <  1:4228702/150‑1 (MQ=255)

                                                                                             gcgcGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGtt                                                         <  2:22960/150‑1 (MQ=255)

                                                                                              cgcgGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCgg                                                                                                                                                   >  2:2274192/1‑59 (MQ=255)

                                                                                              cgcgGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTg                                                        <  1:1073516/150‑1 (MQ=255)

                                                                                                cgGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGt                                                      <  1:1972554/150‑1 (MQ=255)

                                                                                                 ggTGTCCACGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGGTCGGCGAGTGAGACCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTc                                                     <  2:3302236/150‑1 (MQ=255)

                                                                                                          gtgCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCtt                                            <  1:4550217/150‑1 (MQ=255)

                                                                                                             ctgctgcCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTc                                         >  2:5391951/1‑150 (MQ=255)

                                                                                                               gctgcCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCa                                       >  1:884913/1‑150 (MQ=255)

                                                                                                                   ccAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTg                                   >  1:5307895/1‑150 (MQ=255)

                                                                                                                   ccAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTg                                   <  1:4191823/150‑1 (MQ=255)

                                                                                                                     aGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTggg                                 >  2:3989962/1‑150 (MQ=255)

                                                                                                                     aGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTGTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTggg                                 >  2:4188119/1‑150 (MQ=255)

                                                                                                                      ggCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGt                                <  2:2719485/150‑1 (MQ=255)

                                                                                                                         ttCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGGGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTg                             >  1:1985078/1‑150 (MQ=255)

                                                                                                                         ttCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTg                             >  1:1885490/1‑150 (MQ=255)

                                                                                                                         ttCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTg                             >  2:873915/1‑150 (MQ=255)

                                                                                                                          tCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTgg                            <  2:1876291/150‑1 (MQ=255)

                                                                                                                            ggCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTg                                                        >  1:3213682/1‑120 (MQ=255)

                                                                                                                            ggCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCg                          >  1:260313/1‑150 (MQ=255)

                                                                                                                             gcgcTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCgg                         <  1:1648595/150‑1 (MQ=255)

                                                                                                                              cgcTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCggg                        <  1:4817411/150‑1 (MQ=255)

                                                                                                                                 tGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCtt                     >  1:3736553/1‑150 (MQ=255)

                                                                                                                                   gCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTtgc                   <  1:3628450/150‑1 (MQ=255)

                                                                                                                                       gTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTgctgct               <  2:3177416/150‑1 (MQ=255)

                                                                                                                                       gTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTgctgct               >  1:1098756/1‑150 (MQ=255)

                                                                                                                                       actCGCCGAACTGCTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTtgct                  >  1:4698753/3‑147 (MQ=255)

                                                                                                                                       actCGCCGAACTGCTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTgctgct               >  1:5391542/3‑150 (MQ=255)

                                                                                                                                             cGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTCCa         <  1:4950174/150‑1 (MQ=255)

                                                                                                                                             cGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACACAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTCCa         <  1:761182/150‑1 (MQ=255)

                                                                                                                                               aGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCTGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTCCAGt       <  1:1762469/150‑1 (MQ=255)

                                                                                                                                                 cAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTCCAGTTc     <  1:3039900/150‑1 (MQ=255)

                                                                                                                                                 cAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTCCAGTTc     <  2:5307895/150‑1 (MQ=255)

                                                                                                                                                 cAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTCCAGTTc     <  2:5307991/150‑1 (MQ=255)

                                                                                                                                                    ttCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGt                                                          <  2:1642204/94‑1 (MQ=255)

                                                                                                                                                    ttCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCgg                         >  1:91505/1‑127 (MQ=255)

                                                                                                                                                    ttCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCgg                         <  2:91505/127‑1 (MQ=255)

                                                                                                                                                    ttCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTCCAGTTCGTc  <  1:3240562/150‑1 (MQ=255)

                                                                                                                                                    ttCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTCCAGTTCGTc  >  2:2984897/1‑150 (MQ=255)

                                                                                                                                                    ttCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGATAGTCGAGGGCGGCCAGTTTCTGCTGTAACTGATGCTGGCATTGAACGCGCGACTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTCCAGTTCGTc  >  1:3509611/1‑150 (MQ=255)

TGGCCTTGAGCTCGCTGGCCAGTTGCACGCCTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTCCAGTTCGTC  >  NC_002947/2299356‑2299653

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