BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,791,400	T→C	Q61Q (CAA→CAG)	btuD ←	ATP‑binding component of vitamin B12 transport system

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,791,400	0	T	C	100.0%	629.9 / NA	197	Q61Q (CAA→CAG)	btuD	ATP‑binding component of vitamin B12 transport system
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (91/106); total (91/106)

Only 100 of 204 total aligned reads displayed.

GCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGTTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCA  >  NC_000913/1791309‑1791488

gCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGtt                                                                                         <  2:3153062/93‑1 (MQ=255)

 ccAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTg                                                                                        <  1:1333449/93‑1 (MQ=255)

   aGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTggc                                                                                      >  2:100379/1‑91 (MQ=255)

   aGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGcc                                                                                      >  2:5579291/1‑93 (MQ=255)

   aGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGcc                                                                                      >  1:2730827/1‑93 (MQ=255)

     aCCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGcccc                                                                                    >  2:2407567/1‑93 (MQ=255)

     aCCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGcccc                                                                                    <  1:2888730/93‑1 (MQ=255)

     aCCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGcccc                                                                                    <  1:1533954/93‑1 (MQ=255)

       cGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCcgc                                                                                  >  2:2908026/1‑93 (MQ=255)

          cGTTGCAACCGGCGGCGTCTGCTGTTGTGAAAGATAGGCACGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGaa                                                                               >  1:2049359/1‑93 (MQ=255)

          cGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGaa                                                                               >  1:3688141/1‑93 (MQ=255)

           gTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAAc                                                                              <  2:2808270/93‑1 (MQ=255)

              gCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGa                                                                           >  2:2286633/1‑93 (MQ=255)

               cAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGaa                                                                          >  1:1631643/1‑93 (MQ=255)

               cAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGaa                                                                          >  1:2449014/1‑93 (MQ=255)

                aaaCGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAAt                                                                         <  1:3819044/93‑1 (MQ=255)

                aaaCGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAAt                                                                         >  1:2851525/1‑93 (MQ=255)

                 aaCGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATg                                                                        <  1:3453436/93‑1 (MQ=255)

                  aCGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGc                                                                       >  2:3995134/1‑93 (MQ=255)

                  aCGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGc                                                                       >  2:4856870/1‑93 (MQ=255)

                   cggcggCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCt                                                                      <  1:4453238/93‑1 (MQ=255)

                    ggcggcGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCtt                                                                     >  2:68597/1‑93 (MQ=255)

                     gcggcgTCTGCTGTTGTGAAAGATAGGCGCTATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTc                                                                    <  1:901244/93‑1 (MQ=255)

                     gcggcgTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTc                                                                    >  1:958953/1‑93 (MQ=255)

                       ggcgTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTccc                                                                  >  2:5204356/1‑93 (MQ=255)

                       gacgTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTccc                                                                  <  2:5755735/91‑1 (MQ=255)

                        gcgTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCt                                                                 <  1:5317613/93‑1 (MQ=255)

                        gcgTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCt                                                                 <  1:1418911/93‑1 (MQ=255)

                        gcgTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCt                                                                 <  2:4355401/93‑1 (MQ=255)

                          gTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTa                                                               >  2:5484251/1‑93 (MQ=255)

                          gTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCCTCCCTTa                                                               >  1:6315398/1‑93 (MQ=255)

                           tCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTAc                                                              >  2:2699750/1‑93 (MQ=255)

                           tCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTAc                                                              <  1:1601377/93‑1 (MQ=255)

                           tCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTAc                                                              >  1:5023714/1‑93 (MQ=255)

                             tgctgTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCg                                                            >  2:5151454/1‑93 (MQ=255)

                             tgctgTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCg                                                            <  1:598170/93‑1 (MQ=255)

                              gctgTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGc                                                           <  2:704842/93‑1 (MQ=255)

                                  ttgtGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGt                                                       >  1:2593086/1‑93 (MQ=255)

                                   tgtgAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTc                                                      <  1:466889/93‑1 (MQ=255)

                                    gtgAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGCACTGAATGCTTCCCTTACCGCTTGTCa                                                     >  1:5484764/1‑93 (MQ=255)

                                     tgAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCAt                                                    <  2:3360430/93‑1 (MQ=255)

                                     tgAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCAt                                                    >  1:755741/1‑93 (MQ=255)

                                        aaGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTcc                                                 <  2:4656944/93‑1 (MQ=255)

                                        aaGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTcc                                                 >  2:2991505/1‑93 (MQ=255)

                                        aaGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTcc                                                 <  1:6418134/93‑1 (MQ=255)

                                        aaGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTcc                                                 <  1:67420/93‑1 (MQ=255)

                                         aGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCg                                                <  2:5852598/93‑1 (MQ=255)

                                         aGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCg                                                >  2:6084722/1‑93 (MQ=255)

                                          gATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCgg                                               <  1:3377242/93‑1 (MQ=255)

                                            tAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGcc                                             <  2:5082465/93‑1 (MQ=255)

                                              ggCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCAt                                           <  2:1503476/93‑1 (MQ=255)

                                              ggCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGGCATTCCGGCCAt                                           >  1:5346359/1‑93 (MQ=255)

                                               gcgcGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCAtt                                          <  2:3850688/93‑1 (MQ=255)

                                               gcgcGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCAtt                                          <  1:4882130/93‑1 (MQ=255)

                                                cgcgATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTc                                         <  2:3700125/93‑1 (MQ=255)

                                                cgcgATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTc                                         <  2:6506725/93‑1 (MQ=255)

                                                 gcgATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTcg                                        <  2:5174944/93‑1 (MQ=255)

                                                 gcgATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTcg                                        >  2:6375252/1‑93 (MQ=255)

                                                  cgATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTcgc                                       <  1:2434180/93‑1 (MQ=255)

                                                  cgATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTcgc                                       <  2:2465018/93‑1 (MQ=255)

                                                   gATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTcgcg                                      >  1:1505555/1‑93 (MQ=255)

                                                     tGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGcc                                    >  1:936707/1‑93 (MQ=255)

                                                     tGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGcc                                    >  1:404292/1‑93 (MQ=255)

                                                        aGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTTGCGCCAGt                                 >  2:272244/1‑93 (MQ=255)

                                                        aGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGt                                 >  1:6635960/1‑93 (MQ=255)

                                                         gcgcGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGtc                                >  1:5140343/1‑92 (MQ=255)

                                                          cgcgAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTaa                               <  2:3059135/93‑1 (MQ=255)

                                                          cgcgAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTaa                               >  2:3208334/1‑93 (MQ=255)

                                                             gAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGt                            >  2:4872991/1‑93 (MQ=255)

                                                             gAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGt                            <  2:393541/93‑1 (MQ=255)

                                                              aGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTa                           >  1:2070898/1‑93 (MQ=255)

                                                              aGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCAGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTa                           <  1:6448066/93‑1 (MQ=255)

                                                               gTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTAc                          <  2:1846342/93‑1 (MQ=255)

                                                                ttttGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTAct                         >  1:3845741/1‑93 (MQ=255)

                                                                 tttGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTActc                        >  1:1361840/1‑93 (MQ=255)

                                                                 tttGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTActc                        <  1:1272582/93‑1 (MQ=255)

                                                                 tttGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTActc                        <  1:3252178/93‑1 (MQ=255)

                                                                 tttGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTActc                        <  1:3084862/93‑1 (MQ=255)

                                                                  ttgttgCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTActct                       <  2:6138070/93‑1 (MQ=255)

                                                                  ttgttgCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTActct                       <  2:1363642/93‑1 (MQ=255)

                                                                    gttgCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTa                     <  1:2923567/93‑1 (MQ=255)

                                                                     ttgCGGGCCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTAc                    <  1:5159352/93‑1 (MQ=255)

                                                                      tgCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTAcc                   <  1:3912918/93‑1 (MQ=255)

                                                                       gCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTAATCTTAccc                  <  1:6564649/93‑1 (MQ=255)

                                                                         ggACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCcgc                >  1:4364853/1‑93 (MQ=255)

                                                                         ggACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCcgc                >  1:2194843/1‑93 (MQ=255)

                                                                          gACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCcgcg               <  1:4152849/93‑1 (MQ=255)

                                                                          gACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACAcgcg               >  1:6065180/1‑93 (MQ=255)

                                                                            ccATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGcc             >  2:3759614/1‑93 (MQ=255)

                                                                            ccATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGcc             >  2:1272725/1‑93 (MQ=255)

                                                                             cATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCa            >  1:1059350/1‑93 (MQ=255)

                                                                              aTGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCAt           >  2:4137325/1‑93 (MQ=255)

                                                                              aTGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCAt           <  2:2524048/93‑1 (MQ=255)

                                                                              aTGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCAt           <  1:352699/93‑1 (MQ=255)

                                                                              aTGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCAt           <  1:1164641/93‑1 (MQ=255)

                                                                               tGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCAtt          >  1:1865582/1‑93 (MQ=255)

                                                                                gCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTc         <  1:5703583/93‑1 (MQ=255)

                                                                                 cTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCg        >  1:6010525/1‑93 (MQ=255)

                                                                                    ccAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGGACTCTTACCCGCGCCATTCCGcc     >  2:4095969/1‑93 (MQ=255)

                                                                                       gTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGccccc  >  1:1848344/1‑92 (MQ=255)

GCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGTTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCA  >  NC_000913/1791309‑1791488

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