BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	304,068	A→G	V136A (GTC→GCC)	yagV ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	304,068	0	A	G	99.3%	461.3 / NA	152	V136A (GTC→GCC)	yagV	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/0); major base G (64/87); minor base T (0/1); total (64/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 = 9.99e-01

Only 100 of 159 total aligned reads displayed.

GTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTCTCACCAGCCTG  >  NC_000913/303979‑304156

gTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGGCCt                                                                                       >  2:1307982/1‑93 (MQ=255)

gTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGGCCt                                                                                       <  1:505727/93‑1 (MQ=255)

 tactacCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTTGGGGCCTc                                                                                      >  2:2091953/1‑93 (MQ=255)

 tactacCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGGCCTc                                                                                      <  1:3067062/93‑1 (MQ=255)

 tactacCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGGCCTc                                                                                      >  1:427218/1‑93 (MQ=255)

 tactacCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGGCCTc                                                                                      >  2:4316763/1‑93 (MQ=255)

   ctacCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGGCCTcgc                                                                                    <  1:355744/93‑1 (MQ=255)

   ctacCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGGCCTcgc                                                                                    <  1:2713525/93‑1 (MQ=255)

   ctacCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGGCCTcgc                                                                                    <  1:4029483/93‑1 (MQ=255)

    tacCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGGCCTcgcg                                                                                   >  2:431937/1‑93 (MQ=255)

     acCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGGCCTCGCGa                                                                                  >  1:1629582/1‑93 (MQ=255)

     acCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGGCCTCGCGa                                                                                  <  2:3635759/93‑1 (MQ=255)

       cAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGGCCTCGCGaaa                                                                                <  1:1966849/93‑1 (MQ=255)

         gAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGGCCTCGCGAAATg                                                                              >  1:666012/1‑93 (MQ=255)

           aTGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGGCCTCGCGAAATgag                                                                            >  2:228282/1‑93 (MQ=255)

            tGGTATCCATCCCCACCGCCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGGCCTCGCGAAATgaga                                                                           <  2:1798299/93‑1 (MQ=255)

            tGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGGCCTCGCGAAATgaga                                                                           >  1:3787235/1‑93 (MQ=255)

            tGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGGCCTCGCGAAATgaga                                                                           >  2:442305/1‑93 (MQ=255)

              gTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGGCCTCGCGAAATGAGAcc                                                                         <  1:4374889/93‑1 (MQ=255)

               tATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGGCCTCGCGAAATGAGAccc                                                                        >  2:609801/1‑93 (MQ=255)

                aTCCATCACCACCACCGGCTCCGTTCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGGCCTCCCGAAATGAGACCCg                                                                       >  2:2955936/1‑93 (MQ=255)

                aTCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGGCCTCGCGAAATGAGACCCg                                                                       >  2:2546040/1‑93 (MQ=255)

                 tCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGGCCTCGCGAAATGAGACCCgg                                                                      >  1:2750772/1‑93 (MQ=255)

                 tCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGGCCTCGCGAAATGAGACCCgg                                                                      <  1:3460589/93‑1 (MQ=255)

                    aTCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGGCCTCGCGAAATGAGACCCGgtag                                                                   >  2:3664957/1‑93 (MQ=255)

                     tCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGGCCTCGCGAAATGAGACCCGgtaga                                                                  >  1:586168/1‑92 (MQ=255)

                       accaccaccGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGGCCTCGCGAAATGAGACCCGGtagtag                                                                >  1:2265401/1‑93 (MQ=255)

                       accaccaccGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGGCCTCGCGAAATGAGACCCGGtagtag                                                                >  1:3918012/1‑93 (MQ=255)

                        ccaccaccGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGGCCTCGCGAAATGAGACCCGGTAGTAgc                                                               <  2:1466423/93‑1 (MQ=255)

                         caccaccGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGGCCTCGCGAAATGAGACCCGGTAGTAgcg                                                              <  2:2553771/93‑1 (MQ=255)

                          accaccGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGGCCTCGCGAAATGAGACCCGGTAGTAgcgc                                                             <  2:4459612/93‑1 (MQ=255)

                          accaccGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGGCCTCGCGAAATGAGACCCGGTAGTAgcgc                                                             >  2:691428/1‑93 (MQ=255)

                            caccGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGCGTGGGGGCCTCGCGAAATGCGACCCGGTAGTAGCGCTc                                                           >  1:644665/1‑93 (MQ=255)

                            caccGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGGCCTCGCGAAATGAGACCCGGTAGTAGCGCTc                                                           >  2:3251972/1‑93 (MQ=255)

                              ccGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGGCCTCGCGAAATGAGACCCGGTAGTAGCGCTcgc                                                         >  2:2719587/1‑93 (MQ=255)

                              ccGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGGCCTCGCGAAATGAGACCCGGTAGTAGCGCTcgc                                                         <  1:4667883/93‑1 (MQ=255)

                              ccGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGGCCTCGCGAAATGAGACCCGGTAGTAGCGCTcgc                                                         >  2:4197563/1‑93 (MQ=255)

                              ccGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGGCCTCGCGAAATGAGACCCGGTAGTAGCGCTcgc                                                         >  1:1721046/1‑93 (MQ=255)

                                ggCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGGCCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCgg                                                       >  1:3427428/1‑93 (MQ=255)

                                 gCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGGCCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGt                                                      >  1:4136964/1‑93 (MQ=255)

                                 gCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGGCCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGt                                                      >  2:4173638/1‑93 (MQ=255)

                                 gCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGGCCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGt                                                      <  1:1131826/93‑1 (MQ=255)

                                 gCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGGCCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGt                                                      >  1:1785397/1‑93 (MQ=255)

                                  cTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGGCCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGtt                                                     >  2:2774272/1‑93 (MQ=255)

                                   tCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGGCCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTa                                                    >  1:2423410/1‑93 (MQ=255)

                                   tCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGGCCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTa                                                    <  2:3136297/93‑1 (MQ=255)

                                    ccGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGGCCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTAt                                                   >  1:3999984/1‑93 (MQ=255)

                                    ccGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGGCCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTAt                                                   <  1:2998905/93‑1 (MQ=255)

                                     cGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGGCCGCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATc                                                  <  2:2993934/93‑1 (MQ=255)

                                      gTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGGCCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATcc                                                 <  1:2392118/93‑1 (MQ=255)

                                      gTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGGCCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATcc                                                 <  2:1395896/93‑1 (MQ=255)

                                       tGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGGCCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATccc                                                >  1:4538320/1‑93 (MQ=255)

                                       tGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGGCCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATccc                                                >  1:985942/1‑93 (MQ=255)

                                       tGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGCGGGGGCCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATccc                                                >  2:3100007/1‑93 (MQ=255)

                                        gCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGGCCTCGCGAAATGAGACCCGGGAGTAGCGCTCGCGGTTATCCCg                                               >  2:1187323/1‑93 (MQ=255)

                                          tGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGGCCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTg                                             <  1:4290445/93‑1 (MQ=255)

                                            aCCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGGCCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGa                                           <  2:1475957/93‑1 (MQ=255)

                                            aCCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGGCCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGa                                           >  1:1701117/1‑93 (MQ=255)

                                             ccACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGGCCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGAc                                          <  2:3390880/93‑1 (MQ=255)

                                             ccACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGGCCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGAc                                          >  1:2336028/1‑93 (MQ=255)

                                             ccACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGGCCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGAc                                          <  2:4600343/93‑1 (MQ=255)

                                               aagccgccGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGGCCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCa                                        >  2:2960111/3‑93 (MQ=255)

                                               aCGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGGCCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCa                                        >  1:2802526/1‑93 (MQ=255)

                                                cgccgcCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGGCCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCAt                                       <  2:4446692/93‑1 (MQ=255)

                                                cgccgcCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGGCCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCAt                                       <  1:1672186/93‑1 (MQ=255)

                                                cgccgcCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGGCCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCAt                                       <  2:814833/93‑1 (MQ=255)

                                                   cgccgGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGGCCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGAt                                    >  2:1292385/1‑93 (MQ=255)

                                                   cgccgGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGGCCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGAt                                    <  1:4562143/93‑1 (MQ=255)

                                                    gccgGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGGCCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATa                                   >  2:3024058/1‑93 (MQ=255)

                                                    gccgGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGGCCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATa                                   <  1:3458723/93‑1 (MQ=255)

                                                    gccgGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGGCCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATa                                   <  2:3271512/93‑1 (MQ=255)

                                                    gccgGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGGCCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATa                                   <  2:1053312/93‑1 (MQ=255)

                                                       ggTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGGCCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTa                                >  2:1922759/1‑93 (MQ=255)

                                                        gTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGGCCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTaa                               >  2:3364151/1‑93 (MQ=255)

                                                           gggCTGCGTCTTGTATGGTTACGAGTGGGGGCCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAAt                            <  2:683820/93‑1 (MQ=255)

                                                           gggCTGCGTCTTGTATGGTTACGAGTGGGGGCCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAAt                            <  2:3416552/93‑1 (MQ=255)

                                                               tGCGTCTTGTATGGTTACGAGTGGGGGCCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTaa                        <  2:4673576/93‑1 (MQ=255)

                                                                gCGTCTTGTATGGTTACGAGTGGGGGCCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTaaa                       >  1:2821086/1‑93 (MQ=255)

                                                                 cGTCTTGTATGGTTACGAGTGGGGGCCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTaaaa                      <  1:147953/93‑1 (MQ=255)

                                                                  gTCTTGTATGGTTACGAGTGGGGGCCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAAt                     <  1:4512047/93‑1 (MQ=255)

                                                                  gTCTTGTATGGTTACGAGTGGGGGCCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAAt                     <  1:2556908/93‑1 (MQ=255)

                                                                  gTCTTGTATGGTTACGAGTGGGGGCCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAAt                     <  2:4140027/93‑1 (MQ=255)

                                                                  gTCTTGTATGGTTACGAGTGGGGGCCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAAt                     <  2:3207171/93‑1 (MQ=255)

                                                                    cTTGTATGGTTACGAGTGGGGGCCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATAc                   >  2:2428885/1‑93 (MQ=255)

                                                                    cTTGTATGGTTACGAGTGGGGGCCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATAc                   >  1:154354/1‑93 (MQ=255)

                                                                    cTTGTATGGTTACGAGTGGGGGCCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATAc                   <  1:4487349/93‑1 (MQ=255)

                                                                    cTTGTATGGTTACGAGTGGGGGCCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATAc                   <  2:1620485/93‑1 (MQ=255)

                                                                    cTTGTATGGTTACGAGTGGGGGCCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATAc                   <  1:2974472/93‑1 (MQ=255)

                                                                      tGTATGGTTACGAGTGGGGGCCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTc                 <  2:4403721/93‑1 (MQ=255)

                                                                      tGTATGGTTACGAGTGGGGGCCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTc                 <  1:1998156/93‑1 (MQ=255)

                                                                        tATGGTTACGAGTGGGGGCCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGc               <  2:1781133/93‑1 (MQ=255)

                                                                          tGGTTACGAGTGGGGGCCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGctc             >  2:4437158/1‑93 (MQ=255)

                                                                             gtACGAGTGGGGGCCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTCTCa          <  2:1522099/92‑1 (MQ=255)

                                                                              tACGAGTGGGGGCCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTCTCAc         <  2:1260256/93‑1 (MQ=255)

                                                                              tACGAGTGGGGGCCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTCTCAc         >  1:1603469/1‑93 (MQ=255)

                                                                              tACGAGTGGGGGCCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTCTCAc         <  1:411158/93‑1 (MQ=255)

                                                                              tACGAGTGGGGGCCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTCTCAc         <  2:3470924/93‑1 (MQ=255)

                                                                                   gTGGGGGCCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTCTCACCAGcc    <  1:2198409/93‑1 (MQ=255)

                                                                                   gTGGGGGCCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTCTCACCAGcc    >  2:3519574/1‑93 (MQ=255)

                                                                                     ggggGCCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTCTCACCAGCCTg  >  2:2908455/1‑93 (MQ=255)

GTACTACCAGAATGGTATCCATCACCACCACCGGCTCCGTGCTGACCACGCCGCCGGTTGGGCTGCGTCTTGTATGGTTACGAGTGGGGACCTCGCGAAATGAGACCCGGTAGTAGCGCTCGCGGTTATCCCGTGGACCATGATAGTAAAATTTAAAATACTCGCTCTCACCAGCCTG  >  NC_000913/303979‑304156

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