BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,615,859	T→G	100%	T403P (ACC→CCC)	yjjI ←	DUF3029 family protein, putative glycine radical enzyme

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,615,859	0	T	G	100.0%	512.2 / ‑6.7	147	T403P (ACC→CCC)	yjjI	DUF3029 family protein, putative glycine radical enzyme
Reads supporting (aligned to +/- strand): ref base T (0/0); major base G (70/76); minor base A (1/0); total (71/76)

Only 100 of 148 total aligned reads displayed.

GCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCG  >  NC_000913/4615762‑4615952

gCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGa                                                                                             <  2:3787787/100‑1 (MQ=255)

  cTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATc                                                                                           >  1:3586819/1‑100 (MQ=255)

    gCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCgg                                                                                         >  2:1706121/1‑100 (MQ=255)

    gCGAGTTCCGTATGATGGTTTCGTCGAGCGTCGGAATGTCGCTTAGGCCGGGATAATGGGGAGGATTATGCGTGGCGAGAGTTTGGATGGGGGGGATCgg                                                                                         >  2:3477527/1‑100 (MQ=255)

     cGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGGCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGa                                                                                        <  2:72887/100‑1 (MQ=255)

      gggtgccGTTTGACGGCTTCGTCGCGCCGCAGTATGTCGCTGATGCCGGAGTAACAATAAGCATGAAGCGCGGCGACAGTTTACAGATGGGGGATCGGan                                                                                       <  2:1004705/94‑2 (MQ=255)

       ggTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGggatcggatc                                                                                      <  1:1788343/100‑1 (MQ=255)

        gTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCt                                                                                     >  1:665704/1‑100 (MQ=255)

        gTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCt                                                                                     >  2:1923763/1‑100 (MQ=255)

         ttACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTg                                                                                    >  2:2905817/1‑100 (MQ=255)

           acattttagggTTGCGTCGAGCGTCAGAATGTCGCTGATGAGGAAATAATAATAAGCATGATGCGGGGTGACAGTTTGTAGATGGGGGATGGAATCTGGc                                                                                  <  2:2656086/91‑1 (MQ=255)

           aCGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGGTCTGGc                                                                                  <  2:115578/100‑1 (MQ=255)

           aCGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGc                                                                                  <  2:2753353/100‑1 (MQ=255)

           aCGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGc                                                                                  <  1:2599265/100‑1 (MQ=255)

            cGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCt                                                                                 >  2:2837681/1‑100 (MQ=255)

            cGTTTGATGGTTTCGGCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCt                                                                                 <  1:3947820/100‑1 (MQ=255)

             gTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTc                                                                                >  1:2641313/1‑100 (MQ=255)

             gTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGACGGAATAATAATAAGCATGATGCGGGGCGACCGGTTGCAGAGAGGGGATCGGCTCTg                                                                                    >  2:2561649/1‑96 (MQ=255)

              tttGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCa                                                                               >  2:3332211/1‑100 (MQ=255)

                 gATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCg                                                                            >  1:4169142/1‑100 (MQ=255)

                 gATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCg                                                                            >  1:3928212/1‑100 (MQ=255)

                  aTGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGc                                                                           <  1:3154346/100‑1 (MQ=255)

                   tGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGcc                                                                          <  1:3416098/100‑1 (MQ=255)

                   tGGTGTCTTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGc                                                                           <  2:3466955/99‑1 (MQ=255)

                     gTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCAt                                                                        <  2:1373287/100‑1 (MQ=255)

                       ttCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATAc                                                                      <  2:4103008/100‑1 (MQ=255)

                          gtcgAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGc                                                                   <  2:1372512/100‑1 (MQ=255)

                          gtcgAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGc                                                                   <  2:1338809/100‑1 (MQ=255)

                           tcgAGCGTCAGAATGTCGCTGATGCCGGAATAATAAAAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCa                                                                  <  1:2933618/100‑1 (MQ=255)

                              aGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAAc                                                               <  1:1375554/100‑1 (MQ=255)

                               gCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAAcg                                                              <  1:664327/100‑1 (MQ=255)

                               gCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAAcg                                                              <  1:226869/100‑1 (MQ=255)

                                cGTCAGACTGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAAcgc                                                             >  2:63651/1‑100 (MQ=255)

                                cGTCAGAATGTCGCTGATGCAGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATAGGATCTGGCTCATCGCCATACGGCAAAagc                                                             >  2:3718314/1‑100 (MQ=255)

                                 gTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAAcgcg                                                            >  1:825914/1‑100 (MQ=255)

                                 gTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAAcgcg                                                            <  2:670509/100‑1 (MQ=255)

                                  tCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAAcgcgc                                                           >  1:1700518/1‑100 (MQ=255)

                                  tCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAAcgcgc                                                           >  1:748364/1‑100 (MQ=255)

                                  tCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGAGATCGGATCTGGCTCATCGCCATACGGCAAAcgcgc                                                           >  2:2297173/1‑100 (MQ=255)

                                   cAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAAcgcgcg                                                          >  2:970165/1‑100 (MQ=255)

                                   cAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAAcgcgcg                                                          <  2:661036/100‑1 (MQ=255)

                                   cAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAAcgcgcg                                                          <  1:3925915/100‑1 (MQ=255)

                                   cAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAAcgcgcg                                                          >  1:3707018/1‑100 (MQ=255)

                                     gAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGcc                                                        <  1:2866462/100‑1 (MQ=255)

                                         gTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGc                                                    >  1:2433836/1‑100 (MQ=255)

                                          tCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCg                                                   <  1:150390/100‑1 (MQ=255)

                                           cGCTGATGCCGGAATAATAATAAGCATGAAGCGGGTCGACTGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCgt                                                  >  2:1221293/1‑100 (MQ=255)

                                            gCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCgtg                                                 <  1:3382262/100‑1 (MQ=255)

                                             cTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCgtgg                                                <  1:2579690/100‑1 (MQ=255)

                                              tGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCgtgg                                                <  2:1604510/99‑1 (MQ=255)

                                                aTGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGc                                             <  1:466986/100‑1 (MQ=255)

                                                 tGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGcc                                            >  2:2162086/1‑100 (MQ=255)

                                                 tGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATAGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGcc                                            >  2:958646/1‑100 (MQ=255)

                                                  gCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCg                                           >  1:1086170/1‑100 (MQ=255)

                                                  gCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCg                                           <  1:1633116/100‑1 (MQ=255)

                                                  gCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCg                                           <  1:2855822/100‑1 (MQ=255)

                                                  gCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCg                                           >  1:3081209/1‑100 (MQ=255)

                                                   ccGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGa                                          <  1:2036421/100‑1 (MQ=255)

                                                    cGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCg                                           <  2:2015513/98‑1 (MQ=255)

                                                    cGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGat                                         <  2:1647345/100‑1 (MQ=255)

                                                         taataataaGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCgg                                    <  2:3512464/100‑1 (MQ=255)

                                                          aataataaGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGa                                   >  1:1095347/1‑100 (MQ=255)

                                                          aataataaGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGa                                   >  1:1095616/1‑100 (MQ=255)

                                                           ataataaGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGTTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGaa                                  >  2:1601158/1‑100 (MQ=255)

                                                           ataataaGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGaa                                  >  2:1526170/1‑100 (MQ=255)

                                                           ataataaGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGaa                                  >  2:2657860/1‑100 (MQ=255)

                                                           ataataaGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCAGACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGaa                                  >  2:2260784/1‑100 (MQ=255)

                                                            taataaGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAAc                                 <  2:1237292/100‑1 (MQ=255)

                                                             aataaGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACt                                <  1:2764913/100‑1 (MQ=255)

                                                             aataaGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACt                                >  2:4032906/1‑100 (MQ=255)

                                                               taaGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGTAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGa                              <  2:2229693/100‑1 (MQ=255)

                                                                 aGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGa                              >  2:52064/1‑98 (MQ=255)

                                                                  gCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATcc                           <  2:23620/100‑1 (MQ=255)

                                                                   cATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATccc                          >  1:4001023/1‑100 (MQ=255)

                                                                   cATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACCGATccc                          >  2:3254641/1‑100 (MQ=255)

                                                                     tgatgCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCg                        <  2:484729/100‑1 (MQ=255)

                                                                      gatgCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGa                       >  1:3538765/1‑100 (MQ=255)

                                                                        tgCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGAc                      <  1:3705745/99‑1 (MQ=255)

                                                                        tgCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACt                     >  1:222171/1‑100 (MQ=255)

                                                                         gCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACtg                    <  1:970169/100‑1 (MQ=255)

                                                                         gCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACtg                    <  1:970218/100‑1 (MQ=255)

                                                                         gCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACtg                    <  1:3364638/100‑1 (MQ=255)

                                                                           ggggCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGCTATCGGAACTGATCCCCGACtgtg                  >  1:1788856/1‑100 (MQ=255)

                                                                             ggCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCg                <  1:1432494/100‑1 (MQ=255)

                                                                                gACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCgtg              <  2:1401825/99‑1 (MQ=255)

                                                                                 aCAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTa            <  1:716082/100‑1 (MQ=255)

                                                                                 aCAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTa            <  2:735773/100‑1 (MQ=255)

                                                                                  cAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTaa           <  2:3357067/100‑1 (MQ=255)

                                                                                  cAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTaa           >  1:588343/1‑100 (MQ=255)

                                                                                    gTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACa         >  2:2260111/1‑100 (MQ=255)

                                                                                     tttGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACAt        <  1:2686041/100‑1 (MQ=255)

                                                                                     tttGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACAt        >  2:2524966/1‑100 (MQ=255)

                                                                                     tttGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACAt        >  2:3578782/1‑100 (MQ=255)

                                                                                      ttGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATg       <  2:2799259/100‑1 (MQ=255)

                                                                                      ttGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATg       >  1:840238/1‑100 (MQ=255)

                                                                                        gCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGgc     >  2:470150/1‑100 (MQ=255)

                                                                                         cAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGgct    >  2:1518012/1‑99 (MQ=255)

                                                                                           gATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGgcgc   >  2:1506958/1‑99 (MQ=255)

                                                                                           gATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGcgcg  >  1:53717/1‑100 (MQ=255)

                                                                                           gATGGGGGATCGAATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGcgcg  <  1:1180632/100‑1 (MQ=255)

GCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCG  >  NC_000913/4615762‑4615952

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