BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,561,348	A→G	C49C (TGT→TGC)	b4563 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,561,348	0	A	G	100.0%	601.4 / NA	181	C49C (TGT→TGC)	b4563	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (80/101); total (80/101)

Only 100 of 195 total aligned reads displayed.

ATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCAACGC  >  NC_000913/4561258‑4561437

 tgtgCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAg                                                                                        <  1:3902682/93‑1 (MQ=255)

 tgtgCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAg                                                                                        <  1:3682176/93‑1 (MQ=255)

  gtgCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGTCTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAgg                                                                                       >  2:4527047/1‑93 (MQ=255)

  gtgCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAgg                                                                                       <  2:1543515/93‑1 (MQ=255)

  gtgCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAgg                                                                                       >  1:1558760/1‑93 (MQ=255)

  gtgCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAgg                                                                                       <  2:136826/93‑1 (MQ=255)

   tgCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGt                                                                                      <  1:408644/93‑1 (MQ=255)

    gCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTc                                                                                     <  1:3599631/93‑1 (MQ=255)

     cAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCa                                                                                    >  2:1211133/1‑93 (MQ=255)

     cAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCa                                                                                    >  1:4611982/1‑93 (MQ=255)

      aaTGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAg                                                                                   >  2:844232/1‑93 (MQ=255)

      aaTGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAg                                                                                   <  2:3598425/93‑1 (MQ=255)

         gCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCa                                                                                <  1:3270097/93‑1 (MQ=255)

              gcCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGc                                                                           >  2:897512/1‑93 (MQ=255)

              gcCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGc                                                                           >  2:1459356/1‑93 (MQ=255)

              cccAGACACCATCCTGCTGGCAATAAATTGACACAGTCAGCCCAGGGACCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGc                                                                           <  1:2612181/92‑1 (MQ=255)

               ccAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCa                                                                          >  2:6370750/1‑93 (MQ=255)

                cAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCAt                                                                         >  1:6103637/1‑93 (MQ=255)

                 aGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATa                                                                        >  1:1735556/1‑93 (MQ=255)

                 aGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATa                                                                        >  1:5062351/1‑93 (MQ=255)

                  gACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATaa                                                                       >  2:5133870/1‑93 (MQ=255)

                    cacCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAAcc                                                                     <  1:6407242/93‑1 (MQ=255)

                    cacCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAAcc                                                                     >  2:4956023/1‑93 (MQ=255)

                     cccATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAAccc                                                                    <  1:5842556/92‑1 (MQ=255)

                       cATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGc                                                                  <  1:2341585/93‑1 (MQ=255)

                        aTCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCa                                                                 >  2:3760124/1‑93 (MQ=255)

                        aTCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCa                                                                 >  2:4749751/1‑93 (MQ=255)

                          ccTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAgc                                                               <  1:1866188/93‑1 (MQ=255)

                            tgctgGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAgccg                                                             <  2:2559523/93‑1 (MQ=255)

                            cgctgGCAATAAGTTGACTCAGTCGGCCAAGCGCTCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAgccg                                                             <  2:2559522/92‑1 (MQ=255)

                              ctgGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAgccgcc                                                           >  2:5676737/1‑93 (MQ=255)

                              ctgGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAgccgcc                                                           <  1:4763137/93‑1 (MQ=255)

                              ctgGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAgccgcc                                                           <  1:2982691/93‑1 (MQ=255)

                              ctgGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAgccgca                                                           <  1:2264970/93‑2 (MQ=255)

                               tgGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGccgccg                                                          >  1:1357604/1‑93 (MQ=255)

                               tgGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGccgccg                                                          <  1:5263436/93‑1 (MQ=255)

                               tgGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGccgccg                                                          >  1:5956350/1‑93 (MQ=255)

                                ggCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCcgccgc                                                         >  2:1909231/1‑93 (MQ=255)

                                ggCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCcgccgc                                                         >  1:4153808/1‑93 (MQ=255)

                                  cAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGt                                                       >  1:215982/1‑93 (MQ=255)

                                   aaTAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTa                                                      >  2:4042073/1‑93 (MQ=255)

                                    aTAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTaa                                                     >  2:4388484/1‑93 (MQ=255)

                                     tAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTaaa                                                    <  2:5275397/93‑1 (MQ=255)

                                      aaaTTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAAt                                                   <  2:5910063/93‑1 (MQ=255)

                                       aaTTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAACTc                                                  >  2:393896/1‑93 (MQ=255)

                                       aaTTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATc                                                  >  1:2361367/1‑93 (MQ=255)

                                        aTTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCCAATGCATAACCCGCAGCCGCCGCGTAAATCa                                                 >  1:4533635/1‑93 (MQ=255)

                                         ttGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCAc                                                <  2:692462/93‑1 (MQ=255)

                                         ttGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCAc                                                <  2:2842877/93‑1 (MQ=255)

                                              tcagtcagCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGa                                           <  1:3696409/93‑1 (MQ=255)

                                               cagtcagCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGaa                                          <  2:6179349/93‑1 (MQ=255)

                                               cagtcagCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGaa                                          <  1:2732495/93‑1 (MQ=255)

                                               cagtcagCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGaa                                          <  2:3224266/93‑1 (MQ=255)

                                                agtcagCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAat                                         >  2:3069583/1‑93 (MQ=255)

                                                 gtcagCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAata                                        >  2:2484101/1‑93 (MQ=255)

                                                 gtcagCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAata                                        <  2:815136/93‑1 (MQ=255)

                                                  tcagCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAatat                                       <  1:2516899/93‑1 (MQ=255)

                                                   cagCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATc                                      <  2:5352845/93‑1 (MQ=255)

                                                   cagCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATc                                      >  1:1201022/1‑93 (MQ=255)

                                                   cagCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATc                                      <  2:3414060/93‑1 (MQ=255)

                                                   cagCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATc                                      <  1:6162244/93‑1 (MQ=255)

                                                    agCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCg                                     <  2:4222393/93‑1 (MQ=255)

                                                        cAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGa                                 >  1:1184954/1‑93 (MQ=255)

                                                         aGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGAc                                >  1:308673/1‑93 (MQ=255)

                                                         aGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGAc                                <  2:303660/93‑1 (MQ=255)

                                                         aGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGAc                                <  1:978385/93‑1 (MQ=255)

                                                         aGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGAc                                >  1:392883/1‑93 (MQ=255)

                                                         aGGGATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGAc                                <  2:3322622/93‑1 (MQ=255)

                                                          gggATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGAcg                               >  2:3285488/1‑93 (MQ=255)

                                                          gggATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGAcg                               >  1:1605908/1‑93 (MQ=255)

                                                          gggATCACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGAcg                               >  1:1117389/1‑93 (MQ=255)

                                                               cACCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCAGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGaaa                          >  2:1714231/1‑93 (MQ=255)

                                                                aCCCATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGaaaa                         >  1:1691131/1‑93 (MQ=255)

                                                                 cccATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAg                        >  1:3410699/1‑93 (MQ=255)

                                                                 cccATTGTATTTCCCCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGGAAATCACGGGGAATATCGCTGACGCGAAAAg                        <  1:280310/93‑1 (MQ=255)

                                                                  ccATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGt                       <  2:4843360/93‑1 (MQ=255)

                                                                  ccATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGt                       <  2:3501924/93‑1 (MQ=255)

                                                                  ccATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGt                       <  2:5504113/93‑1 (MQ=255)

                                                                   cATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTc                      >  2:5152553/1‑93 (MQ=255)

                                                                   cATTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTc                      >  1:1043115/1‑93 (MQ=255)

                                                                    aTTGTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCa                     >  1:2001633/1‑93 (MQ=255)

                                                                       gTATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGc                  >  2:3579592/1‑93 (MQ=255)

                                                                        tATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGcc                 <  1:1132048/93‑1 (MQ=255)

                                                                        tATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGcc                 <  2:1636606/93‑1 (MQ=255)

                                                                        tATTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGcc                 <  2:6241216/93‑1 (MQ=255)

                                                                         aTTTCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCAGACGCGAAAAGTCATGccc                <  1:6314766/93‑1 (MQ=255)

                                                                           ttCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGc              <  2:382530/93‑1 (MQ=255)

                                                                           ttCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGc              <  1:5596036/93‑1 (MQ=255)

                                                                            tCCTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGct             >  2:2221986/1‑93 (MQ=255)

                                                                              cTCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGctgc           <  2:2535163/93‑1 (MQ=255)

                                                                               tCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGctgcc          <  2:3313425/93‑1 (MQ=255)

                                                                               tCGCCGGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGctgcc          <  1:4910523/93‑1 (MQ=255)

                                                                                    ggggggGCAGGTCAGCCAAATGCATACCCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCaa     <  2:2963134/93‑1 (MQ=255)

                                                                                    ggggggGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCaa     <  1:1643510/93‑1 (MQ=255)

                                                                                    ggggggGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCaa     >  1:496099/1‑93 (MQ=255)

                                                                                     ggggggCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCAAc    <  2:4497768/93‑1 (MQ=255)

                                                                                     ggggggCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCAAc    <  1:4431089/93‑1 (MQ=255)

                                                                                     ggggggCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCAAc    >  1:3519343/1‑93 (MQ=255)

                                                                                      ggggGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCAACg   <  2:5062306/93‑1 (MQ=255)

                                                                                       ggggCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCAACGc  <  1:3247592/93‑1 (MQ=255)

ATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCAACGC  >  NC_000913/4561258‑4561437

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