BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,054,441	T→C	Y267H (TAC→CAC)	shiA →	shikimate:H(+) symporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,054,441	0	T	C	100.0%	1601.6 / NA	381	Y267H (TAC→CAC)	shiA	shikimate:H(+) symporter
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (180/201); total (180/201)

Only 100 of 386 total aligned reads displayed.

CCGCGAAAAAACGCATCCCGGTTATCGAAGCGCTGTTACGACATCCCGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGTACATCGTTACTGCCTTTGCACTTAATTATTCAACCCAGAATATGGGGCTACCGCGCGAACTTTTCCTTAATATTGGTTTGCTGGTAGGTGGATTAAG  >  NC_000913/2054343‑2054538

ccGCGAAAAAACGCATCCCGGTTATCGAAGCGCTGTTACGACATCCCGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGCAc                                                                                                 <  2:3933299/101‑1 (MQ=255)

ccGCGAAAAAACGCATCCCGGTTATCGAAGCGCTGTTACGACATCCCGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGCAc                                                                                                 >  2:3934010/1‑101 (MQ=255)

  gcgAAAAAACGCATCCCGGTTATCGAAGCGCTGTTACGACATCCCGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGCACAt                                                                                               >  1:3934018/1‑101 (MQ=255)

  gcgAAAAAACGCATCCCGGTTATCGAAGCGCTGTTACGACATCCCGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGCACAt                                                                                               >  1:3934015/1‑101 (MQ=255)

   cgAAAAAACGCATCCCGGTTATCGAAGCGCTGTTACGACATCCCGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGCACATc                                                                                              <  1:3933191/101‑1 (MQ=255)

   cgAAAAAACGCATCCCGGTTATCGAAGCGCTGTTACGACATCCCGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGCACATc                                                                                              >  2:3934022/1‑101 (MQ=255)

    gAAAAAACGCATCCCGGTTATCGAAGCGCTGTTACGACATCCCGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGCACATCg                                                                                             <  1:3933512/101‑1 (MQ=255)

     aaaaaaCGCATCCCGGTTATCGAAGCGCTGTTACGACATCCCGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGCACATc                                                                                              <  2:3933639/99‑1 (MQ=255)

      aaaaaCGCATCCCGGTTATCGAAGCGCTGTTACGACATCCCGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGCACATCGtt                                                                                           >  1:3934029/1‑101 (MQ=255)

      aaaaaCGCATCCCGGTTATCGAAGCGCTGTTACGACATCCCGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGCACATCGtt                                                                                           <  1:3933143/101‑1 (MQ=255)

      aaaaaCGCATCCCGGTTATCGAAGCGCTGTTACGACATCCCGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGCACATCGtt                                                                                           >  2:3934028/1‑101 (MQ=255)

         aaCGCATCCCGGTTATCGAAGCGCTGTTACGACATCCCGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGCACATCGTTACt                                                                                        <  2:3933374/101‑1 (MQ=255)

          aCGCATCCCGGTTATCGAAGCGCTGTTACGACATCCCGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGCACATCGTTACTg                                                                                       <  1:3933350/101‑1 (MQ=255)

           cGCATCCCGGTTATCGAAGCGCTGTTACGACATCCCGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGCACATCGTTACTGc                                                                                      <  2:3933411/101‑1 (MQ=255)

              aTCCCGGTTATCGAAGCGCTGTTACGACATCCCGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGCACATCGTTACTGCCtt                                                                                   >  1:3934035/1‑101 (MQ=255)

                cccGGTTATCGAAGCGCTGTTACGACATCCCGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGCACATCGTTACTGCCTTTg                                                                                 <  2:3933547/101‑1 (MQ=255)

                 ccGGTTATCGAAGCGCTGTTACGACATCCCGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGCACATCGTTACTGCCTTTGc                                                                                >  2:3934039/1‑101 (MQ=255)

                 ccGGTTATCGAAGCGCTGTTACGACATCCCGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGCACATCGTTACTGCCTTTGc                                                                                <  1:3933442/101‑1 (MQ=255)

                  cGGTTATCGAAGCGCTGTTACGACATCCCGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGCACATCGTTACTGCCTTTGCa                                                                               >  1:3934042/1‑101 (MQ=255)

                     ttATCGAAGCGCTGTTACGACATCCCGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGCACATCGTTACTGCCTTTGCACtt                                                                            >  2:3934047/1‑101 (MQ=255)

                      tATCGAAGCGCTGTTACGACATCCCGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGCACATCGTTACTGCCTTTGCACTTa                                                                           <  1:3933471/101‑1 (MQ=255)

                      tATCGAAGCGCTGTTACGACATCCCGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGCACATCGTTACTGCCTTTGCACTTa                                                                           <  2:3933620/101‑1 (MQ=255)

                       aTCGAAGCGCTGTTACGACATCCCGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGCACATCGTTACTGCCTTTGCACTTaa                                                                          <  1:3933687/101‑1 (MQ=255)

                       aTCGAAGCGCTGTTACGACATCCCGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGCACATCGTTACTGCCTTTGCACTTaa                                                                          >  1:3934051/1‑101 (MQ=255)

                       aTCGAAGCGCTGTTACGACATCCCGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGCACATCGTTACTGCCTTTGCACTTaa                                                                          >  2:3934052/1‑101 (MQ=255)

                          gAAGCGCTGTTACGACATCCCGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGCACATCGTTACTGCCTTTGCACTTAatta                                                                       >  1:3934059/1‑101 (MQ=255)

                          gAAGCGCTGTTACGACATCCCGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGCACATCGTTACTGCCTTTGCACTTAatta                                                                       >  2:3934061/1‑101 (MQ=255)

                           aaGCGCTGTTACGACATCCCGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGCACATCGTTACTGCCTTTGCACTTAattat                                                                      >  2:3934065/1‑101 (MQ=255)

                              cgcTGTTACGACATCCCGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGCACATCGTTACTGCCTTTGCACt                                                                             <  2:3933447/91‑1 (MQ=255)

                               gcTGTTACGACATCCCGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGCACATCGTTACTGCCTTTGCACTTAATTATTCaa                                                                  >  1:3934073/1‑101 (MQ=255)

                               gcTGTTACGACATCCCGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGCACATCGTTACTGCCTTTGCACTTAATTATTCaa                                                                  >  2:3934071/1‑101 (MQ=255)

                                 tGTTACGACATCCCGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGCACATCGTTACTGCCTTTGCACTTAATTATTCAAcc                                                                >  2:3934078/1‑101 (MQ=255)

                                  gTTACGACATCCCGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGCACATCGTTACTGCCTTTGCACTTAATTATTCAAccc                                                               >  1:3934079/1‑101 (MQ=255)

                                    tACGACATCCCGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGCACATCGTTACTGCCTTTGCACTTAATTATTCAACCCAg                                                             <  1:3933592/101‑1 (MQ=255)

                                     aCGACATCCCGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGCACATCGTTACTGCCTTTGCACTTAATTATTCAACCCAGa                                                            <  1:3933264/101‑1 (MQ=255)

                                     aCGACATCCCGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGCACATCGTTACTGCCTTTGCACTTAATTATTCAACCCAGa                                                            <  2:3933149/101‑1 (MQ=255)

                                      cGACATCCCGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGCACATCGTTACTGCCTTTGCACTTAATTATTCAACCCAGaa                                                           >  2:3934084/1‑101 (MQ=255)

                                        aCATCCCGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGCACATCGTTACTGCCTTTGCACTTAATTATTCAACCCAGAata                                                         >  1:3934087/1‑101 (MQ=255)

                                         cATCCCGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGCACATCGTTACTGCCTTTGCACTTAATTATTCAACCCAGAatat                                                        <  2:3933666/101‑1 (MQ=255)

                                          aTCCCGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGCACATCGTTACTGCCTTTGCACTTAATTATTCAACCCAGAATATg                                                       >  1:3934094/1‑101 (MQ=255)

                                          aTCCCGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGCACATCGTTACTGCCTTTGCACTTAATTATTCAACCCAGAATATg                                                       >  2:3934093/1‑101 (MQ=255)

                                           tCCCGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGCACATCGTTACTGCCTTTGCACTTAATTATTCAACCCAGAATATgg                                                      <  2:3933472/101‑1 (MQ=255)

                                           tCCCGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGCACATCGTTACTGCCTTTGCACTTAATTATTCAACCCAGAATATgg                                                      <  2:3933525/101‑1 (MQ=255)

                                            cccGGTGCTTTCCTGAATATTATTGCGCTACGACTGTGCGAATTGCTGACGATGCACGTCGTTACTGCCTTTGCACTTAATTATTCAACCCAGAATATggg                                                     <  2:3933342/101‑1 (MQ=255)

                                            cccGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGCACATCGTTACTGCCTTTGCACTTAATTATTCAACCCAGAATATggg                                                     <  1:3933536/101‑1 (MQ=255)

                                             ccGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGCACATCGTTACTGCCTTTGCACTTAATTATTCAACCCAGAat                                                          <  2:3933439/95‑1 (MQ=255)

                                              cGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGCACATCGTTACTGCCTTTGCACTTAATTATTCAACCCAGAATATGGGGc                                                   >  2:3934099/1‑101 (MQ=255)

                                              cGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGCACATCGTTACTGCCTTTGCACTTAATTATTCAACCCAGAATATGGGGc                                                   <  2:3933719/101‑1 (MQ=255)

                                                gTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGCACATCGTTACTGCCTTTGCACTTAATTATTCAAc                                                                 <  1:3933370/85‑1 (MQ=255)

                                                 tGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGCACATCGTTACTGCCTTTGCACTTAATTATTCAACCCAGAATATGGGGCTAc                                                <  1:3933609/101‑1 (MQ=255)

                                                 tGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGCACATCGTTACTGCCTTTGCACTTAATTATTCAACCCAGAATATGGGGCTAc                                                >  1:3934103/1‑101 (MQ=255)

                                                 tGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGCACATCGTTACTGCCTTTGCACTTAATTATTCAACCCAGAATATGGGGCTAc                                                <  2:3933430/101‑1 (MQ=255)

                                                  gCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGCACATCGTTACTGCCTTTGCACTTAATTATTCAACCCAGAATATGGGGCTAcc                                               <  1:3933464/101‑1 (MQ=255)

                                                   cTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGCACATCGTTACTGCCTTTGCACTTAATTATTCAACCCAGAATCTGGGGCTACcg                                              <  1:3933803/101‑1 (MQ=255)

                                                    tttCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGCACATCGTTACTGCCTTTGCACTTAATTATTCAACCCAGAATATGGGGCTACcgc                                             >  1:3934108/1‑101 (MQ=255)

                                                      tCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGCACATCGTTACTGCCTTTGCACTTAATTATTCAACCCAGAATATGGGGCTACcgcgc                                           >  2:3934111/1‑101 (MQ=255)

                                                      tCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGCACATCGTTACTGCCTTTGCACTTAATTATTCAACCCAGAATATGGGGCTACcgcgc                                           <  1:3933557/101‑1 (MQ=255)

                                                         tGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGCACATCGTTACTGCCTTTGCACTTAATTATTCAACCCAGAATATGGGGCTACCGCGCGaa                                        <  2:3933694/101‑1 (MQ=255)

                                                          gAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGCACATCGTTACTGCCTTTGCACTTAATTATTCAACCCAGAATATGGGGCTACCGCGCGAAc                                       >  2:3934117/1‑101 (MQ=255)

                                                          gAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGCACATCGTTACTGCCTTTGCACTTAATTATTCAACCCAGAATATGGGGCTACCGCGCGAAc                                       >  1:3934120/1‑101 (MQ=255)

                                                           aaGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGCACATCGTTACTGCCTTTGCACTTAATTATTCAACCCAGAATATGGGGCTACCGCGCGAACt                                      >  2:3934122/1‑101 (MQ=255)

                                                            aGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGCACATCGTTACTGCCTTTGCACTTAATTATTCAACCCAGAATATGGGGCTACCGCGCGAACtt                                     <  1:3933723/101‑1 (MQ=255)

                                                             gATTATTGCGCTACGACTGTGCGAATTGCTGACGATGCACATCGTTACTGCCTTTGCACTTAATTATTCAACCCAGAATATGGGGCTACCGCGCGAACttt                                    >  1:3934129/1‑101 (MQ=255)

                                                               ttattGCGCTACGACTGTGCGAATTGCTGACGATGCACATCGTTACTGCCTTTGCACTTAATTATTCAACCCAGAATATGGGGCTACCGCGCGAACTTTTc                                  <  2:3933582/101‑1 (MQ=255)

                                                                tattGCGCTACGACTGTGCGAATTGCTGACGATGCACATCGTTACTGCCTTTGCACTTAATTATTCAACCCAGAATATGGGGCTACCGCGCGAACTTTTcc                                 >  1:3934132/1‑101 (MQ=255)

                                                                  ttGCGCTACGACTGTGCGAATTGCTGACGATGCACATCGTTACTGCCTTTGCACTTAATTATTCAACCCAGAATATGGGGCTACCGCGCGAACTTTTCCtt                               >  2:3934137/1‑101 (MQ=255)

                                                                  ttGCGCTACGACTGTGCGAATTGCTGACGATGCACATCGTTACTGCCTTTGCACTTAATTATTCAACCCAGAATATGGGGCTACCGCGCGAACTTTTCCtt                               >  1:3934136/1‑101 (MQ=255)

                                                                   tGCGCTACGACTGTGCGAATTGCTGACGATGCACATCGTTACTGCCTTTGCACTTAATTATTCAACCCAGAATATGGGGCTACCGCGCGAACTTTTCCTTa                              <  1:3933415/101‑1 (MQ=255)

                                                                   tGCGCTACGACTGTGCGAATTGCTGACGATGCACATCGTTACTGCCTTTGCACTTAATTATTCAACCCAGAATATGGGGCTACCGCGCGAACTTTTCCTTa                              <  1:3933931/101‑1 (MQ=255)

                                                                    gcgcTACGACTGTGCGAATTGCTGACGATGCACATCGTTACTGCCTTTGCACTTAATTATTCAACCCAGAATATGGGGCTACCGCGCGAACTTTTCCTTaa                             <  1:3933711/101‑1 (MQ=255)

                                                                     cgcTACGACTGTGCGAATTGCTGACGATGCACATCGTTACTGCCTTTGCACTTAATTATTCAACCCAGAATATGGGGCTACCGCGCGAACTTTTCCTTAat                            >  2:3934140/1‑101 (MQ=255)

                                                                      gcTACGACTGTGCGAATTGCTGACGATGCACATCGTTACTGCCTTTGCACTTAATTATTCAACCCAGAATATGGGGCTACCGCGCGAACTTTTCCTTAata                           <  2:3933492/101‑1 (MQ=255)

                                                                      gcTACGACTGTGCGAATTGCTGACGATGCACATCGTTACTGCCTTTGCACTTAATTATTCAACCCAGAATATGGGGCTACCGCGCGAACTTTTCCTTAata                           <  1:3933552/101‑1 (MQ=255)

                                                                        tACGACTGTGCGAATTGCTGACGATGCACATCGTTACTGCCTTTGCACTTAATTATTCAACCCAGAATATGGGGCTACCGCGCGAACTTTTCCTTAATAtt                         >  1:3934147/1‑101 (MQ=255)

                                                                         aCGACTGTGCGAATTGCTGACGATGCACATCGTTACTGCCTTTGCACTTAATTATTCAACCCAGAATATGGGGCTACCGCGCGAACTTTTCCTTAATATTg                        <  2:3933598/101‑1 (MQ=255)

                                                                         aCGACTGTGCGAATTGCTGACGATGCACATCGTTACTGCCTTTGCACTTAATTATTCAACCCAGAATATGGGGCTACCGCGCGAACTTTTCCTTAATATTg                        >  2:3934148/1‑101 (MQ=255)

                                                                           gACTGTGCGAATTGCTGACGATGCACATCGTTACTGCCTTTGCACTTAATTATTCAACCCAGAATATGGGGCTACCGCGCGAACTTTTCCTTAATATTGGt                      <  1:3933766/101‑1 (MQ=255)

                                                                            aCTGTGCGAATTGCTGACGATGCACATCGTTACTGCCTTTGCACTTAATTATTCAACCCAGAATATGGGGCTACCGCGCGAACTTTTCCTTAATAtt                         <  2:3933507/97‑1 (MQ=255)

                                                                            aCTGTGCGAATTGCTGACGATGCACATCGTTACTGCCTTTGCACTTAATTATTCAACCCAGAATATGGGGCTACCGCGCGAACTTTTCCTTAATATTGGtt                     <  1:3933423/101‑1 (MQ=255)

                                                                              tgtgCGAATTGCTGACGATGCACATCGTTACTGCCTTTGCACTTAATTATTCAACCCAGAATATGGGGCTACCGCGCGAACTTTTCCTTAATATTGGTTTg                   >  2:3934155/1‑101 (MQ=255)

                                                                              tgtgCGAATTGCTGACGATGCACATCGTTACTGCCTTTGCACTTAATTATTCAACCCAGAATATGGGGCTACCGCGCGAACTTTTCCTTAATATTGGTTTg                   >  1:3934157/1‑101 (MQ=255)

                                                                                tgCGAATTGCTGACGATGCACATCGTTACTGCCTTTGCACTTAATTATTCAACCCAGAATATGGGGCTACCGCGCGAACTTTTCCTTAATATTGGTTTGCt                 <  1:3933572/101‑1 (MQ=255)

                                                                                tgCGAATTGCTGACGATGCACATCGTTACTGCCTTTGCACTTAATTATTCAACCCAGAATATGGGGCTACCGCGCGAACTTTTCCTTAATATTGGTTTGCt                 >  2:3934161/1‑101 (MQ=255)

                                                                                  cGAATTGCTGACGATGCACATCGTTACTGCCTTTGCACTTAATTATTCAACCCAGAATATGGGGCTACCGCGCGAACtttt                                   <  1:3933705/81‑1 (MQ=255)

                                                                                  cGAATTGCTGACGATGCACATCGTTACTGCCTTTGCACTTAATTATTCAACCCAGAATATGGGGCTACCGCGCGAACTTTTCCTTAATATTGGTTTGCTgg               >  2:3934166/1‑101 (MQ=255)

                                                                                   gAATTGCTGACGATGCACATCGTTACTGCCTTTGCACTTAATTATTCAACCCAGAATATGGGGCTACCGCGCGAACTTTTCCTTAATATTGGTTTGCTGGt              >  1:3934171/1‑101 (MQ=255)

                                                                                    aaTTGCTGACGATGCACATCGTTACTGCCTTTGCACTTAATTATTCAACCCAGAATATGGGGCTACCGCGCGAACTTTTCCTTAATATTGGTTTGCTGGTa             <  2:3933836/101‑1 (MQ=255)

                                                                                    aaTTGCTGACGATGCACATCGTTACTGCCTTTGCACTTAATTATTCAACCCAGAATATGGGGCTACCGCGCGAACTTTTCCTTAATATTGGTTTGCTGGTa             <  2:3933793/101‑1 (MQ=255)

                                                                                    aaTTGCTGACGATGCACATCGTTACTGCCTTTGCACTTAATTATTCAACCCAGAATATGGGGCTACCGCGCGAACTTTTCCTTAATATTGGTTTGCTGGTa             >  2:3934173/1‑101 (MQ=255)

                                                                                    aaTTGCTGACGATGCACATCGTTACTGCCTTTGCACTTAATTATTCAACCCAGAATATGGGGCTACCGCGCGAACTTTTCCTTAATATTGGTTTGCTGGTa             <  1:3933838/101‑1 (MQ=255)

                                                                                     aTTGCTGACGATGCACATCGTTACTGCCTTTGCACTTAATTATTCAACCCAGAATATGGGGCTACCGCGCGAACTTTTCCTTAATATTGGTTTGCTGGTAg            <  1:3933654/101‑1 (MQ=255)

                                                                                      ttGCTGACGATGCACATCGTTACTGCCTTTGCACTTAATTATTCAACCCAGAATATGGGGCTACCGCGCGAACTTTTCCTTAATATTGGTTTGCTGGTAgg           <  2:3933734/101‑1 (MQ=255)

                                                                                      ttGCTGACGATGCACATCGTTACTGCCTTTGCACTTAATTATTCAACCCAGAATATGGGGCTACCGCGCGAACTTTTCCTTAATATTGGTTTGCTGGTAgg           <  1:3933627/101‑1 (MQ=255)

                                                                                       tGCTGACGATGCACATCGTTACTGCCTTTGCACTTAATTATTCAACCCAGAATATGGGGCTACCGCGCGAACTTTTCCTTAATATTGGTTTGCTGGTAGGt          >  1:3934178/1‑101 (MQ=255)

                                                                                       tGCTGACGATGCACATCGTTACTGCCTTTGCACTTAATTATTCAACCCAGAATATGGGGCTACCGCGCGAACTTTTCCTTAATATTGGTTTGCTGGTAGGt          <  1:3933498/101‑1 (MQ=255)

                                                                                          tGACGATGCACATCGTTACTGCCTTTGCACTTAATTATTCAACCCAGAATATGGGGCTACCGCGCGAACTTTTCCTTAATATTGGTTTGCTGGTAGGTGGa       <  2:3933537/101‑1 (MQ=255)

                                                                                           gACGATGCACATCGTTACTGCCTTTGCACTTAATTATTCAACCCAGAATATGGGGCTACCGCGCGAACTTTTCCTTAATATTGGTTTGCTGGTAGGTGGAt      >  1:3934183/1‑101 (MQ=255)

                                                                                            aCGATGCACATCGTTACTGCCTTTGCACTTAATTATTCAACCCAGAATATGGGGCTACCGCGCGAACTTTTCCTTAATATTGGTTTGCTGGTAGGTGGAtt     >  2:3934185/1‑101 (MQ=255)

                                                                                            aCGATGCACATCGTTACTGCCTTTGCACTTAATTATTCAACCCAGAATATGGGGCTACCGCGCGAACTTTTCCTTAATATTGGTTTGCTGGTAGGTGGAtt     <  2:3933756/101‑1 (MQ=255)

                                                                                               aTGCACATCGTTACTGCCTTTGCACTTAATTATTCAACCCAGAATATGGGGCTACCGCGCGAACTTTTCCTTAATATTGGTTTGCTGGTAGGTGGATTAAg  >  1:3934191/1‑101 (MQ=255)

CCGCGAAAAAACGCATCCCGGTTATCGAAGCGCTGTTACGACATCCCGGTGCTTTCCTGAAGATTATTGCGCTACGACTGTGCGAATTGCTGACGATGTACATCGTTACTGCCTTTGCACTTAATTATTCAACCCAGAATATGGGGCTACCGCGCGAACTTTTCCTTAATATTGGTTTGCTGGTAGGTGGATTAAG  >  NC_000913/2054343‑2054538

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