BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,112,176	T→C	H65R (CAC→CGC)	b2972 ←	putative peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,112,176	0	T	C	100.0%	990.8 / NA	316	H65R (CAC→CGC)	b2972	putative peptidase
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (149/167); total (149/167)

Only 100 of 330 total aligned reads displayed.

TGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGTGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTCCGCCATTTGTTGGCGCAGCATGATTGGG  >  NC_000913/3112086‑3112265

tGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAgcgc                                                                                         <  2:5060428/93‑1 (MQ=255)

tGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAgcgc                                                                                         >  2:780112/1‑93 (MQ=255)

tGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAgcgc                                                                                         >  2:578490/1‑93 (MQ=255)

 gACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGcgcg                                                                                        >  1:3897028/1‑93 (MQ=255)

    aGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAAc                                                                                     >  2:3890729/1‑93 (MQ=255)

    aGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAAc                                                                                     <  1:2987005/93‑1 (MQ=255)

    aGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAAc                                                                                     >  1:1640453/1‑93 (MQ=255)

     gTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAAcg                                                                                    <  2:5181524/93‑1 (MQ=255)

     gTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTCACAATGTGGGCAGCGCGAAcg                                                                                    <  2:3327985/93‑1 (MQ=255)

      acgcgGCATCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAAcgc                                                                                   >  2:6892895/2‑93 (MQ=255)

       cgcgGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAAcgcg                                                                                  <  2:6099582/93‑1 (MQ=255)

        gcgGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGGCTGCTGACAATGTGGGCAGCGCGAAAGCgg                                                                                 >  2:1467035/1‑93 (MQ=255)

         cgGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGt                                                                                <  1:6335826/93‑1 (MQ=255)

         cgGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGt                                                                                >  1:6825015/1‑93 (MQ=255)

         cgGCAGCCCCCTTTGAGCATCCGCCAGGAGAACAGGGGAATATTGTCCACTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGt                                                                                <  1:7293705/93‑1 (MQ=255)

           gCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAg                                                                              >  2:858814/1‑93 (MQ=255)

            cAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAgc                                                                             <  2:5525242/93‑1 (MQ=255)

              gcgcCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAgcgc                                                                           >  1:2568442/1‑93 (MQ=255)

              gcgcCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAgcgc                                                                           <  1:7208448/93‑1 (MQ=255)

               cgcCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCa                                                                          >  2:3589602/1‑93 (MQ=255)

                gcCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCaa                                                                         <  1:7363363/93‑1 (MQ=255)

                gcCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCaa                                                                         >  1:5251623/1‑93 (MQ=255)

                  ccTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAgg                                                                       <  2:3193555/93‑1 (MQ=255)

                   cTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGc                                                                      >  2:5616046/1‑93 (MQ=255)

                   cTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGc                                                                      >  1:8753/1‑93 (MQ=255)

                    tttGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCt                                                                     >  1:2065126/1‑93 (MQ=255)

                    tttGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCt                                                                     <  2:2475528/93‑1 (MQ=255)

                      tGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTa                                                                   <  1:4642054/93‑1 (MQ=255)

                        aGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTAtt                                                                 <  2:4663829/93‑1 (MQ=255)

                        aGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTAtt                                                                 <  2:6971351/93‑1 (MQ=255)

                        aGCATCAGCCAGGAGAACAGCGGAATATTGTAACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTAtt                                                                 <  2:2551793/93‑1 (MQ=255)

                         gCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTAttt                                                                >  1:5798529/1‑93 (MQ=255)

                         gCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTAttt                                                                >  1:3695673/1‑93 (MQ=255)

                         gCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTAttt                                                                <  1:5446457/93‑1 (MQ=255)

                          cATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTAtttt                                                               <  2:4448507/93‑1 (MQ=255)

                           aTCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTAttttt                                                              >  1:4220479/1‑93 (MQ=255)

                             cAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGa                                                            <  1:4078821/93‑1 (MQ=255)

                             cAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGa                                                            <  1:1750711/93‑1 (MQ=255)

                              aGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGAc                                                           <  2:2782406/93‑1 (MQ=255)

                              aGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGAc                                                           <  2:2928044/93‑1 (MQ=255)

                              aGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGAc                                                           <  1:2677725/93‑1 (MQ=255)

                              aGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGAc                                                           <  1:6974634/93‑1 (MQ=255)

                              aGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGAc                                                           <  1:100015/93‑1 (MQ=255)

                               gCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACt                                                          >  2:6569750/1‑93 (MQ=255)

                               gCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACt                                                          <  1:1065595/93‑1 (MQ=255)

                                 cAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTgc                                                        >  2:6833532/1‑93 (MQ=255)

                                 cAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTgc                                                        <  2:1915804/93‑1 (MQ=255)

                                 cAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTgc                                                        <  2:1841857/93‑1 (MQ=255)

                                  aGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTgcg                                                       >  1:2379999/1‑93 (MQ=255)

                                  aGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTgcg                                                       >  2:3458067/1‑93 (MQ=255)

                                   ggAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTgcgc                                                      >  1:5131102/1‑93 (MQ=255)

                                        aCAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCg                                                 <  2:7067444/93‑1 (MQ=255)

                                         cAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGa                                                >  1:3855389/1‑93 (MQ=255)

                                           gCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGAct                                              >  1:11797/1‑92 (MQ=255)

                                           gCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACa                                              >  2:2279990/1‑93 (MQ=255)

                                           gCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACa                                              >  2:3741542/1‑93 (MQ=255)

                                             ggAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACAtt                                            <  1:4924683/93‑1 (MQ=255)

                                              gAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACAttt                                           >  1:981768/1‑93 (MQ=255)

                                              gAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACAttt                                           <  2:6259254/93‑1 (MQ=255)

                                              gAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTAGTTTTGACTGCGCACTCGACAttt                                           >  2:4306477/1‑93 (MQ=255)

                                               aaTATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTc                                          <  2:6657551/93‑1 (MQ=255)

                                                 tatTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCAc                                        >  2:4532667/1‑93 (MQ=255)

                                                 tatTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCAc                                        >  1:1942019/1‑93 (MQ=255)

                                                 tatTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCAc                                        >  2:3064332/1‑93 (MQ=255)

                                                 tatTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCAc                                        >  2:7359274/1‑93 (MQ=255)

                                                  atTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCAcc                                       <  1:4779992/93‑1 (MQ=255)

                                                   ttGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCg                                      >  2:7272811/1‑93 (MQ=255)

                                                    tGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCATTCGACATTTCACCGt                                     <  1:7709154/93‑1 (MQ=255)

                                                         cGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAAc                                <  2:5346878/93‑1 (MQ=255)

                                                          gTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACt                               >  2:4591382/1‑93 (MQ=255)

                                                           tATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTc                              >  2:4226898/1‑93 (MQ=255)

                                                            aTCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTcc                             <  1:286116/93‑1 (MQ=255)

                                                               cGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTCCGcc                          >  1:617979/1‑93 (MQ=255)

                                                                ggATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTccacca                         >  2:6433356/1‑93 (MQ=255)

                                                                 gATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTCCGCCAt                        <  2:5393741/93‑1 (MQ=255)

                                                                    ggTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTCCGCCATttg                     <  2:1069976/93‑1 (MQ=255)

                                                                     gTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTCCGCCATttgt                    <  1:5490906/93‑1 (MQ=255)

                                                                      tCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTCCGCCATttgtt                   <  2:3657179/93‑1 (MQ=255)

                                                                      tCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTCCGCCATttgtt                   >  1:2175213/1‑93 (MQ=255)

                                                                      tCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTCCGCCATttgtt                   >  1:4270587/1‑93 (MQ=255)

                                                                      tCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTCCGCCATttgtt                   <  1:6736916/93‑1 (MQ=255)

                                                                      tCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTCCGCCATttgtt                   >  1:3548150/1‑93 (MQ=255)

                                                                      tCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTCCGCCATttgtt                   <  1:7484792/93‑1 (MQ=255)

                                                                         gctgACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTCCGCCATTTGTTGgc                >  1:808390/1‑93 (MQ=255)

                                                                          ctgACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTCCGCCATTTGTTGgcg               >  1:3303153/1‑93 (MQ=255)

                                                                          ctgACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTCCGCCATTTGTTGgcg               >  2:1636745/1‑93 (MQ=255)

                                                                            gACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTCCGCCATTTGTTGGCgca             <  2:4731053/93‑1 (MQ=255)

                                                                               agtgtgGGCAGCGCGAGCGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTCCGCCATTTGTTGGCgcagca          <  1:6705232/91‑1 (MQ=255)

                                                                               aaTGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTCCGCCATTTGTTGGCgcagca          <  2:6132183/93‑1 (MQ=255)

                                                                                aTGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTCCGCCATTTGTTGGCGCAGCAt         >  2:5882971/1‑93 (MQ=255)

                                                                                aTGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTCCGCCATTTGTTGGCGCAGCAt         <  1:5858492/93‑1 (MQ=255)

                                                                                aTGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTCCGCCATTTGTTGGCGCAGCAt         <  2:925811/93‑1 (MQ=255)

                                                                                 tgtgGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTCCGCCATTTGTTGGCGCAGCATg        >  2:4119872/1‑93 (MQ=255)

                                                                                  gtgGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTCCGCCATTTGTTGGCGCAGCATGa       <  1:5699159/93‑1 (MQ=255)

                                                                                  gtgGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTCCGCCATTTGTTGGCGCAGCATGa       <  2:7188959/93‑1 (MQ=255)

                                                                                    gggCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTCCGCCATTTGTTGGCGCAGCATGAtt     <  1:1307262/93‑1 (MQ=255)

                                                                                     ggCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTCCGCCATTTGTTGGCGCAGCATGATTg    <  1:5354041/93‑1 (MQ=255)

                                                                                      gCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTCCGCCATTTGTTGGCGCAGCATGATTgg   <  2:5570840/93‑1 (MQ=255)

                                                                                      gCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTCCGCCATTTGTTGGCGCAGCATGATTgg   >  1:5997484/1‑93 (MQ=255)

                                                                                       cAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTCCGCCATTTGTTGGCGCAGCATGATTggg  <  1:5629959/93‑1 (MQ=255)

TGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGTGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTCCGCCATTTGTTGGCGCAGCATGATTGGG  >  NC_000913/3112086‑3112265

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