BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	332,019	T→C	V142A (GTG→GCG)	yahA →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	332,019	0	T	C	99.4%	535.8 / ‑5.0	168	V142A (GTG→GCG)	yahA	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (80/87); total (80/88)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.83e-01

Only 100 of 177 total aligned reads displayed.

CCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAA  >  NC_000913/331929‑332108

ccATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGgcgc                                                                                         <  1:3313540/93‑1 (MQ=255)

ccATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGgcgc                                                                                         >  1:755550/1‑93 (MQ=255)

 cATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCt                                                                                        >  2:840919/1‑93 (MQ=255)

 cATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCt                                                                                        >  2:4437452/1‑93 (MQ=255)

 cATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCt                                                                                        >  1:4179761/1‑93 (MQ=255)

 cATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCt                                                                                        >  1:574458/1‑93 (MQ=255)

 cATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCt                                                                                        >  2:2486180/1‑93 (MQ=255)

 cATCAGTCTGGAGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCt                                                                                        >  2:193901/1‑93 (MQ=255)

  aTCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCtt                                                                                       >  1:3755073/1‑93 (MQ=255)

  aTCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCtt                                                                                       <  2:1629624/93‑1 (MQ=255)

  aTCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCtt                                                                                       <  2:21667/93‑1 (MQ=255)

  aTCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCtt                                                                                       >  1:140473/1‑93 (MQ=255)

  aTCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCtt                                                                                       >  2:2938213/1‑93 (MQ=255)

   tCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTg                                                                                      >  1:4638431/1‑93 (MQ=255)

        cTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGc                                                                                 <  1:2123673/93‑1 (MQ=255)

         tGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCt                                                                                <  1:3658873/93‑1 (MQ=255)

          ggCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTg                                                                               >  2:273214/1‑93 (MQ=255)

          ggCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTg                                                                               >  1:4140994/1‑93 (MQ=255)

             gTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGa                                                                            >  2:2434017/1‑93 (MQ=255)

             gTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGa                                                                            >  1:4648969/1‑93 (MQ=255)

                ggAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACa                                                                         >  2:213866/1‑93 (MQ=255)

                ggAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACa                                                                         <  1:4147230/93‑1 (MQ=255)

                  aaaaCCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATc                                                                       <  2:346532/93‑1 (MQ=255)

                    aaCCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCca                                                                     >  2:2367725/1‑93 (MQ=255)

                      ccACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCcaca                                                                   >  1:2527483/1‑93 (MQ=255)

                           aaTTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACgg                                                              <  1:2008455/93‑1 (MQ=255)

                           aaTTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTGGTCCGCTGGGAACATCCACAAACgg                                                              <  2:1046495/93‑1 (MQ=255)

                            aTTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACTGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACggg                                                             <  2:428849/93‑1 (MQ=255)

                            aTTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACggg                                                             >  2:575684/1‑93 (MQ=255)

                            aTTCAAACCGTGGATCCAACCGGTTTTCCGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACggg                                                             <  2:571455/93‑1 (MQ=255)

                             ttCAAACCGTGGATCCACCCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGa                                                            <  1:1769530/93‑1 (MQ=255)

                             ttCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAATGGGa                                                            >  1:3567428/1‑93 (MQ=255)

                             ttCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGa                                                            <  1:977582/93‑1 (MQ=255)

                             ttCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGa                                                            >  2:3182440/1‑93 (MQ=255)

                             ttCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGa                                                            >  1:1900379/1‑93 (MQ=255)

                               cAAACCGTGGATCCAACCGGTTTTCGGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAAt                                                          <  1:4656090/93‑1 (MQ=255)

                                aaaCCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAAtt                                                         <  2:1900379/93‑1 (MQ=255)

                                  aCCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTAt                                                       <  2:1166270/93‑1 (MQ=255)

                                    cGTGGATCACACCTGTTTTGTGCGCGGAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATcc                                                     <  1:2664813/93‑1 (MQ=255)

                                     gTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATccc                                                    >  1:3743701/1‑93 (MQ=255)

                                     gTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATccc                                                    >  2:816854/1‑93 (MQ=255)

                                      tGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCa                                                   <  2:1749543/93‑1 (MQ=255)

                                      tGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCa                                                   >  2:1968667/1‑93 (MQ=255)

                                       ggATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCAc                                                  <  1:3272842/93‑1 (MQ=255)

                                        gATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCAcc                                                 <  2:2214602/93‑1 (MQ=255)

                                        gATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCAcc                                                 <  1:1256397/93‑1 (MQ=255)

                                        gATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCAcc                                                 <  1:4490498/93‑1 (MQ=255)

                                        gATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCAcc                                                 >  1:452121/1‑93 (MQ=255)

                                         aTCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCg                                                <  1:1867000/93‑1 (MQ=255)

                                            cAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGAt                                             <  1:1169806/93‑1 (MQ=255)

                                             caCCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATc                                            <  1:404976/92‑1 (MQ=255)

                                             aaCCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATc                                            >  1:908811/1‑93 (MQ=255)

                                               ccGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAg                                          >  2:4573739/1‑93 (MQ=255)

                                               ccGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAg                                          >  1:1008170/1‑93 (MQ=255)

                                                cGGTTTGCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGt                                         <  2:1345974/93‑1 (MQ=255)

                                                cGGGTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTACGCTGGGAACATCCACAAACGGGAATTATCCCACCAGATCAGt                                         >  2:4653081/1‑93 (MQ=255)

                                                 ggTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGtt                                        >  2:3547052/1‑93 (MQ=255)

                                                 ggTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGtt                                        >  1:179499/1‑93 (MQ=255)

                                                 cggtttCTGCGCGCACACTTGCGTACTGACGGGCTGTGAGGCGCTTGGCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGTATCAGtt                                        <  1:3909279/90‑1 (MQ=255)

                                                  gTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGttt                                       <  2:3827672/93‑1 (MQ=255)

                                                  gTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGttt                                       >  1:2906633/1‑93 (MQ=255)

                                                  gTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGttt                                       <  1:3922644/93‑1 (MQ=255)

                                                  gTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGttt                                       >  1:221799/1‑93 (MQ=255)

                                                   ttttCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTa                                      <  2:1851556/93‑1 (MQ=255)

                                                    tttCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTAt                                     <  1:3417087/93‑1 (MQ=255)

                                                     ttCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTAtt                                    >  2:714870/1‑93 (MQ=255)

                                                          cgcgcAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCctct                               <  1:1326712/93‑1 (MQ=255)

                                                           gcgcAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTg                              >  1:822164/1‑93 (MQ=255)

                                                            cgcAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTgg                             <  1:3989797/93‑1 (MQ=255)

                                                             gcAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGc                            <  1:925065/93‑1 (MQ=255)

                                                               aGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGGACATCCACAAACGGGAATTATCCCACCGGATCAGATTATTCCTCTGGCgg                          >  2:3386318/1‑93 (MQ=255)

                                                               aGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCgg                          >  1:353707/1‑93 (MQ=255)

                                                               aGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCgg                          >  2:382430/1‑93 (MQ=255)

                                                                gACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGa                         >  2:4431224/1‑93 (MQ=255)

                                                                gACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGa                         <  1:4388205/93‑1 (MQ=255)

                                                                 aCTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAg                        <  1:647511/93‑1 (MQ=255)

                                                                 aCTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAg                        <  1:647510/93‑1 (MQ=255)

                                                                 aCTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCCGGGAACATCCACAAACGGGAATTATCCCACCGGATCCGTTTATTCCTCTGgaggag                        >  2:1210098/1‑93 (MQ=255)

                                                                   tGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTc                      <  2:2481247/93‑1 (MQ=255)

                                                                   tGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTc                      <  1:286529/93‑1 (MQ=255)

                                                                    ggCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCa                     >  2:395742/1‑93 (MQ=255)

                                                                    ggCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCa                     <  2:185389/93‑1 (MQ=255)

                                                                     gCGTAATGACGGGTTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCAt                    <  2:1336048/93‑1 (MQ=255)

                                                                       gTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATcc                  <  1:1037744/93‑1 (MQ=255)

                                                                       gTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATcc                  <  2:3456046/93‑1 (MQ=255)

                                                                        tACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCg                 >  2:1596561/1‑93 (MQ=255)

                                                                           tGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTc              >  2:1081859/1‑93 (MQ=255)

                                                                            gACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCt             <  2:3368326/93‑1 (MQ=255)

                                                                            gACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCt             <  2:2294204/93‑1 (MQ=255)

                                                                                ggCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTAtt         <  2:3167393/93‑1 (MQ=255)

                                                                                ggCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTAtt         <  2:3899684/93‑1 (MQ=255)

                                                                                ggCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTAtt         >  1:3783248/1‑93 (MQ=255)

                                                                                 gCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTg        >  1:2929246/1‑93 (MQ=255)

                                                                                   tgtgAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTc      <  2:2843308/93‑1 (MQ=255)

                                                                                   tgtgAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTc      >  2:2637828/1‑93 (MQ=255)

                                                                                    gtgAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCa     <  1:2356922/93‑1 (MQ=255)

                                                                                    gtgAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCa     >  2:1684440/1‑93 (MQ=255)

                                                                                     tgAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCAt    <  1:3589492/93‑1 (MQ=255)

                                                                                     tgAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGACGAGTCATCCGGTCTTATTGTCAt    >  2:2622114/1‑93 (MQ=255)

                                                                                       aGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATaa  >  2:3277598/1‑93 (MQ=255)

CCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAA  >  NC_000913/331929‑332108

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