BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,736,520	A→G	Q218R (CAG→CGG)	purR →	PurR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,736,520	0	A	G	99.5%	667.5 / NA	205	Q218R (CAG→CGG)	purR	PurR
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (100/104); total (101/104)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 215 total aligned reads displayed.

TGGAACGTAACACCGGCGCAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCAGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGG  >  NC_000913/1736430‑1736607

tGGAACGTAACACCGGCGCAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCggg                                                                                       <  1:2820690/93‑1 (MQ=255)

tGGAACGTAACACCGGCGCAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCggg                                                                                       >  2:2843627/1‑93 (MQ=255)

tGGAACGTAACACCGGCGCAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCggg                                                                                       >  1:6099521/1‑93 (MQ=255)

   aaCGTAACACCGGCGCAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTg                                                                                    <  1:356012/93‑1 (MQ=255)

     cGTAACACCGGCGCAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGAc                                                                                  <  2:410183/93‑1 (MQ=255)

     cGTAACACCGGCGCAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGAc                                                                                  <  2:4350701/93‑1 (MQ=255)

       tAACACCGGCGCAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTTGATTGTGCGGGGTGACtt                                                                                <  1:5831731/93‑1 (MQ=255)

       tAACACCGGCGCAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGAGtt                                                                                >  2:4778449/1‑93 (MQ=255)

        aaCACCGGCGCAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACttt                                                                               <  1:4517230/93‑1 (MQ=255)

         acacCGGCGCAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTg                                                                              >  2:6440069/1‑93 (MQ=255)

           acCGGCGCAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGaa                                                                            <  2:205842/93‑1 (MQ=255)

            ccGGCGCAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGGGACTTTGAAc                                                                           >  1:3018909/1‑93 (MQ=255)

             cGGCGCAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAAcc                                                                          <  2:4111711/93‑1 (MQ=255)

                    ggCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATGAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATcc                                                                   <  2:1822767/93‑1 (MQ=255)

                    ggCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATcc                                                                   >  1:2242561/1‑93 (MQ=255)

                    ggCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATcc                                                                   <  1:5067042/93‑1 (MQ=255)

                     gccgccTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCg                                                                  >  2:2474267/1‑93 (MQ=255)

                     gccgccTTGCAGGTATTATGAAGGCGATGGAAGAAGCGATGAGCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCg                                                                  <  2:2748930/93‑1 (MQ=255)

                      ccgccTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCgg                                                                 <  1:2835616/93‑1 (MQ=255)

                      ccgccTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCgg                                                                 >  1:187359/1‑93 (MQ=255)

                      ccgccTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCgg                                                                 <  1:6624625/93‑1 (MQ=255)

                           ttGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATc                                                            <  1:2611995/93‑1 (MQ=255)

                           ttGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATc                                                            <  1:5948363/93‑1 (MQ=255)

                           ttGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATc                                                            >  1:1620863/1‑93 (MQ=255)

                           ttGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATc                                                            >  1:6357229/1‑93 (MQ=255)

                            tGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATcg                                                           <  2:6112966/93‑1 (MQ=255)

                             gCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATcgc                                                          >  2:6042529/1‑93 (MQ=255)

                             gCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATcgc                                                          <  2:6609285/93‑1 (MQ=255)

                                 gTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCa                                                      <  1:3210534/93‑1 (MQ=255)

                                  ttttATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCAt                                                     >  1:1745226/1‑93 (MQ=255)

                                   tttATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATg                                                    >  2:3572435/1‑93 (MQ=255)

                                   tttATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATg                                                    >  2:6325656/1‑93 (MQ=255)

                                    ttATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATgc                                                   >  2:1058356/1‑93 (MQ=255)

                                    ttATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATgc                                                   >  1:4236268/1‑93 (MQ=255)

                                    ttATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATgc                                                   >  1:6255629/1‑93 (MQ=255)

                                     tATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATgca                                                  <  2:2273121/93‑1 (MQ=255)

                                       tGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATgcagc                                                <  1:68842/93‑1 (MQ=255)

                                       tGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATgcagc                                                >  1:5421037/1‑93 (MQ=255)

                                        gAAGGCGTGGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGCCTTTGAACCTGAATCCGGTTATCGCGCCCTgcagca                                               >  1:3538648/1‑93 (MQ=255)

                                         aaGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCaa                                              >  1:5772859/1‑93 (MQ=255)

                                          aGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCaaa                                             >  1:6016567/1‑93 (MQ=255)

                                          aGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCaaa                                             >  2:553444/1‑93 (MQ=255)

                                          aGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCaaa                                             >  2:423541/1‑93 (MQ=255)

                                          aGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCCGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGGCATGGAGCaaa                                             <  1:2069211/93‑1 (MQ=255)

                                           ggCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAAt                                            >  1:1435875/1‑93 (MQ=255)

                                           ggCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAAt                                            <  1:4871919/93‑1 (MQ=255)

                                             cGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATcc                                          <  1:5870506/93‑1 (MQ=255)

                                              gATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCt                                         <  1:4161469/93‑1 (MQ=255)

                                               aTGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTg                                        >  2:3224301/1‑93 (MQ=255)

                                                tGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGt                                       >  2:162784/1‑93 (MQ=255)

                                                tGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTAGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGt                                       <  2:1541152/93‑1 (MQ=255)

                                                 ggAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAGCCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTc                                      <  2:3285481/93‑1 (MQ=255)

                                                  gaagaaGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCg                                     >  1:3255209/1‑93 (MQ=255)

                                                   aagaagCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGc                                    <  1:2543253/93‑1 (MQ=255)

                                                   aagaagCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGc                                    <  2:1101448/93‑1 (MQ=255)

                                                   aagaagCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGc                                    <  1:3305980/93‑1 (MQ=255)

                                                      aagCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGc                                 <  1:4323964/93‑1 (MQ=255)

                                                      aagCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGc                                 <  1:1059857/93‑1 (MQ=255)

                                                       agCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGcc                                <  2:1581756/93‑1 (MQ=255)

                                                       agCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGcc                                >  1:117815/1‑93 (MQ=255)

                                                       agCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGcc                                >  1:61424/1‑93 (MQ=255)

                                                       agCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGcc                                >  2:3861244/1‑93 (MQ=255)

                                                         cGATGATCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGc                              <  1:551898/93‑1 (MQ=255)

                                                            tgatCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATc                           >  2:257777/1‑93 (MQ=255)

                                                             gatCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCg                          >  1:5689461/1‑93 (MQ=255)

                                                              atCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGc                         >  1:922913/1‑93 (MQ=255)

                                                               tCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGcc                        >  2:775184/1‑93 (MQ=255)

                                                               tCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGcc                        >  2:4914223/1‑93 (MQ=255)

                                                               tCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGcc                        >  2:4676113/1‑93 (MQ=255)

                                                               tCAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGcc                        >  2:2133674/1‑93 (MQ=255)

                                                                cAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGccc                       >  2:4397843/1‑93 (MQ=255)

                                                                cAAGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGccc                       >  2:5778096/1‑93 (MQ=255)

                                                                 aaGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCt                      <  2:1223005/93‑1 (MQ=255)

                                                                  aGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTa                     >  1:1289387/1‑93 (MQ=255)

                                                                  aGGTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTa                     >  1:6720378/1‑93 (MQ=255)

                                                                  aGGTGCCGGAAAGCTGGATTGTGCGGGGCGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAGATCCTGTCGCAGCCGCATCGCCCTa                     >  2:6479165/1‑93 (MQ=255)

                                                                  aGGTGCCGGAAAGCTGGATTGTGCAGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTa                     >  2:2998418/1‑93 (MQ=255)

                                                                   ggTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTAc                    <  2:2628160/93‑1 (MQ=255)

                                                                   ggTGCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTAc                    <  1:3754759/93‑1 (MQ=255)

                                                                      gCCGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGc                 <  2:2133114/93‑1 (MQ=255)

                                                                       ccGGAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGcc                <  1:6042302/93‑1 (MQ=255)

                                                                          gAAAGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGtc             <  2:5938478/93‑1 (MQ=255)

                                                                           aaaGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGtct            <  2:152890/93‑1 (MQ=255)

                                                                           aaaGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGtct            <  1:2304683/93‑1 (MQ=255)

                                                                            aaGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGtctt           <  1:1966564/93‑1 (MQ=255)

                                                                             aGCTGGATTGTGCGGTGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGtcttc          <  2:6696810/93‑1 (MQ=255)

                                                                             aGCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGtcttc          <  2:3475624/93‑1 (MQ=255)

                                                                              gCTGGATTGTGCGGGGTGACTTTGACCCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTATCGCAGCCGCATCGCCCTACTGCCGtcttct         <  1:1671363/93‑1 (MQ=255)

                                                                              gCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGtcttct         >  1:3049439/1‑93 (MQ=255)

                                                                              gCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGtcttct         <  2:2819372/93‑1 (MQ=255)

                                                                              gCTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGtcttct         >  2:3713284/1‑93 (MQ=255)

                                                                               cTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCtg        <  1:2330243/93‑1 (MQ=255)

                                                                               cTGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCtg        <  1:4060274/93‑1 (MQ=255)

                                                                                tGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCtgt       >  2:1998044/1‑93 (MQ=255)

                                                                                tGGATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCtgt       <  2:341841/93‑1 (MQ=255)

                                                                                  gATTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTgtgg     >  2:172854/1‑93 (MQ=255)

                                                                                   aTTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTgtggt    >  1:3915430/1‑93 (MQ=255)

                                                                                   aTTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTgtggt    >  2:6215834/1‑93 (MQ=255)

                                                                                   aTTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTgtggt    <  1:1933512/93‑1 (MQ=255)

                                                                                   aTTGTGCGGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTgtggg    >  2:5273425/1‑92 (MQ=255)

TGGAACGTAACACCGGCGCAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGATTGTGCAGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTCGCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGG  >  NC_000913/1736430‑1736607

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