BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,263,722	G→C	6.1%	L336V (CTG→GTG)	PP_3734 ←	ABC transporter permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,263,722	0	G	C	6.1%	486.0 / 14.6	147	L336V (CTG→GTG)	PP_3734	ABC transporter permease
Reads supporting (aligned to +/- strand): ref base G (67/71); new base C (3/6); total (70/77)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 151 total aligned reads displayed.

CAACCATGGCCCCTGCGCCAGTCGCCAGGCCTGGGGTTGCGCCAGCGCCAGCAAAGGCAACCGGGCAGCGCGCATTACACTCTCCAAGCCGGCCAGCAATGCACCGAGGACACTCACCAGCACACCCAGCAGCCACCAGCCCGCCGGCAGGCTCAACGTGCCGGCGACCTGGGCGCCATACAGGCCACGCAGGCTGGCGGCCACGTCGGGCAGCAGCCACGCGGCCAATGCGTAACCACTGGCCACCCCGGCCAGCCCACCCAGTGCCGCGAACAAGCCCAGCTCCAGCACCAGGGCG  >  NC_002947/4263573‑4263870

cAACCATGGCCCCTGCGCCAGTCGCCAGGCCTGGGGTTGCGCCAGCGCCAGCAAAGGCAACCGGGCAGCGCGCATTACACTCTCCAAGCCGGCCAGCAATGCACCGAGGACACTCACCAGCACACCCAGCAGCCACCAGCCCGCCGGCAg                                                                                                                                                      >  2:5096900/1‑150 (MQ=255)

 aaCCATGGCCCCTGCGCCAGTCGCCAGGCCTGGGGTTGCGCCAGCGCCAGCAAAGGCAACCGGGCAGCGCGCATTACACTCTCCAAGCCGGCCAGCAATGCACCGAGGACACTCACCAGCACACCCAGCAGCCACCAGCCCGCCGGCAgg                                                                                                                                                     >  1:142981/1‑150 (MQ=255)

 aaCCATGGCCCCTGCGCCAGTCGCCAGGCCTGGGGTTGCGCCAGCGCCAGCAAAGGCAACCGGGCAGCGCGCATTACACTCTCCAAGCCGGCCAGCAATGCACCGAGGACACTCACCAGCACACCCAGCAGCCACCAGCCCGCCGGCAgg                                                                                                                                                     >  1:3499344/1‑150 (MQ=255)

   ccATGGCCCCCGCGCCAGTCGCCAGGCCTGGGGTTGCGCCAGCGCCAGCAAAGGCAACCGGGCCGCGCGCATTACACTCTCCAAGCCGGCCAGCAATGCACCGAGGACACTCACCAGCACACCCAGCAGCCACCAGCCCGCCGGCAGGCt                                                                                                                                                   >  1:3408423/1‑150 (MQ=255)

     aTGGCCCCTGCGCCAGTCGCCAGGCCTGGGGTTGCGCCAGCGCCAGCAAAGGCAACCGGGCAGCGCGCATTACACTCTCCAAGCCGGCCAGCAATGCACCGAGGACACTCACCAGCACACCCAGCAGCCACCAGCCCGCCGGCAGGCTCa                                                                                                                                                 >  1:3894167/1‑150 (MQ=255)

       ggCCCCTGCGCCAGTCGCCAGGCCTGGGGTTGCGCCAGCGCCAGCAAAGGCAACCGGGCAGCGCGCATTACACTCTCCAAGCCGGCCAGCAATGCACCGAGGACACTCACCAGCACACCCAGCAGCCACCAGCCCGCCGGCAGGCTCAAc                                                                                                                                               <  1:2619197/150‑1 (MQ=255)

        gCCCCTGCGCCAGTCGCCAGGCCTGGGGTTGCGCCAGCGCCAGCAAAGGCAACCGGGCAGCGCGCATTACACTCTCCAAGCCGGCCAGCAATGCACCGAGGACACTCACCAGCACACCCAGCAGCCACCAGCCCGCCGGCAGGCTCAACg                                                                                                                                              >  1:219312/1‑150 (MQ=255)

           ccTGCGCCAGTCGCCAGGCCTGGGGTTGCGCCAGCGCCAGCAAAGGCAACCGGGCCGCGCGCATTACACTCTCCAAGCCGGCCAGCAATGCACCGAGGACACTCACCAGCACACCCAGCAGCCACCAGCCCGCCGGCAGGCTCAACGTGc                                                                                                                                           >  2:917769/1‑150 (MQ=255)

                gcCAGTCGCCAGGCCTGGGGTTGCGCCAGCGCCAGCAAAGGCAACCGGGCCGCGCGCATTACACTCTCCAAGCCGGCCAGCAATGCACCGAGGACACTCACCAGCACACCCAGCAGCCACCAGCCCGCCGGCAGGCTCAACGTGCCGGCg                                                                                                                                      >  2:2064060/1‑150 (MQ=255)

                 ccAGTCGCCAGGCCTGGGGTTGCGCCAGCGCCAGCAAAGGCAACCGGGCAGCGCGCATTACACTCTCCAAGCCGGCCAGCAATGCACCGAGGACACTCACCAGCACACCCAGCAGCCACCAGCCCGCCGGCACGTTGAGCCTGCCGGCgg                                                                                                                                     <  2:2828253/150‑2 (MQ=255)

                 ccAGTCGCCAGGCCTGGGGTTGCGCCAGCGCCAGCAAAGGCAACCGGGCAGCGCGCATTACACTCTCCAAGCCGGCCAGCAATGCACCGAGGACACTCACCAGCACACCCAGCAGCCACCAGCCCGCCGGCACGTTGAGCCTGCCGGCgg                                                                                                                                     <  2:2828320/150‑2 (MQ=255)

                  cAGTCGCCAGGCCTGGGGTTGCGCCAGCGCCAGCAAAGGCAACCGGGCAGCGCGCATTACACTCTCCAAGCCGGCCAGCAATGCACCGAGGACACTCACCAGCACACCCAGCAGCCACCAGCCCGCCGGCACGTTGAGCCTGCCGGCggg                                                                                                                                    <  1:1485236/150‑3 (MQ=255)

                     tCGCCAGGCCTGGGGTTGCGCCAGCGCCAGCAAAGGCAACCGGGCAGCGCGCATTACACTCTCCAAGCCGGCCAGCAATGCACCGAGGACACTCACCAGCACACCCAGCAGCCACCAGCCCGCCGGCAGGCTCAACGTGCCGGCGACCTg                                                                                                                                 <  2:3499344/150‑1 (MQ=255)

                       gCCAGGCCTGGGGTTGCGCCAGCGCCAGCAAAGGCAACCGGGCAGCGCGCATTACACTCTCCAAGCCGGCCAGCAATGCACCGAGGACACTCACCAGCACACCCAGCAGCCACCAGCCCGCCGGCACGTTGAGCCTGCCGGCGGGCTggt                                                                                                                               <  1:4493297/150‑2 (MQ=255)

                        ccAGGCCTGGGGTTGCGCCAGCGCCAGCAAAGGCAACCGGGCAGCGCGCATTACACTCTCCAAGCCGGCCAGCAATGCACCGAGGACACTCACCAGCACACCCAGCAGCCACCAGCCCGCCGGCAGGCTCAACGTGCCGGCGACCTGGgc                                                                                                                              <  1:2495923/150‑1 (MQ=255)

                         cAGGCCTGGGGTTGCGCCAGCGCCAGCAAAGGCAACCGGGCAGCGCGCATTACACTCTCCAAGCCGGCCAGCACTGCACCGAGGACACTCACCAGCACACCCAGCAGCCACCAGCCCGCCGGCAGGCTCAACGTGCCGGCGACCTGGgcg                                                                                                                             >  2:3797106/1‑150 (MQ=255)

                         cAGGCCTGGGGTTGCGCCAGCGCCAGCAAAGGCAACCGGGCAGCGCGCATTACACTCTCCAAGCCGGCCAGCAATGCACCGAGGACACTCACCAGCACACCCAGCAGCCACCAGCCCGCCGGCAGGCTCAACGTGCCGGCGACCTGGgcg                                                                                                                             >  2:4095073/1‑150 (MQ=255)

                           ggCCTGGGGTTGCGCCAGCGCCAGCAAAGGCAACCGGGCAGCGCGCATTACACTCTCCAAGCCGGCCAGCAATGCACCGAGGACACTCACCAGCACACCCAGCAGCCACCAGCCCGCCGGCAGGCTCAACGTGCCGGCGACCTGGGCGcc                                                                                                                           <  1:1460723/150‑1 (MQ=255)

                              cTGGGGTTGCGCCAGCGCCAGCAAAGGCAACCGGGCAGCGCGCATTACACTCTCCAAGCCGGCCAGCAATGCACCGAGGACACTCACCAGCACACCCAGCAGCCACCAGCCCGCCGGCAGGCTCAACGTGCCGGCGACCTGGGCGCCATa                                                                                                                        <  1:3790695/150‑1 (MQ=255)

                              cTGGGGTTGCGCCAGCGCCAGCAAAGGCAACCGGGCAGCGCGCATTACACTCTCCAAGCCGGCCAGCAATGCACCGAGGACACTCACCAGCACACCCAGCAGCCACCAGCCCGCCGGCAGGCTCAACGTGCCGGCGACCTGGGCGCCATa                                                                                                                        <  1:4358334/150‑1 (MQ=255)

                              cTGGGGTTGCGCCAGCGCCAGCAAAGGCAACCGGGCAGCGCGCATTACACTCTCCAAGCCGGCCAGCAATGCACCGAGGACACTCACCAGCACACCCAGCAGCCACCAGCCCGCCGGCAGGCTCAACGTGCCGGCGACCTGGGCGCCATa                                                                                                                        <  1:3790667/150‑1 (MQ=255)

                                ggggTTGCGCCAGCGCCAGCAAAGGCAACCGGGCAGCGCGCATTACACTCTCCAAGCCGGCCAGCAATGCACCGAGGACACTCACCAGCACACCCAGCAGCCACCAGCCCGCCGGCAGGCTCAACGTGCCGGCGACCTGGGCGCCATACa                                                                                                                      <  1:3598778/150‑1 (MQ=255)

                                ggggTTGCGCCAGCGCCAGCAAAGGCAACCGGGCAGCGCGCATTACACTCTCCAAGCCGGCCAGCAATGCACCGAGGACACTCACCAGCACACCCAGCAGCCACCAGCCCGCCGGCAGGCTCAACGTGCCGGCGACCTGGGCGCCATACa                                                                                                                      <  1:3666641/150‑1 (MQ=255)

                                 gggTTGCGCCAGCGCCAGCAAAGGCAACCGGGCAGCGCGCATTACACTCTCCAAGCCGGCCAGCAATGCACCGAGGACACTCACCAGCACACCCAGCAGCCACCAGCCCGCCGGCAGGCTCAACGTGCCGGCGACCTGGGCGCCATACAg                                                                                                                     <  2:4592314/150‑1 (MQ=255)

                                  ggTAGCGCCAGCGCCAGCAAAGGCAACCGGGCAGCGCGCATTACACTCTCCAAGCCGGCCAGCAATGCACCGAGGACACTCACCAGCACACCCAGCAGCCACCAGCCCGCCGGCAGGCTCAACGTGCCGGCGACCTGGGCGCCATACAgg                                                                                                                    <  2:3650117/150‑1 (MQ=255)

                                        gccagcgccagcAAAGGCAACCGGGCAGCGCGCATTACACTCTCCAAGCCGGCCAGCAATGCACCGAGGACACTCACCAGCACACCCAGCAGCCACCAGCCCGCCGGCAGGCTCAACGTGCCGGCGACCTGGGCGCCATACAGGCCACGc                                                                                                              >  1:2898954/1‑150 (MQ=255)

                                        gccagcgccagcAAAGGCAACCGGGCAGCGCGCATTACACTCTCCAAGCCGGCCAGCAATGCACCGAGGACACTCACCAGCACACCCAGCAGCCACCAGCCCGCCGGCAGGCTCAACGTGCCGGCGACCTGGGCGCCATACAGGCCACGc                                                                                                              >  1:3895145/1‑150 (MQ=255)

                                         ccagcgccagcAAAGGCAACCGGGCAGCGCGCATTACACTCTCCAAGCCGGCCAGCAATGCACCGAGGACACTCACCAGCACACCCAGCAGCCACCAGCCCGCCGGCAGGCTCAACGTGCCGGCGACCTGGGCGCCATACAGGCCACGCa                                                                                                             <  2:1739711/150‑1 (MQ=255)

                                             cgccagcAAAGGCAACCGGGCAGCGCGCATTACACTCTCCAAGCCGGCCAGCAATGCACCGAGGACACTCACCAGCACACCCAGCAGCCACCAGCCCGCCGGCAGGCTCAACGTGCCGGCGACCTGGGCGCCATACAGGCCACGCAGGCt                                                                                                         <  1:250556/150‑1 (MQ=255)

                                              gccagcAAAGGCAACCGGGCAGCGCGCATTACACTCTCCAAGCCGGCCAGCAATGCACCGAGGACACTCACCAGCACACCCAGCAGCCACCAGCCCGCCGGCAGGCTCAACGTGCCGGCGACCTGGGCGCCATACAGGCCACGCAGGCTg                                                                                                        <  2:1468674/150‑1 (MQ=255)

                                                  gcAAAGGCAACCGGGCCGCGCGCATTACACTCTCCAAGCCGGCCAGCAATGCACCGAGGACACTCACCAGCACACCCAGCAGCCACCAGCCCGCCGGCAGGCTCAACGTGCCGGCGACCTGGGCGCCATACAGGCCACGCAGGCTggcgg                                                                                                    >  1:215562/1‑150 (MQ=255)

                                                         cAACCGGGCAGCGCGCATTACACTCTCCAAGCCGGCCAGCAATGCACCGAGGACACTCACCAGCACACCCAGCAGCCACCAGCCCGCCGGCCGGCTCGACGTGCCGGCGACCTGGGCGCCATCCAGGCCACGCCGGCTGGCGGCCACggc                                                                                             >  2:3700876/1‑148 (MQ=255)

                                                         cAACCGGGCAGCGCGCATTACACTCTCCAAGCCGGCCAGCAATGCACCGAGGACACTCACCAGCACACCCAGCAGCCACCAGCCCGCCGGCAGGCTCAACGTGCCGGCGACCTGGGCGCCATACAGGCCACGCAGGCTGGCGGCCACGTc                                                                                             <  2:3554871/150‑1 (MQ=255)

                                                         cAACCGGGCAGCGCGCATTACACTCTCCAAGCCGGCCAGCAATGCACCGAGGACACTCACCAGCACACCCAGCAGCCACCAGCCCGCCGGCAGGCTCAACGTGCCGGCGACCTGGGCGCCATACAGGCCACGCAGGCTGGCGGCCACGTc                                                                                             <  2:4946969/150‑1 (MQ=255)

                                                         cAACCGGGCAGCGCGCATTACACTCTCCAAGCCGGCCAGCAATGCACCGAGGACACTCACCAGCACACCCAGCAGCCACCAGCCCGCCGGCAGGCTCAACGTGCCGGCGACCTGGGCGCCATACAGGCCACGCAGGCTGGCGGCCACGTc                                                                                             <  2:1729052/150‑1 (MQ=255)

                                                         cAACCGGGCAGCGCGCATTACACTCTCCAAGCCGGCCAGCAATGCACCGAGGACACTCACCAGCACACCCAGCAGCCACCAGCCCGCCGGCAGGCTCAACGTGCCGGCGACCTGGGCGCCATACAGGCCACGCAGGCTGGCGGCCACGTc                                                                                             <  1:1608725/150‑1 (MQ=255)

                                                          aaCCGGGCCGCGCGCATTACACTCTCCAAGCCGGCCAGCAATGCACCGAGGACACTCACCAGCACACCCAGCAGCCACCAGCCCGCCGGCAGGCTCAACGTGCCGGCGACCTGGGCGCCATACAGGCCACGCAGGCTGGCGGCCACGTCg                                                                                            >  1:4602136/1‑150 (MQ=255)

                                                          aaCCGGGCAGCGCGCATTACACTCTCCAAGCCGGCCAGCAATGCACCGAGGACACTCACCAGCACACCCAGCAGCCACCAGCCCGCCGGCAGGCTCAACGTGCCGGCGACCTGGGCGCCATACAGGCCACGCAGGCTGGCGGCCACGTCg                                                                                            <  1:4790126/150‑1 (MQ=255)

                                                            ccGGGCAGCGCGCATTACACTCTCCAAGCCGGCCAGCAATGCACCGAGGACACTCACCAGCACACCCAGCAGCCACCAGCCCGCCGGCAGGCTCAACGTGCCGGCGACCTGGGCGCCATACAGGCCACGCAGGCTGGCGGCCACGTCggg                                                                                          <  2:1498613/150‑1 (MQ=255)

                                                            ccGGGCAGCGCGCATTACACTCTCCAAGCCGGCCAGCAATGCACCGAGGACACTCACCAGCACACCCAGCAGCCACCAGCCCGCCGGCAGGCTCAACGTGCCGGCGACCTGGGCGCCATACAGGCCACGCAGGCTGGCGGCCACGTCggg                                                                                          <  2:1498457/150‑1 (MQ=255)

                                                            ccGGGCAGCGCGCATTACACTCTCCAAGCCGGCCAGCAATGCACCGAGGACACTCACCAGCACACCCAGCAGCCACCAGCCCGCCGGCAGGCTCAACGTGCCGGCGACCTGGGCGCCATACAGGCCACGCAGGCTGGCGGCCACGTCggg                                                                                          <  1:1418992/150‑1 (MQ=255)

                                                                 cAGCGCGCATTACACTCTCCAAGCCGGCCAGCAATGCACCGAGGACACTCACCAGCACACCCAGCAGCCACCAGCCCGCCGGCAGGCTCAACGTGCCGGCGACCTGGGCGCCATACAGGCCACGCAGGCTGGCGGCCACGTCGGgcagca                                                                                     <  2:3048606/150‑1 (MQ=255)

                                                                 cAGCGCGCATTACACTCTCCAAGCCGGCCAGCAATGCACCGAGGACACTCACCAGCACACCCAGCAGCCACCAGCCCGCCGGCAGGCTCAACGTGCCGGCGACCTGGGCGCCATACAGGCCACGCAGGCTGGCGGCCACGTCGGgcagca                                                                                     <  2:3048537/150‑1 (MQ=255)

                                                                 cAGCGCGCATTACACTCTCCAAGCCGGCCAGCAATGCACCGAGGACACTCACCAGCACACCCAGCAGCCACCAGCCCGCCGGCAGGCTCAACGTGCCGGCGACCTGGGCGCCATACAGGCCACGCAGGCTGGCGGCCACGTCGGgcagca                                                                                     <  1:1671059/150‑1 (MQ=255)

                                                                 cAGCGCGCATTACACTCTCCAAGCCGGCCAGCAATGCACCGAGGACACTCACCAGCACACCCAGCAGCCACCAGCCCGCCGGCAGGCTCAACGTGCCGGCGACCTGGGCGCCATACACGCCACGCAGGCTGGCGGCCACGTCGGgcagca                                                                                     <  2:3672105/150‑1 (MQ=255)

                                                                    cgcgcATTACACTCTCCAAGCCGGCCAGCAATGCACCGAGGACACTCACCAGCACACCCAGCAGCCACCAGCCCGCCGGCAGGCTCAACGTGCCGGCGACCTGGGCGCCATACAGGCCACGCAGGCTGGCGGCCACGTCGGGCAGCAGcc                                                                                  >  2:2243485/1‑150 (MQ=255)

                                                                     gcgcATTACACTCTCCAAGCCGGCCAGCAATGCACCGAGGACACTCACCAGCACACCCAGCAGCCACCAGCCCGCCGGCAGGCTCAACGTGCCGGCGACCTGGGCGCCATACAGGCCACGCAGGCTGGCGGCCACGTCGGGCAGCAGCCa                                                                                 <  1:999078/150‑1 (MQ=255)

                                                                       gcATTACACTCTCCAAGCCGGCCAGCAATGCACCGAGGACACTCACCAGCACACCCAGCAGCCACCAGCCCGCCGGCAGGCTCAACGTGCCGGCGACCTGGGCGCCATACAGGCCACGCAGGCTGGCGGCCACGTCGGGCAGCAGCCAcg                                                                               >  2:4712835/1‑150 (MQ=255)

                                                                           tACACTCTCCAAGCCGGCCAGCAATGCACCGAGGACACTCACCAGCACACCCAGCAGCCACCAGCCCGCCGGCAGGCTCAACGTGCCGGCGACCTGGGCGCCATACAGGCCACGCAGGCTGGCGGCCACGTCGGGCAGCAGCCACGCGGc                                                                           <  1:509538/150‑1 (MQ=255)

                                                                            acacTCTCCAAGCCGGCCAGCAATGCACCGAGGACACTCACCAGCACACCCAGCAGCCACCAGCCCGCCGGCAGGCTCAACGTGCCGGCGACCTGGGCGCCATACAGGCCACGCAGGCTGGCGGCCACGTCGGGCAGCAGCCACGCGGcc                                                                          <  1:3152698/150‑1 (MQ=255)

                                                                             cacTCTCCAAGCCGGCCAGCAATGCACCGAGGACACTCACCAGCACACCCAGCAGCCACCAGCCCGCCGGCAGGCTCAACGTGCCGGCGACCTGGGCGCCATACAGGCCACGCAGGCTGGCGGCCACGTCGGGCAGCAGCCACGCGGCCa                                                                         >  1:4501625/1‑150 (MQ=255)

                                                                              acTCTCCAAGCCGGCCAGCAATGCACCGAGGACACTCACCAGCACACCCAGCAGCCACCAGCCCGCCGGCAGGCTCAACGTg                                                                                                                                            <  1:2662146/82‑1 (MQ=255)

                                                                              acTCTCCAAGCCGGCCAGCAATGCACCGAGGACACTCACCAGCACACCCAGCAGCCACCAGCCCGCCGGCAGGCTCAACGTg                                                                                                                                            >  2:2662146/1‑82 (MQ=255)

                                                                                tctcCAAGCCGGCCAGCAATGCACCGAGGACACTCACCAGCACACCCAGCAGCCACCAGCCCGCCGGCAGGCTCAACGTGACGGCGACCTGGGCGCCATACAGGCCACGCAGGCTGGCGGCCACGTCGGGCAGCAGCCACGCGGCCAATg                                                                      >  2:4827799/1‑150 (MQ=255)

                                                                                 ctcCACGCCGGCCAGCAATGCACCGAGGACACTCACCAGCACACCCAGCAGCCACCAGCCCGCCGGCAGGCTCAACGTGCCGGCGACCTGGGCGCCATACAGGCCACGCAGGCTGGCGGCCACGTCGGGCAGCAGCCACGCGGCCAATGc                                                                     <  2:1334391/150‑1 (MQ=255)

                                                                                  tcCAAGCCGGCCAGCAATGCACCGAGGACACTCACCAGCACACCCAGCAGCCACCAGCCCGCCGGCAGGCTCAACGTGCCGGCGACCTGGGCGCCATACAGGCCACGCAGGCTGGCGGCCACGTCGGGCAGCAGCCACGCGGCCAATGCg                                                                    >  1:185946/1‑150 (MQ=255)

                                                                                  tcCAAGCCGGCCAGCAATGCACCGAGGACACTCACCAGCACACCCAGCAGCCACCAGCCCGCCGGCAGGCTCAACGTGCCGGCGACCTGGGCGCCATACAGGCCACGCAGGCTGGCGGCCACGTCGGGCAGCAGCCACGCGGCCAATGCg                                                                    >  1:3614511/1‑150 (MQ=255)

                                                                                  tcCAAGCCGGCCAGCAATGCACCGAGGACACTCACCAGCACACCCAGCAGCCACCAGCCCGCCGGCAGGCTCAACGTGCCGGCGACCTGGGCGCCATACAGGCCACGCAGGCTGGCGGCCACGTCGGGCAGCAGCCACGCGGCCAATGCg                                                                    >  1:2612562/1‑150 (MQ=255)

                                                                                       gCCGGCCAGCAATGCACCGAGGACACTCACCAGCACACCCAGCAGCCACCAGCCCGCCGGCAGGCTCAACGTGCCGGCGACCTGGGCGCCATACAGGCCACGCAGGCTGGCGGCCACGTCGGGCAGCAGCCACGCGGCCAATGCGTAAcc                                                               >  1:558017/1‑150 (MQ=255)

                                                                                         cGGCCAGCAATGCACCGAGGACACTCACCAGCACACCCAGCAGCCACCAGCCCGCCGGCAGGCTCAACGTGCCGGCGACCTGGGCGCCATACAGGCCACGCAGGCTGGCGGCCACGTCGGGCAGCAGCCACGCGGCCAATGCGTAACCAc                                                             <  1:1074310/150‑1 (MQ=255)

                                                                                          ggCCAGCAATGCACCGAGGACACTCACCAGCACACCCAGCAGCCACCAGCCCGCCGGCAGGCTCAACGTGCCGGCGACCTGGGCGCCATACAGGCCACGCAGGCTGGCGGCCACGTCGGGCAGCAGCCACGCGGCCAATGCGTAACCACt                                                            >  1:2422792/1‑150 (MQ=255)

                                                                                           gCCAGCAATGCACCGAGGACACTCACCAGCACACCCAGCAGCCACCAGCCCGCCGGCAGGCTCAACGTGCCGGCGACCTGGGCGCCATACAGGCCACGCAGGCTGGCGGCCACGTCGGGCAGCAGCCACGCGGCCAATGCGTAACCACTg                                                           <  2:1453829/150‑1 (MQ=255)

                                                                                            ccAGCAATGCACCGAGGACACTCACCAGCACACCCAGCAGCCACCAGCCCGCCGGCAGGCTCAACGTGCCGGCGACCTGGGCGCCATACAGGCCACGCAGGCTGGCGGCCACGTCGGGCAGCAGCCACGCGGCCAATGCGTAACCACTgg                                                          <  2:3919459/150‑1 (MQ=255)

                                                                                             cAGCAATGCACCGAGGACACTCACCAGCACACCCAGCAGCCACCAGCCCGCCGGCAGGCTCAACGTGCCGGCGACCTGGGCGCCATACAGGCCACGCAGGCTGGCGGCCACGTCGGGCAGCAGCCACGCGGCCAATGCGTAACCACTGGc                                                         >  2:1382191/1‑150 (MQ=255)

                                                                                              aGCAATGCACCGAGGACACTCACCAGCACACCCAGCAGCCACCAGCCCGCCGGCAGGCTCAACGTGCCGGCGACCTGGGCGCCATACAGGCCACGCAGGCTGGCGGCCACGTCGGGCAGCAGCCACGCGGCCAATGCGTAACCACTGGcc                                                        >  1:2822311/1‑150 (MQ=255)

                                                                                              aGCAATGCACCGAGGACACTCACCAGCACACCCAGCAGCCACCAGCCCGCCGGCAGGCTCAACGTGCCGGCGACCTGGGCGCCATACAGGCCACGCAGGCTGGCGGCCACGGCGGGCAGCAGCCACGCGGCCAATGCGTAACCACTGGcc                                                        >  2:1110610/1‑150 (MQ=255)

                                                                                                  aTGCACCGAGGACACTCACCAGCACACCCAGCAGCCACCAGCCCGCCGGCAGGCTCAACGTGCCGGCGACCTGGGCGCCATACAGGCCACGCAGGCTGGCGGCCACGTCGGGCAGCAGCCACGCGGCCAATGCGTAACCACTGGCCAccc                                                    <  1:1641474/150‑1 (MQ=255)

                                                                                                      aCCGAGGACACTCACCAGCACACCCAGCAGCCACCAGCCCGCCGGCAGGCTCAACGTGCCGGCGACCTGGGCGCCATACAGGCCACGCAg                                                                                                            >  1:3286625/1‑90 (MQ=255)

                                                                                                      aCCGAGGACACTCACCAGCACACCCAGCAGCCACCAGCCCGCCGGCAGGCTCAACGTGCCGGCGACCTGGGCGCCATACAGGCCACGCAg                                                                                                            <  2:3286625/90‑1 (MQ=255)

                                                                                                      aCCGAGGACACTCACCAGCACACCCAGCAGCCACCAGCCCGCCGGCAGGCTCAACGTGCCGGCGACCTGGGCGCCATACAGGCCACGCAGGCTGGCGGCCACGTCGGGCAGCAGCCACGCGGCCAATGCGTAACCACTGGCCACCCCGGc                                                <  2:4281236/150‑1 (MQ=255)

                                                                                                      aCCGAGGACACTCACCAGCACACCCAGCAGCCACCAGCCCGCCGGCAGGCTCAACGTGCCGGCGACCTGGGCGCCATACAGGCCACGCAGGCTGGCGGCCACGTCGGGCAGCAGCCACGCGGCCAATGCGTAACCACTGGCCACCCCGGc                                                <  2:3090034/150‑1 (MQ=255)

                                                                                                      aCCGAGGACACTCACCAGCACACCCAGCAGCCACCAGCCCGCCGGCAGGCTCAACGTGCCGGCGACCTGGGCGCCATACAGGCCACGCAGGCTGGCGGCCACGTCGGGCAGCAGCCACGCGGCCAATGCGTAACCACTGGCCACCCCGGc                                                >  1:4930688/1‑150 (MQ=255)

                                                                                                        cGAGGACACTCACCAGCACACCCAGCAGCCACCAGCCCGCCGGCAGGCTCAACGTGCCGGCGACCTGGGCGCCATACAGGCCACGCAGGCTGGCGGCCACGTCGGGCAGCAGCCACGCGGCCAATGCGTAACCACTGGCCACCCCGGCCa                                              <  2:407184/150‑1 (MQ=255)

                                                                                                           ggACACTCACCAGCACACCCAGCAGCCACCAGCCCGCCGGCAGGCTCAACGTGCCGGCGACCTGGGCGCCATACAGGCCACGCAGGCTGGCGGCCACGTCGGGCAGCAGCCACGCGGCCAATGCGTAACCACTGGCCACCCCGGCCAGcc                                           >  1:125173/1‑150 (MQ=255)

                                                                                                               acTCACCAGCACACCCAGCAGCCACCAGCCCGCCGGCAGGCTCAACGTGCCGGCGACCTGGGCGCCATACAGGc                                                                                                                   <  1:3681556/74‑1 (MQ=255)

                                                                                                                 tCACCAGCACACCCAGCAGCCACCAGCCCGCCGGCAGGCTCAACGTGCCGGCGACCTGGGCGCCATACAGGCCACGCAGGCTGGCGGCCACGTCGGGCAGCAGCCACGCGGCCAATGCGTAACCACTGGCCACCCCGGCCAGcccaccca                                     >  1:1888688/1‑150 (MQ=255)

                                                                                                                   aCCAGCACACCCAGCAGCCACCAGCCCGCCGGCAGGCTCAACGTGCCGGCGACCTGGGCGCCATACAGGCCACGCAGGCTGGCGGCCACGTCGGGCAGCAGCCACGCGGCCAATGCGTAACCACTGGCCACCCCGGCCAGCCCACCCAGt                                   >  2:5022146/1‑150 (MQ=255)

                                                                                                                    ccAGCACACCCAGCAGCCACCAGCCCGCCGGCAGGCTCAACGTGCCGGCGACCTGGGCGCCATACAGGCCACGCAGGCTGGCGGCCACGTCGGGCAGCAGCCACGCGGCCAATGCGTAACCACTGGCCACCCCGGCCAGCCCACCCAGTg                                  <  2:2724676/150‑1 (MQ=255)

                                                                                                                       gCACACCCAGCAGCCACCAGCCCGCCGGCAGGCTCAACGTGCCGGCGACCTGGGCGCCATACAGGCCACGCAGGCTGGCGGCCACGTCGGGCAGCAGCCACGCGGCCAATGCGTAACCACTGGCCACCCCGGCCAGCCCACCCAGTGCcg                               >  2:547212/1‑150 (MQ=255)

                                                                                                                        cacaCCCAGCAGCCACCAGCCCGCCGGCAGGCTCAACGTGCCGGCGACCTGGGCGCCATACAGGCCACGCAGGCTGGCGGCCACGTCGGGCAGCAGCCACGCGGCCAATGCGTAACCACTGGCCACCCCGGCCAGCCCACCCAGTGCcgc                              <  2:4743219/150‑1 (MQ=255)

                                                                                                                        cacaCCCAGCAGCCACCAGCCCGCCGGCAGGCTCAACGTGCCGGCGACCTGGGCGCCATACAGGCCACGCAGGCTGGCGGCCACGTCGGGCAGCAGCCACGCGGCCAATGCGTAACCACTGGCCACCCCGGCCAGCCCACCCAGTGCcgc                              <  2:1302537/150‑1 (MQ=255)

                                                                                                                        cacaCCCAGCAGCCACCAGCCCGCCGGCAGGCTCAACGTGCCGGCGACCTGGGCGCCATACAGGCCACGCAGGCTGGCGGCCACGTCGGGCAGCAGCCACGCGGCCAATGCGTAACCACTGGCCACCCCGGCCAGCCCACCCAGTGCcgc                              <  2:334869/150‑1 (MQ=255)

                                                                                                                           acCCAGCAGCCACCAGCCCGCCGGCAGGCTCAACGTGCCGGCGACCTGGGCGCCATACAGGCCACGCAGGCTGGCGGCCACGTCGGGCAGCAGCCACGCGGCCAATGCGTAACCACTGGCCACCCCGGCCAGCCCACCCAGTGCCGCGaa                           >  2:4185257/1‑150 (MQ=255)

                                                                                                                              cagcagCCACCAGCCCGCCGGCAGGCTCAACGTGCCGGCGACCTGGGCGCCATACAGGCCACGCAGGCTGGCGGCCACGTCGGGCAGCAGCCACGCGGCCAATGCGTAACCACTGGCCACCCCGGCCAGCCCACCCAGTGCCGCGAACaa                        >  2:3337646/1‑150 (MQ=255)

                                                                                                                                gcagcCACCAGCCCGCCGGCAGGCTCAACGTGCCGGCGACCTGGGCGCCATACAGGCCACGCAGGCTGGCGGCCACGTCGGGCAGCAGCCACGCGGCCAATGCGTAACCACTGGCCACCCCGGCCAGCCCACCCAGTGCCGCGAACAAGc                      <  1:3364569/150‑1 (MQ=255)

                                                                                                                                 cagcCACCAGCCCGCCGGCAGGCTCAACGTGCCGGCGACCTGGGCGCCATACAGGCCACGCAGGCTGGCGGCCACGTCGGGCAGCAGCCACGCGGCCAATGCGTAACCACTGGCCACCCCGGCCAGCCCACCCAGTGCCGCGAACAAGcc                     <  2:2255375/150‑1 (MQ=255)

                                                                                                                                 cagcCACCAGCCCGCCGGCAGGCTCAACGTGCCGGCGACCTGGGCGCCATACAGGCCACGCAGGCTGGCGGCCACGTCGGGCAGCAGCCACGCGGCCAATGCGTAACCACTGGCCACCCCGGCCAGCCCACCCAGTGCCGCGAACAAGcc                     <  1:412657/150‑1 (MQ=255)

                                                                                                                                  agcCACCAGCCCGCCGGCAGGCTCAACGTGCCGGCGACCTGGGCGCCATACAGGCCACGCAGGCTGGCGGCCACGTCGGGCAGCAGCCACGCGGCCAATGCGTAACCACTGGCCACCCCGGCCAGCCCACTCAGTGCCGCGAACAAGccc                    >  1:2412176/1‑150 (MQ=255)

                                                                                                                                      accaGCCCGCCGGCAGGCTCAACGTGCCGGCGACCTGGGCGCCATACAGGCCACGCAGGCTGGCGGCCACGTCGGGCAGCAGCCACGCGGCCAATGCGTAACCACTGGCCACCCCGGCCAGCCCACCCAGTGCCGCGAACAAGCCCAGCt                <  2:5051112/150‑1 (MQ=255)

                                                                                                                                       ccaggtcgccgGCACGTTGAGCCTGCCGGCGGGCTGGGCGCCATACAGGCCACGCAGGCTGGCGGCCACGTCGGGCAGCAGCCACGCGGCCAATGCGTAACCACTGGCCACCCCGGCCAGCCCACCCAGTGCCGCGAACAAGccc                    <  1:4938065/139‑1 (MQ=38)

                                                                                                                                       ccaGCCCGCCGGCAGGCTCAACGTGCCGGCGACCTGGGCGCCATACAGGCCACGCAGGCTGGCGGCCACGTCGGGCAGCAGCCACGCGGCCAATGCGTAACCACTGGCCACCCCGGCCAGCCCACCCAGTGCCGCGAACAAGCCCAGCTc               <  2:1999991/150‑1 (MQ=255)

                                                                                                                                         aGCCCGCCGGCAGGCTCAACGTGCCGGCGACCTGGGCGCCATACAGGCCACGCAGGCTGGCGGCCACGTCGGGCAGCAGCCACGCGGCCAATGCGTAACCACTGGCCACCCCGGCCAGCCCACCCAGTGCCGCGAACAAGCCCAGCTCCa             <  2:4930688/150‑1 (MQ=255)

                                                                                                                                         aGCCCGCCGGCAGGCTCAACGTGCCGGCGACCTGGGCGCCATACAGGCCACGCAGGCTGGCGGCCACGTCGGGCAGCAGCCACGCGGCCAATGCGTAACCACTGGCCACCCCGGCCAGCCCACCCAGTGCCGCGAACAAGCCCAGCTCCa             <  2:818963/150‑1 (MQ=255)

                                                                                                                                         aGCCCGCCGGCAGGCTCAACGTGCCGGCGACCTGGGCGCCATACAGGCCACGCAGGCTGGCGGCCACGTCGGGCAGCAGCCACGCGGCCAATGCGTAACCACTGGCCACCCCGGCCAGCCCACCCAGTGCCGCGAACAAGCCCAGCTCCa             <  2:818822/150‑1 (MQ=255)

                                                                                                                                          ggccgccgGCAGGCTCAACGTGCCGGCGACCTGGGCGCCATACAGGCCACGCAGGCTGGCGGCCACGTCGGGCGGCAGCCACGCGGCCAATGCGTAACCACTGGCCACCCCGGCCAGCCCACCCAGTGCCGCGAACAAGCCCAGCTCCAg            >  2:5071676/3‑150 (MQ=255)

                                                                                                                                          gCCCGCCGGCAGGCTCAACGTGCCGGCGACCTGGGCGCCATACAGGCCACGCAGGCTGGCGGCCACGTCGGGCAGCAGCCACGCGGCCAATGCGTAACCACTGGCCACCCCGGCCAGCCCACCCAGTGCCGCGAACAAGCCCAGCTCCAg            >  2:4749283/1‑150 (MQ=255)

                                                                                                                                            tcgccgGCACGTTGAGCCTGCCGGCGACCTGGGCGCCATACAGGCCACGCAGGCTGGCGGCCACGTCGGGCAGCAGCCACGCGGCCAATGCGTAACCACTGGCCACCCCGGCCAGCCCACCCAGTGCCGCGAACAAGCCCAGCTCCAGCa          >  2:3905373/2‑150 (MQ=255)

                                                                                                                                                cgGCAGGCTCAACGTGCCGGCGACCTGGGCGCCATACAGGCCACGCAGGCTGGCGGCCACGTCGGGCAGCAGCCACGCGGCCAATGCGTAACCACTGGCCACCCCGGCCAGCCCACCCAGTGCCGCGAACAAGCCCAGCTCCAGCACCAg      >  2:2541178/1‑150 (MQ=255)

                                                                                                                                                  gCAGGCTCAACGTGCCGGCGACCTGGGCGCCATACAGGCCACGCAGGCTGGCGGCCACGTCGGGCAGCAGCCACGCGGCCAATGCGTAACCACTGGCCACCCCGGCCAGCCCACCCAGTGCCGCGAACAAGCCCAGCTCCAGCACCAggg    >  1:1773983/1‑150 (MQ=255)

                                                                                                                                                   cAGGCTCAACGTGCCGGCGACCTGGGCGCCATACAGGCCACGCAGGCTGGCGGCCACGTCGGGCAGCAGCCACGCGGCCAATGCGTAACCACTGGCCACCCCGGCCAGCCCACCCAGTGCCGCGAACAAGCCCAGCTCCAGCACCAGGgc   <  2:1955147/150‑1 (MQ=255)

CAACCATGGCCCCTGCGCCAGTCGCCAGGCCTGGGGTTGCGCCAGCGCCAGCAAAGGCAACCGGGCAGCGCGCATTACACTCTCCAAGCCGGCCAGCAATGCACCGAGGACACTCACCAGCACACCCAGCAGCCACCAGCCCGCCGGCAGGCTCAACGTGCCGGCGACCTGGGCGCCATACAGGCCACGCAGGCTGGCGGCCACGTCGGGCAGCAGCCACGCGGCCAATGCGTAACCACTGGCCACCCCGGCCAGCCCACCCAGTGCCGCGAACAAGCCCAGCTCCAGCACCAGGGCG  >  NC_002947/4263573‑4263870

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