BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	311,727	G→C	11.5%	D217E (GAC→GAG)	PP_0256 ←	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	311,727	0	G	C	11.5%	560.6 / 51.8	183	D217E (GAC→GAG)	PP_0256	putative oxidoreductase
Reads supporting (aligned to +/- strand): ref base G (73/89); new base C (6/15); total (79/104)
Fisher's exact test for biased strand distribution p-value = 1.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

Only 100 of 194 total aligned reads displayed.

AAGCGCTCCAGACCACGCTCCTTGAGCGGGTTGATGCACACCACCTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCA‑GACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGC  >  NC_002947/311578‑311876

acccgcTCCAGACCACGCTCCTTGAGCGGGTTGATGCACACCACCTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCA‑GACAAAAATGGCGTCGGCGTGCTCGAAgtcg                                                                                                                                                        >  2:3163121/4‑150 (MQ=255)

 aGCGCTCCAGACCACGCTCCTTGAGCGGGTTGATGCACACCACCTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCA‑GACAAAAATGGCGTCGGCGTGCTCGAAgtcgt                                                                                                                                                       <  2:2757710/150‑1 (MQ=255)

      tCCAGACCACGCTCCTTGAGCGGGTTGATGCACACCACCTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCCGGCCCCA‑GACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAg                                                                                                                                                  >  1:3147046/1‑150 (MQ=255)

           aCCACGCTCCTTGAGCGGGTTGATGCACACCACCTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCA‑GACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑Ac                                                                                                                                            >  1:3603045/1‑150 (MQ=255)

             cACGCTCCTTGAGCGGGTTGATGCACACCACCTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCA‑GACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTg                                                                                                                                          >  1:1208242/1‑150 (MQ=255)

              aCGCTCCTTGAGCGGGTTGATGCACACCACCTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCA‑GACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGt                                                                                                                                         >  2:4568972/1‑150 (MQ=255)

              aCGCTCCTTGAGCGGGTTGATGCACACCACCTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCA‑GACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGt                                                                                                                                         <  1:651250/150‑1 (MQ=255)

               cGCTCCTTGAGCGGGTTGATGCACACCACCTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCA‑GACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTa                                                                                                                                        >  1:1338981/1‑150 (MQ=255)

                       gAGCGGGTTGATGCACACCACCTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCA‑GACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGcc                                                                                                                                <  2:4894692/150‑1 (MQ=255)

                          cGGGTTGATGCACACCACCTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCA‑GACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGAc                                                                                                                             <  1:5400436/150‑1 (MQ=255)

                          cGGGTTGATGCACACCACCTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCA‑GACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGAc                                                                                                                             <  1:2887925/150‑1 (MQ=255)

                            ggTTGATGCACACCACCTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCA‑GACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGAccc                                                                                                                           >  2:2399458/1‑150 (MQ=255)

                                 aTGCACACCACCTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCA‑GACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACg                                                                                                                      >  1:4252207/1‑150 (MQ=255)

                                 aTGCACACCACCTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCA‑GACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACg                                                                                                                      >  2:1142252/1‑150 (MQ=255)

                                  tGCACACCACCTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCA‑GACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGc                                                                                                                     >  2:537889/1‑150 (MQ=255)

                                    cacaCCACCTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCA‑GACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCtt                                                                                                                   >  2:5285238/1‑150 (MQ=255)

                                    cacaCCACCTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCA‑GACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCtt                                                                                                                   <  1:2618999/150‑1 (MQ=255)

                                     gccccaccTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCA‑GACCAAAATGGCGTCGGCGTGCTCGAAGTCGTCGACGGTG‑ACTGTACCTTTGCCGACCCCGACGCttt                                                                                                                  >  1:3666457/4‑150 (MQ=255)

                                     acacCACCTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCA‑GACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCttt                                                                                                                  >  1:1092420/1‑150 (MQ=255)

                                     acacCACCTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCA‑GACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCttt                                                                                                                  >  1:1042701/1‑150 (MQ=255)

                                       accaccTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCA‑GACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGc                                                                                                                <  2:2306805/150‑1 (MQ=255)

                                       accaccTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCA‑GACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGc                                                                                                                <  2:2306839/150‑1 (MQ=255)

                                        ccaccTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCA‑GACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGcc                                                                                                               <  1:591406/150‑1 (MQ=255)

                                            cTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCA‑GACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCaa                                                                                                           >  2:2033693/1‑150 (MQ=255)

                                            cTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCA‑GACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCaa                                                                                                           <  1:903978/150‑1 (MQ=255)

                                             tGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCA‑GACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAac                                                                                                          >  2:694469/1‑150 (MQ=255)

                                             tGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCA‑GACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAac                                                                                                          <  2:3829754/150‑1 (MQ=255)

                                                  gcCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCA‑GACCAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCCacgccac                                                                                                     >  2:4342832/1‑150 (MQ=255)

                                                  gcCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCA‑GACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAacgccac                                                                                                     >  2:1879979/1‑150 (MQ=255)

                                                    cACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCA‑GACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAacgccacgc                                                                                                   >  2:5347513/1‑150 (MQ=255)

                                                     aCGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCA‑GACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAacgccacgcc                                                                                                  >  1:3475558/1‑150 (MQ=255)

                                                         tttACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCA‑GACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTg                                                                                              <  2:1284917/150‑1 (MQ=255)

                                                         tttACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCA‑GACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTg                                                                                              <  1:5181672/150‑1 (MQ=255)

                                                         tttACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCA‑GACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTg                                                                                              >  1:420534/1‑150 (MQ=255)

                                                         tttACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCA‑GACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTg                                                                                              >  1:4862958/1‑150 (MQ=255)

                                                           tACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCA‑GACAAAAATGGCGTCGGCGTGCTCGAAGTAGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGc                                                                                            <  2:3829389/150‑1 (MQ=255)

                                                               gCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCA‑GACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCa                                                                                        >  2:2858569/1‑150 (MQ=255)

                                                                         gTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCA‑GACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGt                                                                              >  1:412278/1‑150 (MQ=255)

                                                                           gggTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCA‑GACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTc                                                                            <  2:860692/150‑1 (MQ=255)

                                                                             gTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCA‑GACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGa                                                                          <  1:1185687/150‑1 (MQ=255)

                                                                             gTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCA‑GACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGa                                                                          <  2:2790975/150‑1 (MQ=255)

                                                                              cccaacATGCGGGGATGATTGGTGCCAGGATTTTGCCCCAC‑GACGAAAATCGCGTCAGCCTGCTCGAAATCCTC‑AAGCAGAACAGTTCCCTTGCCAACCCCTACCGATTGTTTGAGCGCCACACCGCTGGCCTCATGGCACATGTTCGaa                                                                         <  2:4885896/145‑2 (MQ=255)

                                                                                gAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCA‑GACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCa                                                                       >  1:4515812/1‑150 (MQ=255)

                                                                                 aacATGCGGGGATGATTGGTGCCAGGATTTTGCCCCAC‑GACGAAAATCGCGTCAGCCTGCTCGAAATCCTC‑AAGCAGAACAGTTCCCTTGCCAACCCCTACCGATTGTTTGAGCGCCACACCGCTGGCCTCATGGCACATGTTCGAACAg                                                                      >  2:3384367/3‑150 (MQ=255)

                                                                                 aGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCA‑GACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGcc                                                                                                               <  2:3589351/109‑1 (MQ=255)

                                                                                 aGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCA‑GACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGcc                                                                                                               >  1:3589351/1‑109 (MQ=255)

                                                                                  gCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCA‑GACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGt                                                                     <  1:5361846/150‑1 (MQ=255)

                                                                                  acATGCGGGGATGATTGGTGCCAGGATTTTGCCCCAC‑GACGAAAATCGCGTCAGCCTGCTCGAAATCCTC‑AAGCAGAACAGGTCCCTTGCCAACCCCTACCGATTGTTTGAGCGCCACACCGCTGGCCTCATGGCACATGTTCGAACAGt                                                                     <  2:4781816/149‑1 (MQ=255)

                                                                                   cATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCA‑GACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTc                                                                    <  2:361604/150‑1 (MQ=255)

                                                                                    aTGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCA‑GACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCa                                                                   <  2:1680453/150‑1 (MQ=255)

                                                                                     tGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCA‑GACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAg                                                                  <  1:2803231/150‑1 (MQ=255)

                                                                                        gcgGGTGATTGGTGCCGGGGTTCTGGCCCCA‑GACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGAc                                                                                                                             <  1:2809268/88‑1 (MQ=255)

                                                                                        gcgGGTGATTGGTGCCGGGGTTCTGGCCCCA‑GACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGa                                                               >  1:5320372/1‑150 (MQ=255)

                                                                                         cgGGTGATTGGTGCCGGGGTTCTGGCCCCA‑GACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGaa                                                              <  2:1102996/150‑1 (MQ=255)

                                                                                          ggatGATTGGTGCCAGGATTTTGCCCCAC‑GACGAAAATCGCGTCAGCCTGCTCGAAATCCTC‑AAGCAGAACAGTTCCCTTGCCAACCCCTACCGATTGTTTGAGCGCCACACCGCTGGCCTCATGGCACATGTTCGAACAGTCAGGGaaa                                                             <  1:2939829/147‑2 (MQ=255)

                                                                                           ggTGATTGGTGCCGGGGTTCTGGCCCCA‑GACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATgg                                                                                     <  1:2296565/125‑1 (MQ=255)

                                                                                           ggTGATTGGTGCCGGGGTTCTGGCCCCA‑GACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAgt                                                            >  2:4015479/1‑150 (MQ=255)

                                                                                              gATTGGTGCCGGGGTTCTGGCCCCA‑GACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAgttgt                                                         <  2:1913405/150‑1 (MQ=255)

                                                                                               aTTGGTGCCAGGATTTTGCCCCAC‑GACGAAAATCGCGTCAGCCTGCTCGAAATCCTC‑AAGCAGAACAGTTCCCTTGCCAACCCCTACCGATTGTTTGAGCGCCACACCGCTGGCCTCATGGCACATGTTCGAACAGTCAGGGAAattgtt                                                        <  2:2295193/150‑1 (MQ=255)

                                                                                                ttGGTGCCGGGGTTCTGGCCCCA‑GACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGttgttg                                                       <  1:4602288/150‑1 (MQ=255)

                                                                                                ttGGTGCCAGGATTTTGCCCCAC‑GACGAAAATCGCGTCAGCCTGCTCGAAATCCTC‑AAGCAGAACAGTTCCCTTGCCAACCCCTACCGATTGTTTGAGCGCCACACCGCTGGCCTCATGGCACATGTTCGAACAGTCAGGGAAattgttc                                                       <  1:2149738/150‑2 (MQ=255)

                                                                                                 tGGTGCCGGGGTTCTGGCCCCA‑GACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTc                                                                            >  1:187281/1‑128 (MQ=255)

                                                                                                 tGGTGCCGGGGTTCTGGCCCCA‑GACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTc                                                                            <  2:187281/128‑1 (MQ=255)

                                                                                                 tGGTGCCAGGATTTTGCCCCAC‑GACGAAAATCGCGTCAGCCTGCTCGAAATCCTC‑AAGCAGAACAGTTCCCTTGCCAACCCCTACCGATTGTTTGAGCGCCACACCGCTGGCCTCATGGCACATGTTCGAACAGTCAGGGAAattgttcg                                                      >  2:1251424/1‑148 (MQ=255)

                                                                                                 tGGTGCCAGGATTTTGCCCCAC‑GACGAAAATCGCGTCAGCCTGCTCGAAATCCTC‑AAGCAGAACAGTTCCCTTGCCAACCCCTACCGATTGTTTGAGCGCCACACCGCTGGCCTCATGGCACATGTTCGAACAGTCAGGGAAattgttcg                                                      <  1:1847329/150‑3 (MQ=255)

                                                                                                  ggTGCCAGGATTTTGCCCCAC‑GACGAAAATCGCGTCAGCCTGCTCGAAATCCTC‑AAGCAGAACAGTTCCCTTGCCAACCCCTACCGATTGTTTGAGCGCCACACCGCTGGCCTCATGGCACATGTTCGAACAGTCAGGGAAATTGTTCGt                                                     <  1:3384367/150‑1 (MQ=255)

                                                                                                    tGCCGGGGTTCTGGCCCCA‑GACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGc                                                   <  2:2686471/150‑1 (MQ=255)

                                                                                                     gCCGGGGTTCTGGCCCCA‑GACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGcc                                                  <  2:5064580/150‑1 (MQ=255)

                                                                                                     gCCGGGGTTCCGGCCCCA‑GACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGcc                                                  <  2:3147046/150‑1 (MQ=255)

                                                                                                       cGGGGTTCTGGCCCCA‑GACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGt                                                <  1:2676809/150‑1 (MQ=255)

                                                                                                         gggTTCTGGCCCCA‑GACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAg                                              >  1:1333434/1‑150 (MQ=255)

                                                                                                          ggTTCTGGCCCCA‑GACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAgg                                             >  2:3872207/1‑150 (MQ=255)

                                                                                                              tttgCCCCACGACGAAAATCGCGTCAGCCTGCTCGAAATCCTC‑AAGCAGAACAGTTCCCTTGCCAACCCCTACCGATTGTTTGAGCGCCACACCGCTGGCCTCATGGCACATGTTCGAACAGTCAGGGAAATTGTTCGTGCCATAAgcgc                                          <  1:402085/147‑1 (MQ=255)

                                                                                                               tGGCCCCA‑GACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGTCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGgcgcgc                                        >  2:3096652/1‑150 (MQ=255)

                                                                                                                ggCCCCA‑GACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCa                                       >  1:3503488/1‑150 (MQ=255)

                                                                                                                 gCCCCA‑GACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCAc                                      <  1:694469/150‑1 (MQ=255)

                                                                                                                  ccccA‑GACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCAcgccgc                                                                                                <  2:3796016/91‑1 (MQ=255)

                                                                                                                  ccccA‑GACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCAcgccgc                                                                                                >  1:3796016/1‑91 (MQ=255)

                                                                                                                   cccA‑GACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCCAAGGTG‑ACTGTACATTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGACTCATGGCACATGATCGAGCAGTCAGGGAAGATGATGGAGCCGTAGGCGCGCACGa                                    >  1:436899/1‑150 (MQ=255)

                                                                                                                    ccacgACGAAAATCGCGTCAGCCTGCTCGAAATCCTC‑AAGCAGAACAGTTCCCTTGCCAACCCCTACCGATTGTTTGAGCGCCACACCGCTGGCCTCATGGCACATGTTCGAACAGTCAGGGAAATTGTTCGTGCCATAAGCGCGCACGa                                    >  1:4986358/5‑150 (MQ=255)

                                                                                                                     cA‑GACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAAc                                  <  2:1494475/150‑1 (MQ=255)

                                                                                                                         aCAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGt                               >  1:2511125/1‑150 (MQ=255)

                                                                                                                          cAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGtt                              >  2:1315288/1‑150 (MQ=255)

                                                                                                                          cAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGtt                              <  2:2308561/150‑1 (MQ=255)

                                                                                                                            aaaaTGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTgg                            <  2:4394451/150‑1 (MQ=255)

                                                                                                                             aaaTGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGt                           <  2:523996/150‑1 (MQ=255)

                                                                                                                              aaTGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTa                          <  1:3178297/150‑1 (MQ=255)

                                                                                                                                tGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACa                        <  2:2120215/150‑1 (MQ=255)

                                                                                                                                 ggCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAg                       <  1:2721429/150‑1 (MQ=255)

                                                                                                                                        gCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTAcgcc                <  1:3731415/150‑1 (MQ=255)

                                                                                                                                         cGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTAcgccg               >  1:3212009/1‑150 (MQ=255)

                                                                                                                                           tGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCt             <  1:2033693/150‑1 (MQ=255)

                                                                                                                                             cTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTc           <  2:3456951/150‑1 (MQ=255)

                                                                                                                                              tCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCg          >  2:1109451/1‑150 (MQ=255)

                                                                                                                                              tCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCg          <  2:3241060/150‑1 (MQ=255)

                                                                                                                                                  aGTCGTCGTAGGTG‑TCTTTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCGGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGc      <  1:1538650/150‑1 (MQ=255)

                                                                                                                                                   gtcgtcGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCt     <  1:4420022/150‑1 (MQ=255)

                                                                                                                                                    tcgtcgAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTg    >  2:1045586/1‑150 (MQ=255)

                                                                                                                                                    tcgtcgAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTg    <  1:3271511/150‑1 (MQ=255)

                                                                                                                                                      gtcgAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGc  >  1:4744636/1‑150 (MQ=255)

AAGCGCTCCAGACCACGCTCCTTGAGCGGGTTGATGCACACCACCTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCA‑GACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGC  >  NC_002947/311578‑311876

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