BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,817,916	T→G	10.6%	intergenic (‑56/+143)	argQ ← / ← argZ	tRNA‑Arg/tRNA‑Arg

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,817,916	0	T	G	10.6%	577.0 / 57.0	198	intergenic (‑56/+143)	argQ/argZ	tRNA‑Arg/tRNA‑Arg
Reads supporting (aligned to +/- strand): ref base T (82/94); new base G (10/11); total (93/105)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01

Only 100 of 202 total aligned reads displayed.

GAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGC  >  NC_000913/2817767‑2818064

gAAGGTTACGGAGAATATGGTGCGTCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAt                                                                                                                                                      <  1:3698631/150‑1 (MQ=255)

 aaGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAtt                                                                                                                                                     <  2:1682277/150‑1 (MQ=255)

   ggTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTaa                                                                                                                                                   <  2:842518/150‑1 (MQ=255)

   ggTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTaa                                                                                                                                                   <  2:4058513/150‑1 (MQ=255)

    gTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTACCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAag                                                                                                                                                  >  1:2983109/1‑150 (MQ=255)

          gagaATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGG                                                                                                                                            <  2:487408‑M1/150‑1 (MQ=255)

            gaATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTG                                                                                                                                          <  2:1892820‑M1/150‑1 (MQ=255)

            gaATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTAGTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTG                                                                                                                                          <  1:1056661‑M1/150‑1 (MQ=255)

             aaTATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGC                                                                                                                                         <  1:2619274‑M1/150‑1 (MQ=255)

             aaTATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGC                                                                                                                                         <  2:3372699‑M1/150‑1 (MQ=255)

                  ggTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCGGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCG                                                                                                                                    >  1:971496‑M1/1‑150 (MQ=255)

                      cATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAG                                                                                                                                <  1:3585551‑M1/150‑1 (MQ=255)

                          cGGGAGGATTCGCACCTCCGCCCGCTCGGTTCGTAGCCGTGTACTCGATCCCGCTGAGCCACGGATTCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATT                                                                                                                            <  2:971496‑M1/150‑1 (MQ=255)

                          cGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATT                                                                                                                            >  2:2428084‑M1/1‑150 (MQ=255)

                                    cGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTg                                                                                                                  <  1:3884327/150‑1 (MQ=255)

                                     gAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGc                                                                                                                 >  2:276209/1‑150 (MQ=255)

                                     gAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTATACTGCTGATATTGATACCGCTCCTAAAGCGAAATCCAAAACGAGATGCTGCATCCCCCAGCATAACTCGCCtttg                                                                                                                 >  1:2450252/1‑148 (MQ=14)

                                        ccTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGc                                                                                                              <  1:2339442/150‑1 (MQ=255)

                                          tCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTc                                                                                                            <  1:1209452/150‑1 (MQ=255)

                                            cGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGc                                                                                                          >  2:302690/1‑150 (MQ=255)

                                            cGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGc                                                                                                          <  1:1321592/150‑1 (MQ=255)

                                              aCCGCTCGGTTCGTAGCCGCGTTCTCTATGCCGCGGAGATTCGGATGCATCGGGTCAACTAAGTCTACTGCTGCTATTGAGACCGCTTCTAACGCCATTTCAATTACGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGccc                                                                                                        <  2:4111477/150‑1 (MQ=17)

                                                  cTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCg                                                                                                    <  1:3152627/150‑1 (MQ=255)

                                                  cTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCg                                                                                                    <  2:521420/150‑1 (MQ=255)

                                                   tCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAAGTTATTATAATGATGATATTTATACCGATACTAAAGCCATATAAATTAGGATATGGTGAATGATGTAGGAAAACTCGCAATGCGGATCGccaaaagg                                                                                                   >  1:176528/1‑144 (MQ=9)

                                                   tCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGATGACATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCgg                                                                                                   >  1:3232009/1‑150 (MQ=255)

                                                     ggTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTCCTGATATTTATACCGCTACTAAAGCCATATCAATTACGAGATGGTGCATCAGGGAGGATTACTCGCGTTTCGTCTCGCCCGTCgcgc                                                                                                 >  2:3471111/1‑150 (MQ=17)

                                                      gTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGcc                                                                                                >  1:1210616/1‑150 (MQ=255)

                                                      gTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGcc                                                                                                <  1:2766564/150‑1 (MQ=255)

                                                      gTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCGACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGcc                                                                                                <  1:3758598/150‑1 (MQ=255)

                                                        tCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTAATCTACTGCTGATATTGATACCTCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGt                                                                                              <  1:2880412/150‑1 (MQ=255)

                                                            aGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACGGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTcgaacctccgaccgctcggttcgtagccgagtac                                                                                          >  2:2861777‑M1/1‑116 (MQ=255)

                                                               cGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTaaa                                                                                       <  2:3612262/150‑1 (MQ=255)

                                                               cGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTaaa                                                                                       <  1:3568126/150‑1 (MQ=255)

                                                                gAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTAATTGTACTGCTGATATTGATACCGCTACTAAAGCCATAGCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCGTGCGGCTCGCCCTTCGGGCCGTTGCTAAAg                                                                                      >  1:1240176/1‑150 (MQ=255)

                                                                    aCTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAAc                                                                                  <  1:2013929/150‑1 (MQ=255)

                                                                     ctctATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACg                                                                                 <  1:688835/150‑1 (MQ=255)

                                                                     ctctATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACg                                                                                 <  1:3364634/150‑1 (MQ=255)

                                                                     ctctATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACg                                                                                 >  2:3223075/1‑150 (MQ=255)

                                                                      tctATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGt                                                                                >  1:4144909/1‑150 (MQ=255)

                                                                      tctATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACGAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGt                                                                                >  1:1068981/1‑150 (MQ=255)

                                                                        tATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACACCCCGTGCGGCTCGCCGTTCGGGCCGTTTCTACAGCACCGTTa                                                                              >  2:933304/1‑150 (MQ=17)

                                                                          tCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTatc                                                                            >  1:2101811/1‑150 (MQ=255)

                                                                             aGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTatcatc                                                                         >  2:1836674/1‑150 (MQ=255)

                                                                              gCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATcc                                                                        <  2:4206417/150‑1 (MQ=255)

                                                                              gCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATcc                                                                        <  1:629648/150‑1 (MQ=255)

                                                                                 gAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTa                                                                     <  2:550752/150‑1 (MQ=255)

                                                                                    cTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTg                                                                  >  2:3555173/1‑150 (MQ=255)

                                                                                       cGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCt                                                               <  1:788731/150‑1 (MQ=255)

                                                                                         gATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGc                                                             <  1:2135436/150‑1 (MQ=255)

                                                                                          aTGCATAGGGAATACTTATTCTTCTGCTGATATTGATACAGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGTAACGTTATCATCCCTAGTGCCTGCg                                                            <  1:408453/150‑1 (MQ=18)

                                                                                              aTCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTa                                                        >  1:3915232/1‑150 (MQ=255)

                                                                                                  ggAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTc                                                    >  1:4140807/1‑150 (MQ=255)

                                                                                                    aaaaCTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATTGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCg                                                  <  2:2798641/150‑1 (MQ=255)

                                                                                                      aaCTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCt                                                <  1:2413650/150‑1 (MQ=255)

                                                                                                       aCTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTa                                               >  1:4241350/1‑150 (MQ=255)

                                                                                                         ttATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTcgaacctccgaccgctcggttcgtagccgagtactctatccagctgagctacggatgcaaatggcggtgaggcggggat                                             <  1:2520924‑M1/150‑80 (MQ=255)

                                                                                                            ttttactGCAGATATTGCTACCGCTACTAACGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTcgaacctccgaccgctcggttcgtagccgagtactctatccagctgagctacggatgcaaatggcggtgaggcggggattcg                                          <  2:78690‑M1/147‑83 (MQ=255)

                                                                                                               tACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCg                                       >  1:1844656/1‑150 (MQ=255)

                                                                                                               tACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCg                                       >  1:1401709/1‑150 (MQ=255)

                                                                                                                aCTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGa                                      <  1:4273911/150‑1 (MQ=255)

                                                                                                                 ctgctgATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGAt                                     <  1:153327/150‑1 (MQ=255)

                                                                                                                  tgctgATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATa                                    >  1:1162208/1‑150 (MQ=255)

                                                                                                                   gctgATTTTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTATAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATAc                                   <  1:212643/150‑1 (MQ=255)

                                                                                                                     tgATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTcgaacctccgaccgctcggttcgtagccgagtactctatccagctgagctacggatgcaaatggcggtgaggcggggattcgaaccccgga                                 >  2:1267507‑M1/1‑59 (MQ=255)

                                                                                                                       atatTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTc                               >  2:2501726/1‑150 (MQ=255)

                                                                                                                       atatTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTcgaacctccgaccgctcggttcgtagccgagtactctatccagctgagctacggatgcaaatggcggtgaggcggggattcgaacccc                                    >  2:1174167‑M1/1‑57 (MQ=255)

                                                                                                                       atatTGAGACCGCTACTATAGGCATATCTATTTAGATATGGTGCATCCGGGTGGAGTACTCTCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTc                               <  1:2318082/150‑1 (MQ=17)

                                                                                                                         atTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTcgaacctccgaccgctcggttcgtagccgagtactctatccagctgagctacggatgcaaatggcggtgaggcggggattcgaaccccgga                                 >  1:1816532‑M1/1‑55 (MQ=255)

                                                                                                                          ttGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCAtt                            >  2:1715402/1‑150 (MQ=255)

                                                                                                                            gATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGc                          >  1:3920236/1‑150 (MQ=255)

                                                                                                                             aTACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGc                                                                                                          <  2:2734695/69‑1 (MQ=40)

                                                                                                                             aTACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGc                                                                                                          >  1:2734695/1‑69 (MQ=40)

                                                                                                                               aCCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCgg                                                                                                               >  2:4282324/1‑62 (MQ=40)

                                                                                                                               aCCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTcgaacctccgaccgctcggttcgtagccgagtactctatccagctgagctacggatgcaaatggcggtgaggcggggattcgaaccccggatgcagctttt                       <  2:2148593‑M1/150‑102 (MQ=255)

                                                                                                                                ccGCTACTCAAGCCATATCTTTTAAGTGTTTGTGCATCCTGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGaa                      <  1:1660675/150‑1 (MQ=17)

                                                                                                                                ccGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGaa                      <  2:2101811/150‑1 (MQ=255)

                                                                                                                                ccGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGaa                      <  1:525713/150‑1 (MQ=255)

                                                                                                                                 cGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCgg                                                                                                   <  2:2279382/72‑1 (MQ=39)

                                                                                                                                 cGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCgg                                                                                                   >  1:2279382/1‑72 (MQ=39)

                                                                                                                                 cGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAg                     >  1:148369/1‑150 (MQ=255)

                                                                                                                                  gCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTcgaacctccgaccgctcggttcgtagccgagtactctatccagctgagctacggatgcaaatggaggtgagcaggggattagaaccccggatgccgcttatgac                    >  2:3148905‑M1/1‑46 (MQ=255)

                                                                                                                                      ctaAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCATGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAAAGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATAACTTCATTGCTGAAGATTTc                >  1:2571901/1‑150 (MQ=255)

                                                                                                                                       taAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCg               <  1:806550/150‑1 (MQ=255)

                                                                                                                                         aaGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGga             <  1:3984592/150‑1 (MQ=255)

                                                                                                                                         aaGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACACGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGATATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGACACCTTCATTGCTGAAGATTTCGga             >  2:2396287/1‑150 (MQ=255)

                                                                                                                                         aaGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTcgaacctccgaccgctcggttcgtagccgagtactctatccagctgagctacggatgcaaatggcg                                                          >  2:1484916‑M1/1‑39 (MQ=255)

                                                                                                                                           gCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGgaga           <  1:1937464/150‑1 (MQ=255)

                                                                                                                                           gCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGgaga           <  1:1148549/150‑1 (MQ=255)

                                                                                                                                            ccATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGaa          >  2:2359641/1‑150 (MQ=255)

                                                                                                                                             cATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAat         >  1:3468344/1‑150 (MQ=255)

                                                                                                                                                atCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATg      <  2:2172853/150‑1 (MQ=255)

                                                                                                                                                atCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATg      >  2:633453/1‑150 (MQ=255)

                                                                                                                                                atCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATg      <  2:2981479/150‑1 (MQ=255)

                                                                                                                                                 tCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATgg     >  2:413208/1‑150 (MQ=255)

                                                                                                                                                   aaTTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTg   >  1:4183595/1‑150 (MQ=255)

                                                                                                                                                   aaTTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTg   >  1:2057325/1‑150 (MQ=255)

                                                                                                                                                   aaTTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCATGCGGTTAATTACGCTAACTCTTCGAAACCTACAATGCTGAAGATTTAGGAGAAAATGGTg   >  1:1780824/1‑150 (MQ=18)

                                                                                                                                                    aTTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGc  >  2:1366204/1‑150 (MQ=255)

                                                                                                                                                    aTTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGc  >  1:1284895/1‑150 (MQ=255)

GAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGC  >  NC_000913/2817767‑2818064

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