BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,817,916	T→G	11.1%	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	11.1%	392.4 / 39.9	136	intergenic (‑56/+143)	argQ/argZ	tRNA‑Arg/tRNA‑Arg
Reads supporting (aligned to +/- strand): ref base T (37/84); new base G (7/8); total (44/92)
Fisher's exact test for biased strand distribution p-value = 2.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.16e-01

Only 100 of 139 total aligned reads displayed.

GAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGC  >  NC_000913/2817767‑2818064

          gagaATGTGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGG                                                                                                                                            <  2:2469267‑M1/150‑1 (MQ=255)

          gagaATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAATACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTATGAGATGG                                                                                                                                            >  1:2782180‑M1/1‑150 (MQ=255)

          gagaATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGG                                                                                                                                            >  2:520431‑M1/1‑150 (MQ=255)

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

           agaATATGGTGCATCCTGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGT                                                                                                                                           <  2:3021651‑M1/150‑1 (MQ=255)

           agaATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGT                                                                                                                                           <  2:298270‑M1/150‑1 (MQ=255)

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

                atGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATC                                                                                                                                      <  2:2575064‑M1/150‑1 (MQ=255)

                atGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATC                                                                                                                                      <  2:2575134‑M1/150‑1 (MQ=255)

                 tGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCC                                                                                                                                     >  1:392479‑M1/1‑150 (MQ=255)

                 tGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCC                                                                                                                                     >  2:1222456‑M1/1‑150 (MQ=255)

                   gTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGG                                                                                                                                   <  1:1865541‑M1/150‑1 (MQ=255)

                    tGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGG                                                                                                                                  >  2:1274952‑M1/1‑150 (MQ=255)

                        tCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGA                                                                                                                              <  1:2153641‑M1/150‑1 (MQ=255)

                             gaggaTTAGAACCTCCTACCGCTCGGTACGTAGACGAGTACTCTATCCACCTGAGATACGTATTCATCTCGAAAACTTATTCTTCTTCTGATATTGATAAATCTACTAAAGACATATCAATTAAGAGATGGTGCATCCGGGAGGATTACt                                                                                                                         <  2:843734/150‑1 (MQ=14)

                                gaTTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGc                                                                                                                      <  1:1438285/150‑1 (MQ=255)

                                gaTTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGc                                                                                                                      <  1:1438366/150‑1 (MQ=255)

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

                                      aaCCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCg                                                                                                                <  1:2032675/150‑1 (MQ=255)

                                       aCCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCgg                                                                                                               <  1:2982372/150‑1 (MQ=255)

                                       aCCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCgg                                                                                                               <  1:659710/150‑1 (MQ=255)

                                           ccGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCg                                                                                                           <  2:480440/150‑1 (MQ=255)

                                             gACCGCTCGTTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATGGTGATTACTTATTCTACTGCTGCTATTGATACCGCTACTAACGACATATCATTTATGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGcc                                                                                                         <  2:1639135/150‑1 (MQ=17)

                                                 gCTCGGTTCGTAGCCGAGTACTCGATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTc                                                                                                     <  1:867796/150‑1 (MQ=255)

                                                     ggTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGc                                                                                                 >  2:1866227/1‑150 (MQ=255)

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

                                                       ttCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCg                                                                                               >  1:1854816/1‑150 (MQ=255)

                                                       ttCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGAGGCTCGCCCTTCGGGCCg                                                                                               >  1:376074/1‑150 (MQ=255)

                                                         cGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAATAGTTAGTCTAGTGCTGATATTGCTTCCGCTCCTAACGCCATATCAATTATGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGtt                                                                                             <  1:2140364/150‑1 (MQ=17)

                                                         cGTAGCCGAGTACTCTATCCAGCTGAGCGACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGCTACCGCTACTAACGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGtt                                                                                             <  1:2140311/150‑1 (MQ=255)

                                                          ttAGCCGTGGACTCGATCGAGCTGACCTACTCATTCATCGTTACAACATATTCTACATCTTTTATTGAGACATCTTCTAAGGACATATCTATTATGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTg                                                                                            <  1:875475/149‑1 (MQ=9)

                                                          gTAGCCGAGTACTCTTTCCTGCTGAGATACGGATGCATCGGTAAAACTTATTCTACTTCTGATATTGATACCTCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTg                                                                                            <  2:3054051/150‑1 (MQ=17)

                                                          gTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTg                                                                                            <  1:897111/150‑1 (MQ=255)

                                                          gTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTg                                                                                            <  1:1225981/150‑1 (MQ=255)

                                                           tAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAATGCCATAGCAATTAAGGGATGGTGCATCCGGGAGGATTACTCGCCGTGCGGCTCGCCATGCGGGCCGTTga                                                                                           >  2:2676095/1‑149 (MQ=17)

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

                                                               cGAGGACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTACTCTACTGCTGATATTGATTCCGCTACTCTAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTaaa                                                                                       <  1:2476906/150‑1 (MQ=17)

                                                                 aGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATGCCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGc                                                                                     <  2:714925/150‑1 (MQ=255)

                                                                 aGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGc                                                                                     <  2:714798/150‑1 (MQ=255)

                                                                  gTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCg                                                                                                    >  2:2035105/1‑134 (MQ=255)

                                                                  gTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCg                                                                                                    <  1:2035105/134‑1 (MQ=255)

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

                                                                      tctATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGCGCCGTTGCTAAAGCAACGt                                                                                >  1:3091322/1‑150 (MQ=255)

                                                                             aGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTatcatc                                                                         <  1:2999377/150‑1 (MQ=255)

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

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

                                                                                           tGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCgg                                                           <  2:2866624/150‑1 (MQ=255)

                                                                                            gCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGt                                                          <  1:1958210/150‑1 (MQ=255)

                                                                                            gCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGt                                                          <  1:1958222/150‑1 (MQ=255)

                                                                                                   gAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTcc                                                   <  1:513195/150‑1 (MQ=255)

                                                                                                        cTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTAGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTAc                                              >  2:2303502/1‑150 (MQ=255)

                                                                                                         ttATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAG                                                                                                                                <  2:1545614‑M1/67‑1 (MQ=255)

                                                                                                           aTTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTcgaacctccgaccgctcggttcgtagccgagtactctatccagctgagctacggatgcaaatggcggtgaggcggggattc                                           <  2:1489635‑M1/150‑82 (MQ=255)

                                                                                                             tttacTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTcgaacctccgaccgctcggttcgtagccgagtactctatccagctgagctacggatgcaaatggcggtgaggcggggattcga                                         >  1:3215612‑M1/3‑67 (MQ=255)

                                                                                                              cTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTc                                        <  1:1222456/150‑1 (MQ=255)

                                                                                                               tACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTcgaacctccgaccgctcggttcgtagccgagtactctatccagctgagctacggatgcaaatggcggataggcggggataagaac                                       >  2:104717‑M1/1‑65 (MQ=255)

                                                                                                                  tgctgATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTcgaacctccgaccgctcggttcgtagccgagtactctatccagctgagctacggatgcaaatggcggtgaggcggggattcgaacccc                                    <  1:12931‑M1/150‑89 (MQ=255)

                                                                                                                     tgATATTGATACCGCTACTAAAGCCATAGCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCt                                 <  2:693773/150‑1 (MQ=255)

                                                                                                                      gATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTcgaacctccgaccgctcggttcgtagccgagtactctatccagctgagctacggatgcaaatggcggtgaggcggggattcgaaccccggat                                >  2:661666‑M1/1‑58 (MQ=255)

                                                                                                                      gATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTcgaacctccgaccgctcggttcgtagccgagtactctatccagctgagctacggatgcaaatggcggtgaggcggggattcgaaccccggat                                >  2:1249797‑M1/1‑58 (MQ=255)

                                                                                                                      gATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTcgaacctccgaccgctcggttcgtagccgagtactctatccagctgagctacggatgcaaatggaggagaggcggggattcgaacaactgat                                >  2:569015‑M1/1‑58 (MQ=255)

                                                                                                                         atTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCAt                             >  2:564745/1‑150 (MQ=255)

                                                                                                                             atcccGCTACTAAAGCCATATCTTTTAAGTGTTGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTATAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCt                         <  1:293766/147‑1 (MQ=17)

                                                                                                                             aTACCGCTACTAAAGCCATATCTATTAAGAGATTGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCt                         <  1:293773/150‑1 (MQ=255)

                                                                                                                             aTACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTcgaacctccgaccgctcggttcgtagccgagtactctatccagctgagctacggatgcaaatggcggtgaggcggggattcgaaccccggatgcagctt                         <  1:360213‑M1/150‑100 (MQ=255)

                                                                                                                              tACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTg                        <  1:2978014/150‑1 (MQ=255)

                                                                                                                              gaCCGCTTCTAAAGAAATATCTTTTATGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTg                        <  1:116732/149‑1 (MQ=17)

                                                                                                                                ccGTTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCAGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGaa                      >  2:1866474/1‑150 (MQ=255)

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

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

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

                                                                                                                                 cGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCAGGCGCCTCGCCCTACTGGCCGTTGCTAAATCAACGATATCATCCCTAGTGCATGCGGGTAATTCCGCTACCTCATCGATCCCTTCATTGCTGAAg                     >  1:1471238/1‑150 (MQ=17)

                                                                                                                                  gCTACTACAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTcgaacctccgaccgctcggttcgtagccgagtactctatccagctgagctacggatgcaaatggcggtgaggcggggattcgaaccccggatgcagcttttgac                    <  1:1249797‑M1/150‑105 (MQ=255)

                                                                                                                                  gCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGa                    >  2:351817/1‑150 (MQ=255)

                                                                                                                                  gCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGa                    <  2:871101/150‑1 (MQ=255)

                                                                                                                                  gCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTcgaacctccgaccgctcggttcgtagccgagtactctatccagctgagctacggatgcaaatggcggtgaggcggggattcgaaccccggatgcagcttttgac                    <  2:3215612‑M1/150‑105 (MQ=255)

                                                                                                                                  gCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTcgaacctccgaccgctcggttcgtagccgagtactctatccagctgagctacggatgcaaatggcggtgaggcggggattcgaaccccggatgcagcttttgac                    <  1:661666‑M1/150‑105 (MQ=255)

                                                                                                                                    tactacCGCCATCTCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGAtt                  <  1:1389427/150‑1 (MQ=255)

                                                                                                                                    tactaAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGAtt                  <  1:814417/150‑1 (MQ=255)

                                                                                                                                     actaAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGAttt                 <  2:2782180/150‑1 (MQ=255)

                                                                                                                                      ctaAAGCGATATCATTTTATAGATTGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTc                <  2:1012185/150‑1 (MQ=17)

                                                                                                                                       taAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCg                                                                                               <  1:1300953/70‑1 (MQ=40)

                                                                                                                                       taAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCg                                                                                               >  2:1300953/1‑70 (MQ=40)

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

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

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

                                                                                                                                           gCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGccc                                                                                                        <  1:1239394/57‑1 (MQ=255)

                                                                                                                                            ccATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCgg                                                                                                   <  1:1483563/61‑1 (MQ=40)

                                                                                                                                            ccATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCgg                                                                                                   >  2:1483563/1‑61 (MQ=40)

                                                                                                                                              atatCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAata        >  2:532978/1‑150 (MQ=255)

                                                                                                                                              atatCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAata        >  2:776710/1‑150 (MQ=255)

                                                                                                                                               tatCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAAc                                                                                  >  1:591793/1‑75 (MQ=34)

                                                                                                                                               tatCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAAc                                                                                  <  2:591793/75‑1 (MQ=34)

                                                                                                                                               tatCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAAc                                                                                  <  2:446394/75‑1 (MQ=34)

                                                                                                                                                atCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAATCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATg      <  1:2406714/150‑1 (MQ=255)

                                                                                                                                                 tCAAGTAAGAGATGGTGCATCCGGGAGGATTcgaacctccgaccgctcggttcgtagccgagtactctatccagctgagctacggatgcaaatggcggtgaggcggggattcgaaccccggatgcagcttttgaccgcatactcccttag     >  2:165432‑M1/1‑31 (MQ=255)

                                                                                                                                                 tCAAGTAAGAGATGGTGCATCCGGGAGGATTcgaacctccgaccgctcggttcgtagccgagtactctatccagctgagctacggatgcaaatggcggtgaggcggggagtcgaaccccggatgcagcttttgaccgcatactcccttag     >  2:2482584‑M1/1‑31 (MQ=255)

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

                                                                                                                                                     ttAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTaaa                                                                                       >  2:515191/1‑64 (MQ=31)

                                                                                                                                                     ttAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTaaa                                                                                       <  1:515191/64‑1 (MQ=31)

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