BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	481,228	C→A	8.9%	intergenic (‑520/+26)	tomB ← / ← acrB	Hha toxicity attenuator; conjugation‑related protein/multidrug efflux system protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	481,228	0	C	A	8.9%	451.8 / 24.6	135	intergenic (‑520/+26)	tomB/acrB	Hha toxicity attenuator; conjugation‑related protein/multidrug efflux system protein
Reads supporting (aligned to +/- strand): ref base C (67/56); new base A (8/4); total (75/60)
Fisher's exact test for biased strand distribution p-value = 5.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.40e-01

Only 100 of 137 total aligned reads displayed.

GCCAGATGAAATTGGTGACACTAATACCAGGATTGCTCTGAATATGACGTAATAACCGAGGAATGAATAAAGAATTACCGCAAATAATTAAGAATAGCCTCTAAATGATTATGGATTATTAAATCTATTAGCGCTACTTAATATAAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGCGTTCTGCGCGCCGGAACCAGCACCAGTACTGATAACCAGCGGCATAACGCCGAGGATAAACGCCAGCGAGGTCATCAGGAT  >  NC_000913/480955‑481481

gCCAGATGAAATTGGTGACACTAATACCAGGATTGCTCTGAATATGACGTAATAACCGAGGAATGAATAAAGAACTACCGCAAATAATTAAGAATAGCCTCTAAATGATTATGGATTATTAAATCTATTAGCGCTACTTAATATAAATTAACTAACATTCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTAcg                                                                                                                                                                                                                                                              <  2:248605/275‑1 (MQ=255)

                                   ctctGAATATGACGTAATAACCGAGGAATGAATAAAGAATTACCGCAAATAATTAAGAATAGCCTCTAAATGATTATGGATTATTAAATCTATTAGCGCTACTTAATATAAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTTCTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACTGTTTATCAAtgatgaa                                                                                                                                                                                                                           >  4:16440/1‑274 (MQ=255)

                                                                           tACCGCAAATAATTAAGAATAGCCTCTAAATGATTATGGATTATTAAATCTATTAGCGCTACTTAATATAAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTtgctgc                                                                                                                                                                                    >  2:286230/1‑274 (MQ=255)

                                                                                                           aTTATGGATTATTAAATCTATTAGCGCTACTTAATATAAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCt                                                                                                                                                                                   >  4:59810/1‑243 (MQ=255)

                                                                                                               tGGATTATTAAATCTATTAGCGCTACTTAATATAAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAg                                                                                                                                                                               <  2:59810/243‑1 (MQ=255)

                                                                                                                              aTTAGCGCTACTTAATATAAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCCTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTAcg                                                                                                                                                                                                                                                              >  8:1276866/1‑149 (MQ=255)

                                                                                                                                    gcTACTTAATATAAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCt                                                                                                                                                                                                                                                        <  7:406975/149‑1 (MQ=255)

                                                                                                                                    gcTACTTAATATAAATTAACTAACAATCAGTATTTTTTATAAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCt                                                                                                                                                                                                                                                        <  7:1140892/149‑1 (MQ=255)

                                                                                                                                      tACTTAATATAAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTa                                                                                                                                                                                                                                                      <  8:733550/149‑1 (MQ=255)

                                                                                                                                       aCTTAATATAAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTAGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAg                                                                                                                                                                                                                                                     <  8:574740/149‑1 (MQ=255)

                                                                                                                                        cTTAATATAAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGt                                                                                                                                                                                                                                                    >  6:1268345/1‑149 (MQ=255)

                                                                                                                                        cTTAATATAAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGt                                                                                                                                                                                                                                                    <  8:1179808/149‑1 (MQ=255)

                                                                                                                                        cTTAATATAAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGt                                                                                                                                                                                                                                                    <  7:1269121/149‑1 (MQ=255)

                                                                                                                                         ttAATATAAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTg                                                                                                                                                                                                                                                   >  6:780613/1‑149 (MQ=255)

                                                                                                                                         ttAATATAAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCATAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTg                                                                                                                                                                                                                                                   >  7:787926/1‑149 (MQ=255)

                                                                                                                                           aaTATAAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGAt                                                                                                                                                                                                                                                 <  5:229168/149‑1 (MQ=255)

                                                                                                                                           aaTATAAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGAt                                                                                                                                                                                                                                                 <  5:273415/149‑1 (MQ=255)

                                                                                                                                             tataAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTa                                                                                                                                                                                                                                               >  7:1213737/1‑149 (MQ=255)

                                                                                                                                              ataAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTac                                                                                                                                                                                                                                              >  5:191621/1‑149 (MQ=255)

                                                                                                                                               taAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTaca                                                                                                                                                                                                                                             <  8:242795/149‑1 (MQ=255)

                                                                                                                                                  aTTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGCTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGt                                                                                                                                                                                                                                          <  6:895622/149‑1 (MQ=255)

                                                                                                                                                    taactaacAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCTGCTTACGCGGCCTTAGTGATTACACGtag                                                                                                                                                                                                                                        <  7:931398/149‑3 (MQ=255)

                                                                                                                                                     aactaacAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGt                                                                                                                                                                                                                                       <  8:1008388/149‑1 (MQ=255)

                                                                                                                                                       ctaacAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGa                                                                                                                                                                                                                                                  >  7:1419127/1‑136 (MQ=255)

                                                                                                                                                        tatcAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTGGTATGAGATCCTGAGTTGGTGGTTCAATTACGCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGtt                                                                                                                                                                                                                                         <  5:7637/141‑1 (MQ=255)

                                                                                                                                                        taacAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCGTAAGCGGCCt                                                                                                                                                                                                                                                        <  6:319009/129‑1 (MQ=255)

                                                                                                                                                          acAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCaa                                                                                                                                                                                                                                  <  7:968108/149‑1 (MQ=255)

                                                                                                                                                          acAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTCAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCaa                                                                                                                                                                                                                                  >  7:963948/1‑149 (MQ=255)

                                                                                                                                                            aaTCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGAt                                                                                                                                                                                                                                                 >  8:1116389/1‑132 (MQ=255)

                                                                                                                                                                 gTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATgatgat                                                                                                                                                                                                                           <  7:1443862/149‑1 (MQ=255)

                                                                                                                                                                  tATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATc                                                                                                                                                                                                                          >  5:680350/1‑149 (MQ=255)

                                                                                                                                                                    ttttttATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGa                                                                                                                                                                                                                        <  6:680350/149‑1 (MQ=255)

                                                                                                                                                                         tATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTa                                                                                                                                                                                                                   <  5:890841/149‑1 (MQ=255)

                                                                                                                                                                            gAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGc                                                                                   >  2:117992/1‑274 (MQ=255)

                                                                                                                                                                             aaTTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGc                                                                                                                                                                                                               >  7:1359762/1‑149 (MQ=255)

                                                                                                                                                                               ttttATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGTGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCtg                                                                                                                                                                                                             >  8:479660/1‑149 (MQ=255)

                                                                                                                                                                               ttttATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGAt                                                                                                                                                                                                                                                 >  8:885878/1‑113 (MQ=255)

                                                                                                                                                                               ttttATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCtg                                                                                                                                                                                                             >  5:599070/1‑149 (MQ=255)

                                                                                                                                                                               ttttATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCtg                                                                                                                                                                                                             >  7:376144/1‑149 (MQ=255)

                                                                                                                                                                                      cGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGa                                                                                                                                                                                                      <  7:1276866/149‑1 (MQ=255)

                                                                                                                                                                                      cGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGa                                                                                                                                                                                                      <  8:1359762/149‑1 (MQ=255)

                                                                                                                                                                                          gTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATc                                                                                                                                                                                                  <  7:642673/149‑1 (MQ=255)

                                                                                                                                                                                           ttAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCt                                                                                                                                                                                                 <  8:251092/149‑1 (MQ=255)

                                                                                                                                                                                             aaTACTGGTTTTTGTTGGGGATCCTGAGTGGGTGGTTCAATTCTCCCTTAATGTTCGTAGGTTTTGCATAAAAAAGGCCGCTTACGCTGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGCGTTCTGCGCGCCGGAAcc                                                                 <  4:286230/275‑1 (MQ=255)

                                                                                                                                                                                             aaTACTGGTTTTGGTATCAGATCCTGAGTTGGTGGTTCAATTACTCTTTAAAGTTCGTAGGTTATTCATAAAAAAGGCCGCTTACGTGGCCTTAGTGATTTCACGTTGTATCAATGCTGATCTACAGTATGGCTGTGCTCGATATCCTc                                                                                                                                                                                               <  8:787926/149‑1 (MQ=255)

                                                                                                                                                                                              aTACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCa                                                                                                                                                                                              <  8:1042579/149‑1 (MQ=255)

                                                                                                                                                                                                    gTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTg                                                                                                                                                                                        >  8:878826/1‑149 (MQ=255)

                                                                                                                                                                                                    gTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATAc                                                                                                                                                     >  4:305207/1‑184 (MQ=255)

                                                                                                                                                                                                    gTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGCCACCATCCCGCCCATTACCCCGGTACCTACTGACTTCTGCGCGCCGGAACCAGCAccg                                                          >  4:305022/1‑274 (MQ=255)

                                                                                                                                                                                                           tATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTaa                                                                                                                                                                                 >  7:1036124/1‑149 (MQ=255)

                                                                                                                                                                                                             tGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAg                                                                                                                                                                               <  8:1036124/149‑1 (MQ=255)

                                                                                                                                                                                                                 aTCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATTGCTGTGCTCGATATCTTCATTCTTGc                                                                                                                                                                                       >  5:174869/1‑137 (MQ=255)

                                                                                                                                                                                                                 aTCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGc                                                                                                                                                                                       <  6:174869/137‑1 (MQ=255)

                                                                                                                                                                                                                 aTCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAgcggc                                                                                                                                                                           >  8:334821/1‑149 (MQ=255)

                                                                                                                                                                                                                        gTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAac                                                                                                                                                                    >  6:1010165/1‑149 (MQ=255)

                                                                                                                                                                                                                            gtggtTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAaccacc                                                                                                                                                                <  7:1434850/149‑1 (MQ=255)

                                                                                                                                                                                                                                 tCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGaa                                                                                                                                                                     <  6:592993/139‑1 (MQ=255)

                                                                                                                                                                                                                                 tCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGaa                                                                                                                                                                     >  5:592993/1‑139 (MQ=255)

                                                                                                                                                                                                                                 tCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACaaa                                                                                                                                                           >  6:543751/1‑149 (MQ=255)

                                                                                                                                                                                                                                     ttACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGTGAACCACCACAAAGAAt                                                                                                                                                       >  8:1194480/1‑149 (MQ=255)

                                                                                                                                                                                                                                       aCTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATAc                                                                                                                                                     <  6:599070/149‑1 (MQ=255)

                                                                                                                                                                                                                                         tCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATAcc                                                                                                                                                    <  6:242491/148‑1 (MQ=255)

                                                                                                                                                                                                                                         tCCTTAATGTTCGTAGGTTACGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCg                                                                                                                                                   <  6:376783/149‑1 (MQ=255)

                                                                                                                                                                                                                                           cTTAGTGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTTTCAATGATGATCGACAGTATGGCTGTCCTCGATATCTTCATTCTTg                                                                                                                                                                                        <  6:942398/110‑1 (MQ=255)

                                                                                                                                                                                                                                           cTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATc                                                                                                                                                                                                  >  7:496938/1‑100 (MQ=255)

                                                                                                                                                                                                                                           cTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATc                                                                                                                                                                                                  <  8:496938/100‑1 (MQ=255)

                                                                                                                                                                                                                                           cTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGa                                                                                                                                                 >  7:306661/1‑149 (MQ=255)

                                                                                                                                                                                                                                           cTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGa                                                                                                                                                 >  7:1028917/1‑149 (MQ=255)

                                                                                                                                                                                                                                           cTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGa                                                                                                                                                 >  5:1354991/1‑149 (MQ=255)

                                                                                                                                                                                                                                           cTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGa                                                                                                                                                 >  5:628117/1‑149 (MQ=255)

                                                                                                                                                                                                                                            ttAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGaa                                                                                                                                                <  7:878826/149‑1 (MQ=255)

                                                                                                                                                                                                                                              aaTGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACg                                                                                                                                              >  8:977668/1‑149 (MQ=255)

                                                                                                                                                                                                                                               aTGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACga                                                                                                                                             >  7:1045738/1‑149 (MQ=255)

                                                                                                                                                                                                                                                tGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACgaa                                                                                                                                            >  8:990181/1‑149 (MQ=255)

                                                                                                                                                                                                                                                  ttCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACgaaga                                                                                                                                          <  8:796815/149‑1 (MQ=255)

                                                                                                                                                                                                                                                   tCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACgaagaa                                                                                                                                         <  5:1241892/149‑1 (MQ=255)

                                                                                                                                                                                                                                                      tAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGAt                                                                                                                                      >  6:969593/1‑149 (MQ=255)

                                                                                                                                                                                                                                                         gTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGc                                                                                                                                   <  7:795415/149‑1 (MQ=255)

                                                                                                                                                                                                                                                         gTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGc                                                                                                                                   <  7:1012493/149‑1 (MQ=255)

                                                                                                                                                                                                                                                          ttATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGcc                                                                                                                                  >  6:1432257/1‑149 (MQ=255)

                                                                                                                                                                                                                                                            aTGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGTAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGCGTTCTGCGCGCCGGAACCCGCACCAGTACTGATAACCAGCGGCATCACCCCCAGGATAAACCCCAGGCACGTCATCCGGaa  >  4:87376/1‑274 (MQ=255)

                                                                                                                                                                                                                                                             tGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGCGTTCTGCGCGCCGGAACCAGCACCAGTACTGATAACCAGCGGCATAACGCCGAGGATAAACGCCAGCGGGGTCATCAGGAt  >  2:8370/1‑274 (MQ=255)

                                                                                                                                                                                                                                                                atcactaaGGCCGCGTAAGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTAcc                                                                                                                            >  5:569361/7‑149 (MQ=255)

                                                                                                                                                                                                                                                                aTAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATAc                                                                                                                                                     <  6:328768/124‑1 (MQ=255)

                                                                                                                                                                                                                                                                aTAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATAc                                                                                                                                                     >  5:328768/1‑124 (MQ=255)

                                                                                                                                                                                                                                                                 tAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATAcc                                                                                                                                                    <  6:418524/124‑1 (MQ=255)

                                                                                                                                                                                                                                                                 tAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATAcc                                                                                                                                                    >  5:418524/1‑124 (MQ=255)

                                                                                                                                                                                                                                                                    aaaaGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTg                                                                                                                        >  5:265190/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                     taaGGCCGCGTAAGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGa                                                                                                                                                                                                      >  8:612776/2‑70 (MQ=39)

                                                                                                                                                                                                                                                                     taaGGCCGCGTAAGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGc                                                                                                                       >  8:684833/2‑149 (MQ=255)

                                                                                                                                                                                                                                                                     aaaGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGc                                                                                                                       <  7:1055030/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                      aaGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCg                                                                                                                      >  5:1309170/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                      aaGGCCGCGTAAGCCGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGTTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGACGAGGATTGCCAGTACCGTTGCg                                                                                                                      >  6:815901/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                       aGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAGGCGGCGGCGAACCACCACaaa                                                                                                                                                           >  7:289743/1‑111 (MQ=255)

                                                                                                                                                                                                                                                                       aGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCgg                                                                                                                     <  5:32626/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                       aGGCCGCGTAAGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCgg                                                                                                                     >  6:922565/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                        ggCCGCTTACGCTGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCTATCTCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGCAACGAAGAAGATCGCCAGTACCGTTGCGGt                                                                                                                    >  5:814330/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                        ggCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACAACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGt                                                                                                                    <  6:113200/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                          ccGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCa                                                                                                                  >  8:2299/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                                aCGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATcc                                                                                                            >  7:841252/1‑149 (MQ=255)

GCCAGATGAAATTGGTGACACTAATACCAGGATTGCTCTGAATATGACGTAATAACCGAGGAATGAATAAAGAATTACCGCAAATAATTAAGAATAGCCTCTAAATGATTATGGATTATTAAATCTATTAGCGCTACTTAATATAAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTACACCGGTACCTACTGCGTTCTGCGCGCCGGAACCAGCACCAGTACTGATAACCAGCGGCATAACGCCGAGGATAAACGCCAGCGAGGTCATCAGGAT  >  NC_000913/480955‑481481

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