BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,818,149	T→C	6.6%	intergenic (‑14/+49)	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,149	0	T	C	6.6%	316.5 / 9.9	112	intergenic (‑14/+49)	argZ/argY	tRNA‑Arg/tRNA‑Arg
Reads supporting (aligned to +/- strand): ref base T (49/48); new base C (5/2); total (62/50)
Fisher's exact test for biased strand distribution p-value = 4.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 118 total aligned reads displayed.

GCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGA‑‑TATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAACTTACTTTACTGCAG  >  NC_000913/2818004‑2818298

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)

  ggTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTAct                                                                                                                                                   >  1:1198651/1‑150 (MQ=255)

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

   gTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAACTTACTTTActg                                                                                                                                                  >  2:1638753/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:1573245/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:372144/150‑1 (MQ=255)

               aCCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGA‑‑TATTGATa                                                                                                                                    <  2:87585/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‑‑TATTGATACTGCTACTAAAGc                                                                                                                       >  1:1285374/1‑150 (MQ=16)

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

                             ttCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGA‑‑TATTGATACCGCTACTAAAGcc                                                                                                                      >  1:1935751/1‑150 (MQ=255)

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

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

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

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

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

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

                                   gCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGATACGGATGCATCGGGAAATTTATTTTACTGCTGA‑‑TATTGATACCGCTACTAAAGCCATATCa                                                                                                                >  2:707668/1‑150 (MQ=255)

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

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

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

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

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

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

                                           ttacGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATATAATTGACGGCTTA‑‑TATTTATACCGCAACTAAAGCCATAACAATTAagag                                                                                                        >  2:2069736/4‑150 (MQ=11)

                                              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:1815662/150‑1 (MQ=40)

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

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

                                                gagaATATGGTGCATCCGGGAGGATTCGAACCCCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGA‑‑TATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGt                                                                                                   <  2:434079/150‑1 (MQ=255)

                                                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)

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

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

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

                                                                  ggaggaTTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGA‑‑TATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGa                                                                                 <  2:1258788/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                                                                                                           >  2:1496216/1‑117 (MQ=255)

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

                                                                          agttgatCCGACCGGTGTGTTCGTAGCCTAGTACACAACCCCGCTCAGCTACGCTTGAATGGGCAAATTTATTTTACTGCTGA‑‑TATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGa                                                                         <  1:1575626/144‑1 (MQ=14)

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

                                                                                       tgtgggTTGGTAGCCGAGTAGTCTATCCAGCTGAGCTAAGGTTGCATGGGGAAATTTATTTTACTGCTTA‑‑TTTTGATACCGCTTCTTATGCCATATCAAGTAAGATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGt                                                            <  2:1175604/146‑1 (MQ=18)

                                                                                       gCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGA‑‑TATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGt                                                            >  2:2338554/1‑150 (MQ=255)

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

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

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

                                                                                                 tAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGA‑‑TATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTAc                                                  >  2:1704001/1‑150 (MQ=255)

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

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

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

                                                                                                            tctATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGA‑‑TATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCt                                       >  1:1677523/1‑150 (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)

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

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

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

                                                                                                                           tACGGATGCATCGGGAAATTTATTTTACTGCTGA‑‑TATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCa                                          <  2:1326258/132‑1 (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)

                                                                                                                                         gAAATTTATTTTACTGCTGA‑‑TATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAacttactt          >  1:2301432/1‑150 (MQ=255)

                                                                                                                                              ttatttTACTGCTGA‑‑TATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTg                                                                                                  <  1:743991/57‑1 (MQ=255)

                                                                                                                                              ttatttTACTGCTGA‑‑TATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTg                                                                                                  <  1:2254038/57‑1 (MQ=255)

                                                                                                                                                 ttttACTGCTGA‑‑TATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAACTTACTTTACTGCAg  >  2:27324/1‑150 (MQ=255)

                                                                                                                                                 ttttACTGCTGA‑‑TATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAACTTACTTTACTGCAg  >  2:2324856/1‑150 (MQ=255)

GCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGA‑‑TATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAACTTACTTTACTGCAG  >  NC_000913/2818004‑2818298

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