BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,551,460	T→C	A72A (GCT→GCC)	hemF →	coproporphyrinogen III oxidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,551,460	0	T	C	100.0%	764.8 / NA	224	A72A (GCT→GCC)	hemF	coproporphyrinogen III oxidase
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (104/120); total (104/120)

Only 100 of 238 total aligned reads displayed.

CGGCGGGCGTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGT  >  NC_000913/2551370‑2551549

 ggcgggcgTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTcc                                                                                        <  2:3438501/93‑1 (MQ=255)

  gcgggcgTAGTCGGGTGTTGCGTAATGGGGGGGTTTTTGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGGGAGGCGATGCCCGCCTccc                                                                                       >  1:2313900/1‑92 (MQ=255)

    gggcgTAGTCGGGTGTTGCGTAATGGTGGTGGTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGcc                                                                                     >  1:6463966/1‑93 (MQ=255)

     ggcgTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCa                                                                                    <  2:790045/93‑1 (MQ=255)

     ggcgTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCa                                                                                    >  2:2956850/1‑93 (MQ=255)

      gcgTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCAc                                                                                   <  2:5600758/93‑1 (MQ=255)

      gcgTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCAc                                                                                   <  2:5552448/93‑1 (MQ=255)

      gcgTAGTCGGGTGTTGCGTAATGGTGGGGCTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACCGTGAGGCGATGCCTGCCTCCGCCAc                                                                                   >  1:5773798/1‑93 (MQ=255)

        gTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCg                                                                                 >  2:76327/1‑93 (MQ=255)

         tAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGc                                                                                <  1:1162136/93‑1 (MQ=255)

          aGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCt                                                                               >  2:1408436/1‑93 (MQ=255)

          aGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCt                                                                               >  1:3032992/1‑93 (MQ=255)

          aGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCt                                                                               <  2:3612439/93‑1 (MQ=255)

           gTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCCCGGTGAGGCGATGCCTGCCTCCGCCACCGCTc                                                                              >  2:2878233/1‑93 (MQ=255)

           gTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTc                                                                              >  1:5997861/1‑93 (MQ=255)

           gTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTc                                                                              >  2:4733684/1‑93 (MQ=255)

             cGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCAt                                                                            <  1:1781478/93‑1 (MQ=255)

             cGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCAt                                                                            <  1:295732/93‑1 (MQ=255)

              gggTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATc                                                                           <  2:4144679/93‑1 (MQ=255)

              gggTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATc                                                                           >  2:428982/1‑93 (MQ=255)

              gggTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATc                                                                           >  1:5260654/1‑93 (MQ=255)

                 tgtTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGcc                                                                        <  2:6095155/93‑1 (MQ=255)

                  gtTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGccc                                                                       >  2:108310/1‑93 (MQ=255)

                       gTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAAc                                                                  <  2:1684297/93‑1 (MQ=255)

                        tAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACt                                                                 <  1:6278413/93‑1 (MQ=255)

                         aaTGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACtt                                                                <  2:2132947/93‑1 (MQ=255)

                         aaTGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACtt                                                                <  2:1862147/93‑1 (MQ=255)

                            ggtggtGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGcc                                                             <  1:975785/93‑1 (MQ=255)

                            ggtggtGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGcc                                                             >  1:1810772/1‑93 (MQ=255)

                             gtggtgTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCg                                                            <  2:3806113/93‑1 (MQ=255)

                             gtggtgTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCg                                                            >  1:184854/1‑93 (MQ=255)

                              tggtgTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCgg                                                           <  1:1059101/93‑1 (MQ=255)

                                gtgtTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGgc                                                         <  1:5655090/93‑1 (MQ=255)

                                  gtTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCGCGGAACTTGCCGGgcgc                                                       >  1:1375384/1‑93 (MQ=255)

                                  gtTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGgcgc                                                       >  2:2048757/1‑93 (MQ=255)

                                     ttCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGt                                                    >  2:2639231/1‑93 (MQ=255)

                                     ttCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGt                                                    <  1:758809/93‑1 (MQ=255)

                                      tCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGtt                                                   <  1:5909156/93‑1 (MQ=255)

                                       cGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGttt                                                  <  2:5199843/93‑1 (MQ=255)

                                         aaCAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCg                                                <  1:1946616/93‑1 (MQ=255)

                                          aCAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGa                                               >  2:6284579/1‑93 (MQ=255)

                                          aCAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGa                                               <  1:1226209/93‑1 (MQ=255)

                                          aCAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGa                                               >  2:5152406/1‑93 (MQ=255)

                                           caggcaggCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGACTGCCTCCGCCACCGCTCATCTCCCCGAACTTGCCGGGCGCAGTTTCGAg                                              >  2:2187290/1‑93 (MQ=255)

                                            aggcaggcGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAgg                                             >  2:5421160/1‑93 (MQ=255)

                                                aggcGTCAACTTTTCGCATGTCCATGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGACGGGCGCAGATTCGAGGCGAt                                         >  2:5896572/1‑93 (MQ=255)

                                                aggcGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGAt                                         >  2:4430845/1‑93 (MQ=255)

                                                aggcGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGAt                                         >  1:4983605/1‑93 (MQ=255)

                                                aggcGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGAt                                         >  1:4477427/1‑93 (MQ=255)

                                                 ggcGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATg                                        <  1:4056052/93‑1 (MQ=255)

                                                 ggcGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATg                                        >  2:1893494/1‑93 (MQ=255)

                                                 ggcGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATg                                        >  2:3725617/1‑93 (MQ=255)

                                                 gccGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATg                                        <  2:136415/91‑1 (MQ=255)

                                                   ggTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATggg                                      <  2:4012198/92‑1 (MQ=255)

                                                   cGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATggg                                      <  1:2759365/93‑1 (MQ=255)

                                                   cGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATggg                                      <  1:5029964/93‑1 (MQ=255)

                                                   cGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATggg                                      <  2:3891372/93‑1 (MQ=255)

                                                    gTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGc                                     >  2:3299565/1‑93 (MQ=255)

                                                    gTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGc                                     >  2:2788925/1‑93 (MQ=255)

                                                    gTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTACCGGGCGCAGTTTCGAGGCGATgcgg                                     >  2:2497809/1‑90 (MQ=255)

                                                      cAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGt                                   >  1:1001366/1‑93 (MQ=255)

                                                      cAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGt                                   <  1:3670885/93‑1 (MQ=255)

                                                      cAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGt                                   <  1:5347314/93‑1 (MQ=255)

                                                       aaCTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGtt                                  >  1:2139079/1‑93 (MQ=255)

                                                       aaCTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGtt                                  >  1:4891359/1‑93 (MQ=255)

                                                       aaCTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGCGCGCAGTTTAGCGGCGGTGGGCGtt                                  >  1:2246492/1‑93 (MQ=255)

                                                        aCTTTTCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGGTTCGAGGCGATGGGCGttt                                 >  2:4397517/1‑93 (MQ=255)

                                                         cTTTTCGCATGTCCACGATGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTc                                <  2:2375382/93‑1 (MQ=255)

                                                            ttCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACt                             <  2:6068933/93‑1 (MQ=255)

                                                            ttCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACt                             <  2:1144395/93‑1 (MQ=255)

                                                            ttCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACt                             <  2:3000500/93‑1 (MQ=255)

                                                            ttCGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACt                             <  2:2695183/93‑1 (MQ=255)

                                                              cGCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTgg                           <  2:4235186/93‑1 (MQ=255)

                                                               gCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGt                          >  2:1455893/1‑93 (MQ=255)

                                                               gCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGt                          >  1:5950845/1‑93 (MQ=255)

                                                               gCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGt                          <  2:4810795/93‑1 (MQ=255)

                                                               gCATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGt                          <  1:780948/93‑1 (MQ=255)

                                                                cATGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGtg                         >  1:330075/1‑92 (MQ=255)

                                                                 aTGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAg                        <  1:3858361/93‑1 (MQ=255)

                                                                  tGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGt                       >  2:6289169/1‑93 (MQ=255)

                                                                  tGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGt                       >  1:492157/1‑93 (MQ=255)

                                                                  tGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGt                       >  1:4746882/1‑93 (MQ=255)

                                                                  tGTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGt                       >  2:5279811/1‑93 (MQ=255)

                                                                   gTCCACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTg                      <  1:2200181/93‑1 (MQ=255)

                                                                     ccACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCa                    >  1:3577561/1‑93 (MQ=255)

                                                                     ccACGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCa                    >  2:1317907/1‑93 (MQ=255)

                                                                       aCGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATc                  <  2:4907708/93‑1 (MQ=255)

                                                                        cGGTGAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATcc                 <  1:2694862/93‑1 (MQ=255)

                                                                            gAGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCAt             <  1:2802720/93‑1 (MQ=255)

                                                                             aGGCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATa            >  1:3939838/1‑93 (MQ=255)

                                                                               gCGATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAAc          <  1:3828195/93‑1 (MQ=255)

                                                                                 gATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAAccc        >  1:4771340/1‑93 (MQ=255)

                                                                                 gATGCCTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAAccc        <  1:2517797/93‑1 (MQ=255)

                                                                                   tgcctgccTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGt      <  2:6194331/93‑1 (MQ=255)

                                                                                    gcctgcctCCGCCACCGCTCATCGCCCTGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTa     <  1:4302654/93‑1 (MQ=255)

                                                                                     ccTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATTGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTAt    <  2:4650851/93‑1 (MQ=255)

                                                                                     ccTGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTAt    <  2:1220531/93‑1 (MQ=255)

                                                                                       tGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGt  <  2:5994270/93‑1 (MQ=255)

                                                                                       tGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGt  >  1:1095887/1‑93 (MQ=255)

                                                                                       tGCCTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGATTCACTGGTAGTGCATCCGCATAACCCGTATGt  >  1:6076110/1‑93 (MQ=255)

CGGCGGGCGTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGT  >  NC_000913/2551370‑2551549

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