BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,066,366	G→C	5.8%	V153L (GTG→CTG)	PP_2678 →	hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,066,366	0	G	C	5.8%	326.8 / 9.5	104	V153L (GTG→CTG)	PP_2678	hydrolase
Reads supporting (aligned to +/- strand): ref base G (51/47); new base C (2/4); total (53/51)
Fisher's exact test for biased strand distribution p-value = 4.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.79e-01

Only 100 of 110 total aligned reads displayed.

TGGGCAACCAGGCCTTTGCCGATGTGCCGATCGGTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCGGGGTGCTGTTCGCCGGCGACCTGGTGTTCTACCAACGCGCCCTGA  >  NC_002947/3066217‑3066513

tGGGCAACCAGGCCTTTGCCGATGTGCCGATCGGTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCggg                                                                                                                                                     >  1:755034/1‑150 (MQ=255)

tGGGCAACCAGGCCTTTGCCGATGTGCCGATCGGTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCggg                                                                                                                                                     >  2:653713/1‑150 (MQ=255)

tGGGCAACCAGGCCTTTGCCGATGTGCCGATCGGTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCggg                                                                                                                                                     >  2:2034045/1‑150 (MQ=255)

   gCAACCAGGCCTTTGCCGATGTGCCGATCGGTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGtgc                                                                                                                                                  <  2:1134646/150‑1 (MQ=255)

   gCAACCAGGCCTTTGCCGATGTGCCGATCGGTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGtgc                                                                                                                                                  >  2:288737/1‑150 (MQ=255)

      aCCAGGCCTTTGCCGATGTGCCGATCGGTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGtgctgc                                                                                                                                               >  1:1754224/1‑150 (MQ=255)

       ccAGGCCTTTGCCGATGTGCCGATCGGTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCa                                                                                                                                              >  1:112819/1‑150 (MQ=255)

       ccAGGCCTTTGCCGATGTGCCGATCGGTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCa                                                                                                                                              <  1:1777673/150‑1 (MQ=255)

        cAGGCCTTTGCCGATGTGCCGATCGGTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGCTgcagcac                                                                                                                                             >  2:2060940/1‑149 (MQ=255)

            ccTTTGCCGATGTGCCGATCGGTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCgg                                                                                                                                         <  2:359596/150‑1 (MQ=255)

            ccTTTGCCGATGTGCCGATCGGTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGCTGCAGCACCCgg                                                                                                                                         <  2:649256/150‑1 (MQ=255)

            ccTTTGCCGATGTGCCGATCGGTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGCTGCAGCACCCgg                                                                                                                                         <  1:1324146/150‑1 (MQ=255)

             cTTTGCCGATGTGCCGATCGGTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCggg                                                                                                                                        >  1:1044509/1‑150 (MQ=255)

              tttGCCGATGTGCCGATCGGTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGTTGCCGACCCGGGTGCTGCAGCCGGGc                                                                                                                                       <  1:321529/150‑1 (MQ=255)

              tttGCCGATGTGCCGATCGGTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGc                                                                                                                                       >  1:764920/1‑150 (MQ=255)

               ttGCCGATGTGCCGATCGGTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGCTGCAGCACCCgggtc                                                                                                                                      <  1:2060940/150‑3 (MQ=255)

                  ccGATGTGCCGATCGGTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGc                                                                                                                                   <  1:2169349/150‑1 (MQ=255)

                   cGATGTGCCGATCGGTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCGGGTCGTGCa                                                                                                                                  >  1:1922281/1‑150 (MQ=255)

                   cGATGTGCCGATCGGTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCa                                                                                                                                  <  2:1124900/150‑1 (MQ=255)

                    gATGTGCCGATCGGTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCAt                                                                                                                                 >  1:399109/1‑150 (MQ=255)

                      tgtgCCGATCGGTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGa                                                                                                                               <  1:907909/150‑1 (MQ=255)

                          ccGATCGGTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAggtg                                                                                                                           >  2:934989/1‑150 (MQ=255)

                           cGATCGGTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAggtgg                                                                                                                          <  1:13934/150‑1 (MQ=255)

                            gATCGGTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGACGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAggtggg                                                                                                                         <  2:399109/150‑1 (MQ=255)

                             aTCGGTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAggtgggt                                                                                                                        >  1:588050/1‑150 (MQ=255)

                             aTCGGTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAggtgggt                                                                                                                        <  2:2218077/150‑1 (MQ=255)

                             aTCGGTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAggtgggg                                                                                                                        >  2:118838/1‑149 (MQ=255)

                                  tGCCCTGGCTGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCa                                                                                                                   >  2:2133941/1‑150 (MQ=255)

                                     ccTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACACTTGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCg                                                                                                                >  1:623574/1‑150 (MQ=255)

                                      cTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGc                                                                                                               >  1:1045071/1‑150 (MQ=255)

                                          ccGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGa                                                                                                           <  1:446523/150‑1 (MQ=255)

                                           cGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAg                                                                                                          >  2:1448159/1‑150 (MQ=255)

                                           cGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAg                                                                                                          <  2:216452/150‑1 (MQ=255)

                                           cGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAg                                                                                                          >  1:53901/1‑150 (MQ=255)

                                               aCCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCTCTGCCTGAGGTTg                                                                                                      >  1:491198/1‑150 (MQ=255)

                                               aCCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTg                                                                                                      <  2:2284491/150‑1 (MQ=255)

                                               aCCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTg                                                                                                      >  1:492730/1‑150 (MQ=255)

                                                  cgcgACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTgctg                                                                                                   <  1:2133941/150‑1 (MQ=255)

                                                  cgcgACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGCTGCAGCACCCgggtc                                                                                                                                      >  2:91445/1‑113 (MQ=255)

                                                   gcgACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTgctgg                                                                                                  <  1:890516/150‑1 (MQ=255)

                                                     gACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGGGGGCACCGCCTGAGGTTgctgggg                                                                                                >  1:783206/1‑150 (MQ=255)

                                                         tgctgcGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGtt                                                                                                       <  2:1410447/139‑1 (MQ=255)

                                                         tgctgcGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGtt                                                                                                       >  1:1410447/1‑139 (MQ=255)

                                                           ctgcGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGGGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTgggc                                                                                          >  2:1491353/1‑150 (MQ=255)

                                                            tgcGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTgggcg                                                                                         <  1:288737/150‑1 (MQ=255)

                                                             gcgcCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTgggcgg                                                                                        <  2:612276/150‑1 (MQ=255)

                                                                ccAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGAGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCa                                                                                     >  2:1726211/1‑150 (MQ=255)

                                                                 cAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCAt                                                                                    >  1:1940784/1‑150 (MQ=255)

                                                                     aaGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACgg                                                                                <  2:1998724/150‑1 (MQ=255)

                                                                       ggCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGt                                                                              >  2:2184992/1‑150 (MQ=255)

                                                                          gATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCCAGAAAGCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATAAGTTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGcc                                                                           >  1:2275161/1‑150 (MQ=255)

                                                                             gCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGAc                                                                        <  1:892952/150‑1 (MQ=255)

                                                                              ccATGGCCGAGAACATGTATCGCCTGGTCGGCTACTGGATGCGCGGTACCTAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGACTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGAcc                                                                       >  1:1177486/1‑150 (MQ=255)

                                                                               cATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGAGGGCATACGGGTGCCGACCt                                                                      >  2:436217/1‑150 (MQ=255)

                                                                               cATGGCCGAGAACATGTATCGCCTGGTCGGCGAATGGATGCGCGGTAACGAGGTGGGGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCt                                                                      >  1:46947/1‑150 (MQ=255)

                                                                                 tGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTgg                                                                    >  1:2285463/1‑150 (MQ=255)

                                                                                  ggCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGc                                                                   <  1:2184992/150‑1 (MQ=255)

                                                                                     cGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGAt                                                                >  2:1563809/1‑150 (MQ=255)

                                                                                     cGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGAt                                                                >  2:1696654/1‑150 (MQ=255)

                                                                                      gagaACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAggtgggtg                                                                                                                       <  2:2260143/94‑1 (MQ=255)

                                                                                      gagaACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAggtgggtg                                                                                                                       >  1:2260143/1‑94 (MQ=255)

                                                                                         aaCATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTCCGGCTGGCCGGGCATACGGGTGCCGACCCGGCGATTTTc                                                            <  2:394972/150‑1 (MQ=255)

                                                                                           cATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGa                                                          <  1:228275/150‑1 (MQ=255)

                                                                                           cATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGa                                                          >  1:1702645/1‑150 (MQ=255)

                                                                                           cATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGa                                                          <  1:934989/150‑1 (MQ=255)

                                                                                           cATGTATCGCCTGGTCGGCGACTCGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTCTCGa                                                          >  2:1184770/1‑150 (MQ=255)

                                                                                               tATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGaa                                                      >  1:1435315/1‑150 (MQ=255)

                                                                                                 tCGCCTGGTTGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGaaaa                                                    >  2:412/1‑150 (MQ=255)

                                                                                                    ccTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGAc                                                 <  2:112819/150‑1 (MQ=255)

                                                                                                        gTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGGGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCggg                                             <  1:1491353/150‑1 (MQ=255)

                                                                                                        gTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCggg                                             <  2:1657803/150‑1 (MQ=255)

                                                                                                         tCGGCGCCTGGATGCGCGGTACCGCGGTGGTGTTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCgggg                                            <  2:622760/150‑1 (MQ=255)

                                                                                                         tCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGc                                                                                                                                   >  1:2048942/1‑63 (MQ=255)

                                                                                                         tCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGc                                                                                                                                   <  2:2048942/63‑1 (MQ=255)

                                                                                                         tCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCgggg                                            <  2:764920/150‑1 (MQ=255)

                                                                                                          cGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCGGGGt                                           >  2:2083343/1‑150 (MQ=255)

                                                                                                             cGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCGGGGTGCt                                        <  2:2153798/150‑1 (MQ=255)

                                                                                                              gACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGAGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCGGGGTGCTg                                       <  1:1726211/150‑1 (MQ=255)

                                                                                                                   gATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCGGGGTGCTGTTCGc                                  >  1:1791221/1‑150 (MQ=255)

                                                                                                                    aTGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCGGGGTGCTGTTCGcc                                 >  1:1097344/1‑150 (MQ=255)

                                                                                                                      gcgcGGTTCCGGGGTGGTGCTGCCTACCCGAGTGCGCCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCGGGGTGCTGTTCGCCgg                               <  2:1097344/150‑1 (MQ=255)

                                                                                                                      gcgcGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCGGGGTGCTGTTCGCCgg                               >  2:1289822/1‑150 (MQ=255)

                                                                                                                       cgcgGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCGGGGTGCTGTTCGCCGGc                              <  2:755034/150‑1 (MQ=255)

                                                                                                                        gcgGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCGGGGTGCTGTTCGCCGGCg                             >  2:1620025/1‑150 (MQ=255)

                                                                                                                        gcgGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCGGGGTGCTGTTCGCCGGCg                             >  2:1575084/1‑150 (MQ=255)

                                                                                                                             aCCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCGGGGTGCTGTTCGCCGGCGACCTg                        <  2:1044509/150‑1 (MQ=255)

                                                                                                                               cGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCGGGGTGCTGTTCGCCGGCGACCTGgt                      >  2:2115351/1‑150 (MQ=255)

                                                                                                                               cGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCGGGGTGCTGTTCGCCGGCGACCTGgt                      <  1:1620025/150‑1 (MQ=255)

                                                                                                                               cGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGCGCATGAGGTGGGTTGGCACCGCCTGAGGTTGATGGGGTTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCGGGGTGCTGTTCGCCGGCGACCTGgt                      <  1:758060/150‑1 (MQ=255)

                                                                                                                                      gtgCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCGGGGTGCTGTTCGCCGGCGACCTGGTGTTCTAc               >  1:1192896/1‑150 (MQ=255)

                                                                                                                                         ctgcCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCGGGGTGCTGTTCGCCGGCGACCTGGTGTTCTACCaa            <  1:398038/150‑1 (MQ=255)

                                                                                                                                            ctgACCCGGGTGCGGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCGGGGTGCTGTTCGCCGGCGACCTGGTGTTCTACCAAcgc         <  2:541379/148‑1 (MQ=255)

                                                                                                                                            ccGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCGGGGTGCTGTTCGCCGGCGACCTGGTGTTCTACCAAcgc         >  1:575213/1‑150 (MQ=255)

                                                                                                                                               aCCNGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCGGGGTGCTGTTCGCCGGCGACCTGGTGTTCTACCAACGCGcc      >  2:1508988/1‑150 (MQ=255)

                                                                                                                                                cccTGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGATCTGGCGATTTTCGATGAAAAGACCGGGGTGCTGTTCGCCGGCGACCTGGTGTTCTACCAACGCGccc     <  2:1045071/150‑1 (MQ=255)

                                                                                                                                                 ccGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCGGGGTGCTGTTCGCCGGCGACCTGGTGTTCTACCAACGCGCCCt    <  1:1696654/150‑1 (MQ=255)

                                                                                                                                                  cGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCGGGGTGCTGTTCGCCGGCGACCTGGTGTTCTACCAACGCGCCCTg   >  1:1957072/1‑150 (MQ=255)

                                                                                                                                                   gggTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCGGGGTGCTGTTCGCCGGCGACCTGGTGTTCTACCAACGCGCCCTGa  >  2:2326036/1‑150 (MQ=255)

                                                                                                                                                    ggTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCg                                               <  2:546071/104‑1 (MQ=255)

                                                                                                                                                    ggTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCg                                               >  1:546071/1‑104 (MQ=255)

TGGGCAACCAGGCCTTTGCCGATGTGCCGATCGGTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCGGGGTGCTGTTCGCCGGCGACCTGGTGTTCTACCAACGCGCCCTGA  >  NC_002947/3066217‑3066513

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