BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	311,761	C→G	10.5%	S206T (AGC→ACC) ‡	PP_0256 ←	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	311,761	0	C	G	10.5%	627.6 / 51.0	200	S206T (AGC→ACC) ‡	PP_0256	putative oxidoreductase
Reads supporting (aligned to +/- strand): ref base C (75/104); new base G (6/15); total (81/119)
Fisher's exact test for biased strand distribution p-value = 3.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

Only 100 of 203 total aligned reads displayed.

TGCACACCACCTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCA‑GACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGCACGTCCCGAGGTATAGAATTCGGCCTGGTCAGG  >  NC_002947/311612‑311909

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

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

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

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

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

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

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

         ccTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCA‑GACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCa                                                                                                                                             >  1:5357177/1‑150 (MQ=255)

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

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

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

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

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

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

                 ccACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCA‑GACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAacgccacg                                                                                                                                     >  2:830706/1‑150 (MQ=255)

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

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

                      cTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCA‑GACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCt                                                                                                                                <  2:3756026/150‑1 (MQ=255)

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

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

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

                              cATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCA‑GACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCAt                                                                                                                        <  2:496085/150‑1 (MQ=255)

                              cATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCA‑GACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCAt                                                                                                                        >  1:284610/1‑150 (MQ=255)

                                tCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCA‑GACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATgg                                                                                                                      >  1:4394451/1‑150 (MQ=255)

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

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

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

                                           gTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCA‑GACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGa                                                                                                           <  1:1185856/150‑1 (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)

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

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

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

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

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

                                                          atGATTGGTGCCAGGATTTTGCCCCAC‑GACGAAAATCGCGTCAGCCTGCTCGAAATCCTC‑AAGCAGAACAGTTCCCTTGCCAACCCCTACCGATTGTTTGAGCGCCACACCGCTGGCCTCATGGCACATGTTCGAACAGTCAGGGAAatt                                                                                            <  2:1112865/149‑1 (MQ=255)

                                                          atGATTGGTGCCAGGATTTTGCCCCAC‑GACGAAAATCGCGTCAGCCTGCTCGAAATCCTC‑AAGCAGAACAGTTCCCTTGCCAACCCCTACCGATTGTTTGAGCGCCACACCGCTGGCCTCATGGCACATGTTCGAACAGTCAGGGAAatt                                                                                            <  1:2798487/149‑1 (MQ=255)

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

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

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

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

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

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

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

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

                                                                  tGCCGGGGTTCTGGCCCCA‑GACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGc                                                                                    <  2:2686471/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:1096766/1‑150 (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)

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

                                                                                  ccA‑GACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGaa                                                                    >  2:4420022/1‑150 (MQ=255)

                                                                                    acgACGAAAATCGCGTCAGCCTGCTCGAAATCCTC‑AAGCAGAACAGTTCCCTTGCCAACCCCTACCGATTGTTTGAGCGCCACTCCGCTGGCCTCATGGCACATGTTCGAACAGTCAGGGAAATTGTTCGTGCCATAAGCGCGCACGaaa                                                                   >  1:3044811/3‑149 (MQ=255)

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

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

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

                                                                                          aaaaTGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTgg                                                             <  2:1092420/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)

                                                                                                 cGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGt                                                      <  1:4568972/150‑1 (MQ=255)

                                                                                                  gTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTa                                                     >  1:4706220/1‑150 (MQ=255)

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

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

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

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

                                                                                                              gAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGtt                                         >  1:2081504/1‑150 (MQ=255)

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

                                                                                                                 gtcgtcGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCt                                      >  2:2783412/1‑150 (MQ=255)

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

                                                                                                                      cgAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGCAc                                 <  1:475047/150‑1 (MQ=255)

                                                                                                                         aGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGCACGTc                              >  1:1823669/1‑150 (MQ=255)

                                                                                                                          ggTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGCACGTcc                             >  2:4347772/1‑150 (MQ=255)

                                                                                                                           atG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGCACGTccc                            <  1:3788868/149‑1 (MQ=255)

                                                                                                                             g‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGCACGTCCCGa                          <  2:1714605/150‑1 (MQ=255)

                                                                                                                               aCTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGCACGTCCCGAg                         <  1:2170066/150‑1 (MQ=255)

                                                                                                                               aCTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGCACGTCCCGAg                         <  1:2031952/150‑1 (MQ=255)

                                                                                                                               aCTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGCACGTCCCGAg                         <  2:3212009/150‑1 (MQ=255)

                                                                                                                               aCTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTAATGGCTCATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGCACGTCCCgtg                         <  1:5406883/150‑3 (MQ=255)

                                                                                                                                   tACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGCACGTCCCGa                          >  1:654267/1‑145 (MQ=255)

                                                                                                                                    aCCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGCACGTCCCGAGGtata                    >  2:647013/1‑150 (MQ=255)

                                                                                                                                        ttGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGCACGTCCCGAGGTATAGAAt                >  2:1405064/1‑150 (MQ=255)

                                                                                                                                            cGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGCACGTCCCGAGGTATAGAATTCgg            <  1:3953223/150‑1 (MQ=255)

                                                                                                                                             gACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCGCATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGCACGTCCCGAGGTATAGAATTCGGc           >  2:1799777/1‑150 (MQ=255)

                                                                                                                                             gACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGCACGTCCCGAGGTATAGAATTCGGc           <  1:537552/150‑1 (MQ=255)

                                                                                                                                                 ccGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGc                                   <  2:4212444/122‑1 (MQ=255)

                                                                                                                                                 ccGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGc                                   >  1:4212444/1‑122 (MQ=255)

                                                                                                                                                  cGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGCACGTCCCGAGGTATAGAATTCGGCCTGGt      >  2:3338493/1‑150 (MQ=255)

                                                                                                                                                  cGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGCACGTCCCGAGGTATAGAATTCGGCCTGGt      >  1:4364285/1‑150 (MQ=255)

                                                                                                                                                  cGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGCACGTCCCGAGGTATAGAATTCGGCCTGGt      >  2:1691195/1‑150 (MQ=255)

                                                                                                                                                   gACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGCAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGCACGTCCCGAGGTATAGAATTCGGCCTGGTc     >  1:2244633/1‑150 (MQ=255)

                                                                                                                                                     cGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGCACGTCCCGAGGTATAGAATTCGGCCTGGTCAg   <  1:5285238/150‑1 (MQ=255)

                                                                                                                                                      gCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGCACGTCCCGAGGTATAGAATTCGGCCTGGTCAgg  <  1:4342832/150‑1 (MQ=255)

TGCACACCACCTGGGCGCCACGCTTTACCGCATCACGCAGTGGGTCGAGCATGCGCGGGTGATTGGTGCCGGGGTTCTGGCCCCA‑GACAAAAATGGCGTCGGCGTGCTCGAAGTCGTCGAAGGTG‑ACTGTACCTTTGCCGACCCCGACGCTTTGCCCCAACGCCACGCCGCTGGCCTCATGGCACATGTTCGAGCAGTCAGGGAAGTTGTTGGTGCCGTAGGCGCGCACGAACAGTTGGTACAGGTACGCCGCTTCGTTGCTGGCACGTCCCGAGGTATAGAATTCGGCCTGGTCAGG  >  NC_002947/311612‑311909

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