BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,589,059	A→G	T489A (ACT→GCT)	acrD →	multidrug efflux pump RND permease AcrD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,589,059	0	A	G	100.0%	907.0 / NA	220	T489A (ACT→GCT)	acrD	multidrug efflux pump RND permease AcrD
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (89/131); total (89/131)

Only 100 of 220 total aligned reads displayed.

TGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCG  >  NC_000913/2588961‑2589156

tGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCGCt                                                                                                 <  2:3351143/101‑1 (MQ=255)

 ggCCTTCTTCGGCGGCACCACCGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCGCTc                                                                                                <  1:3351163/101‑1 (MQ=255)

  gCCTTCTTCGGCGGCACCACCGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCGCTcc                                                                                               <  1:3351262/101‑1 (MQ=255)

  gCCTTCTTCGGCGGCACCACCGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCGCTcc                                                                                               <  1:3351296/101‑1 (MQ=255)

      tcttcGGCGGCACCACCGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCGCTCCGGct                                                                                           <  1:3351281/101‑1 (MQ=255)

        ttcGGCGGCACCACCGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCGCTCCGGctct                                                                                         >  1:3351610/1‑101 (MQ=255)

        ttcGGCGGCACCACCGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCGCTCCGGctct                                                                                         >  1:3351612/1‑101 (MQ=255)

        ttcGGCGGCACCACCGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCGCTCCGGctct                                                                                         >  2:3351609/1‑101 (MQ=255)

           ggcggcACCACCGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCGCTCCGGCTCtgtg                                                                                      >  2:3351616/1‑101 (MQ=255)

            gcggcACCACCGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCGCTCCGGCTCtgtgt                                                                                     <  1:3351213/101‑1 (MQ=255)

                caccacCGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCca                                                                                 >  1:3351618/1‑101 (MQ=255)

                 accaccGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCcac                                                                                >  2:3351620/1‑101 (MQ=255)

                  ccaccGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCcaca                                                                               <  2:3351368/101‑1 (MQ=255)

                  ccaccGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCcac                                                                                <  2:3351493/100‑1 (MQ=255)

                   caccGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCCacac                                                                              <  2:3351443/101‑1 (MQ=255)

                   caccGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCCacac                                                                              <  1:3351474/101‑1 (MQ=255)

                     ccGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCCACacta                                                                            <  2:3351184/101‑1 (MQ=255)

                         tGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCCACACTACTTa                                                                        >  2:3351623/1‑101 (MQ=255)

                          gCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCCACACTACTTaa                                                                       <  2:3351346/101‑1 (MQ=255)

                            cATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCCACACTACTTAAGc                                                                     >  1:3351625/1‑101 (MQ=255)

                             aTCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCCACACTACTTAAGcc                                                                    <  2:3351357/101‑1 (MQ=255)

                             aTCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCCACACTACTTAAGcc                                                                    <  1:3351424/101‑1 (MQ=255)

                              tCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCCACACTACTTAAGCCa                                                                   <  1:3351196/101‑1 (MQ=255)

                              tCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCCACACTACTTAAGCCa                                                                   >  2:3351627/1‑101 (MQ=255)

                               cTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCCACACTACTTAAGCCAc                                                                  <  2:3351261/101‑1 (MQ=255)

                               cTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCCACACTACTTAAGCCAc                                                                  >  1:3351628/1‑101 (MQ=255)

                               cTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCCACACTACTTAAGCCAc                                                                  <  1:3351286/101‑1 (MQ=255)

                               cTATCGCCAGTTCTATATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCCACACTACTTAAg                                                                      <  1:3351147/97‑1 (MQ=255)

                                tATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCCACACTACTTAAGCCACt                                                                 <  2:3351399/101‑1 (MQ=255)

                                tATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCCACACTACTTAAGCCACt                                                                 <  1:3351343/101‑1 (MQ=255)

                                 aTCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTg                                                                <  2:3351204/101‑1 (MQ=255)

                                    gCCAGTTCTCTATTACCATTGTTGCGGCGATGGGGCTGTCAGTACTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGaaa                                                             >  2:3351631/1‑101 (MQ=255)

                                     ccAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGaaaa                                                            >  2:3351633/1‑101 (MQ=255)

                                      cAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGaaaaa                                                           <  1:3351152/101‑1 (MQ=255)

                                       aGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGaaaaaa                                                          >  2:3351638/1‑101 (MQ=255)

                                       aGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGaaaaaa                                                          >  1:3351637/1‑101 (MQ=255)

                                       aGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGaaaaaa                                                          <  1:3351322/101‑1 (MQ=255)

                                       aGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGaaaaaa                                                          <  2:3351528/101‑1 (MQ=255)

                                         ttCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGaaa                                                             >  1:3351642/1‑96 (MQ=255)

                                         ttCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAgg                                                        >  2:3351641/1‑101 (MQ=255)

                                           ctctATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTg                                                      <  2:3351293/101‑1 (MQ=255)

                                            tctATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGa                                                     >  2:3351644/1‑101 (MQ=255)

                                             ctATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAg                                                    <  1:3351300/101‑1 (MQ=255)

                                              tATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGc                                                   <  1:3351409/101‑1 (MQ=255)

                                              tATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGc                                                   >  2:3351647/1‑101 (MQ=255)

                                                ttACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGcat                                                 <  1:3351445/101‑1 (MQ=255)

                                                ttACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGcat                                                 <  1:3351429/101‑1 (MQ=255)

                                                ttACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGcat                                                 <  2:3351380/101‑1 (MQ=255)

                                                  aCCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGcatca                                               <  2:3351396/101‑1 (MQ=255)

                                                   ccATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGcatcat                                              >  2:3351651/1‑101 (MQ=255)

                                                   ccATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGcatcat                                              <  2:3351450/101‑1 (MQ=255)

                                                   ccATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGcatcat                                              >  1:3351652/1‑101 (MQ=255)

                                                    cATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATg                                             <  1:3351508/101‑1 (MQ=255)

                                                     aTTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATgg                                            <  1:3351372/101‑1 (MQ=255)

                                                     aTTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATgg                                            >  1:3351654/1‑101 (MQ=255)

                                                      ttgttgCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATggg                                           <  2:3351430/101‑1 (MQ=255)

                                                      ttgttgCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATggg                                           <  1:3351437/101‑1 (MQ=255)

                                                       tgttgCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGc                                          <  2:3351305/101‑1 (MQ=255)

                                                       tgttgCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGc                                          >  2:3351656/1‑101 (MQ=255)

                                                       tgttgCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGc                                          >  2:3351658/1‑101 (MQ=255)

                                                        gttgCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCa                                         >  2:3351660/1‑101 (MQ=255)

                                                        gttgCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCa                                         <  1:3351275/101‑1 (MQ=255)

                                                         ttgCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCaa                                        <  2:3351287/101‑1 (MQ=255)

                                                          tgCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCaaa                                       <  2:3351475/101‑1 (MQ=255)

                                                            cggcgATGGTGCTGTCAGTACTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCaaaaa                                     <  2:3351596/101‑1 (MQ=255)

                                                             ggcgATGGTGCTGTCAGTACTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAg                                    <  2:3351385/101‑1 (MQ=255)

                                                               cgATGGTGCTGTCAGTACTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGc                                  >  2:3351665/1‑101 (MQ=255)

                                                               cgATGGTGCTGTCAGTACTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGc                                  >  1:3351664/1‑101 (MQ=255)

                                                                gATGGTGCTGTCAGTACTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGc                                  <  1:3351479/100‑1 (MQ=255)

                                                                gATGGTGCTGTCAGTACTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCt                                 >  2:3351667/1‑101 (MQ=255)

                                                                 aTGGTGCTGTCAGTACTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCtt                                <  2:3351317/101‑1 (MQ=255)

                                                                   ggTGCTGTCAGTACTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCtttt                              <  2:3351323/101‑1 (MQ=255)

                                                                    gTGCTGTCAGTACTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAgg                                                        >  2:3351670/1‑74 (MQ=255)

                                                                    gTGCTGTCAGTACTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCttttt                             <  1:3351340/101‑1 (MQ=255)

                                                                    gTGCTGTCAGTACTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCttttt                             <  2:3351463/101‑1 (MQ=255)

                                                                     tGCTGTCAGTACTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATgg                                            <  1:3351390/85‑1 (MQ=255)

                                                                     tGCTGTCAGTACTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCtttttt                            >  1:3351671/1‑101 (MQ=255)

                                                                     tGCTGTCAGTACTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCtttttt                            <  1:3351453/101‑1 (MQ=255)

                                                                      gCTGTCAGTACTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTg                           <  1:3351349/101‑1 (MQ=255)

                                                                      gCTGTCAGTACTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTg                           >  1:3351674/1‑101 (MQ=255)

                                                                       cTGTCAGTACTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGc                          >  2:3351675/1‑101 (MQ=255)

                                                                       cTGTCAGTACTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGc                          <  1:3351521/101‑1 (MQ=255)

                                                                         gTCAGTACTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCt                        >  2:3351678/1‑101 (MQ=255)

                                                                         gTCAGTACTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCt                        <  2:3351329/101‑1 (MQ=255)

                                                                           cAGTACTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTgg                      >  2:3351681/1‑101 (MQ=255)

                                                                            aGTACTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGt                     <  2:3351504/101‑1 (MQ=255)

                                                                             gTACTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGtt                    <  1:3351416/101‑1 (MQ=255)

                                                                             gTACTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGtt                    >  2:3351683/1‑101 (MQ=255)

                                                                               aCTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTa                  >  2:3351686/1‑101 (MQ=255)

                                                                               aCTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTa                  >  1:3351687/1‑101 (MQ=255)

                                                                               aCTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTa                  <  1:3351606/101‑1 (MQ=255)

                                                                                cTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTaa                 <  1:3351234/101‑1 (MQ=255)

                                                                                cTGGTAGCGATGATCCTCGCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTaa                 <  1:3351366/101‑1 (MQ=255)

                                                                                      gCGATGATCCTCGCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGAt           <  2:3351255/101‑1 (MQ=255)

                                                                                      gCGATGATCCTCGCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGAt           >  2:3351689/1‑101 (MQ=255)

                                                                                      gCGATGATCCTCGCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGAt           <  2:3351490/101‑1 (MQ=255)

                                                                                         atgatCCTCGCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCa                                         >  2:3351691/1‑68 (MQ=255)

                                                                                         atgatCCTCGCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGtt        <  1:3351359/101‑1 (MQ=255)

                                                                                             tCCTCGCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAAc    >  1:3351694/1‑101 (MQ=255)

                                                                                             tCCTCGCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAAc    <  1:3351337/101‑1 (MQ=255)

TGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCG  >  NC_000913/2588961‑2589156

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