BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,073,577	T→C	intergenic (‑6/‑162)	yihU ← / → yihV	3‑sulfolactaldehyde reductase/6‑deoxy‑6‑sulfofructose kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,073,577	0	T	C	100.0%	1657.3 / NA	393	intergenic (‑6/‑162)	yihU/yihV	3‑sulfolactaldehyde reductase/6‑deoxy‑6‑sulfofructose kinase
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (173/220); total (173/220)

Only 100 of 410 total aligned reads displayed.

ATCAAAGACGCGAAGTTGATGCCCTTGCTGCAATAAATTGCTCGCCATTGGCGAACCCATTTGTCCTAAACCGATAAACGCGATTGCTGCCATAACCCTCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGTATTTGTAATTTATCGTCAAAAAATTGACAGC  >  NC_000913/4073479‑4073674

aTCAAAGACGCGAAGTTGATGCCCTTGCTGCAATAAATTGCTCGCCATTGGCGAACCCATTTGTCCTAAACCGATAAACGCGATTGCTGCCATAACCccct                                                                                                 >  1:7896709/1‑101 (MQ=255)

   aaaGACGCGAAGTTGATGCCCTTGCTGCAATAAATTGCTCGCCATTGGCGAACCCATTTGTCCTAAACCGATAAACGCGATTGCTGCCATAACCCcctcct                                                                                              >  1:7896714/1‑101 (MQ=255)

    aaGACGCGAAGTTGATGCCCTTGCTGCAATAAATTGCTCGCCATTGGCGAACCCATTTGTCCTAAACCGATAAACGCGATTGCTGCCATAACCCCCTCCTg                                                                                             <  2:7895982/101‑1 (MQ=255)

    aaGACGCGAAGTTGATGCCCTTGCTGCAATAAATTGCTCGCCATTGGCGAACCCATTTGTCCTAAACCGATAAACGCGATTGCTGCCATAACCCCCTCCTg                                                                                             <  1:7896075/101‑1 (MQ=255)

     aGACGCGAAGTTGATGCCCTTGCTGCAATAAATTGCTCGCCATTGGCGAACCCATTTGTCCTAAACCGATAAACGCGATTGCTGCCATAACCCCCTCCTGa                                                                                            <  1:7895871/101‑1 (MQ=255)

       aCGCGAAGTTGATGCCCTTGCTGCAATAAATTGCTCGCCATTGGCGAACCCATTTGTCCTAAACCGATAAACGCGATTGCTGCCATAACCCCCTCCTGAAt                                                                                          >  2:7896720/1‑101 (MQ=255)

        cgcgAAGTTGATGCCCTTGCTGCAATAAATTGCTCGCCATTGGCGAACCCATTTGTCCTAAACCGATAAACGCGATTGCTGCCATAACCCCCTCCTGAATa                                                                                         >  1:7896721/1‑101 (MQ=255)

          cgAAGTTGATGCCCTTGCTGCAATAAATTGCTCGCCATTGGCGAACCCATTTGTCCTAAACCGATAAACGCGATTGCTGCCATAACCCCCTCCTGAATACa                                                                                       >  1:7896725/1‑101 (MQ=255)

             aGTTGATGCCCTTGCTGCAATAAATTGCTCGCCATTGGCGAACCCATTTGTCCTAAACCGATAAACGCGATTGCTGCCATAACCCCCTCCTGAATACAGtt                                                                                    >  1:7896731/1‑101 (MQ=255)

               ttGATGCCCTTGCTGCAATAAATTGCTCGCCATTGGCGAACCCATTTGTCCTAAACCGATAAACGCGATTGCTGCCATAACCCCCTCCTGAATACAGTTAt                                                                                  >  2:7896735/1‑101 (MQ=255)

                 gATGCCCTTGCTGCAATAAATTGCTCGCCATTGGCGAACCCATTTGTCCTAAACCGATAAACGCGATTGCTGCCATAACCCCCTCCTGAATACAGTTATGt                                                                                <  2:7896277/101‑1 (MQ=255)

                  aTGCCCTTGCTGCAATAAATTGCTCGCCATTGGCGAACCCATTTGTCCTAAACCGATAAACGCGATTGCTGCCATAACCCCCTCCTGAATACAGTTATGTc                                                                               >  1:7896739/1‑101 (MQ=255)

                  aTGCCCTTGCTGCAATAAATTGCTCGCCATTGGCGAACCCATTTGTCCTAAACCGATAAACGCGATTGCTGCCATAACCCCCTCCTGAATACAGTTATGTc                                                                               >  2:7896742/1‑101 (MQ=255)

                   tGCCCTTGCTGCAATAAATTGCTCGCCATTGGCGAACCCATTTGTCCTAAACCGATAAACGCGATTGCTGCCATAACCCCCTCCTGAATACAGTTATGTCa                                                                              <  2:7896254/101‑1 (MQ=255)

                         tgctgcAATAAATTGCTCGCCATTGGCGAACCCATTTGTCCTAAACCGATAAACGCGATTGCTGCCATAACCCCCTCCTGAATACAGTTATGTCACttttt                                                                        <  1:7895980/101‑1 (MQ=255)

                          gctgcAATAAATTGCTCGCCATTGGCGAACCCATTTGTCCTAAACCGATAAACGCGATTGCTGCCATAACCCCCTCCTGAATACAGTTATGTCACTTTTTg                                                                       >  1:7896748/1‑101 (MQ=255)

                            tgcAATAAATTGCTCGCCATTGGCGAACCCATTTGTCCTAAACCGATAAACGCGATTGCTGCCATAACCCCCTCCTGAATACAGTTATGTCACTTTTTGTc                                                                     <  1:7896250/101‑1 (MQ=255)

                            tgcAATAAATTGCTCGCCATTGGCGAACCCATTTGTCCTAAACCGATAAACGCGATTGCTGCCATAACCCCCTCCTGAATACAGTTATGTCACTTTTTGTc                                                                     <  1:7896109/101‑1 (MQ=255)

                            tgcAATAAATTGCTCGCCATTGGCGAACCCATTTGTCCTAAACCGATAAACGCGATTGCTGCCATAACCCCCTCCTGAATACAGTTATGTCACTTTTTGTc                                                                     <  1:7896028/101‑1 (MQ=255)

                            tgcAATAAATTGCTCGCCATTGGCGAACCCATTTGTCCTAAACCGATAAACGCGATTGCTGCCATAACCCCCTCCTGAATACAGTTATGTCACTTTTTGTc                                                                     <  1:7895728/101‑1 (MQ=255)

                            tgcAATAAATTGCTCGCCATTGGCGAACCCATTTGTCCTAAACCGATAAACGCGATTGCTGCCATAACCCCCTCCTGAATACAGTTATGTCACTTTTTGTc                                                                     <  2:7896352/101‑1 (MQ=255)

                            tgcAATAAATTGCTCGCCATTGGCGAACCCATTTGTCCTAAACCGATAAACGCGATTGCTGCCATAACCCCCTCCTGAATACAGTTATGTCACTTTTTGTc                                                                     <  1:7896386/101‑1 (MQ=255)

                             gcAATAAATTGCTCGCCATTGGCGAACCCATTTGTCCTAAACCGATAAACGCGATTGCTGCCATAACCCCCTCCTGAATACAGTTATGTCACTTTTTGTCa                                                                    >  1:7896757/1‑101 (MQ=255)

                              cAATAAATTGCTCGCCATTGGCGAACCCATTTGTCCTAAACCGATAAACGCGATTGCTGCCATAACCCCCTCCTGAATACAGTTATGTCACTTTTTGTCAt                                                                   >  2:7896761/1‑101 (MQ=255)

                               aaTAAATTGCTCGCCATTGGCGAACCCATTTGTCCTAAACCGATAAACGCGATTGCTGCCATAACCCCCTCCTGAATACAGTTATGTCACTTTTTGTCAtt                                                                  <  2:7896186/101‑1 (MQ=255)

                                aTAAATTGCTCGCCATTGGCGAACCCATTTGTCCTAAACCGATAAACGCGATTGCTGCCATAACCCCCTCCTGAATACAGTTATGTCACTTTTTGTCAttt                                                                 >  2:7896767/1‑101 (MQ=255)

                                aTAAATTGCTCGCCATTGGCGAACCCATTTGTCCTAAACCGATAAACGCGATTGCTGCCATAACCCCCTCCTGAATACAGTTATGTCACTTTTTGTCAttt                                                                 <  1:7896126/101‑1 (MQ=255)

                                 tAAATTGCTCGCCATTGGCGAACCCATTTGTCCTAAACCGATAAACGCGATTGCTGCCATAACCCCCTCCTGAATACAGTTATGTCACTTTTTGTCATTTa                                                                >  2:7896772/1‑101 (MQ=255)

                                 tAAATTGCTCGCCATTGGCGAACCCATTTGTCCTAAACCGATAAACGCGATTGCTGCCATAACCCCCTCCTGAATACAGTTATGTCACTTTTTGTCATTTa                                                                >  1:7896769/1‑101 (MQ=255)

                                   aaTTGCTCGCCATTGGCGAACCCATTTGTCCTAAACCGATAAACGCGATTGCTGCCATAACCCCCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATg                                                              <  2:7896369/101‑1 (MQ=255)

                                    aTTGCTCGCCATTGGCGAACCCATTTGTCCTAAACCGATAAACGCGATTGCTGCCATAACCCCCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGa                                                             <  2:7896146/101‑1 (MQ=255)

                                    aTTGCTCGCCATTGGCGAACCCATTTGTCCTAAACCGATAAACGCGATTGCTGCCATAACCCCCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGa                                                             <  2:7896317/101‑1 (MQ=255)

                                    aTTGCTCGCCATTGGCGAACCCATTTGTCCTAAACCGATAAACGCGATTGCTGCCATAACCCCCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGa                                                             >  1:7896776/1‑101 (MQ=255)

                                     ttGCTCGCCATTGGCGAACCCATTTGTCCTAAACCGATAAACGCGATTGCTGCCATAACCCCCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGAc                                                            <  1:7895940/101‑1 (MQ=255)

                                        cTCGCCATTGGCGAACCCATTTGTCCTAAACCGATAAACGCGATTGCTGCCATAACCCCCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACAtg                                                         >  1:7896783/1‑101 (MQ=255)

                                        cTCGCCATTGGCGAACCCATTTGTCCTAAACCGATAAACGCGATTGCTGCCATAACCCCCGCCTGAATACAGTTAt                                                                                  <  2:7895949/76‑1 (MQ=255)

                                         tCGCCATTGGCGAACCCATTTGTCCTAAACCGATAAACGCGATTGCTGCCATAACCCCCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACAtgc                                                        >  2:7896785/1‑101 (MQ=255)

                                           gCCATTGGCGAACCCATTTGTCCTAAACCGATAAACGCGATTGCTGCCATAACCCCCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACAtgctt                                                      >  1:7896789/1‑101 (MQ=255)

                                             cATTGGCGAACCCATTTGTCCTAAACCGATAAACGCGATTGCTGCCATAACCCCCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACAtgctttg                                                    >  2:7896794/1‑101 (MQ=255)

                                             cATTGGCGAACCCATTTGTCCTAAACCGATAAACGCGATTGCTGCCATAACCCCCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACAtgctttg                                                    <  2:7896012/101‑1 (MQ=255)

                                              aTTGGCGAACCCATTTGTCCTAAACCGATAAACGCGATTGCTGCCATAACCCCCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACAtgctttgc                                                   >  2:7896799/1‑101 (MQ=255)

                                               ttGGCGAACCCATTTGTCCTAAACCGATAAACGCGATTGCTGCCATAACCCCCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACAtgctttgct                                                  >  1:7896800/1‑101 (MQ=255)

                                                tGGCGAACCCATTTGTCCTAAACCGATAAACGCGATTGCTGCCATAACCCCCTCCTGAATAc                                                                                        <  2:7896294/62‑1 (MQ=255)

                                                tGGCGAACCCATTTGTCCTAAACCGATAAACGCGATTGCTGCCATAACCCCCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACAtgctttgctt                                                 <  1:7896590/101‑1 (MQ=255)

                                                tGGCGAACCCATTTGTCCTAAACCGATAAACGCGATTGCTGCCATAACCCCCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACAtgctttgctt                                                 <  2:7896323/101‑1 (MQ=255)

                                                    gAACCCATTTGTCCTAAACCGATAAACGCGATTGCTGCCATAACCCCCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCt                                             <  2:7896023/101‑1 (MQ=255)

                                                    gAACCCATTTGTCCTAAACCGATAAACGCGATTGCTGCCATAACCCCCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCt                                             >  2:7896808/1‑101 (MQ=255)

                                                     aaCCCATTTGTCCTAAACCGATAAACGCGATTGCTGCCATAACCCCCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTg                                            <  2:7896125/101‑1 (MQ=255)

                                                     aaCCCATTTGTCCTAAACCGATAAACGCGATTGCTGCCATAACCCCCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTg                                            >  1:7896809/1‑101 (MQ=255)

                                                      aCCCATTTGTCCTAAACCGATAAACGCGATTGCTGCCATAACCCCCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGt                                           <  2:7896439/101‑1 (MQ=255)

                                                      aCCCATTTGTCCTAAACCGATAAACGCGATTGCTGCCATAACCCCCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGt                                           <  1:7896353/101‑1 (MQ=255)

                                                       cccATTTGTCCTAAACCGATAAACGCGATTGCTGCCATAACCCCCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGtt                                          <  2:7896165/101‑1 (MQ=255)

                                                       cccATTTGTCCTAAACCGATAAACGCGATTGCTGCCATAACCCCCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGtt                                          <  2:7896265/101‑1 (MQ=255)

                                                         cATTTGTCCTAAACCGATAAACGCGATTGCTGCCATAACCCCCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGtttt                                        >  2:7896814/1‑101 (MQ=255)

                                                          aTTTGTCCTAAACCGATAAACGCGATTGCTGCCATAACCCCCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGttttt                                       <  2:7896236/101‑1 (MQ=255)

                                                            ttGTCCTAAACCGATAAACGCGATTGCTGCCATAACCCCCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGa                                     <  2:7896339/101‑1 (MQ=255)

                                                            ttGTCCTAAACCGATAAACGCGATTGCTGCCATAACCCCCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGa                                     >  1:7896822/1‑101 (MQ=255)

                                                            ttGTCCTAAACCGATAAACGCGATTGCTGCCATAACCCCCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGa                                     >  2:7896823/1‑101 (MQ=255)

                                                             tGTCCTAAACCGATAAACGCGATTGCTGCCATAACCCCCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGAt                                    <  1:7896467/101‑1 (MQ=255)

                                                             tGTCCTAAACCGATAAACGCGATTGCTGCCATAACCCCCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGAt                                    <  1:7896432/101‑1 (MQ=255)

                                                              gTCCTAAACCGATAAACGCGATTGCTGCCATAACCCCCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTg                                      >  1:7896831/1‑98 (MQ=255)

                                                              gTCCTAAACCGATAAACGCGATTGCTGCCATAACCCCCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATc                                   >  2:7896830/1‑101 (MQ=255)

                                                              gTCCTAAACCGATAAACGCGATTGCTGCCATAACCCCCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATc                                   <  2:7896106/101‑1 (MQ=255)

                                                              gTCCTAAACCGATAAACGCGATTGCTGCCATAACCCCCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATc                                   <  1:7896229/101‑1 (MQ=255)

                                                               tCCTAAACCGATAAACGCGATTGCTGCCATAACCCCCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATc                                   <  2:7895898/100‑1 (MQ=255)

                                                                ccTAAACCGATAAACGCGATTGCTGCCATAACCCCCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGt                                 <  1:7896154/101‑1 (MQ=255)

                                                                 cTAAACCGATAAACGCGATTGCTGCCATAACCCCCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGTa                                <  2:7896639/101‑1 (MQ=255)

                                                                  tAAACCGATAAACGCGATTGCTGCCATAACCCCCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGTAt                               >  1:7896836/1‑101 (MQ=255)

                                                                   aaaCCGATAAACGCGATTGCTGCCATAACCCCCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGTAtt                              >  2:7896839/1‑101 (MQ=255)

                                                                      ccGATAAACGCGATTGCCGCCATAACCCCCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGTATTTGt                           >  1:7896845/1‑101 (MQ=255)

                                                                       cGATAAACGCGATTGCTGCCATAACCCCCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGTATTTGTa                          >  2:7896847/1‑101 (MQ=255)

                                                                       cGATAAACGCGATTGCTGCCATAACCCCCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGTATTTGTa                          <  1:7896365/101‑1 (MQ=255)

                                                                          tAAACGCGATTGCTGCCATAACCCCCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGTATTTGTAAtt                       <  2:7896399/101‑1 (MQ=255)

                                                                           aaaCGCGATTGCTGCCATAACCCCCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGTATTTGTAAttt                      <  2:7896209/101‑1 (MQ=255)

                                                                           aaaCGCGATTGCTGCCATAACCCCCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGTATTTGTAAttt                      <  1:7896199/101‑1 (MQ=255)

                                                                           aaaCGCGATTGCTGCCATAACCCCCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGTATTTGTAAttt                      >  1:7896852/1‑101 (MQ=255)

                                                                             aCGCGATTGCTGCCATAACCCCCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGTATTTGTAATTTAt                    <  2:7896091/101‑1 (MQ=255)

                                                                             aCGCGATTGCTGCCATAACCCCCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGTATTTGTAATTTAt                    <  1:7896182/101‑1 (MQ=255)

                                                                              cgcgATTGCTGCCATAACCCCCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGTATTTGTAATTTATc                   >  2:7896856/1‑101 (MQ=255)

                                                                               gcgATTGCTGCCATAACCCCCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGTATTTGTAATTTATCg                  <  1:7896622/101‑1 (MQ=255)

                                                                               gcgATTGCTGCCATAACCCCCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGTATTTGTAATTTATCg                  <  1:7896299/101‑1 (MQ=255)

                                                                                  aTTGCTGCCATAACCCCCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGTATTTGTAATTTATCGTCa               >  1:7896862/1‑101 (MQ=255)

                                                                                   ttGCTGCCATAACCCCCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGTATTTGTAATTTATCGTCaa              >  2:7896865/1‑101 (MQ=255)

                                                                                   ttGCTGCCATAACCCCCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGTATTTGTAATTTATCGTCaa              <  2:7896490/101‑1 (MQ=255)

                                                                                     gctgcCATAACCCCCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGTATTTGTAATTTATCGTCaaaa            <  1:7896325/101‑1 (MQ=255)

                                                                                     gctgcCATAACCCCCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGTATTTGTAATTTATCGTCaaaa            <  2:7896393/101‑1 (MQ=255)

                                                                                     gctgcCATAACCCCCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGTATTTGTAATTTATCGTCaaaa            >  1:7896869/1‑101 (MQ=255)

                                                                                      ctgcCATAACCCCCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGTATTTGTAATTTATCGTCaaaaa           <  1:7896489/101‑1 (MQ=255)

                                                                                       tgcCATAACCCCCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGTATTTGTAATTTATCGTCaaaaaa          >  2:7896873/1‑101 (MQ=255)

                                                                                        gcCATAACCCCCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGTATTTGTAATTTATCGTCAAAAAAt         >  1:7896876/1‑101 (MQ=255)

                                                                                         ccATAACCCCCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGTATTTGTAATTTATCGTCAAAAAAtt        <  1:7896285/101‑1 (MQ=255)

                                                                                          cATAACCCCCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGTATTTGTAATTTATCGTCAAAAAATTg       <  2:7896418/101‑1 (MQ=255)

                                                                                          cATAACCCCCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGTATTTGTAATTTATCGTCAAAAAATTg       <  1:7896268/101‑1 (MQ=255)

                                                                                            tAACCCCCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGTATTTGTAATTTATCGTCAAAAAATTGAc     >  2:7896884/1‑101 (MQ=255)

                                                                                            tAACCCCCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGTATTTGTAATTTATCGTCAAAAAATTGAc     >  1:7896882/1‑101 (MQ=255)

                                                                                             aaCCCCCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGTATTTGTAATTTATCGTCAAAAAATTGACa    <  2:7896222/101‑1 (MQ=255)

                                                                                             aaCCCCCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGTATTTGTAATTTATCGTCAAAAAATTGACa    <  2:7896558/101‑1 (MQ=255)

                                                                                              aCCCCCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGTATTTGTAATTTATCGTCAAAAAATTGACa    >  1:7896889/1‑100 (MQ=255)

                                                                                              aCCCCCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGTATTTGTAATTTATCGTCAAAAAATTGACAg   <  1:7896408/101‑1 (MQ=255)

                                                                                               ccccCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGTATTTGTAATTTATCGTCAAAAAATTGACAGc  <  2:7896753/101‑1 (MQ=255)

ATCAAAGACGCGAAGTTGATGCCCTTGCTGCAATAAATTGCTCGCCATTGGCGAACCCATTTGTCCTAAACCGATAAACGCGATTGCTGCCATAACCCTCTCCTGAATACAGTTATGTCACTTTTTGTCATTTATGACATGCTTTGCTTGTCTGTTTTTGATCGTATTTGTAATTTATCGTCAAAAAATTGACAGC  >  NC_000913/4073479‑4073674

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: `‑`