BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	847,215	A→T	9.9%	intergenic (‑96/+43)	glnP ← / ← glnH	glutamine transporter subunit/glutamine transporter subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	847,215	0	A	T	9.9%	709.8 / 55.4	223	intergenic (‑96/+43)	glnP/glnH	glutamine transporter subunit/glutamine transporter subunit
Reads supporting (aligned to +/- strand): ref base A (107/92); new base T (10/12); total (119/104)
Fisher's exact test for biased strand distribution p-value = 5.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.39e-01

Only 100 of 225 total aligned reads displayed.

ACCAGCGCGACGTGGTTGGCTATCCAACCTCCGAAGGTGCGTGCAAAACCTGCCAGCAATCCGATTACCAGACCGCCTGCCAGACCGAGGACCGAAATCCACAGGGTCATTTTGGCACCTTCAATCAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGC  >  NC_000913/846950‑847482

aCCAGCGCGACGTGGTTGGCTATCCAACCTCCGAAGGTGCGTGCAAAACCTGCCAGCAATCCGATTACCAGACCGCCTGCCAGACCGAGGACCGAAATCCACAGGGTCATTTTGGCACCTTCAATCAGAAGCGGAATGNCCGACCCGATGGCACTCCATGCAAACTGCATATGTTGTTCCGGTTACCGTGGTGTCCATATTCACATTAAGGGTAATTAANGGCCCCGGAAGGGGTTTGAGCCCCTNACCTCCCTCCTCCCCCCCCAGGGGGGAgg                                                                                                                                                                                                                                                                    >  3:62759/1‑275 (MQ=255)

                                                    ccAGCAATCCGATTACCAGACCGCCTGCCAGACCGAGGACCGAAATCCACAGGGTCATTTTGGCACCTTCAATCAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGCGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGcc                                                                                                                                                                                                                                                            >  4:142431/1‑231 (MQ=255)

                                                        cAATCCGATTACCAGACCGCCTGCCAGACCGAGGACCGAAATCCACAGGGTCATTTTGGCACCTTCAATCAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAAGAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTg                                                                                                                                                                                                                                                        <  2:142431/231‑1 (MQ=255)

                                                                                                                       ttCAATCAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAgg                                                                                                                                                                                                                                                                           >  5:1470443/1‑149 (MQ=255)

                                                                                                                        tCAATCAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGttatt                                                                                                                                                                                                                                 >  4:373636/1‑190 (MQ=255)

                                                                                                                           aTCAGAATCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGgag                                                                                                                                                                                                                                                                       >  6:982918/1‑149 (MQ=255)

                                                                                                                           aTCAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCTATGGGgag                                                                                                                                                                                                                                                                       >  7:616379/1‑149 (MQ=255)

                                                                                                                            tCAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGgaga                                                                                                                                                                                                                                                                      <  5:1255065/149‑1 (MQ=255)

                                                                                                                            tCAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTAttt                                                                                                                                                                                                                             <  2:373637/190‑1 (MQ=255)

                                                                                                                            tCAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAACCCTCTCCCCTATGGGgaga                                                                                                                                                                                                                                                                      <  7:1128247/149‑1 (MQ=255)

                                                                                                                             cAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGagag                                                                                                                                                                                                                                                                     <  7:923310/149‑1 (MQ=255)

                                                                                                                             cAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGagag                                                                                                                                                                                                                                                                     <  8:174044/149‑1 (MQ=255)

                                                                                                                              aGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCTATGGGGAGAgg                                                                                                                                                                                                                                                                    <  7:868090/149‑1 (MQ=255)

                                                                                                                              aGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAgg                                                                                                                                                                                                                                                                    <  8:732893/149‑1 (MQ=255)

                                                                                                                               gAAGCGGAATGGCAGGCCATATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGTTGTAAATAATCAAATTCAGGGTAAGTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAAACCTCTCCCCGATGGGGAGAgga                                                                                                                                                                                                                                                                   <  8:1124181/149‑2 (MQ=255)

                                                                                                                                  gCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCTATGGGGAGAggatta                                                                                                                                                                                                                                                                <  6:1084817/149‑5 (MQ=255)

                                                                                                                                    ggAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAACCCTCTCCCCTATGGGGAGAggattagg                                                                                                                                                                                                                                                              <  5:67215/149‑7 (MQ=255)

                                                                                                                                      aaTGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGcc                                                                                                                                                                                                                                                            >  7:780689/1‑149 (MQ=255)

                                                                                                                                       aTGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCa                                                                                                                                                                                                                                                           <  6:936824/149‑1 (MQ=255)

                                                                                                                                            aGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATCGGGGAGAGGGAACTGCCAGTGCg                                                                                                                                                                                                                                                      >  7:257803/1‑149 (MQ=255)

                                                                                                                                            aGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCg                                                                                                                                                                                                                                                      >  6:272973/1‑149 (MQ=255)

                                                                                                                                              gCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGtt                                                                                                                                                                                                                                                    >  7:629903/1‑149 (MQ=255)

                                                                                                                                               ccAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGttt                                                                                                                                                                                                                                                   <  6:861473/149‑1 (MQ=255)

                                                                                                                                                cAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTGATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAATTTCCGTTCTCGCGCAGGGTTTTCAACGCTCCGTTGACTTTGTCACGCAGCTCGTCGCTACCtt                                                                                                                     >  2:95379/1‑274 (MQ=255)

                                                                                                                                                   aTGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACa                                                                                                                                                                                                                                               <  5:459519/149‑1 (MQ=255)

                                                                                                                                                    tGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAg                                                                                                                                                                                                                                              <  5:832645/149‑1 (MQ=255)

                                                                                                                                                     ggCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGTGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAgg                                                                                                                                                                                                                                             >  8:1377191/1‑149 (MQ=255)

                                                                                                                                                          tCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAg                                                                                                                                                                                                                                        >  6:427139/1‑149 (MQ=255)

                                                                                                                                                             aGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGt                                                                                                                                                                                                                                     >  8:856511/1‑149 (MQ=255)

                                                                                                                                                                 aaaCTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCTCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGttatt                                                                                                                                                                                                                                 >  8:790341/1‑149 (MQ=255)

                                                                                                                                                                 aaaCTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTGCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGttatt                                                                                                                                                                                                                                 >  6:894731/1‑149 (MQ=255)

                                                                                                                                                                     tGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTAttt                                                                                                                                                                                                                             >  7:808332/1‑149 (MQ=255)

                                                                                                                                                                      gCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTc                                                                                                                                                                                                                            <  8:983652/149‑1 (MQ=255)

                                                                                                                                                                      gCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTc                                                                                                                                                                                                                            <  8:313216/149‑1 (MQ=255)

                                                                                                                                                                      gCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTc                                                                                                                                                                                                                            <  5:414697/149‑1 (MQ=255)

                                                                                                                                                                      gCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGATTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTc                                                                                                                                                                                                                            <  7:98418/149‑1 (MQ=255)

                                                                                                                                                                       cATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCCTAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCCGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTTTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCTTTGACTTTGTCCCGCAGCCCGTCGCTACCTTTCCGGTACGCAATACCGTATtgct                                                                                             >  2:310615/1‑275 (MQ=255)

                                                                                                                                                                        atatGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGaa                                                                                                                                                                                                                                                                 <  8:1088241/110‑1 (MQ=255)

                                                                                                                                                                        atatGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGaa                                                                                                                                                                                                                                                                 >  7:1088241/1‑110 (MQ=255)

                                                                                                                                                                        atatGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCgg                                                                                                                                                                                                                          >  8:368483/1‑149 (MQ=255)

                                                                                                                                                                          gttgtgttACTGTTAACGTTGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGGTTGCGCCCCTCACCCTAATACTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTAGTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTgcg                                                                                          <  1:62759/271‑1 (MQ=255)

                                                                                                                                                                          atGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGgt                                                                                                                                                                                                                                            >  5:599111/1‑129 (MQ=255)

                                                                                                                                                                          atGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGgt                                                                                                                                                                                                                                            <  6:599111/129‑1 (MQ=255)

                                                                                                                                                                             ttgttCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAg                                                                                                                                                                                                                     >  5:365972/1‑149 (MQ=255)

                                                                                                                                                                                ttCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTAc                                                                                                                                                                                                                  >  8:6756/1‑149 (MQ=255)

                                                                                                                                                                                ttCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTAc                                                                                                                                                                                                                  >  8:5821/1‑149 (MQ=255)

                                                                                                                                                                                ttCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTAc                                                                                                                                                                                                                  >  6:530240/1‑149 (MQ=255)

                                                                                                                                                                                  ccTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCg                                                                                                                                                                                                                <  6:296709/149‑1 (MQ=255)

                                                                                                                                                                                   cTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGt                                                                                                                                                                                                                                                         >  6:768700/1‑107 (MQ=255)

                                                                                                                                                                                    tGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGaa                                                                                                                                                                                                              >  7:486852/1‑149 (MQ=255)

                                                                                                                                                                                    tGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGaa                                                                                                                                                                                                              >  5:895618/1‑149 (MQ=255)

                                                                                                                                                                                          cGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCAttt                                                                                                                                                                                                        >  5:1333570/1‑149 (MQ=255)

                                                                                                                                                                                              gtgtAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTg                                                                                                                                                                                                    >  6:1017113/1‑149 (MQ=255)

                                                                                                                                                                                              gtgtAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTtc                                                                                                                                                                            >  4:190910/1‑173 (MQ=255)

                                                                                                                                                                                               tgtAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGt                                                                                                                                                                                                   >  5:841640/1‑149 (MQ=255)

                                                                                                                                                                                               tgtAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCGCCCCATAGGGGAGAGGGACCTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTACGGTTCAGTACCGAACCATTCTTTGt                                                                                                                                                                                                   <  8:927669/149‑1 (MQ=255)

                                                                                                                                                                                                gtAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTa                                                                                                                                                                                                  <  6:365972/149‑1 (MQ=255)

                                                                                                                                                                                                 tAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGGGCGTTTTACAGTTTTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAg                                                                                                                                                                                                 >  6:223747/1‑149 (MQ=255)

                                                                                                                                                                                                  aaaTAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGa                                                                                                                                                                                                <  8:822093/149‑1 (MQ=255)

                                                                                                                                                                                                       atCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCg                                                                                                                                                                                           >  7:1164411/1‑149 (MQ=255)

                                                                                                                                                                                                       atCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCg                                                                                                                                                                                           >  7:387531/1‑149 (MQ=255)

                                                                                                                                                                                                        tcgaaTTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCTTAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGt                                                                                                                                                                                          <  6:1260446/146‑1 (MQ=255)

                                                                                                                                                                                                        tCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGt                                                                                                                                                                                          >  8:1038962/1‑149 (MQ=255)

                                                                                                                                                                                                        tCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGt                                                                                                                                                                                          <  5:272973/149‑1 (MQ=255)

                                                                                                                                                                                                         cAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGtt                                                                                                                                                                                         <  8:241002/149‑1 (MQ=255)

                                                                                                                                                                                                          aaaTTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTg                                                                                                                                                                                        >  7:50385/1‑149 (MQ=255)

                                                                                                                                                                                                           aaTTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGt                                                                                                                                                                                       <  7:1345257/149‑1 (MQ=255)

                                                                                                                                                                                                             ttCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTaa                                                                                                                                                                                     <  7:101590/149‑1 (MQ=255)

                                                                                                                                                                                                                 gggTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTc                                                                                                                                                                                 <  7:1442098/149‑1 (MQ=255)

                                                                                                                                                                                                                  ggTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAAcgcg                                                                                                                                                       >  4:183649/1‑174 (MQ=255)

                                                                                                                                                                                                                         aatGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTcg                                                                                                                                                                         >  5:214010/1‑149 (MQ=255)

                                                                                                                                                                                                                            ggCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCgcgc                                                                                                                                                                      >  5:1263507/1‑149 (MQ=255)

                                                                                                                                                                                                                                ccGGGGCGGGTTTGCTCCCCTCCCCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTtgatg                                     <  4:95379/274‑1 (MQ=255)

                                                                                                                                                                                                                                 cGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTtgatga                                    <  4:252922/274‑1 (MQ=255)

                                                                                                                                                                                                                                    aGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGtttt                                                                                                                                                              <  7:1185967/149‑1 (MQ=255)

                                                                                                                                                                                                                                       gggTTTGCGCCGCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCaa                                                                                                                                                           >  5:773285/1‑149 (MQ=255)

                                                                                                                                                                                                                                        ggTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATCTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAAc                                                                                                                                                          <  5:1065753/149‑1 (MQ=255)

                                                                                                                                                                                                                                             gcgcCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGcc                                                                                                                                                     >  7:716991/1‑149 (MQ=255)

                                                                                                                                                                                                                                               gcCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGt                                                                                                                                                   <  6:773285/149‑1 (MQ=255)

                                                                                                                                                                                                                                                 cccTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTg                                                                                                                                                 >  6:1473661/1‑149 (MQ=255)

                                                                                                                                                                                                                                                 cccTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTg                                                                                                                                                 >  7:226885/1‑149 (MQ=255)

                                                                                                                                                                                                                                                 cccTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTg                                                                                                                                                 >  7:1384355/1‑149 (MQ=255)

                                                                                                                                                                                                                                                 cccTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGGTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACTCGCCGTTg                                                                                                                                                 >  7:1284514/1‑149 (MQ=255)

                                                                                                                                                                                                                                                   cTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGAc                                                                                                                                               >  5:920750/1‑149 (MQ=255)

                                                                                                                                                                                                                                                   cTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGTTGAAGTCAAGAATGTTTGGAGt                 >  2:246323/1‑275 (MQ=255)

                                                                                                                                                                                                                                                   cTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTTCAGAATGTTTGGCg                  >  2:163628/1‑274 (MQ=255)

                                                                                                                                                                                                                                                         cTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGAttt                                                                                                                                                                                             >  5:763498/1‑97 (MQ=255)

                                                                                                                                                                                                                                                           ttcccTCTCCCCTATGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCAc                                                                                                                                       >  5:1032777/4‑149 (MQ=255)

                                                                                                                                                                                                                                                            aTCCTCT‑CCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCAtttt                                                                                                                                                                                                       <  7:680222/83‑1 (MQ=255)

                                                                                                                                                                                                                                                              ccTCTCCCCTATGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGt                                                                                                                                                                                          <  8:1154338/95‑1 (MQ=255)

                                                                                                                                                                                                                                                              ccTCTCCCCTATGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAAc                                                                                                                                                          <  7:998617/127‑1 (MQ=255)

                                                                                                                                                                                                                                                              ccTCTCCCCTATGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCa                                                                                                                                    >  7:575112/1‑149 (MQ=255)

                                                                                                                                                                                                                                                               ctctCCCCTATGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCa                                                                                                                                        <  5:532349/144‑1 (MQ=255)

                                                                                                                                                                                                                                                               ctctCCCCTATGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAg                                                                                                                                   >  8:9246/1‑149 (MQ=255)

                                                                                                                                                                                                                                                               ctctCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAg                                                                                                                                   >  5:283165/1‑149 (MQ=255)

                                                                                                                                                                                                                                                               ctctCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAg                                                                                                                                   >  8:951020/1‑149 (MQ=255)

                                                                                                                                                                                                                                                               ctctCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAg                                                                                                                                   >  8:46209/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                    cccATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCtcgt                                                                                                                              >  6:133132/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                      cATAGGGGAGAGGGAACTGCCAGTTCTTTTTACAGGTGTAGCGTTATTATTGCGGTTCAGTACCGAACCATTTTTTGTAGATGGCGTTGTAAGTTCCGTTCGCGCGCAGGGTTTTCAACGCGCCGTTGACTTTCGCACGCAGCTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTTCAGCGAGTCACCTa                                                                        <  1:2463/201‑1 (MQ=255)

                                                                                                                                                                                                                                                                      cATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCtcgtcg                                                                                                                            >  6:792136/1‑149 (MQ=255)

ACCAGCGCGACGTGGTTGGCTATCCAACCTCCGAAGGTGCGTGCAAAACCTGCCAGCAATCCGATTACCAGACCGCCTGCCAGACCGAGGACCGAAATCCACAGGGTCATTTTGGCACCTTCAATCAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCTCGTCGCTACCTTTCGGGAACGCAATACCGTATTGCTGCGCTTCCAGAGAGTCACCTACCGCTTTGAACTGACCGTTACCGGCGGTTTTGATGAAGTACAGAATGTTTGGCGTATCGTGCAGAACGGC  >  NC_000913/846950‑847482

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