BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,471,402	G→T	6.2%	R478R (CGG→AGG)	alg8 ←	glycosyltransferase Alg8

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,471,402	0	G	T	6.2%	572.3 / 18.4	172	R478R (CGG→AGG)	alg8	glycosyltransferase Alg8
Reads supporting (aligned to +/- strand): ref base G (75/86); new base T (3/8); total (78/94)
Fisher's exact test for biased strand distribution p-value = 3.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 175 total aligned reads displayed.

GGCGCTGGGCTTCGGACTCATGCACGACATTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATAGAGAATCAGCG  >  NC_002947/1471254‑1471551

ggCGCTGGGCTTCGGACTCATGCACGACATTGACGTTCACGGCGGTCTTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCcgc                                                                                                                                                      >  2:3117639/1‑150 (MQ=255)

  cgcTGGGCTTCGGACTCATGCACGACATTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACTGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGa                                                                                                                                                    <  2:2307137/150‑1 (MQ=255)

     tGGGCTTCGGACTCATGCACGACATTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGa                                                                                                                                                 <  1:4514458/150‑1 (MQ=255)

      gggCTTCGGACTCATGCACGACATTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGAc                                                                                                                                                <  1:3701077/150‑1 (MQ=255)

         cTTCGGACTCATGCACGACATTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAg                                                                                                                                             <  2:4973977/150‑1 (MQ=255)

          ttCGGACTCATGCACGACATTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCTCGCGGACCatg                                                                                                                                            <  1:4430219/150‑3 (MQ=255)

           tCGGACTCATGCACGACATTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCTCGCGGACCatga                                                                                                                                           <  1:4389754/150‑4 (MQ=255)

           tCGGACTCATGCACGACATTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCTCGCGGACCatga                                                                                                                                           <  1:1764455/150‑4 (MQ=255)

           tCGGACTCATGCACGACATTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCTCGCGGACCatga                                                                                                                                           <  1:3006286/150‑4 (MQ=255)

            cGGACTCATGCACGACATTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGgtg                                                                                                                                          >  2:2411336/1‑150 (MQ=255)

              gACTCATGCACGACATTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGtt                                                                                                                                        <  2:952050/150‑1 (MQ=255)

                cTCATGCACGACATTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGa                                                                                                                                      <  1:2426660/150‑1 (MQ=255)

                cTCATGCACGACATTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGa                                                                                                                                      <  2:424785/150‑1 (MQ=255)

                 tCATGCACGACATTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCGGTTAGATCCGCTCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGaa                                                                                                                                     >  1:5083047/1‑150 (MQ=255)

                   aTGCACGACATTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAAcc                                                                                                                                   >  2:272657/1‑150 (MQ=255)

                   aTGCACGACATTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAAcc                                                                                                                                   <  1:166253/150‑1 (MQ=255)

                     gCACGACATTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCAt                                                                                                                                 <  2:5415847/150‑1 (MQ=255)

                       aCGACATTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCAttg                                                                                                                               >  2:2115743/1‑150 (MQ=255)

                       aCGACATTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCAttg                                                                                                                               >  1:5322146/1‑150 (MQ=255)

                         gACATTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCTCGCGGACCatga                                                                                                                                           >  1:2729581/1‑133 (MQ=255)

                         gACATTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCAttgtt                                                                                                                             <  1:2690474/150‑1 (MQ=255)

                            aTTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGaa                                                                                                                          <  1:5385062/150‑1 (MQ=255)

                             ttGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGaaa                                                                                                                         <  1:2200215/150‑1 (MQ=255)

                             ttGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTAGaaa                                                                                                                         <  2:3484370/150‑1 (MQ=255)

                                 cGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTg                                                                                                                     >  1:2546827/1‑150 (MQ=255)

                                  gTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTgg                                                                                                                    >  2:3784841/1‑150 (MQ=255)

                                  gTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTgg                                                                                                                    >  1:695765/1‑150 (MQ=255)

                                     cACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGa                                                                                                                 <  2:1046738/150‑1 (MQ=255)

                                      aCGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAg                                                                                                                <  1:4781021/150‑1 (MQ=255)

                                      aCGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAg                                                                                                                <  2:2293288/150‑1 (MQ=255)

                                      aCGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAg                                                                                                                >  2:248950/1‑150 (MQ=255)

                                        ggcggTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGt                                                                                                              <  2:2850020/150‑1 (MQ=255)

                                                 cATGGTGGTTTGTTCCTGTTAGATCCGATCCGTGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGa                                                                                                     <  1:1952127/150‑1 (MQ=255)

                                                  aTGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAg                                                                                                    <  2:1737532/150‑1 (MQ=255)

                                                    ggtggtTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAg                                                                                                                                             <  1:2997782/107‑1 (MQ=255)

                                                     gtggtTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCCTGGTCCTCGCGGACCatga                                                                                                                                           >  1:627197/1‑105 (MQ=255)

                                                     gtggtTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCCTGGTCCTCGCGGACCatga                                                                                                                                           <  2:627197/108‑4 (MQ=255)

                                                       ggtTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAg                                                                                               >  2:92840/1‑150 (MQ=255)

                                                        gtTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGt                                                                                              <  2:964573/150‑1 (MQ=255)

                                                        gtTTGTTCATGTTAGATCCGATCCTGTTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTTAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGt                                                                                              <  1:937699/150‑1 (MQ=255)

                                                         tttGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTg                                                                                             >  1:3411177/1‑150 (MQ=255)

                                                          ttGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTgg                                                                                            <  2:4642981/150‑1 (MQ=255)

                                                            gTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGc                                                                                          >  1:486189/1‑150 (MQ=255)

                                                             ttCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCt                                                                                         >  2:4604582/1‑150 (MQ=255)

                                                             ttCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCt                                                                                         >  1:3296579/1‑150 (MQ=255)

                                                             ttCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCt                                                                                         <  1:2153392/150‑1 (MQ=255)

                                                               ccTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGTTGa                                                                                       >  2:338063/1‑150 (MQ=255)

                                                                 tGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAGCCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGAcg                                                                                     >  1:4077262/1‑150 (MQ=255)

                                                                 tGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGAcg                                                                                     >  1:4300810/1‑150 (MQ=255)

                                                                 tGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGAcg                                                                                     >  1:1945282/1‑150 (MQ=255)

                                                                   ttAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGa                                                                                                     >  1:1891130/1‑132 (MQ=255)

                                                                   ttAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGa                                                                                                     <  2:1891130/132‑1 (MQ=255)

                                                                   ttAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGAcgcg                                                                                   >  1:505106/1‑150 (MQ=255)

                                                                   ttAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGAcgcg                                                                                   >  2:1758104/1‑150 (MQ=255)

                                                                       atccgatccgGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGgt                                                                                                                                           <  1:3109936/90‑1 (MQ=255)

                                                                       atccgatccgGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCa                                                                               <  1:5218007/150‑1 (MQ=255)

                                                                        tccgatccgGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTggg                                                                                           <  2:4523555/137‑1 (MQ=255)

                                                                        tccgatccgGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTggg                                                                                           >  1:4523555/1‑137 (MQ=255)

                                                                         ccgatccgGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACg                                                                             >  1:1139025/1‑150 (MQ=255)

                                                                          cgatccgGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGa                                                                            <  2:1358677/150‑1 (MQ=255)

                                                                          cgatccgGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGCCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGa                                                                            <  2:3004790/150‑1 (MQ=255)

                                                                           gatccgGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGAc                                                                           >  1:3730215/1‑150 (MQ=255)

                                                                              ccgGGTTGGGCTCACACGACAATTAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGc                                                                        <  1:3632462/150‑1 (MQ=255)

                                                                                gggTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCtt                                                                      <  1:5089390/150‑1 (MQ=255)

                                                                                gggTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCtt                                                                      <  1:4204211/150‑1 (MQ=255)

                                                                                gggTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCtt                                                                      >  2:4220988/1‑150 (MQ=255)

                                                                                  gTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGt                                                                    >  2:1952718/1‑150 (MQ=255)

                                                                                   ttGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTc                                                                   <  2:3730215/150‑1 (MQ=255)

                                                                                       gCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAgg                                                               <  2:1628506/150‑1 (MQ=255)

                                                                                            cacGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGgaa                                                          <  1:4235565/150‑1 (MQ=255)

                                                                                            cacGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGgaa                                                          <  1:1758104/150‑1 (MQ=255)

                                                                                             acGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGgaag                                                         >  1:539958/1‑150 (MQ=255)

                                                                                               gACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGgaagaa                                                       >  1:4413992/1‑150 (MQ=255)

                                                                                                 ccATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAaca                                                     <  2:2988115/150‑1 (MQ=255)

                                                                                                  cATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAacac                                                    <  2:3234926/150‑1 (MQ=255)

                                                                                                   aTGAACCGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACg                                                   <  1:4648715/150‑1 (MQ=255)

                                                                                                     gAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTa                                                 >  2:1152866/1‑150 (MQ=255)

                                                                                                      aaCAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGGGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAg                                                <  2:2826024/150‑1 (MQ=255)

                                                                                                      aaCAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAg                                                <  1:1849556/150‑1 (MQ=255)

                                                                                                          gcacagcaACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCt                                            >  2:752711/1‑150 (MQ=255)

                                                                                                          gcacagcaACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCt                                            <  1:1274513/150‑1 (MQ=255)

                                                                                                           cacagcaACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCtt                                           <  1:832895/150‑1 (MQ=255)

                                                                                                            acagcaACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTtc                                          >  1:4115119/1‑150 (MQ=255)

                                                                                                              agcaACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTtcat                                        <  2:4115119/150‑1 (MQ=255)

                                                                                                                 aaCGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTtcatca                                      <  2:2215059/149‑1 (MQ=255)

                                                                                                                 aaCGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAg                                     <  1:3193724/150‑1 (MQ=255)

                                                                                                                    gAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAgcgc                                  >  1:4117396/1‑150 (MQ=255)

                                                                                                                           cTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCgatg                           >  2:1653725/1‑150 (MQ=255)

                                                                                                                              gCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATc                        <  1:2262908/150‑1 (MQ=255)

                                                                                                                                     agaAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTggtggt                 >  2:4812311/1‑150 (MQ=255)

                                                                                                                                     agaAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGt                 <  1:4871800/150‑1 (MQ=255)

                                                                                                                                      gaAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCgat                            <  1:982988/138‑1 (MQ=255)

                                                                                                                                         ggTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATa             <  1:248950/150‑1 (MQ=255)

                                                                                                                                           tCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATaga           >  2:623877/1‑150 (MQ=255)

                                                                                                                                           tCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATaga           <  1:1572727/150‑1 (MQ=255)

                                                                                                                                             ctGGTCCTCGCGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCt                                                                                         <  1:464367/69‑1 (MQ=255)

                                                                                                                                               ggTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATAGAGAAt       >  2:4437957/1‑150 (MQ=255)

                                                                                                                                                 tCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTGGTAATAGAGAATCa     >  2:5355569/1‑150 (MQ=255)

                                                                                                                                                   cgcgAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATAGAGAATCAGc   >  2:3280493/1‑150 (MQ=255)

                                                                                                                                                    gcgAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATAGAGAATCAGCg  <  2:4392988/150‑1 (MQ=255)

GGCGCTGGGCTTCGGACTCATGCACGACATTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATAGAGAATCAGCG  >  NC_002947/1471254‑1471551

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