BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,576,170	A→G	R246R (CGT→CGC)	mcrB ←	McrB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,576,170	0	A	G	100.0%	659.1 / NA	180	R246R (CGT→CGC)	mcrB	McrB
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (115/65); total (115/65)

Only 100 of 190 total aligned reads displayed.

AAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAA  >  NC_000913/4576080‑4576259

aaaaaTATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCg                                                                                         >  1:5366698/1‑93 (MQ=255)

aaaaaTATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCg                                                                                         >  2:3531202/1‑93 (MQ=255)

 aaaaTATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGa                                                                                        >  1:3908492/1‑93 (MQ=255)

  aaaTATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGAt                                                                                       >  2:5448692/1‑93 (MQ=255)

   aaTATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCTGACGCCATTCGGGCGATa                                                                                      <  2:251729/93‑1 (MQ=255)

   aaTATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATa                                                                                      >  1:3727585/1‑93 (MQ=255)

   aaTATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATa                                                                                      <  1:3600115/93‑1 (MQ=255)

   aaTATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATa                                                                                      <  2:3849138/93‑1 (MQ=255)

    atatACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAg                                                                                     >  1:108050/1‑93 (MQ=255)

    atatACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAg                                                                                     >  1:6088151/1‑93 (MQ=255)

    atatACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAg                                                                                     >  1:5970603/1‑93 (MQ=255)

    atatACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAAATGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAg                                                                                     >  2:4673826/1‑93 (MQ=255)

     tataCTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGc                                                                                    >  2:1436251/1‑93 (MQ=255)

       taCTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGccc                                                                                  <  1:4859094/93‑1 (MQ=255)

         cTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTg                                                                                >  2:4471618/1‑93 (MQ=255)

          tttttCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGt                                                                               >  1:2460431/1‑93 (MQ=255)

           ttttCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGta                                                                              <  1:4266160/93‑1 (MQ=255)

            tttCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGtat                                                                             <  2:2817844/93‑1 (MQ=255)

              tctcTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATaa                                                                           <  1:830539/93‑1 (MQ=255)

              tctcTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATaa                                                                           >  1:6609805/1‑93 (MQ=255)

              tctcTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATaa                                                                           >  1:2286244/1‑93 (MQ=255)

              tctcTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATaa                                                                           >  2:6131946/1‑93 (MQ=255)

                tctGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATaaaa                                                                         >  2:2096967/1‑93 (MQ=255)

                tctGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATaaaa                                                                         >  1:6494560/1‑93 (MQ=255)

                tctGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATaaaa                                                                         >  2:4570397/1‑93 (MQ=255)

                tctGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATaaaa                                                                         >  1:4073085/1‑93 (MQ=255)

                  tGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATc                                                                       >  2:5311016/1‑93 (MQ=255)

                   ggCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATcc                                                                      <  1:3062267/93‑1 (MQ=255)

                    gctgctCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCt                                                                     <  2:4731133/93‑1 (MQ=255)

                    gctgctCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCt                                                                     >  1:5076254/1‑93 (MQ=255)

                    gctgctCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCt                                                                     >  2:6513860/1‑93 (MQ=255)

                    gctgctCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCt                                                                     <  2:818898/93‑1 (MQ=255)

                     ctgctCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTc                                                                    >  1:1055151/1‑93 (MQ=255)

                      tgctCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCa                                                                   <  2:3956252/93‑1 (MQ=255)

                       gctCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGGGAGTGCGCTGTATAAAATCCTcct                                                                  >  1:2837064/1‑91 (MQ=255)

                        ctctTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATa                                                                 >  1:264897/1‑93 (MQ=255)

                          ctTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGc                                                               >  2:2727395/1‑93 (MQ=255)

                          ctTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGc                                                               <  2:3679264/93‑1 (MQ=255)

                          ctTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGc                                                               >  1:5693828/1‑93 (MQ=255)

                             tAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCtat                                                            >  1:2166284/1‑93 (MQ=255)

                             tAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCtat                                                            <  2:4823825/93‑1 (MQ=255)

                              aGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCtata                                                           >  2:1698663/1‑93 (MQ=255)

                               gCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATaa                                                          <  1:345046/93‑1 (MQ=255)

                               gCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATaa                                                          >  2:4129964/1‑93 (MQ=255)

                                  tGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGAt                                                       >  2:4941183/1‑93 (MQ=255)

                                  tGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGAt                                                       <  1:2976434/93‑1 (MQ=255)

                                   gCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGAtt                                                      <  1:6113327/93‑1 (MQ=255)

                                   gCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGAtt                                                      <  1:5299429/93‑1 (MQ=255)

                                    cTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTg                                                     >  1:4087145/1‑93 (MQ=255)

                                     tGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGa                                                    >  1:3289319/1‑93 (MQ=255)

                                       aCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATg                                                  >  1:4960475/1‑93 (MQ=255)

                                       aCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATg                                                  <  2:6215106/93‑1 (MQ=255)

                                          aaaaTTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGaa                                               <  1:6559712/93‑1 (MQ=255)

                                           aaaTTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAAc                                              >  2:1246409/1‑93 (MQ=255)

                                             aTTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTg                                            <  2:422756/93‑1 (MQ=255)

                                             aTTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTg                                            >  2:1388872/1‑93 (MQ=255)

                                              ttGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGa                                           <  2:5218270/93‑1 (MQ=255)

                                              ttGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGa                                           >  1:1989524/1‑93 (MQ=255)

                                               tGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGaa                                          <  2:1446488/93‑1 (MQ=255)

                                                gTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAAc                                         <  1:2235131/93‑1 (MQ=255)

                                                gTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAAc                                         >  2:6326981/1‑93 (MQ=255)

                                                   aaaTATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCat                                      >  1:1305549/1‑93 (MQ=255)

                                                   aaaTATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCat                                      >  1:4514966/1‑93 (MQ=255)

                                                    aaTAGGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCata                                     <  1:4700122/93‑1 (MQ=255)

                                                      tatGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATAtt                                   >  1:2187619/1‑93 (MQ=255)

                                                           cGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGAcgc                              <  2:1925765/93‑1 (MQ=255)

                                                           cGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGAcgc                              >  1:3781715/1‑93 (MQ=255)

                                                           cGTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGAcgc                              >  1:6445337/1‑93 (MQ=255)

                                                            gTCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGAcgcg                             >  1:643838/1‑93 (MQ=255)

                                                             tCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGt                            <  1:279109/93‑1 (MQ=255)

                                                             tCTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGt                            <  2:1880822/93‑1 (MQ=255)

                                                              cTTTACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGtt                           <  2:3058279/93‑1 (MQ=255)

                                                               tttACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTg                          >  2:31676/1‑93 (MQ=255)

                                                                ttACGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGc                         <  1:5122173/93‑1 (MQ=255)

                                                                   cGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGa                      <  1:1475996/93‑1 (MQ=255)

                                                                   cGTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGa                      <  2:3488906/93‑1 (MQ=255)

                                                                    gTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAg                     >  1:4174033/1‑93 (MQ=255)

                                                                    gTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAg                     >  1:366919/1‑93 (MQ=255)

                                                                    gTCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAg                     >  1:5854042/1‑93 (MQ=255)

                                                                     tCGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGc                    <  1:626789/93‑1 (MQ=255)

                                                                      cGGAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGcc                   <  1:3083611/93‑1 (MQ=255)

                                                                       ggAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCt                  <  2:2619070/93‑1 (MQ=255)

                                                                        gAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCtt                 <  2:434726/93‑1 (MQ=255)

                                                                        gAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCtt                 >  1:1247449/1‑93 (MQ=255)

                                                                        gAAGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCtt                 >  1:1824789/1‑93 (MQ=255)

                                                                          aGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCtttt               >  2:2313187/1‑93 (MQ=255)

                                                                          aGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCtttt               >  2:6297274/1‑93 (MQ=255)

                                                                          aGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCtttt               >  1:4814584/1‑93 (MQ=255)

                                                                          aGCCGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCtttt               <  2:5806091/93‑1 (MQ=255)

                                                                            ccGACGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTtct             >  1:6529879/1‑93 (MQ=255)

                                                                               aCGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCt          >  1:1369637/1‑93 (MQ=255)

                                                                               aCGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCt          <  2:5591997/93‑1 (MQ=255)

                                                                               aCGCCATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCt          >  1:6358944/1‑93 (MQ=255)

                                                                                  ccATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTc       >  1:1835674/1‑93 (MQ=255)

                                                                                   cATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCa      <  1:972229/93‑1 (MQ=255)

                                                                                   cATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCa      >  2:707892/1‑93 (MQ=255)

                                                                                   cATTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATGGCCGCGTTGCGGAGCCTTTTCTCCTGTCa      <  2:2873113/93‑1 (MQ=255)

                                                                                    aTTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAg     <  1:2799326/93‑1 (MQ=255)

                                                                                    aTTCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAg     >  1:1930386/1‑93 (MQ=255)

                                                                                      tCGGGCGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCa   <  1:5226735/93‑1 (MQ=255)

AAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAA  >  NC_000913/4576080‑4576259

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