BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,008,469	+G	intergenic (+36/‑197)	b4555 → / → yigM	hypothetical protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,008,466	1	.	G	100.0%	539.9 / NA	145	intergenic (+33/‑200)	b4555/yigM	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (0/0); new base G (78/67); total (78/67)

Only 100 of 147 total aligned reads displayed.

AACTGGCTACAG‑CTGAAGGATTG‑ATTGCGTT‑‑GTGACGGGATGAACGACGCCGAAATGCTGTCGATGGC‑‑GGGAAA‑G‑CTGCATTATGGGCAGTGGCCACCAGCGTCTGAAAGACCTTCATCCCGAGCT‑GAAGTGATTGGTACTAATGCCGACGAC  >  NC_000913/4008399‑4008551

actgggCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCG‑GGGAAAGG‑CTGCATTATggg                                                                     <  1:6583209/89‑1 (MQ=255)

actgggCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCG‑GGGAAAGG‑CTGCATTATggg                                                                     <  2:2992138/89‑1 (MQ=255)

 ctgggCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGc                                                                    >  1:4536471/4‑93 (MQ=255)

 ctgggCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGc                                                                    >  1:2280110/4‑93 (MQ=255)

 ctgggCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGc                                                                    <  2:4408754/90‑1 (MQ=255)

  tgggCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCa                                                                   >  2:713210/3‑93 (MQ=255)

   gggCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAg                                                                  >  2:3677025/2‑93 (MQ=255)

   gggCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAg                                                                  <  2:2523121/92‑1 (MQ=255)

    ggCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGt                                                                 >  2:6397383/1‑93 (MQ=255)

     gCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTg                                                                <  1:2520996/93‑1 (MQ=255)

      cTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTgc                                                               <  2:795572/93‑2 (MQ=255)

       tACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTgcg                                                              >  2:1902341/1‑91 (MQ=255)

       tACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTgcg                                                              >  2:3897677/1‑91 (MQ=255)

       tACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTgcg                                                              >  2:3111292/1‑91 (MQ=255)

           gccTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGGAAG‑G‑CTGGATTATGGGCAGTGCgcacc                                                          >  1:5604860/3‑93 (MQ=255)

           gccTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCgcacc                                                          >  2:5136573/3‑93 (MQ=255)

           gccTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCgcacc                                                          >  2:5720606/3‑93 (MQ=255)

           gccTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCgcacc                                                          >  2:3643001/3‑93 (MQ=255)

             cTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCGCACCAg                                                        <  1:5836693/93‑1 (MQ=255)

             cTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCGCACCAg                                                        <  1:4041548/93‑1 (MQ=255)

              tGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCGCACCAGc                                                       >  1:866450/1‑93 (MQ=255)

              tGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCGCACCAGc                                                       >  1:4922669/1‑93 (MQ=255)

              tGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCGCACCAGc                                                       <  1:5078385/93‑1 (MQ=255)

              tGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCGCACCAGc                                                       >  2:2011403/1‑93 (MQ=255)

               gAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCGCACCAGCg                                                      <  2:410613/93‑1 (MQ=255)

               gAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCGCACCAGCg                                                      >  1:4902761/1‑93 (MQ=255)

               gAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCGCACCAGCg                                                      >  2:3864395/1‑93 (MQ=255)

                 aGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTc                                                    <  2:4593789/93‑1 (MQ=255)

                 aGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTc                                                    <  1:5924210/93‑1 (MQ=255)

                 aGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTc                                                    >  1:1838594/1‑93 (MQ=255)

                  ggATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCt                                                   >  2:4017896/1‑93 (MQ=255)

                               tttggtgACGGGATGAACGACGCCGAAATGCTGTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCAt                                      <  1:1056897/89‑1 (MQ=255)

                               tttggtgACGGGATGAACGACGCCGAAATGCTGTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCAt                                      <  2:5152570/89‑1 (MQ=255)

                               tttggtgACGGGATGAACGACGCCGAAATGCTGTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCAt                                      >  2:6403989/5‑93 (MQ=255)

                                ttggtgACGGGATGAACGACGCCGAAATGCTGTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATc                                     <  2:335385/90‑1 (MQ=255)

                                ttggtgACGGGATGAACGACGCCGAAATGCTGTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATc                                     >  2:175197/4‑93 (MQ=255)

                                ttggtgACGGGATGAACGACGCCGAAATGCTGTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATc                                     <  2:929718/90‑1 (MQ=255)

                                 tggtgACGGGATGAACGACGCCGAAATGCTGTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATcc                                    <  2:2541985/91‑1 (MQ=255)

                                  ggtgACGGGATGAACGACGCCGAAATGCTGTCGATGGCG‑GGGAAAGG‑CTGCATGATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATccc                                   <  2:2522323/92‑1 (MQ=255)

                                    tgACGGGATGAACGACGCCGAAATGCTGTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGa                                 >  2:5674995/1‑93 (MQ=255)

                                    tgACGGGATGAACGACGCCGAAATGCTGTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGa                                 >  1:2185073/1‑93 (MQ=255)

                                    tgACGGGATGAACGACGCCGAAATGCTGTCGATGGCG‑GGGAAAGG‑CTGCATCATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGa                                 >  2:6247931/1‑93 (MQ=255)

                                     gANGGGATGAACGACGCCGAANNGCTGTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAg                                <  2:2435353/93‑1 (MQ=255)

                                      aCGGGATGAACGACGCCGAAATGCTGTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGc                               <  2:1770585/93‑1 (MQ=255)

                                      aCGGGATGAACGACGCCGAAATGCTGTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGc                               >  2:931356/1‑93 (MQ=255)

                                       cGGGATGAACGACGCCGAAATGCTGTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCt                              <  1:737491/93‑1 (MQ=255)

                                         ggATGAACGACGCCGAAATGCTGTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCT‑gg                           >  1:3896910/1‑92 (MQ=255)

                                         ggATGAACGACGCCGAAATGCTGTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCT‑gg                           >  1:5501570/1‑92 (MQ=255)

                                          gATGAACGACGCCGAAATGCTGTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCT‑gga                          >  2:1986235/1‑91 (MQ=255)

                                          gATGAACGACGCCGAAATGCTGTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCT‑gga                          >  2:5522726/1‑91 (MQ=255)

                                          gATGAACGACGCCGAAATGCTGTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCT‑gga                          >  2:3644748/1‑91 (MQ=255)

                                           aTGAACGACGCCGAAATGCTGTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCT‑ggaa                         <  1:302638/93‑4 (MQ=255)

                                           aTGAACGACGCCGAAATGCTGTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCT‑ggaa                         <  2:5301800/93‑4 (MQ=255)

                                            tGAACGACGCCGAAATGCTGTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCT‑ggaag                        >  2:2414163/1‑89 (MQ=255)

                                            tGAACGACGCCGAAATGCTGTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCT‑ggaag                        <  1:1068601/93‑5 (MQ=255)

                                              aaCGACGCCGAAATGCTGTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTg                       >  2:2852590/1‑93 (MQ=255)

                                               acgacgCCGAAATGCTGTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGa                      >  1:1174402/1‑93 (MQ=255)

                                               acgacgCCGAAATGCTGTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGa                      <  2:5431799/93‑1 (MQ=255)

                                                cgacgCCGAAATGCTGTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGAt                     <  1:3177658/93‑1 (MQ=255)

                                                cgacgCCGAAATGCTGTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGAt                     >  2:3203999/1‑93 (MQ=255)

                                                 gacgCCGAAATGCTGTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCGCCCCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGAtt                    <  2:6547914/93‑1 (MQ=255)

                                                 gacgCCGAAATGCTGTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGAtt                    <  1:871684/93‑1 (MQ=255)

                                                 gacgCCGAAATGCTGTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGAtt                    <  1:2244435/93‑1 (MQ=255)

                                                 gacgCCGAAATGCTGTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGAtt                    <  2:3756621/93‑1 (MQ=255)

                                                  acgCCGAAATGCTGTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTg                   <  2:4805245/93‑1 (MQ=255)

                                                  acgCCGAAATGCTGTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTg                   <  2:6052003/93‑1 (MQ=255)

                                                   cgCCGAAATGCTGTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTgg                  <  2:5347832/93‑1 (MQ=255)

                                                    gCCGAAATGCTGTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGt                 >  1:1206236/1‑93 (MQ=255)

                                                     ccGAAATGCTGTCGATGGCG‑GGGAAA‑GTCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTa                >  2:303501/1‑93 (MQ=255)

                                                      cGAAATGCTGTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTAc               <  1:460624/93‑1 (MQ=255)

                                                       gAAATGCTGTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACt              >  1:1199901/1‑93 (MQ=255)

                                                       gAAATGCTGTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACt              >  1:1742587/1‑93 (MQ=255)

                                                       gAAATGCTGTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACt              >  2:608246/1‑93 (MQ=255)

                                                        aaaTGCTGTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTa             >  1:5147812/1‑93 (MQ=255)

                                                         aaTGCTGTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTaa            >  1:535230/1‑93 (MQ=255)

                                                          aTGCTGTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTagt           >  2:1366886/1‑91 (MQ=255)

                                                          aTGCTGTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAAt           <  2:3317426/93‑1 (MQ=255)

                                                          aTGCTGTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAAt           <  2:2542082/93‑1 (MQ=255)

                                                           tGCTGTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATg          <  1:6184572/93‑1 (MQ=255)

                                                            gCTGTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGc         <  1:2924265/93‑1 (MQ=255)

                                                             cTGTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGcc        <  1:3534800/93‑1 (MQ=255)

                                                             cTGTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGcc        >  1:1743840/1‑93 (MQ=255)

                                                             cTGTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGcc        >  2:6133468/1‑93 (MQ=255)

                                                              tGTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcg       >  2:371689/1‑93 (MQ=255)

                                                              tGTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcg       >  1:6401371/1‑93 (MQ=255)

                                                              tGTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcg       >  1:4101900/1‑93 (MQ=255)

                                                              cgTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcg       <  1:506992/92‑1 (MQ=255)

                                                               gTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcga      <  1:4260896/93‑1 (MQ=255)

                                                               gTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcga      <  2:4865780/93‑1 (MQ=255)

                                                               gTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcga      <  1:5663612/93‑1 (MQ=255)

                                                               gTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcga      >  2:1185427/1‑93 (MQ=255)

                                                               gTCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcga      <  1:3422428/93‑1 (MQ=255)

                                                                tCGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcgac     <  2:3028737/93‑1 (MQ=255)

                                                                 cGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcgacg    <  1:3548998/93‑1 (MQ=255)

                                                                 cGATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcgacg    >  1:523193/1‑93 (MQ=255)

                                                                  gATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcgacga   <  2:5635097/93‑1 (MQ=255)

                                                                  gATGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcgacga   >  2:1573723/1‑93 (MQ=255)

                                                                   aTGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCgacgac  <  2:1135102/93‑1 (MQ=255)

                                                                   aTGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCgacgac  >  1:5701826/1‑93 (MQ=255)

                                                                   aTGGCG‑GGGAAAGG‑CTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCgacgac  <  2:841724/93‑1 (MQ=255)

AACTGGCTACAG‑CTGAAGGATTG‑ATTGCGTT‑‑GTGACGGGATGAACGACGCCGAAATGCTGTCGATGGC‑‑GGGAAA‑G‑CTGCATTATGGGCAGTGGCCACCAGCGTCTGAAAGACCTTCATCCCGAGCT‑GAAGTGATTGGTACTAATGCCGACGAC  >  NC_000913/4008399‑4008551

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