BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,257,838	+G	coding (2222/2241 nt)	tdcE ←	putative formate acetyltransferase 3

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,257,837	1	.	G	100.0%	820.2 / NA	234	A741A (GCG→GCC)	tdcE	putative formate acetyltransferase 3
Reads supporting (aligned to +/- strand): ref base . (0/0); new base G (121/113); total (121/113)

Only 100 of 247 total aligned reads displayed.

TTTTTCATACATCCTCCGGCGTCAGAGCGCCTGGGTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCAC‑GCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTT  >  NC_000913/3257745‑3257923

tttttCATACATCCTCCGGCGTCAGAGCGCCTGGGTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCAc                                                                                         <  2:372957/93‑1 (MQ=255)

  tttCATACATCCTCCGGCGTCAGAGCGCCTGGGTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCAc‑gg                                                                                      <  1:3199065/93‑2 (MQ=255)

  tttCATACATCCTCCGGCGTCAGAGCGCCTGGGTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCAc‑gg                                                                                      >  1:310203/1‑92 (MQ=255)

     cATACACCCTCCGGCGTCAGAGCGCCTGGGTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCAc‑ggcgt                                                                                   <  1:2315830/93‑5 (MQ=255)

        aCATCCTCCGGCGTCAGAGCGCCTGGGTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGc                                                                                 >  2:6409377/1‑93 (MQ=255)

        aCATCCTCCGGCGTCAGAGCGCCTGGGTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGc                                                                                 <  1:4654377/93‑1 (MQ=255)

         caccctccGGCGTCAGAGCGCATGGGTAAAGGTACGTGAAATACCATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCCCGGCGTAGcc                                                                                <  1:3117086/90‑1 (MQ=255)

         cATCCTCCGGCGTCAGAGCGCCTGGGTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGcc                                                                                <  1:5724834/93‑1 (MQ=255)

         cATCCTCCGGCGTCAGAGCGCCTGGGTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGcc                                                                                <  1:2201531/93‑1 (MQ=255)

         cATCCTCCGGCGTCAGAGCGCCTGGGTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGcc                                                                                >  1:1352660/1‑93 (MQ=255)

          aTCCTCCGGCGTCAGAGCGCCTGGGTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCa                                                                               >  2:6030191/1‑93 (MQ=255)

          aTCCTCCGGCGTCAGAGCGCCTGGGTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCa                                                                               >  2:4381656/1‑93 (MQ=255)

          aTCCTCCGGCGTCAGAGCGCCTGGGTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCa                                                                               >  2:5656182/1‑93 (MQ=255)

           tcctccGGCGTCAGAGCGCCTGGGTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCag                                                                              >  1:4219754/1‑93 (MQ=255)

           tcctccGGCGTCAGAGCGCCTGGGTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCag                                                                              >  1:5208319/1‑93 (MQ=255)

            cctccGGCGTCAGAGCGCCTGGGTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCaga                                                                             <  2:1959391/93‑1 (MQ=255)

              tccGGCGTCAGAGCGCCTGGGTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAgaga                                                                           <  1:5112816/93‑1 (MQ=255)

                cGGCGTCAGAGCGCCTGGGTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGaca                                                                         >  1:4381376/1‑93 (MQ=255)

                 ggCGTCAGAGCGCCTGGGTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGacac                                                                        <  1:6466485/93‑1 (MQ=255)

                 ggCGTCAGAGCGCCTGGGTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGacac                                                                        <  2:1585230/93‑1 (MQ=255)

                  gCGTCAGAGCGCCTGGGTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGATACACg                                                                       >  2:3231910/1‑93 (MQ=255)

                  gCGTCAGAGCGCCTGGGTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACg                                                                       >  2:5521235/1‑93 (MQ=255)

                  gCGTCAGAGCGCCTGGGTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACg                                                                       >  2:5717669/1‑93 (MQ=255)

                  gCGTCAGAGCGCCTGGGTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACg                                                                       >  1:548737/1‑93 (MQ=255)

                   cGTCAGAGCGCCTGGGTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACgg                                                                      <  1:5338495/93‑1 (MQ=255)

                    gTCAGAGCGCCTGGGTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGa                                                                     >  1:5511778/1‑93 (MQ=255)

                    gTCAGAGCGCCTGGGTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGa                                                                     <  2:3985775/93‑1 (MQ=255)

                    gTCAGAGCGCCTGGGTAAAGGCACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGACGCGCACGGCGTAGCCAGCGACACGGa                                                                     >  1:4982530/1‑93 (MQ=255)

                     tCAGAGCGCCTGGGTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGAt                                                                    >  2:4979307/1‑93 (MQ=255)

                     tCAGAGCGCCTGGGTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGAt                                                                    >  2:1998489/1‑93 (MQ=255)

                       agagCGCCTGGGTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTg                                                                  <  2:3324723/93‑1 (MQ=255)

                       agagCGCCTGGGTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTg                                                                  <  2:891072/93‑1 (MQ=255)

                          gcgcCTGGGTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCa                                                               <  1:4032910/93‑1 (MQ=255)

                          gcgcCTGGGTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCa                                                               <  2:5842383/93‑1 (MQ=255)

                              cTGGGTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGtt                                                           >  2:1805169/1‑93 (MQ=255)

                              cTGGGTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGtt                                                           >  1:5131760/1‑93 (MQ=255)

                              cTGGGTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGtt                                                           <  1:289831/93‑1 (MQ=255)

                              cTGGGTAAAGGTACGTGAAATAACATCCTGTTACTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGtt                                                           >  2:3978647/1‑93 (MQ=255)

                               tGGGTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTa                                                          >  1:6415589/1‑93 (MQ=255)

                                gggTAAAGGTACGTGAAATAACATCCTCTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGGTAg                                                         >  1:101069/1‑93 (MQ=255)

                                  gTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGa                                                       <  2:1323607/93‑1 (MQ=255)

                                    aaaGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGata                                                     <  2:3196981/93‑1 (MQ=255)

                                    aaaGGTAAGTGAAATAACATCCTGTTGCGGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGata                                                     <  2:4680052/93‑1 (MQ=255)

                                     aaGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGatat                                                    >  1:4999762/1‑93 (MQ=255)

                                     aaGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGatat                                                    >  2:881157/1‑93 (MQ=255)

                                      aGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATAtt                                                   <  1:6035896/93‑1 (MQ=255)

                                       ggTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATAttt                                                  >  2:1153652/1‑93 (MQ=255)

                                       ggTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATAttt                                                  >  2:5805844/1‑93 (MQ=255)

                                       ggTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATAttt                                                  >  1:2442464/1‑93 (MQ=255)

                                        gTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATAtttt                                                 >  2:5232300/1‑93 (MQ=255)

                                        gTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATAtttt                                                 >  2:3018260/1‑93 (MQ=255)

                                         tACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATAttttt                                                >  1:1211530/1‑93 (MQ=255)

                                          aCGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTc                                               <  1:4583788/93‑1 (MQ=255)

                                          aCGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTc                                               >  2:3667013/1‑93 (MQ=255)

                                             tGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGGCAGGTTAGGATATTTTTCCgg                                            >  2:4138669/1‑93 (MQ=255)

                                              gAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCggg                                           <  1:4904801/93‑1 (MQ=255)

                                               aaaTAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGt                                          >  2:171535/1‑93 (MQ=255)

                                                aaTAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTg                                         <  1:947292/93‑1 (MQ=255)

                                                  tAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCt                                       <  2:427455/93‑1 (MQ=255)

                                                  tAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCt                                       >  1:1130696/1‑93 (MQ=255)

                                                   aaCATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTc                                      <  2:608993/93‑1 (MQ=255)

                                                   aaCATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTc                                      >  2:5404993/1‑93 (MQ=255)

                                                   aaCATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTc                                      >  1:5224735/1‑93 (MQ=255)

                                                         cTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGc                                >  2:1207944/1‑93 (MQ=255)

                                                          tGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCa                               <  2:339228/93‑1 (MQ=255)

                                                            ttGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATc                             <  1:5941155/93‑1 (MQ=255)

                                                            ttGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATc                             >  2:2878292/1‑93 (MQ=255)

                                                             tGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATcc                            <  1:5660391/93‑1 (MQ=255)

                                                             tGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATcc                            <  1:601208/93‑1 (MQ=255)

                                                             tGCTGTGCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGAATATTTTTCCGGGTGCTCGATGGCATcc                            <  2:1296182/93‑1 (MQ=255)

                                                              gCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCca                           >  2:5082587/1‑93 (MQ=255)

                                                              gCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCca                           >  2:5287172/1‑93 (MQ=255)

                                                               cTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCcag                          >  2:5201111/1‑93 (MQ=255)

                                                               cTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCcag                          >  1:3458602/1‑93 (MQ=255)

                                                                tGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCcagc                         <  1:1561728/93‑1 (MQ=255)

                                                                tGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCcagc                         >  2:1549106/1‑93 (MQ=255)

                                                                tGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCcagc                         <  2:2072252/93‑1 (MQ=255)

                                                                 gtgcACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCcagca                        <  2:1116376/90‑1 (MQ=255)

                                                                 gTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCcagca                        <  2:677971/93‑1 (MQ=255)

                                                                 gTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCcagca                        <  2:2461331/93‑1 (MQ=255)

                                                                  ttCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCcagcag                       <  2:3922515/93‑1 (MQ=255)

                                                                    cACGTGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAgcagca                     <  1:882246/93‑1 (MQ=255)

                                                                      cGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCAtt                   <  1:3954654/93‑1 (MQ=255)

                                                                      cGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCAtt                   <  1:1754988/93‑1 (MQ=255)

                                                                       gggTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTGGATGGCATCCAGCAGCAttt                  <  1:2786828/93‑1 (MQ=255)

                                                                       gggTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCCttt                  >  1:1867757/1‑93 (MQ=255)

                                                                       gcgTCAGTGCGGTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCAttt                  <  1:1378069/91‑1 (MQ=37)

                                                                        ggTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTc                 >  2:4232746/1‑93 (MQ=255)

                                                                         gTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTcc                <  2:3065357/93‑1 (MQ=255)

                                                                         gTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTcc                >  2:1966754/1‑93 (MQ=255)

                                                                          tCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTccc               <  1:2605043/93‑1 (MQ=255)

                                                                          tCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTccc               <  1:1949842/93‑1 (MQ=255)

                                                                            aGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGa             >  1:4475135/1‑93 (MQ=255)

                                                                             gTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGAt            >  1:727034/1‑93 (MQ=255)

                                                                              tGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGAtt           >  2:403949/1‑93 (MQ=255)

                                                                               gCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTc          >  2:2509257/1‑93 (MQ=255)

                                                                                  ttGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCAtt       <  1:5785105/93‑1 (MQ=255)

                                                                                   tGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTa      <  2:3561945/93‑1 (MQ=255)

                                                                                    gAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTac     <  1:3898292/93‑1 (MQ=255)

                                                                                      aGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTacgt   <  1:5453723/93‑1 (MQ=255)

TTTTTCATACATCCTCCGGCGTCAGAGCGCCTGGGTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCAC‑GCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTT  >  NC_000913/3257745‑3257923

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