BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,713,746	C→T	P56P (CCG→CCA)	pdxY ←	pyridoxal kinase 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,713,746	0	C	T	95.5%	415.0 / 2.5	132	P56P (CCG→CCA)	pdxY	pyridoxal kinase 2
Reads supporting (aligned to +/- strand): ref base C (5/1); new base T (80/46); total (85/47)
Fisher's exact test for biased strand distribution p-value = 4.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 139 total aligned reads displayed.

CGATCCCAGATAGCCACTTAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAGGCGGCG  >  NC_000913/1713656‑1713835

cgATCCCAGATAGCCACTTAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGggtgg                                                                                         >  1:3139693/1‑93 (MQ=255)

 gATCCCAGATAGCCACTTAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGggtggc                                                                                        >  2:1955330/1‑93 (MQ=255)

 gATCCCAGATAGCCACTTAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGggtggc                                                                                        >  1:2668214/1‑93 (MQ=255)

  aTCCCAGATAGCCACTTAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCa                                                                                       >  1:1071612/1‑93 (MQ=255)

   tCCCAGATAGCCACTTAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCAt                                                                                      >  2:1504130/1‑93 (MQ=255)

      cAGATAGCCACTTAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCAc                                                                                   >  2:1971689/1‑93 (MQ=255)

       aGATAGCCACTTAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGGTGGGTGGCATCACg                                                                                  >  2:157861/1‑93 (MQ=255)

       aGATAGCCACTTAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACg                                                                                  >  2:1937674/1‑93 (MQ=255)

        gATAGCCACTTAATACGGCATCACAGGTGTGTAATTTATCAATGGTGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGc                                                                                 >  2:2773966/1‑93 (MQ=255)

        gATAGCCACTTAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGc                                                                                 >  2:2668694/1‑93 (MQ=255)

           aGCCACTTAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGc                                                                              <  2:3938064/93‑1 (MQ=255)

           aGCCACTTAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGGTAAATGGCTGGGTGGCATCACGCAGc                                                                              >  2:4207367/1‑93 (MQ=255)

            gCCACTTAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGcc                                                                             >  1:4028985/1‑93 (MQ=255)

            gCCACTTAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGcc                                                                             >  2:1083779/1‑93 (MQ=255)

            gCCACTTAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGcc                                                                             >  2:3383642/1‑93 (MQ=255)

            gCCACTTAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGcc                                                                             >  1:2220717/1‑93 (MQ=255)

             ccACTTAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCa                                                                            <  1:3130294/93‑1 (MQ=255)

               aCTTAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGt                                                                          >  2:3510827/1‑93 (MQ=255)

                cTTAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTc                                                                         <  1:1019476/93‑1 (MQ=255)

                cTTAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGCGGCATCACGCAGACAGTc                                                                         >  1:3829530/1‑93 (MQ=255)

                cTTAATACGGCATCACAGGTATGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTc                                                                         <  2:76159/93‑1 (MQ=255)

                 ttAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTcc                                                                        >  1:2319772/1‑93 (MQ=255)

                  tAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCa                                                                       <  1:3030809/93‑1 (MQ=255)

                  tAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCa                                                                       >  1:1478298/1‑93 (MQ=255)

                   aaTACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCAt                                                                      >  2:270027/1‑93 (MQ=255)

                     tACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCAttt                                                                    >  2:4689287/1‑93 (MQ=255)

                     tACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCAttt                                                                    >  1:810781/1‑93 (MQ=255)

                      aCGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTg                                                                   >  2:4137828/1‑93 (MQ=255)

                        ggCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGcc                                                                 >  1:2428783/1‑93 (MQ=255)

                        ggCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGcc                                                                 <  1:2776361/93‑1 (MQ=255)

                         gCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCg                                                                >  1:4750173/1‑93 (MQ=255)

                           aTCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTa                                                              <  2:4390271/93‑1 (MQ=255)

                           aTCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTa                                                              <  1:1924223/93‑1 (MQ=255)

                            tCACAGGTGTGTAATTTATCAATGGCGGCGATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTAt                                                             >  2:788485/1‑93 (MQ=255)

                            tCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTAt                                                             <  1:2818246/93‑1 (MQ=255)

                            tCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTAt                                                             <  2:2831004/93‑1 (MQ=255)

                             cacaGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTGCATTTGCCGTAtt                                                            >  1:822321/1‑93 (MQ=255)

                             cacaGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTAtt                                                            <  1:1177282/93‑1 (MQ=255)

                             cacaGGAGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTAtt                                                            >  2:3381729/1‑93 (MQ=255)

                              acaGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTg                                                           >  2:1001662/1‑93 (MQ=255)

                               caGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTgg                                                          >  2:2888745/1‑93 (MQ=255)

                                aGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTggg                                                         <  1:1920911/93‑1 (MQ=255)

                                 ggTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGgt                                                        >  1:3057705/1‑93 (MQ=255)

                                  gtgtgtAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGgtg                                                       >  1:4358796/1‑93 (MQ=255)

                                     tgtAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGa                                                    >  2:2908302/1‑93 (MQ=255)

                                     tgtAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGa                                                    >  1:42030/1‑93 (MQ=255)

                                     tgtAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGa                                                    >  2:1508346/1‑93 (MQ=255)

                                      gtAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGAt                                                   <  1:4306218/93‑1 (MQ=255)

                                      gtAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCAGCACGCAGCCAGTCCAGTTGCCGTATCGTGTGCGAt                                                   >  1:3542193/1‑93 (MQ=255)

                                       tAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGAtt                                                  >  2:2345052/1‑93 (MQ=255)

                                               cAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAAtt                                          >  2:2551693/1‑93 (MQ=255)

                                               cAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAAtt                                          <  2:1828020/93‑1 (MQ=255)

                                                aaTGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTg                                         >  2:1344633/1‑93 (MQ=255)

                                                 aTGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAAttta                                        >  1:4070157/1‑91 (MQ=255)

                                                   ggcggcAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGaag                                      >  1:874380/1‑92 (MQ=255)

                                                   ggcggcAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAAc                                      >  2:4075599/1‑93 (MQ=255)

                                                   ggcggcAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAAc                                      <  1:2844144/93‑1 (MQ=255)

                                                   ggcggcAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAAc                                      >  1:649606/1‑93 (MQ=255)

                                                   ggcggcAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAAc                                      <  2:1906257/93‑1 (MQ=255)

                                                    gcggcAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACg                                     <  2:4122943/93‑1 (MQ=255)

                                                    gcggcAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACg                                     <  2:4493281/93‑1 (MQ=255)

                                                     cgtcAATGCCTCGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACgg                                    <  2:2135233/90‑1 (MQ=255)

                                                       gcAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGgtg                                  >  1:1487966/1‑93 (MQ=255)

                                                       gcAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGgtg                                  >  2:2492381/1‑93 (MQ=255)

                                                        cAATGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGgtgt                                 <  2:1260184/93‑1 (MQ=255)

                                                         aaTGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGtt                                >  2:47708/1‑93 (MQ=255)

                                                          aTGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTc                               <  1:3304269/93‑1 (MQ=255)

                                                          aTGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTc                               <  1:4485036/93‑1 (MQ=255)

                                                           tGCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCa                              <  1:2696137/93‑1 (MQ=255)

                                                            gCCTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAg                             >  2:678841/1‑93 (MQ=255)

                                                              cTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCg                           >  1:4064234/1‑93 (MQ=255)

                                                              cTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCg                           >  1:2580775/1‑93 (MQ=255)

                                                              cTTGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCg                           >  2:3132425/1‑93 (MQ=255)

                                                               ttGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCgg                          >  1:3356171/1‑93 (MQ=255)

                                                               ttGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCgg                          >  2:368586/1‑93 (MQ=255)

                                                               ttGCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGACAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCgg                          >  1:3113083/1‑93 (MQ=255)

                                                                 gCACAATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTTAACGGTGTTCAGCGGcc                        >  1:408360/1‑93 (MQ=255)

                                                                  cacaATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCa                       >  1:3860015/1‑93 (MQ=255)

                                                                  cacaATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCa                       <  1:2131445/93‑1 (MQ=255)

                                                                   acaATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAg                      >  2:1013909/1‑93 (MQ=255)

                                                                   acaATTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAg                      >  1:2346876/1‑93 (MQ=255)

                                                                     aaTTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGAc                    <  2:2920051/93‑1 (MQ=255)

                                                                     aaTTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGAc                    >  1:1708339/1‑93 (MQ=255)

                                                                     aaTTTCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGAc                    <  2:3458445/93‑1 (MQ=255)

                                                                         tCGGTTAAATGGCTGTGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGTTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTc                <  2:2672870/93‑1 (MQ=255)

                                                                         tCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTc                >  2:4490362/1‑93 (MQ=255)

                                                                         tCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTc                <  2:3341155/93‑1 (MQ=255)

                                                                         tCGGTTAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTc                <  1:1609675/93‑1 (MQ=255)

                                                                             ttAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCgcgc            >  2:534941/1‑93 (MQ=255)

                                                                             ttAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGAGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCgcgc            >  1:1455623/1‑93 (MQ=255)

                                                                              tAAATGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGcc           >  1:3327643/1‑93 (MQ=255)

                                                                                 aTGGCTGGGTGGCATCACGCAGCCAGTCGATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAg        <  2:2545300/93‑1 (MQ=255)

                                                                                  tGGCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAgg       <  1:4263795/93‑1 (MQ=255)

                                                                                   ggCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAggc      >  1:614027/1‑93 (MQ=255)

                                                                                   ggCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAggc      >  2:3094644/1‑93 (MQ=255)

                                                                                   ggCTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAggc      >  1:4660133/1‑93 (MQ=255)

                                                                                     cTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAggcgg    >  2:3805974/1‑93 (MQ=255)

                                                                                     cTGGGTGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAggcgg    <  2:604608/93‑1 (MQ=255)

                                                                                     cTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAggcgc    >  1:3866410/1‑92 (MQ=255)

                                                                                       gggCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGCGTTCACCGGCCAGACGTTCGCGCCCCGgcggcg  >  2:4006147/1‑93 (MQ=255)

CGATCCCAGATAGCCACTTAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAGGCGGCG  >  NC_000913/1713656‑1713835

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