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	99.4%	510.8 / NA	163	H65R (CAC→CGC)	b2972	putative peptidase
Reads supporting (aligned to +/- strand): ref base T (0/0); major base C (95/67); minor base G (0/1); total (95/68)
Fisher's exact test for biased strand distribution p-value = 4.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 175 total aligned reads displayed.

TGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGTGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTCCGCCATTTGTTGGCGCAGCATGATTGG  >  NC_000913/3112086‑3112264

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

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

  aCAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGa                                                                                      >  2:594912/1‑93 (MQ=255)

   cAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGaa                                                                                     <  1:1928385/93‑1 (MQ=255)

     gTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAAcg                                                                                   >  2:2614184/1‑93 (MQ=255)

     gTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAAcg                                                                                   >  1:3082804/1‑93 (MQ=255)

      tCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAAcgc                                                                                  >  1:2747528/1‑93 (MQ=255)

      tCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAAcgc                                                                                  >  1:4313899/1‑93 (MQ=255)

        gcgGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCgg                                                                                <  1:3540215/93‑1 (MQ=255)

        ccgGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCgg                                                                                <  1:3042780/92‑1 (MQ=255)

           gCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAg                                                                             <  1:4146965/93‑1 (MQ=255)

           gCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAg                                                                             <  1:2182090/93‑1 (MQ=255)

           gCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAg                                                                             <  2:1311369/93‑1 (MQ=255)

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

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

             aGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAgcg                                                                           >  2:2865921/1‑93 (MQ=255)

              gcgcCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAATGCGGTAgcgc                                                                          >  2:747891/1‑93 (MQ=255)

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

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

                 ccttttGTGCATCAGCTAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGGGCGAACGCGGTAGCGCAAg                                                                       <  2:522234/90‑1 (MQ=255)

                 cccTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAg                                                                       <  2:4011613/93‑1 (MQ=255)

                 cccTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAg                                                                       <  1:399159/93‑1 (MQ=255)

                 cccTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAg                                                                       >  2:3285414/1‑93 (MQ=255)

                 cccTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAg                                                                       >  1:16691/1‑93 (MQ=255)

                  ccTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAgg                                                                      >  1:2130135/1‑93 (MQ=255)

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

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

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

                     ttGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTGGCGCAAGGCtt                                                                   <  1:3299123/93‑1 (MQ=255)

                       gAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTAt                                                                 <  1:829803/93‑1 (MQ=255)

                       gAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTAt                                                                 <  2:654863/93‑1 (MQ=255)

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

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

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

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

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

                            tCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTg                                                            >  2:2075879/1‑93 (MQ=255)

                            tCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTg                                                            >  2:3993997/1‑93 (MQ=255)

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

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

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

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

                                     agaACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTATCGCAAGGCTTATTTTTGACTGCGCAc                                                   >  2:3792832/1‑93 (MQ=255)

                                     agaACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCAc                                                   <  2:3576832/93‑1 (MQ=255)

                                     agaACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCAc                                                   <  2:3539761/93‑1 (MQ=255)

                                     agaACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCAc                                                   >  1:2193179/1‑93 (MQ=255)

                                      gaACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACt                                                  <  2:1436787/93‑1 (MQ=255)

                                       aaCAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTc                                                 <  2:2444213/93‑1 (MQ=255)

                                        aCAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCg                                                >  2:1665428/1‑93 (MQ=255)

                                          aGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGAGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGAc                                              >  2:4422199/1‑93 (MQ=255)

                                           gCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGGAGCGCAAGGCTTATTTTTGACTGCGCACTCcaca                                             >  2:3343333/1‑93 (MQ=255)

                                           gCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGAGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACa                                             <  2:488976/93‑1 (MQ=255)

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

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

                                              gAATATTGTCACGGATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACAttt                                          >  1:1802939/1‑93 (MQ=255)

                                               aaTATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTc                                         >  2:3608684/1‑93 (MQ=255)

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

                                               aaTATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTc                                         >  2:2584887/1‑93 (MQ=255)

                                                atatTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCa                                        <  1:3019252/93‑1 (MQ=255)

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

                                                 tatTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCAc                                       <  1:387374/93‑1 (MQ=255)

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

                                                 tatTGTCACGTATCCGGATGGTCTGCTGACAATGTGGCCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCAc                                       <  2:3352860/93‑1 (MQ=255)

                                                   ttGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCg                                     <  2:4335961/93‑1 (MQ=255)

                                                     gTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTg                                   >  1:4176528/1‑93 (MQ=255)

                                                      tCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGa                                  >  2:919421/1‑93 (MQ=255)

                                                      tCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGa                                  >  1:1487400/1‑93 (MQ=255)

                                                       cACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGaa                                 >  2:4052007/1‑93 (MQ=255)

                                                       cACGTATCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGaa                                 >  1:2821974/1‑93 (MQ=255)

                                                            aTCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTcc                            >  2:783822/1‑93 (MQ=255)

                                                             tCCGGATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTCCg                           <  2:3320614/93‑1 (MQ=255)

                                                              ccGGATGGTCTGCTGACAATGTGGGCAGCGTGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTCCGc                          <  2:2306593/93‑1 (MQ=255)

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

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

                                                                ggATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAACATCCGCCa                        >  2:2469938/1‑93 (MQ=255)

                                                                ggATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTCCGCCa                        >  1:4636247/1‑93 (MQ=255)

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

                                                                 gATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTCCGCCAt                       >  1:256982/1‑93 (MQ=255)

                                                                 gATGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTCCGCCAt                       >  1:4241063/1‑93 (MQ=255)

                                                                  aTGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTCCGCCAtt                      <  2:128089/93‑1 (MQ=255)

                                                                  aTGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTCCGCCAtt                      <  2:4249744/93‑1 (MQ=255)

                                                                  aTGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTCCGCCAtt                      <  1:3200908/93‑1 (MQ=255)

                                                                   tGGTCTGCTGACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTCCGCCAttt                     >  1:2288486/1‑93 (MQ=255)

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

                                                                       ctgctgACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTCCGCCATttgttg                 <  2:1651828/93‑1 (MQ=255)

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

                                                                            gACAATGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTCCGCCATTTGTTGGCgca            >  1:1549486/1‑93 (MQ=255)

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

                                                                               aaTGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTCCGCCATTTGTTGGCgcagca         >  1:4195765/1‑93 (MQ=255)

                                                                               aaTGTGGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTCCGCCATTTGTTGGCgcagca         >  1:144834/1‑93 (MQ=255)

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

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

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

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

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

                                                                                  gtgGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTCCGCCATTTGTTGGCGCAGCATGa      >  2:3195561/1‑93 (MQ=255)

                                                                                   tgGGCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTCCGCCATTTGTTGGCGCAGCATGAt     >  1:1232272/1‑93 (MQ=255)

                                                                                    gggCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTCCGCCATTTGTTGGCGCAGCATGAtt    >  1:3826635/1‑93 (MQ=255)

                                                                                     ggCAGCGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCCCTCGACATTTCACCGTGAAACTCCGCCATTTGTTGGCGCAGCATGATTg   <  2:2967510/93‑1 (MQ=255)

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

TGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGTGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTGCGCACTCGACATTTCACCGTGAAACTCCGCCATTTGTTGGCGCAGCATGATTGG  >  NC_000913/3112086‑3112264

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