BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,153,604	G→A	R102Q (CGG→CAG)	pabC →	aminodeoxychorismate lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,153,604	0	G	A	100.0%	730.6 / NA	171	R102Q (CGG→CAG)	pabC	aminodeoxychorismate lyase
Reads supporting (aligned to +/- strand): ref base G (0/0); new base A (41/130); total (41/130)

Only 100 of 174 total aligned reads displayed.

GCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCT  >  NC_000913/1153506‑1153701

   ggCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCAGATtc                                                                                              <  1:1685741/101‑1 (MQ=255)

      agcagAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCAGATtctct                                                                                           <  1:1685715/101‑1 (MQ=255)

      agcagAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACACAGATtctct                                                                                           <  1:1685779/101‑1 (MQ=255)

         agAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCAGTTTCTCTCCg                                                                                        <  1:1685818/101‑1 (MQ=255)

         agAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCAGATTCTCTCCg                                                                                        <  2:1685811/101‑1 (MQ=255)

         agAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCAGATTCTCTCCg                                                                                        <  1:1685778/101‑1 (MQ=255)

          gAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCAGATTCTCTCCGt                                                                                       <  2:1686137/101‑1 (MQ=255)

          gAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCAGATTCTCTCCGt                                                                                       <  2:1685804/101‑1 (MQ=255)

           aaCAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCAGATTCTCTCCGtt                                                                                      <  1:1685541/101‑1 (MQ=255)

           aaCAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCAGATTCTCTCCGtt                                                                                      <  1:1685449/101‑1 (MQ=255)

           aaCAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCAGATTCTCTCCGtt                                                                                      <  1:1685877/101‑1 (MQ=255)

           aaCAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCAGATTCTCTCCGtt                                                                                      <  1:1685805/101‑1 (MQ=255)

            aCAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCAGATTCTCTCCGTTa                                                                                     <  1:1685713/101‑1 (MQ=255)

             cAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCAGATTCTCTCCGTTAc                                                                                    <  2:1685682/101‑1 (MQ=255)

               gCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCAGATTCTCTCCGTTACgg                                                                                  <  1:1685748/101‑1 (MQ=255)

               gCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCAGATTCACTCCGTTACgg                                                                                  <  1:1685735/101‑1 (MQ=255)

                cAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCAGATTCTCTCCGTTACGGc                                                                                 <  1:1685840/101‑1 (MQ=255)

                      ggTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCAGATTCTCTCCGTTACGGCTTATcc                                                                           <  2:1685866/101‑1 (MQ=255)

                         gtgCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCAGATTCTCTCCGTTACGGCTTATCCTGc                                                                        <  1:1685679/101‑1 (MQ=255)

                            cTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCAGATTCTCTCCGTTACGGCTTATCCTGcaca                                                                     <  1:1685691/101‑1 (MQ=255)

                            cTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCAGATTCTCTCCGTTACGGCTTATCCTGcaca                                                                     <  2:1685601/101‑1 (MQ=255)

                                 aGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCAGATTCTCTCCGTTACGGCTTATCCTGCACATTACg                                                                <  1:1685896/101‑1 (MQ=255)

                                 aGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCAGATTCTCTCCGTTACGGCTTATCCTGCACATTACg                                                                <  1:1685808/101‑1 (MQ=255)

                                 aGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCAGATTCTCTCCGTTACGGCTTATCCTGCACATTACg                                                                <  2:1685835/101‑1 (MQ=255)

                                  gTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCAGATTCTCTCCGTTACGGCTTATCCTGca                                                                       <  2:1685722/93‑1 (MQ=255)

                                  gTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCAGATTCTCTCCGTTACGGCTTATCCTGCACATTACGa                                                               <  1:1685663/101‑1 (MQ=255)

                                   tCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCAGATTCTCTCCGTTACGGCTTATCCTGCACATTACGAc                                                              <  2:1686120/101‑1 (MQ=255)

                                   tCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCAGATTCTCTCCGTTACGGCTTATCCTGCACATTACGAc                                                              <  1:1685966/101‑1 (MQ=255)

                                    cGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCAGATTCTCTCCGTTACGGCTTATCCTGCACATTACGAcc                                                             <  1:1685870/101‑1 (MQ=255)

                                    cGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCAGATTCTCTCCGTTACGGCTTATCCTGCACATTACGAcc                                                             <  1:1685858/101‑1 (MQ=255)

                                    cGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCAGATTCTCTCCGTTACGGCTTATCCTGCACATTACGAcc                                                             <  1:1685757/101‑1 (MQ=255)

                                       gATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCAGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGtt                                                          <  2:1685911/101‑1 (MQ=255)

                                        aTCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCAGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGttt                                                         <  2:1685726/101‑1 (MQ=255)

                                           aGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCAGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCg                                                      <  2:1685873/101‑1 (MQ=255)

                                           aGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCAGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCg                                                      <  2:1685951/101‑1 (MQ=255)

                                            gTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCAGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGt                                                     <  2:1685833/101‑1 (MQ=255)

                                            gTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCAGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGt                                                     <  1:1685959/101‑1 (MQ=255)

                                             tCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCAGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTa                                                    <  2:1686022/101‑1 (MQ=255)

                                               gcgGTAGTGGTGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCAGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAAc                                                  >  1:1686240/1‑101 (MQ=255)

                                               gcgGTAGTGGCGGGCGAGTGTACAGCACATTGAACAGCGGACCGGCAACGCAGATTCTCTCCGTTACGGCTTATCCTGCACATTACGa                                                               <  2:1685764/88‑1 (MQ=255)

                                                cgGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCAGATTCTCTCCGTTACGGCTTATCCTGCACAtt                                                                   <  2:1685689/83‑1 (MQ=255)

                                                cgGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCAGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACg                                                 <  2:1686186/101‑1 (MQ=255)

                                                  gTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCAGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAg                                               <  1:1686095/101‑1 (MQ=255)

                                                  gTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCAGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAg                                               <  1:1685947/101‑1 (MQ=255)

                                                    aGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCAGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAggg                                             <  1:1685895/101‑1 (MQ=255)

                                                     gTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCAGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAgggg                                            <  2:1685784/101‑1 (MQ=255)

                                                      tGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCAGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGa                                           <  1:1685750/101‑1 (MQ=255)

                                                       ggcgggcgAGGGTACAGCACATTGAACAGCGGACCGGCAACGCAGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGAt                                          <  1:1685642/101‑1 (MQ=255)

                                                        gcgggcgAGGGTACAGCACATTGAACAGCGGACCGGCAACGCAGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGAtt                                         >  1:1686245/1‑101 (MQ=255)

                                                        gcgggcgAGGGTACAGCACATTGAACAGCGGACCGGCAACGCAGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGAtt                                         >  2:1686244/1‑101 (MQ=255)

                                                        gcgggcgAGGGTACAGCACATTGAACAGCGGACCGGCAACGCAGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGAtt                                         <  2:1685980/101‑1 (MQ=255)

                                                            gcgAGGGTACAGCACATTGAACAGCGGACCGGCAACGCAGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGt                                     >  2:1686248/1‑101 (MQ=255)

                                                            gcgAGGGTACAGCACATTGAACAGCGGACCGGCAACGCAGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGt                                     <  2:1685794/101‑1 (MQ=255)

                                                            gcgAGGGTACAGCACATTGAACAGCGGACCGGCAACGCAGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGt                                     <  2:1685820/101‑1 (MQ=255)

                                                             cgAGGGTACAGCACATTGAACAGCGGACCGGCAACGCAGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGtt                                    >  1:1686249/1‑101 (MQ=255)

                                                              gAGGGTACAGCACATTGAACAGCGGACCGGCAACGCAGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTg                                   >  1:1686251/1‑101 (MQ=255)

                                                              gAGGGTACAGCACATTGAACAGCGGACCGGCAACGCAGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTg                                   >  2:1686253/1‑101 (MQ=255)

                                                              gAGGGTACAGCACATTGAACAGCGGACCGGCAACGCAGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTg                                   >  2:1686252/1‑101 (MQ=255)

                                                               aGGGTACAGCACATTGAACAGCGGACCGGCAACGCAGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTgg                                  <  1:1685914/101‑1 (MQ=255)

                                                               aGGGTACAGCACATTGAACAGCGGACCGGCAACGCAGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTgg                                  <  1:1685823/101‑1 (MQ=255)

                                                                 ggTACAGCACATTGAACAGCGGACCGGCAACGCAGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGgcg                                <  1:1685962/101‑1 (MQ=255)

                                                                   tACAGCACATTGAACAGCGGACCGGCAACGCAGATtctc                                                                                            <  1:1686002/39‑1 (MQ=255)

                                                                   tACAGCACATTGAACAGCGGACCGGCAACGCAGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCt                              >  2:1686256/1‑101 (MQ=255)

                                                                    aCAGCACATTGAACAGCGGACCGGCAACGCAGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTa                             <  2:1686016/101‑1 (MQ=255)

                                                                     cAGCACATTGAACAGCGGACCGGCAACGCAGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTaa                            >  1:1686258/1‑101 (MQ=255)

                                                                      aGCACATTGAACAGCGGACCGGCAACGCAGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAg                           <  1:1685676/101‑1 (MQ=255)

                                                                      aGCACATTGAACAGCGGACCGGCAACGCAGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAg                           <  1:1685871/101‑1 (MQ=255)

                                                                       gCACATTGAACAGCGGACCGGCAACGCAGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGc                          >  2:1686260/1‑101 (MQ=255)

                                                                       gCACATTGAACAGCGGACCGGCAACGCAGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGc                          <  1:1686149/101‑1 (MQ=255)

                                                                        cacaTTGAACAGCGGACCGGCAACGCAGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGcc                         <  1:1686017/101‑1 (MQ=255)

                                                                          caTTGAACAGCGGACCGGCAACGCAGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCg                       <  2:1685660/101‑1 (MQ=255)

                                                                           aTTGAACAGCGGACCGGCAACGCAGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGc                                                       >  2:1686262/1‑68 (MQ=255)

                                                                            ttGAACAGCGGACCGGCAACGCAGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGt                     >  1:1686263/1‑101 (MQ=255)

                                                                             tGAACAGCGGACCGGCAACGCAGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTg                    <  1:1685927/101‑1 (MQ=255)

                                                                              gAACAGCGGACCGGCAACGCAGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGc                   <  1:1685944/101‑1 (MQ=255)

                                                                              gAACAGCGGACCGGCAACGCAGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGc                   >  2:1686264/1‑101 (MQ=255)

                                                                               aaCAGCGGACCGGCAACGCAGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCg                  <  2:1685686/101‑1 (MQ=255)

                                                                               aaCAGCGGACCGGCAACGCAGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCg                  <  2:1685670/101‑1 (MQ=255)

                                                                                 cAGCGGACCGGCAACGCAGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGc                >  1:1686267/1‑101 (MQ=255)

                                                                                 cAGCGGACCGGCAACGCAGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGc                <  1:1685999/101‑1 (MQ=255)

                                                                                  aGCGGACCGGCAACGCAGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCt               <  1:1686168/101‑1 (MQ=255)

                                                                                   gCGGACCGGCAACGCAGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTg              <  1:1685979/101‑1 (MQ=255)

                                                                                    cGGACCGGCAACGCAGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTgg             >  2:1686268/1‑101 (MQ=255)

                                                                                     ggACCGGCAACGCAGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGAt                                          >  2:1686270/1‑71 (MQ=255)

                                                                                     ggACCGGCAACGCAGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTggg            <  2:1685878/101‑1 (MQ=255)

                                                                                     ggACCGGCAACGCAGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTggg            <  1:1685771/101‑1 (MQ=255)

                                                                                      gACCGGCAACGCAGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTgggg           >  1:1686271/1‑101 (MQ=255)

                                                                                       aCCGGCAACGCAGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGgc          <  2:1685930/101‑1 (MQ=255)

                                                                                       aCCGGCAACGCAGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGgc          <  2:1685852/101‑1 (MQ=255)

                                                                                       aCCGGCAACGCAGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGgc          <  1:1685994/101‑1 (MQ=255)

                                                                                        ccGGCAACGCAGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGgcg         <  2:1685881/101‑1 (MQ=255)

                                                                                         cGGCAACGCAGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGgcgc        >  2:1686273/1‑101 (MQ=255)

                                                                                         cGGCAACGCAGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGgcgc        >  2:1686274/1‑101 (MQ=255)

                                                                                          ggCAACGCAGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCg                       >  2:1686276/1‑85 (MQ=255)

                                                                                           gCAACGCAGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCaa      <  1:1685950/101‑1 (MQ=255)

                                                                                            cAACGCAGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAAt     <  2:1686078/101‑1 (MQ=255)

                                                                                             aaCGCAGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTggg            >  2:1686280/1‑93 (MQ=255)

                                                                                             aaCGCAGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATc    >  1:1686279/1‑101 (MQ=255)

                                                                                              aCGCAGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATcc   <  2:1685711/101‑1 (MQ=255)

                                                                                               cgcaGATTCTCTCCGTTACGGCTTATTCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTgggg           >  1:1686283/1‑92 (MQ=255)

GCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCT  >  NC_000913/1153506‑1153701

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