BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,817,916	T→G	8.0%	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	8.0%	761.9 / 50.9	247	intergenic (‑56/+143)	argQ/argZ	tRNA‑Arg/tRNA‑Arg
Reads supporting (aligned to +/- strand): ref base T (104/123); new base G (12/8); total (116/131)
Fisher's exact test for biased strand distribution p-value = 2.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.71e-01

Only 100 of 258 total aligned reads displayed.

GAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCA  >  NC_000913/2817767‑2818065

gAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAt                                                                                                                                                       <  1:1118859/150‑1 (MQ=255)

gAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAt                                                                                                                                                       <  1:826911/150‑1 (MQ=255)

          gagaATATGATGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGG                                                                                                                                             <  2:3585450‑M1/150‑1 (MQ=255)

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

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

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

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

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

                     gCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGA                                                                                                                                  <  1:1268537‑M1/150‑1 (MQ=255)

                     gCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGA                                                                                                                                  <  1:2111232‑M1/150‑1 (MQ=255)

                       aTCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGG                                                                                                                                >  2:2184886‑M1/1‑150 (MQ=255)

                       aTCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGG                                                                                                                                >  2:4698964‑M1/1‑150 (MQ=255)

                           gggAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAt                                                                                                                                                       >  1:2470106/1‑123 (MQ=255)

                           gggAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTa                                                                                                                            >  1:1870273/1‑150 (MQ=255)

                             gaggaTTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACt                                                                                                                          >  1:3726905/1‑150 (MQ=255)

                                     gcaCCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGTCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGc                                                                                                                  <  1:2882592/148‑1 (MQ=255)

                                     gAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGc                                                                                                                  <  1:3063837/150‑1 (MQ=255)

                                     gAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACACGccaggc                                                                                                                  >  1:2512586/1‑146 (MQ=255)

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

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

                                         cTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTCTTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCt                                                                                                              <  2:1951990/150‑1 (MQ=255)

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

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

                                                    cGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACCGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCggg                                                                                                   <  1:3991749/150‑1 (MQ=255)

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

                                                          gTAGCCGAGTACTCTATCCAGCTGAGCTAAGGATTCATCGGGAAAACTTATTCTACTTCTGATATTGATACCGCTACTAAAGCCATATCAATTTAGAGATGGTGCATACGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTg                                                                                             <  1:2369821/150‑1 (MQ=17)

                                                           tAGCCGAGTACTCTATCCAGCTGAGCTACGGATTCATCGGGAATACTTATTCTACTGCTGATATTGATACAGCTACTAAAGCCATATATATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGc                                                                                            <  1:866572/150‑1 (MQ=17)

                                                           tAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGc                                                                                            >  2:3344751/1‑150 (MQ=255)

                                                           gaGCCGAGTACTCTATCCAGCTGCGCTACGGATGCATCGGGAAAACTTATTCTACTGCTTATATTGATACCGATACTTTAGCCATATCAATTAAGCGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGc                                                                                            <  1:1505857/149‑1 (MQ=17)

                                                              ccGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTaa                                                                                         >  1:2765184/1‑150 (MQ=255)

                                                               cGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCGACTGCTGATATTGAGACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTaaa                                                                                        <  1:1812912/150‑1 (MQ=255)

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

                                                                       ctATCCAGCTGAGCTACGGATGCATCGGGAAAAATTATTCTACTGCTGATATTGATACCGCTACTAAACCCATAGCAATTAAGAGATGGTGAATCCGGGAGGATTCCTAGCCGTGCGGCTCGCGCATCGGGCGGTTGCTAAAGCAACGta                                                                                >  2:923876/1‑149 (MQ=17)

                                                                        tATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACACCCACGGACGCTCGCCCTGCGGGACGTGGCGGGAGCAAGGttt                                                                               >  2:1868610/1‑149 (MQ=14)

                                                                             aGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTaaa                                                                                        <  2:3546471/136‑1 (MQ=255)

                                                                             aGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTaaa                                                                                        >  1:3546471/1‑136 (MQ=255)

                                                                              gCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTAGCCTAGCGGCTAGCCCATAGGGACGTTGCTAAAGCAACGTTATCATcc                                                                         >  2:4807025/1‑150 (MQ=18)

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

                                                                                 gAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACACGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTa                                                                      >  2:1445762/1‑150 (MQ=255)

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

                                                                                         gATGCATCGGGAATACTTATTCTACTGCTGATATTGAGACCGCTACTAAAGCCATATCTTTTATGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGc                                                              <  1:2567125/150‑1 (MQ=18)

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

                                                                                         gATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTcgaacctccgaccgctcggttcgtagccgagtactctatccagctgagctacggatgcatcgg                                                              <  2:1746500‑M1/150‑64 (MQ=255)

                                                                                           tGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCgg                                                            >  1:509890/1‑150 (MQ=255)

                                                                                           tGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCgg                                                            >  1:1690882/1‑150 (MQ=255)

                                                                                           tGCATCGGGAAAACTTAGTCGACTGCTGATAATGATACCGCTACTAACGCCATATCAATTAAGAGATGGTGCATCCGGGTGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCgg                                                            <  2:290943/150‑1 (MQ=18)

                                                                                            gCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCTATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGt                                                           <  2:1270349/150‑1 (MQ=255)

                                                                                            gCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGt                                                           >  2:2144712/1‑150 (MQ=255)

                                                                                             cATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGtt                                                          >  2:3074328/1‑150 (MQ=255)

                                                                                             cATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGtt                                                          >  2:845780/1‑150 (MQ=255)

                                                                                             cATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGtt                                                          >  2:1036440/1‑150 (MQ=255)

                                                                                             cATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGtt                                                          <  1:1401795/150‑1 (MQ=255)

                                                                                             cATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGtt                                                          <  1:4068143/150‑1 (MQ=255)

                                                                                               tCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTaa                                                        <  1:703870/150‑1 (MQ=255)

                                                                                               tCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTaa                                                        <  1:418752/150‑1 (MQ=255)

                                                                                                cGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAAt                                                       <  2:2028234/150‑1 (MQ=255)

                                                                                                  ggAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTc                                                     <  2:118068/150‑1 (MQ=255)

                                                                                                    aaaaCTTATTCTACTTCTGATATTGATACCTATACTAAAGCCATATCAATTAAGAGATGTTGCATACGGGAGTATTACTCGCCTTGAGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCg                                                   <  2:2148596/150‑1 (MQ=18)

                                                                                                    aaaaCTTATTCTACTGCTGATATTGATACCGCTACTACAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCg                                                   <  1:4672457/150‑1 (MQ=255)

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

                                                                                                        tttaTTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTcgaacctccgaccgctcggttcgtagccgagtactctatccagcagagctacggatgcaaatggaggagaggcgggga                                               >  2:2491483‑M1/2‑72 (MQ=255)

                                                                                                        cTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTAc                                               >  1:3599575/1‑150 (MQ=255)

                                                                                                         ttATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTAcc                                              >  1:2266539/1‑150 (MQ=255)

                                                                                                          tATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACct                                             >  2:1631225/1‑150 (MQ=255)

                                                                                                            ttCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACctct                                           >  1:4087739/1‑150 (MQ=255)

                                                                                                             tttacTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTcgaacctccgaccgctcggttcgtagccgagtactctatccagctgagctacggatgcaaatggcggtgaggcggggattcga                                          <  2:388928‑M1/148‑84 (MQ=255)

                                                                                                              cTACTGCTGATATTGATACCGCTACTACAGCCATATCAATTACGTGATGGTGCCTCCGGGAGGAGTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTc                                         <  2:4449733/150‑1 (MQ=17)

                                                                                                                aCTGCTTATTTTTATACCGCTACTTATGCGATATCAAGTATGATATGGTGCATCCGGGAGGATTcgaacctccgaccgctcggttcgtagccgagtactctatccagctgagctacggatgcaaatggcggtgaggcggggattcgaacc                                       <  2:2626469‑M1/150‑87 (MQ=255)

                                                                                                                    ctgATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCATAGCGGCTCGCCCTTCGGGCCGTTGCTAAAGAAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATaac                                   >  2:3854439/1‑148 (MQ=255)

                                                                                                                     tgATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCt                                  <  2:3726905/150‑1 (MQ=255)

                                                                                                                         atTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCATCTACCTCTTCGAAACCTTCAt                              >  1:2984802/1‑150 (MQ=255)

                                                                                                                          ttGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACATCTTCGATCCCTTCAtt                             >  1:3921686/1‑150 (MQ=255)

                                                                                                                            gagaCCGCTACCAAACCAATATCTCTTGAGAGTTGGTGCTTCGTGTTGTATTACTATCCCTGCGGCTCGCCCTTCGGGCCGTTGCTTTAGCAACGTTATCATCCCTTGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGc                           <  1:3806733/147‑1 (MQ=12)

                                                                                                                                ccGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGaa                       >  1:3777868/1‑150 (MQ=255)

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

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

                                                                                                                                 cGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTAGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGTAg                      >  2:2416879/1‑150 (MQ=255)

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

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

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

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

                                                                                                                                   ctactaAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTcgaacctccgaccgctcggttcgtagccgtgtactctatccagctgagctacggatgcaaatggcggtgaggcggggattcgaaccccggatgcagcttatgacc                    >  2:3680670‑M1/1‑45 (MQ=255)

                                                                                                                                   ctactaAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTcgaacctccgaccgctcggttcgtagccgagtactctatccagctgagctacggatgcaaatggcggtgaggcggg                                                 <  1:4717698‑M1/121‑77 (MQ=255)

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

                                                                                                                                      ctaAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTc                 <  1:2184886/150‑1 (MQ=255)

                                                                                                                                      ctaAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTc                 <  1:1286773/150‑1 (MQ=255)

                                                                                                                                       taAAGCCATATCAATTTAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTTAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCg                <  2:3120887/150‑1 (MQ=255)

                                                                                                                                       taAAGCCATATCAATTAAGAGATGGTGCCTCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTTTGGCCGTTGCTAACGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCg                <  2:467092/150‑1 (MQ=255)

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

                                                                                                                                         aaGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTcgaacctccgaccgctcggttcgtagccgagtactctatccagctgagctacggatgcaaatggaggggaggcggggataagaacaaatgatgcagctattgaacccataa              >  2:286688‑M1/1‑39 (MQ=255)

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

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

                                                                                                                                                atCACTTAAGAGAGGGTGCCTCCGGGAGGAGTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCACCGTTATCCTCCCTAGTGCGTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATg       <  2:223648/150‑1 (MQ=17)

                                                                                                                                                atCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTACAGCAACGTTCTCATCCCTAGTGCCTGCGGTTAATTCCGCTACCACTTCGATGCCTTCGGTGCTGAAGATTTCGGAGAATATg       >  2:3374800/1‑150 (MQ=18)

                                                                                                                                                atCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTAGGAGAATATg       >  1:4377118/1‑150 (MQ=255)

                                                                                                                                                 tCAATTAAGTGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCt                                                                                                        <  2:3757837/52‑1 (MQ=255)

                                                                                                                                                 tCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCt                                                                                                        <  2:3951737/52‑1 (MQ=35)

                                                                                                                                                 tCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCt                                                                                                        >  1:3757837/1‑52 (MQ=35)

                                                                                                                                                 tCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCt                                                                                                        >  1:3951737/1‑52 (MQ=35)

                                                                                                                                                     ttAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAAGATGGTGCa  >  2:1556429/1‑150 (MQ=255)

GAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCA  >  NC_000913/2817767‑2818065

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