BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,818,145	T→C	6.4%	intergenic (‑10/+53)	argZ ← / ← argY	tRNA‑Arg/tRNA‑Arg

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,818,145	0	T	C	6.4%	325.7 / 9.1	112	intergenic (‑10/+53)	argZ/argY	tRNA‑Arg/tRNA‑Arg
Reads supporting (aligned to +/- strand): ref base T (54/46); new base C (4/3); total (62/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

Only 100 of 116 total aligned reads displayed.

AGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGA‑‑TATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAACTTACTTTAC  >  NC_000913/2817997‑2818293

aGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAatt                                                                                                                                                       <  2:289148/150‑1 (MQ=255)

aGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAatt                                                                                                                                                       <  2:2364097/150‑1 (MQ=255)

aGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAatt                                                                                                                                                       <  2:1549245/150‑1 (MQ=255)

       gCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTa                                                                                                                                                <  2:875588/150‑1 (MQ=255)

        cGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTAc                                                                                                                                               >  1:1008592/1‑150 (MQ=255)

        cGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTAc                                                                                                                                               >  1:1008500/1‑150 (MQ=255)

         ggTTAGTTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTAct                                                                                                                                              >  1:1607561/1‑150 (MQ=255)

         ggTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGACCCACCGACCCGGGGGATCGTAGCCGCGTCGTCTATCGAGATGAGCTCCGGATGGGTCGGGAAATATGTGTTAcg                                                                                                                                              >  2:913986/1‑149 (MQ=14)

         ggTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTAGAACCTCCGACCGCTCGGTTCGTAGCAGAGGACTCTATCCAGCTGAGCTACGGATGCCTCGGGAAATTTAGTTTAct                                                                                                                                              >  1:1198826/1‑150 (MQ=14)

          gTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAACTTACTTTActg                                                                                                                                             >  2:1638753/1‑150 (MQ=255)

          gTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAACTTACTTTActg                                                                                                                                             >  2:1638703/1‑150 (MQ=255)

               ttCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTCCGGATGCATCGGGGACTGTATTGTACTGCTGag                                                                                                                                        >  1:827226/1‑149 (MQ=14)

               ttCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGa‑‑t                                                                                                                                      >  1:1077918/1‑150 (MQ=255)

               ttCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGa‑‑t                                                                                                                                      >  1:1442821/1‑150 (MQ=255)

                  cGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGA‑‑TAtt                                                                                                                                   <  2:1194106/150‑1 (MQ=255)

                   gCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGA‑‑TATTg                                                                                                                                  <  1:490120/150‑1 (MQ=255)

                   gCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGA‑‑TATTg                                                                                                                                  <  2:2256130/150‑1 (MQ=255)

                    cTACCTCTTCGATACCTTCCTTGCTGAAGATGTCGGAGTATATGGTGCATCCGGGAGGATTCGACCCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGA‑‑TATTGa                                                                                                                                 <  2:1397995/150‑1 (MQ=255)

                    cTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGA‑‑TATTGa                                                                                                                                 <  2:1391725/150‑1 (MQ=255)

                    cTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGA‑‑TATTGa                                                                                                                                 <  2:1573248/150‑1 (MQ=255)

                    cTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGA‑‑TATGGa                                                                                                                                 >  1:1642339/1‑150 (MQ=255)

                    cTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTAATTTACTGCTGA‑‑TATTGa                                                                                                                                 >  1:924145/1‑150 (MQ=255)

                    cTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTAGAACCTCAGACCGCACGGGTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGAGGCATCGGGAAATTTAGTTTACTGGTGA‑‑TATTGa                                                                                                                                 >  1:156395/1‑150 (MQ=14)

                     tACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGA‑‑TATTGAt                                                                                                                                <  2:1819793/150‑1 (MQ=255)

                      aCCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGA‑‑TATTGATa                                                                                                                               <  2:87585/150‑1 (MQ=255)

                      aCCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGA‑‑TATTGATa                                                                                                                               <  2:372144/150‑1 (MQ=255)

                      aCCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGA‑‑TATTGATa                                                                                                                               <  2:2235337/150‑1 (MQ=255)

                      aCCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAACTTACTTTACTGCAGA‑‑ttttttga                                                                                                                               <  1:2210182/150‑5 (MQ=14)

                         tctTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATGCAGCTGAGCTACGGATGCATCTGGATATTTATTTTAATGCTGA‑‑TATTGATAACg                                                                                                                            >  1:120965/1‑150 (MQ=14)

                           ttCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGAACTCAAACCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGA‑‑TATTGATACCGct                                                                                                                          >  1:1354409/1‑150 (MQ=255)

                                tACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAACTTACTTTACTGCAGATTTTTTGATACCGctact                                                                                                                       >  1:693066/1‑150 (MQ=11)

                                   cTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGA‑‑TATTGATACCGCTACTAAAGc                                                                                                                  >  2:2397472/1‑150 (MQ=18)

                                    ttCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGA‑‑TATTTATACAGCTACTAAAGcc                                                                                                                 >  1:122226/1‑150 (MQ=255)

                                    ttCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGA‑‑TATTGATACCGCTACTAAAGcc                                                                                                                 <  1:1493068/150‑1 (MQ=255)

                                    ttCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGA‑‑TATTGATACCGCTACTAAAGcc                                                                                                                 <  2:1362030/150‑1 (MQ=255)

                                    ttCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGA‑‑TATTGATACCGCTACTAAAGcc                                                                                                                 <  1:1736396/150‑1 (MQ=255)

                                         tGCTGAAGATTTCGGAGAATCTGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGA‑‑TATTGATACCGCTACTAAAGCCATATc                                                                                                            <  2:1589820/150‑1 (MQ=255)

                                          gCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGA‑‑TATTGATACCGCTACTAAAGCCATATCa                                                                                                           >  2:1265846/1‑150 (MQ=255)

                                          gCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGA‑‑TATTGATACCGCTACTAAAGCCATATCa                                                                                                           >  2:1241668/1‑150 (MQ=255)

                                           cTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGA‑‑TATTGATACCGCTACTAAAGCCATATCaa                                                                                                          >  2:2287270/1‑150 (MQ=255)

                                           cTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAGTTACTTTACTGCTGA‑‑TATTGATACCGCTACTAAAGCCATATCaa                                                                                                          >  2:164079/1‑150 (MQ=255)

                                           cTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAACTTACTTTACTGCAGA‑‑ttttt                                                                                                                                  <  1:482201/126‑2 (MQ=11)

                                           cTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAACTTACTTTACTGCAGA‑‑ttttt                                                                                                                                  >  2:482201/1‑125 (MQ=11)

                                                 tttacGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTTGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGA‑‑TATTGATACCGCTACTAAAGCCATATCAAGTAaga                                                                                                    <  1:970837/146‑1 (MQ=255)

                                                 gttacGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGA‑‑TATTGATACTGCTACTAAAGCCATATCAAGTAaga                                                                                                    <  2:1285374/146‑1 (MQ=255)

                                                 gttacGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGA‑‑TATTGATACCGCTACTAAAGCCATATCAAGTAaga                                                                                                    <  1:2397472/146‑1 (MQ=38)

                                                     cGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTAATTTACTGCTGA‑‑TATTGATACCGCTACTAAAGCCAGATAAAGAAAGAGATg                                                                                                >  1:1711474/1‑150 (MQ=255)

                                                     cGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAATTAATTTACTGCTGA‑‑TATTGATAACGCTCCTCAAGCCCAATCAAGAGAGAGATg                                                                                                >  1:849633/1‑150 (MQ=11)

                                                      ggAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGA‑‑TATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATgg                                                                                               <  1:1142079/150‑1 (MQ=40)

                                                       gagaATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGA‑‑TATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGt                                                                                              <  2:937765/150‑1 (MQ=40)

                                                       gagaATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGA‑‑TATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGt                                                                                              <  2:1815662/150‑1 (MQ=40)

                                                       gagaATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAATTATTTTACTGCATA‑‑TATTGATAACGCTAATAAATACATATCAAATAAGAGATgtt                                                                                              >  2:2324888/1‑148 (MQ=11)

                                                       gagaATATGGTGCATCCGGGAGGACTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGA‑‑TATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGt                                                                                              <  2:680522/150‑1 (MQ=255)

                                                        agaATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGA‑‑TATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTg                                                                                             <  2:220833/150‑1 (MQ=40)

                                                        agaATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTGTACTGCTGA‑‑TATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTg                                                                                             <  2:832101/150‑1 (MQ=38)

                                                         gaATATGGTGCATCCGGGAGGATTCGAACCTACGACCGCTCGGTTCGTAGCCGAGAACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGA‑‑TATTGATACCGCTAGTAAAGCCATATCTAGTGGGAGAGGGTGc                                                                                            >  1:448565/1‑150 (MQ=14)

                                                                gTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGA‑‑TATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCggg                                                                                     <  1:2069736/150‑1 (MQ=40)

                                                                gTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGA‑‑TATTGATAACGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCggg                                                                                     >  2:740191/1‑150 (MQ=255)

                                                                gTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCAACGGGAAATTTATTTAAATGCTGA‑‑TATTGATAACGCTACTCAAGCCATATCAAGTAAGAGATGGTGCATCCggg                                                                                     >  2:25099/1‑150 (MQ=14)

                                                                 tGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGA‑‑TATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCgg                                                                                      <  1:2052677/148‑1 (MQ=40)

                                                                 tGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGA‑‑TATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCgg                                                                                      >  2:2052677/1‑148 (MQ=40)

                                                                         ggaggaTTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGA‑‑TATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGa                                                                            <  2:1258788/150‑1 (MQ=255)

                                                                          gaggaTTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGA‑‑TATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGaa                                                                           <  1:206143/150‑1 (MQ=255)

                                                                           aggaTTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGA‑‑TATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAAc                                                                          >  1:2190368/1‑150 (MQ=255)

                                                                           aggaTTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGA‑‑TATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAAc                                                                          <  2:1165598/150‑1 (MQ=255)

                                                                             gaTTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGA‑‑TATTGATACCGCTACTAAAGCCATATCAAGTACGAGATGGTGCATCCGGGAGGAGTAGaagcg                                                                        >  1:630009/1‑147 (MQ=255)

                                                                              tttCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTCTCCAGCTGAGCTACGGATGCTTCGGGATATTTATTTTACTGCTGA‑‑TATTGCTACCGCTACTTTAGCCATATCAAGTATGAGATGGTGCATCCGGGAGGATTCGAACCTc                                                                       <  1:235090/149‑1 (MQ=14)

                                                                                tCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGA‑‑TATTGATACCGCTACTAAAGCCATATCAAGTaa                                                                                                      <  1:1496216/117‑1 (MQ=255)

                                                                                tCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGA‑‑TATTGATACCGCTACTAAAGCCATATCAAGTaa                                                                                                      >  2:1496216/1‑117 (MQ=255)

                                                                                 cGAACCTCCGACCGCCCGGTTCGTAGCCGTGTACTCTATCCCGCTGCGCTACGGATGCATCGGGAAATTTATTTTACTGCTGA‑‑TATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGa                                                                    <  2:575765/150‑1 (MQ=255)

                                                                                   aaCCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGA‑‑TATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGAcc                                                                  >  1:788792/1‑150 (MQ=255)

                                                                                   aaCCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGA‑‑TATTGATACCGCTACTAAAGCCATATCAAGAAAGAGATGGTGAATCCGGGAGGATTCGAAACTCCGAcc                                                                  >  1:1272655/1‑150 (MQ=255)

                                                                                               cTCGGTTCGTAGCCGCGTTCTCCATCCAGCTGTGATACGGATGCATCGGGAAAGTTATTTTACTGCTGA‑‑TCTTGATACCGCTACTACAGCCATATCTAGTATGATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTa                                                      <  1:1194272/150‑1 (MQ=37)

                                                                                                tCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGA‑‑TATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTac                                                     >  1:1175604/1‑149 (MQ=255)

                                                                                                 cGGTTCGGAGCCGAGTTCTCCATCCAGCTGCGCATCGGCTGCATCGGTAAACTTATTTTACTGCTGA‑‑TATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGc                                                    <  2:296696/150‑1 (MQ=18)

                                                                                                       gTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGA‑‑TATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTa                                              <  1:2243468/150‑1 (MQ=255)

                                                                                                         aGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGA‑‑TATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTAct                                            <  1:2338554/150‑1 (MQ=255)

                                                                                                              aGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGA‑‑TATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATc                                       <  1:28200/150‑1 (MQ=255)

                                                                                                                  ctctATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGA‑‑TATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGc                                   >  2:1448977/1‑150 (MQ=255)

                                                                                                                   tctATCCAGCTGCGCTACGGATGCATCGGGAAATTTATTTTACTGCTGA‑‑TATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCt                                  <  2:1808054/150‑1 (MQ=255)

                                                                                                                   tctATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGA‑‑TATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCt                                  >  2:435840/1‑150 (MQ=255)

                                                                                                                   tctATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGA‑‑TATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCt                                  >  1:1677523/1‑150 (MQ=255)

                                                                                                                   tctATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGA‑‑TATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCt                                  <  1:1984076/150‑1 (MQ=255)

                                                                                                                   tctATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGA‑‑TATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCACGGTTCGTAGCCGAGTACTCTATCCAGCt                                  >  1:2409511/1‑150 (MQ=255)

                                                                                                                    ctATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGA‑‑TATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCTGAGCCGAGTACTCTATCCACCTg                                 >  2:1898269/1‑150 (MQ=255)

                                                                                                                    ctATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGA‑‑TATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTg                                 >  2:2127412/1‑150 (MQ=255)

                                                                                                                     tATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGA‑‑TATTGATACCGCTACTACAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGa                                <  1:707668/150‑1 (MQ=255)

                                                                                                                     tATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGA‑‑TATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGa                                <  1:1241668/150‑1 (MQ=255)

                                                                                                                     tATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGA‑‑TATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGa                                <  1:1265846/150‑1 (MQ=255)

                                                                                                                          aGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGA‑‑TATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTAc                           >  1:686932/1‑150 (MQ=255)

                                                                                                                              gAGCTACGTATGTATCGGTAAATTTATTTTACTGCTGA‑‑TATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGAt                       <  2:1834691/150‑1 (MQ=255)

                                                                                                                               aGCTACGGATGCATCGGGAAATTTATTTTACTGCTGA‑‑TATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATg                      >  2:452341/1‑150 (MQ=255)

                                                                                                                               aGCTACGGATGCATCGGGAAATTTATTTTACTGCTGA‑‑TATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATg                      >  2:2296251/1‑150 (MQ=255)

                                                                                                                                  tACGGATGCATCGGGAAATTTATTTTACTGCTGA‑‑TATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCa                                     >  1:1326258/1‑132 (MQ=255)

                                                                                                                                  tACGGATGCATCGGGAAATTTATTTTACTGCTGA‑‑TATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCAt                   >  1:101422/1‑150 (MQ=255)

                                                                                                                                  tACGGATGCATCGGGAAATTTATTTTACTGCTGA‑‑TATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGATGACCTACGGATGCaa                   >  1:1735555/1‑149 (MQ=255)

                                                                                                                                     ggATGCATCGGGAAATTTATTTTACTGCTGA‑‑TATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGGTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCgg                >  2:1594131/1‑150 (MQ=255)

                                                                                                                                     ggATGCATCGGGAAATTTATTTTACTGCTGA‑‑TATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCgg                >  2:502517/1‑150 (MQ=255)

                                                                                                                                     ggATGCATCGGGAAATTTATTTTACTGCTGA‑‑TATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCgg                >  2:1311046/1‑150 (MQ=255)

                                                                                                                                                   atttatttTACTGCTGA‑‑TATTGATACCGCTACTACAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAACttactttac  <  2:1270714/150‑1 (MQ=255)

AGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGA‑‑TATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAACTTACTTTAC  >  NC_000913/2817997‑2818293

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