BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	961,713	G→C	5.9%	H235Q (CAC→CAG)	selD ←	selenide, water dikinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	961,713	0	G	C	5.9%	839.4 / 29.9	248	H235Q (CAC→CAG)	selD	selenide, water dikinase
Reads supporting (aligned to +/- strand): ref base G (92/141); new base C (7/8); total (99/149)
Fisher's exact test for biased strand distribution p-value = 5.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 253 total aligned reads displayed.

CTGATCTTGTGCCCGTAGCTGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCCTTCTTTTCGGCAGT  >  NC_002947/961564‑961862

cTGATCTTGTGCCCGTAGCTGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAgg                                                                                                                                                       <  2:2526603/150‑1 (MQ=255)

  gATCTTGTGCCCGTAGCTGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTg                                                                                                                                                     >  2:3694549/1‑150 (MQ=255)

     cTTGTGCCCGTAGCTGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGcc                                                                                                                                                  <  2:536099/150‑1 (MQ=255)

      ttGTGCCCGTAGCTGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCg                                                                                                                                                 >  1:3789430/1‑150 (MQ=255)

         tgCCCGTAGCTGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCACCTGGTCGAGCTGGCCGaag                                                                                                                                              <  2:5413724/150‑3 (MQ=255)

           cccGTAGCTGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAg                                                                                                                                            <  2:1491778/150‑1 (MQ=255)

            ccGTAGCTGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAgg                                                                                                                                           >  1:3024100/1‑150 (MQ=255)

             cGTAGCTGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGc                                                                                                                                          >  2:932413/1‑150 (MQ=255)

             cGTAGCTGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGc                                                                                                                                          <  1:1322027/150‑1 (MQ=255)

             cGTAGCTGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGc                                                                                                                                          <  2:458141/150‑1 (MQ=255)

                 gCTGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGaa                                                                                                                                      <  1:5457480/150‑1 (MQ=255)

                  cTGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAg                                                                                                                                     >  1:652263/1‑150 (MQ=255)

                   tGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGc                                                                                                                                    >  2:4015560/1‑150 (MQ=255)

                     tCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGGCCACCTGGTCGAGCTGGc                                                                                                                                          <  1:4411447/142‑1 (MQ=255)

                     tCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGGCCACCTGGTCGAGCTGGc                                                                                                                                          >  2:4410941/1‑142 (MQ=255)

                       gTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCTGGTGGCCGAGCAGGCCGAAGCCGGt                                                                                                                                <  1:2117074/150‑1 (MQ=255)

                       gTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGt                                                                                                                                >  1:1361017/1‑150 (MQ=255)

                       ctgcatGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGt                                                                                                                                <  2:722448/145‑1 (MQ=255)

                        tAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGc                                                                                                                                          <  2:3489189/139‑1 (MQ=255)

                        tAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTg                                                                                                                               <  2:4143764/150‑1 (MQ=255)

                            tGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGt                                                                                                                           <  2:1928846/150‑1 (MQ=255)

                              cgcAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCg                                                                                                                         >  1:1149424/1‑150 (MQ=255)

                              cgcAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCg                                                                                                                         >  2:3601296/1‑150 (MQ=255)

                               gcAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCgg                                                                                                                        >  2:1526845/1‑150 (MQ=255)

                               gcAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCgg                                                                                                                        >  2:2764416/1‑150 (MQ=255)

                                cAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGt                                                                                                                       <  2:227581/150‑1 (MQ=255)

                                    gTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATg                                                                                                                   <  1:3399321/150‑1 (MQ=255)

                                    gTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATg                                                                                                                   <  1:3228291/150‑1 (MQ=255)

                                    gTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATg                                                                                                                   <  1:1966321/150‑1 (MQ=255)

                                       ccACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGcc                                                                                                                <  1:2457223/150‑1 (MQ=255)

                                        cACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGATACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCCGGCCGAAGCCGGTGACGTCGGTCATGGCCt                                                                                                               >  1:1030919/1‑150 (MQ=255)

                                         aCCGGGGATGCAGCCTTCGGCCAGGTACTGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCCGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCtt                                                                                                              >  1:4963755/1‑150 (MQ=255)

                                            ggggATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGAc                                                                                                           <  2:4846840/150‑1 (MQ=255)

                                             gggATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGAcc                                                                                                          <  2:5516820/150‑1 (MQ=255)

                                                aTGCAGCATTCGGCCAGGTAATGGTCGACACTAGTCAGGCGTGGCACGGCGGAGTATTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCATGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCg                                                                                                       <  1:353482/150‑1 (MQ=255)

                                                  gCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTc                                                                                                     <  2:2941037/150‑1 (MQ=255)

                                                   cAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTcc                                                                                                    <  2:4756868/150‑1 (MQ=255)

                                                    aGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCa                                                                                                   <  1:2273030/150‑1 (MQ=255)

                                                    aGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCa                                                                                                   >  2:4715498/1‑150 (MQ=255)

                                                    aGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCa                                                                                                   <  1:2345239/150‑1 (MQ=255)

                                                      ccTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGt                                                                                                 >  1:5038529/1‑150 (MQ=255)

                                                        ttCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGttt                                                                                               <  2:853808/150‑1 (MQ=255)

                                                        ttCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCACCTGGTCGAGCTGGCCGaa                                                                                                                                               <  2:1192592/102‑2 (MQ=38)

                                                          cGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGc                                                                                             <  1:1263181/150‑1 (MQ=255)

                                                           ggCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCGGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGcc                                                                                            <  1:1500960/150‑1 (MQ=255)

                                                              cAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGaa                                                                                         >  1:2695961/1‑150 (MQ=255)

                                                               aGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAg                                                                                        >  2:2318691/1‑150 (MQ=255)

                                                                ggTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGc                                                                                       >  2:4276693/1‑150 (MQ=255)

                                                                  tAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCgg                                                                                     <  1:4505763/150‑1 (MQ=255)

                                                                   aaTGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGc                                                                                    <  2:3640571/150‑1 (MQ=255)

                                                                   aaTGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGAGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGc                                                                                    >  2:281506/1‑150 (MQ=255)

                                                                         cGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCgg                                                                              <  1:4274957/150‑1 (MQ=255)

                                                                            cacTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCg                                                                           >  2:2228666/1‑150 (MQ=255)

                                                                            cacTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGACTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCg                                                                           >  2:4917513/1‑150 (MQ=255)

                                                                                 ggcaggcGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGc                                                                            >  2:3092191/1‑144 (MQ=255)

                                                                                   caggcGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAg                                                                    <  1:4842413/150‑1 (MQ=255)

                                                                                    aggcGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAgg                                                                   <  1:5523803/150‑1 (MQ=255)

                                                                                      gcGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAATCGGCTGCCGGGCGTGTTGAGGGt                                                                 <  2:2083651/150‑1 (MQ=255)

                                                                                      gcGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGt                                                                 <  1:2774981/150‑1 (MQ=255)

                                                                                      gcGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGt                                                                 <  2:780959/150‑1 (MQ=255)

                                                                                      gcGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGt                                                                 <  1:3504226/150‑1 (MQ=255)

                                                                                      gcGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGt                                                                 >  2:2024475/1‑150 (MQ=255)

                                                                                       cGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTg                                                                <  1:2016897/150‑1 (MQ=255)

                                                                                        gTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCACGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGc                                                               >  1:4402444/1‑150 (MQ=255)

                                                                                         tGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCATGACCCCGTCCAGTTTGCCCAAGCGGCTGCCCGGCGGGTTGAGGGTGca                                                              >  1:2828493/1‑150 (MQ=255)

                                                                                                  ggcgTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTc                                                     <  2:4295592/150‑1 (MQ=255)

                                                                                                  ggcgTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTc                                                     <  1:3694549/150‑1 (MQ=255)

                                                                                                    cgTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCACGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAATCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTcgc                                                   >  1:5064313/1‑150 (MQ=255)

                                                                                                      tAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGa                                                 <  1:4593686/150‑1 (MQ=255)

                                                                                                      gcagggtggTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGa                                                 <  2:3598525/143‑1 (MQ=255)

                                                                                                        gTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCGCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGGGTTGTGGGTGCACATCCAGTCGCGAGc                                               >  1:1399149/1‑150 (MQ=255)

                                                                                                             ggTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGc                                          <  2:5105361/150‑1 (MQ=255)

                                                                                                              gTGCGCGGTGTGGCCGCTGCCTTCGGCCAGCTCGACCTGGTGGCCGAGCAGGCCAAAGCCGGTGAGGTCGGTGATGGCCTTGACCCCGTCCTGTTTGCCGTAGCGGCTGCCGGGCGTGTTGAGGGTGGACATCCAGTCGCGAGCCAgtgc                                         <  2:4268031/150‑4 (MQ=255)

                                                                                                               ggcgcGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCt                                        <  2:176771/149‑1 (MQ=255)

                                                                                                                gcgcGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCtt                                       <  2:2513235/150‑1 (MQ=255)

                                                                                                                 cgcgGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTg                                      <  2:3187137/150‑1 (MQ=255)

                                                                                                                 cgcgGTGAGGCCGCTGCCATCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTg                                      <  1:251809/150‑1 (MQ=255)

                                                                                                                  gcgGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTgg                                     <  2:1331037/150‑1 (MQ=255)

                                                                                                                       gAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTg                                <  1:4053904/150‑1 (MQ=255)

                                                                                                                         ggCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACAACCAGTCGCGAGCCAGGCCTTGGTCCTGca                              >  1:3479882/1‑150 (MQ=255)

                                                                                                                           ccGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGGGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGcaca                            >  2:624841/1‑150 (MQ=255)

                                                                                                                               tGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATc                                                          <  1:5322730/116‑1 (MQ=255)

                                                                                                                                    tCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTg                   <  1:1580281/150‑1 (MQ=255)

                                                                                                                                    tCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTg                   <  2:5209302/150‑1 (MQ=255)

                                                                                                                                      ggCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGc                 <  1:1071301/150‑1 (MQ=255)

                                                                                                                                      ggCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGc                 <  2:2728627/150‑1 (MQ=255)

                                                                                                                                       gCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCgg                                                                                     <  2:5548140/81‑1 (MQ=255)

                                                                                                                                       gCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCgg                                                                                     >  1:5548140/1‑81 (MQ=255)

                                                                                                                                       gCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGcc                >  2:3097431/1‑150 (MQ=255)

                                                                                                                                       gCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGcc                <  2:4464774/150‑1 (MQ=255)

                                                                                                                                         cAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCctt              >  1:2037682/1‑150 (MQ=255)

                                                                                                                                          tgCTCGGCCACCTGGTCGAGCTGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCcttc             >  2:3041604/2‑150 (MQ=255)

                                                                                                                                          aGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCcttc             <  1:3621294/150‑1 (MQ=255)

                                                                                                                                            cTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGTCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCcttctt           <  1:799703/150‑1 (MQ=255)

                                                                                                                                             tCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCCTTCttt          <  1:4006830/150‑1 (MQ=255)

                                                                                                                                             tCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCCTTCttt          >  2:2143502/1‑150 (MQ=255)

                                                                                                                                               gACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCCTTCTTTTc        >  1:2223920/1‑150 (MQ=255)

                                                                                                                                                aCCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCCTTCTTTTCg       <  1:3083574/150‑1 (MQ=255)

                                                                                                                                                    ggTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCCTTCTTTTCGGCAg   >  1:943066/1‑150 (MQ=255)

                                                                                                                                                     gTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCCTTCTTTTCGGCAgt  >  2:3820202/1‑150 (MQ=255)

CTGATCTTGTGCCCGTAGCTGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCCTTCTTTTCGGCAGT  >  NC_002947/961564‑961862

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