BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,853,551	T→C	Q16R (CAG→CGG)	yidG ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,853,551	0	T	C	100.0%	788.5 / NA	227	Q16R (CAG→CGG)	yidG	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (113/114); total (113/114)

Only 100 of 238 total aligned reads displayed.

GCCTGGTGCCAGTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATAT  >  NC_000913/3853461‑3853640

gCCTGGTGCCAGTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGc                                                                                         <  2:3161854/93‑1 (MQ=255)

 ccTGGTGCCAGTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCa                                                                                        >  1:3606628/1‑93 (MQ=255)

  cTGGTGCCAGTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAg                                                                                       <  2:2388950/93‑1 (MQ=255)

  cTGGTGCCAGTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAg                                                                                       <  2:3628312/93‑1 (MQ=255)

  cTGGTGCCAGTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAg                                                                                       <  1:3437467/93‑1 (MQ=255)

   tGGTGCCAGTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCCGCCGCAGt                                                                                      >  1:5000126/1‑93 (MQ=255)

     gTGCCAGTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTcc                                                                                    <  1:1769419/93‑1 (MQ=255)

      tGCCAGTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCt                                                                                   >  2:5342919/1‑93 (MQ=255)

      tGCCAGTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCt                                                                                   <  2:31791/93‑1 (MQ=255)

       gCCAGTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTg                                                                                  >  2:4912404/1‑93 (MQ=255)

       gCCAGTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTg                                                                                  <  1:630584/93‑1 (MQ=255)

        ccAGTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTgg                                                                                 >  2:4476177/1‑93 (MQ=255)

        ccAGTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTgg                                                                                 >  2:1240735/1‑93 (MQ=255)

         cAGTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTggg                                                                                >  1:3368016/1‑93 (MQ=255)

         cAGTTGTGTTTGATAGCCAACGCCATCAACGCGCCATAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTggg                                                                                <  2:3982397/93‑1 (MQ=255)

          aattGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGt                                                                               >  1:6373842/3‑93 (MQ=255)

          aGTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGt                                                                               >  2:2729917/1‑93 (MQ=255)

           gTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTc                                                                              >  1:561328/1‑93 (MQ=255)

           gTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTc                                                                              >  2:2479323/1‑93 (MQ=255)

            ttGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCg                                                                             <  2:2830361/93‑1 (MQ=255)

            ttGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCg                                                                             >  1:5150517/1‑93 (MQ=255)

            ttGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCg                                                                             <  1:744008/93‑1 (MQ=255)

             tgtgTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCgg                                                                            <  2:5132064/93‑1 (MQ=255)

                 tttGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAAt                                                                        <  2:1507782/93‑1 (MQ=255)

                 tttGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAAt                                                                        >  2:3304574/1‑93 (MQ=255)

                  ttGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATg                                                                       >  2:4185374/1‑93 (MQ=255)

                     aTAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGt                                                                    <  1:611471/93‑1 (MQ=255)

                     aTAGCCAACGCCAGCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGt                                                                    <  1:1916476/93‑1 (MQ=255)

                      tAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTc                                                                   <  2:136068/93‑1 (MQ=255)

                      tAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTc                                                                   <  1:5431404/93‑1 (MQ=255)

                      tAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTc                                                                   >  1:6175220/1‑93 (MQ=255)

                       aGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCt                                                                  <  1:6462732/93‑1 (MQ=255)

                       aGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCt                                                                  <  1:2424428/93‑1 (MQ=255)

                        gCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTg                                                                 >  1:2359389/1‑93 (MQ=255)

                        gCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTg                                                                 <  2:1682359/93‑1 (MQ=255)

                           caCGCCATCAACGCGCCGTAGCCCAGCATGGTAGGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCt                                                              <  2:545426/92‑1 (MQ=255)

                           aaCGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCt                                                              >  1:4060569/1‑93 (MQ=255)

                            aCGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCtt                                                             <  2:2143321/93‑1 (MQ=255)

                            aCGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCtt                                                             >  2:732413/1‑93 (MQ=255)

                                cATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCg                                                         >  2:4392934/1‑93 (MQ=255)

                                  tCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGc                                                       >  2:796602/1‑93 (MQ=255)

                                   cAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCt                                                      <  1:2230468/93‑1 (MQ=255)

                                    aaCGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTa                                                     <  2:3615928/93‑1 (MQ=255)

                                    aaCGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTa                                                     >  2:6084827/1‑93 (MQ=255)

                                      cgcgCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATc                                                   <  2:2626480/93‑1 (MQ=255)

                                      cgcgCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATc                                                   <  1:4267635/93‑1 (MQ=255)

                                      cgcgCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATc                                                   <  1:5750158/93‑1 (MQ=255)

                                       gcgtCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATcc                                                  <  2:5319814/93‑1 (MQ=255)

                                        cgcCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCg                                                 >  1:2777844/1‑93 (MQ=255)

                                         gcCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCgg                                                <  2:6109749/93‑1 (MQ=255)

                                         gcCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCgg                                                <  1:1648123/93‑1 (MQ=255)

                                          ccGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGc                                               >  2:5964480/1‑93 (MQ=255)

                                          ccGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGc                                               <  2:3235680/93‑1 (MQ=255)

                                           cGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCa                                              <  2:2205305/93‑1 (MQ=255)

                                            gTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCAt                                             >  2:5500094/1‑93 (MQ=255)

                                             tAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATa                                            >  2:5653813/1‑93 (MQ=255)

                                             tAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATa                                            >  2:605217/1‑93 (MQ=255)

                                             tAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATa                                            >  2:423947/1‑93 (MQ=255)

                                              aGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATac                                           >  1:1063086/1‑93 (MQ=255)

                                              aGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATac                                           <  2:6132508/93‑1 (MQ=255)

                                               gCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATaca                                          >  2:134387/1‑93 (MQ=255)

                                                cccAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATacaa                                         <  2:5417140/93‑1 (MQ=255)

                                                  cAGCATGGTACGAAACCACCCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATacaaca                                       <  2:466982/93‑1 (MQ=255)

                                                     cATGGTCCGAGACCACGCCAGTGATTTCCGCTCCGTCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCa                                    <  1:5322013/93‑1 (MQ=255)

                                                      aTGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAg                                   >  2:4577517/1‑93 (MQ=255)

                                                      aTGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAg                                   >  1:5971779/1‑93 (MQ=255)

                                                      aTGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAg                                   <  1:580466/93‑1 (MQ=255)

                                                        ggTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTc                                 <  1:4118288/93‑1 (MQ=255)

                                                         gTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTcc                                >  1:3835253/1‑93 (MQ=255)

                                                          tACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTccc                               <  1:1188357/93‑1 (MQ=255)

                                                           aCGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCa                              <  1:5838132/93‑1 (MQ=255)

                                                           aCGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCa                              <  1:2766372/93‑1 (MQ=255)

                                                           aCGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCa                              <  1:5036171/93‑1 (MQ=255)

                                                            cGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCAt                             <  2:3943052/93‑1 (MQ=255)

                                                              aaaCCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATaa                           >  2:2542761/1‑93 (MQ=255)

                                                              aaaCCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATaa                           >  1:1254364/1‑93 (MQ=255)

                                                              aaaCCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATaa                           >  1:3506563/1‑93 (MQ=255)

                                                                aCCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACg                         <  1:5183062/93‑1 (MQ=255)

                                                                aCCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACg                         >  2:4938197/1‑93 (MQ=255)

                                                                aCCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACg                         >  2:3000983/1‑93 (MQ=255)

                                                                  cACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGAt                       >  2:1048156/1‑93 (MQ=255)

                                                                  cACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCAGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGAt                       <  2:5023222/93‑1 (MQ=255)

                                                                   aCGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATc                      <  1:982094/93‑1 (MQ=255)

                                                                    cGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCa                     >  1:429798/1‑93 (MQ=255)

                                                                     gCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCAc                    <  1:5138309/93‑1 (MQ=255)

                                                                      ccAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCAcc                   <  2:3559498/93‑1 (MQ=255)

                                                                       cAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCg                  >  1:4546567/1‑93 (MQ=255)

                                                                         gTGATGTGCGCTGCGGCCGCAGTCCTGGTTCGGCAATGCGTCTGGCTTTGCGGCTATCCGTCATACAACACCAGTCCCATAACGATCACCGcc                <  1:2927374/93‑1 (MQ=255)

                                                                           gATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCAc              <  1:4751226/93‑1 (MQ=255)

                                                                           gATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCAc              <  1:3657655/93‑1 (MQ=255)

                                                                              gtgCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGAc           <  2:2882202/93‑1 (MQ=255)

                                                                               tgCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGAcc          <  1:1514325/93‑1 (MQ=255)

                                                                                   cTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCa      >  2:1969664/1‑93 (MQ=255)

                                                                                    tCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCaa     >  1:451150/1‑93 (MQ=255)

                                                                                    tCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCaa     <  2:3452596/93‑1 (MQ=255)

                                                                                    tCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCaa     >  2:4752986/1‑93 (MQ=255)

                                                                                    tCCGGCCCCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCaa     >  1:5896959/1‑93 (MQ=255)

                                                                                     ccGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAat    <  2:1671255/93‑1 (MQ=255)

                                                                                       ggCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAatat  <  1:5486541/93‑1 (MQ=255)

                                                                                       ggCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAatat  >  1:2110039/1‑93 (MQ=255)

GCCTGGTGCCAGTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATAT  >  NC_000913/3853461‑3853640

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