BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,206,988	A→G	T442A (ACC→GCC)	ygjE →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,206,988	0	A	G	82.8%	506.8 / 99.5	238	T442A (ACC→GCC)	ygjE	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (15/26); new base G (79/118); total (94/144)
Fisher's exact test for biased strand distribution p-value = 7.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

Only 100 of 244 total aligned reads displayed.

CGCTGGCGATGCCGGAAATCCCGCTGCCGGTGTTCTGCCTGATGGTTGGTGCGGCAATTGGTCTGGGGAGCATTCTTACACCATACGCCACCGGCCCCAGTCCGATTTACTACGGCAGTGGTTATCTGCCAACGGCGGATTACTGGCGACTGGGGGCGATTTTTGGGCTGATATTCCTCGTA  >  NC_000913/3206899‑3207080

cgcTGGCGATGCCGGAAATCCCGCTGCCGGTGTTCTGCCTGATGGTTGGTGCGGCAATTGGTCTGGGGAGCATTCTTACACCATACGccgccg                                                                                           <  1:1179581/93‑1 (MQ=255)

  cTGGCGATGCCGGAAATCCCGCTGCCGGTGTTCTGCCTGATGGTTGGTGCGGCAATTGGTCTGGGGAGCATTCTTACACCATACGCCGCCGGc                                                                                         >  2:3675011/1‑93 (MQ=255)

   tGGCGATGCCGGAAATCCCGCTGCCGGTGTTCTGCCTGATGGTTGGTGCGGCAATTGGTCTGGGGAGCATTCTTACACCATACGCCGCCGGcc                                                                                        >  1:4487451/1‑93 (MQ=255)

   tGGCGATGCCGGAAATCCCGCTGCCGGTGTTCTGCCTGATGGTTGGTGCGGCAATTGGTCTGGGGAGCATTCTTACACCATACGCCGCCGGcc                                                                                        >  2:4677975/1‑93 (MQ=255)

    ggCGATGCCGGAAATCCCGCTGCCGGTGTTCTGCCTGATGGTTGGTGCGGCAATTGGTCTGGGGAGCATTCTTACACCATACGCCGCCGGccc                                                                                       <  1:6738532/93‑1 (MQ=255)

    ggCGATGCCGGAAATCCCGCTGCCGGTGTTCTGCCTGATGGTTGGTGCGGCAATTGGTCTGGGGAGCATTCTTACACCATACGCCGCCGGccc                                                                                       >  2:5970510/1‑93 (MQ=255)

    ggCGATGCCGGAAATCCCGCTGCCGGTGTTCTGCCTGATGGTTGGTGCGGCAATTGGTCTGGGGAGCATTCTTACACCATACGCCGCCGGccc                                                                                       <  2:1268801/93‑1 (MQ=255)

      cGATGCCGGAAATCCCGCTGCCGGTGTTCTGCCTGATGGTTGGTGCGGCAATTGGTCTGGGGAGCATTCTTACACCATACGCCGCCGGCCCCa                                                                                     <  1:3598688/93‑1 (MQ=255)

       gATGCCGGAAATCCCGCTGCCGGTGTTCTGCCTGATGGTTGGTGCGGCAATTGGTCTGGGGAGCATTCTTACACCATACGCCGCCGGCCCCAg                                                                                    <  1:1477346/93‑1 (MQ=255)

       gATGCCGGAAATCCCGCTGCCGGTGTTCTGCCTGATGGTTGGTGCGGCAATTGGTCTGGGGAGCATTCTTACACCATACGCCACCGGCCCCAg                                                                                    >  1:2841333/1‑93 (MQ=255)

         tGCCGGAAATCCCGCTGCCGGTGTTCTGCCTGATGGTTGGTGCGGCAATTGGTCTGGGGAGCATTCTTACACCATACGCCGCCGGCCCCAGTc                                                                                  <  1:4741864/93‑1 (MQ=255)

         tGCCGGAAATCCCGCCGCCGGTGTTCTGCCTGATGGTTGGTGCGGCAATTGGTCTGGGGAGCATTCTTACACCATACGCCGCCGGCCCCAGTc                                                                                  >  1:1036405/1‑93 (MQ=255)

          gCCGGAAACCCCGCTGCCGGTTTTCTGCCTGATGGTTGGTGCGGCAATTGGTCTGGGGAGCATTCTTCCACCATACGCCGCCGGCCCCAGTcc                                                                                 <  1:6165559/93‑1 (MQ=255)

           ccGGAAATCCCGCTGCCGGTGTTCTGCCTGATGGTTGGTGCGGCAATTGGTCTGGGGAGCATTCTTACACCATACGCCGCCGGCCCCAGTCCg                                                                                >  2:2007692/1‑93 (MQ=255)

            cGGAAATCCCGCTGCCGGTGTTCTGCCTGATGGTTGGTGCGGCAATTGGTCTGGGGAGCATTCTTACACCATACGCCGCCGGCCCCAGTCCGa                                                                               >  1:6555855/1‑93 (MQ=255)

                  tCCCGCTGCCGGTGTTCTGCCTGATGGTTGGTGCGGCAATTGGTCTGGGGAGCATTCTTACACCATACGCCGCCGGCCCCAGTCCGATTtact                                                                         >  1:6925267/1‑93 (MQ=255)

                  tCCCGCTGCCGGTGTTCTGCCTGATGGTTGGTGCGGCAATTGGTCTGGGGAGCATTCTTACACCATACGCCGCCGGCCCCAGTCCGATTtact                                                                         >  1:5623754/1‑93 (MQ=255)

                   cccGCTGCCGGTGTTCTGCCTGATGGTTGGTGCGGCAATTGGTCTGGGGAGCATTCTTACCCCATACGCCACCGGCCCCAGTCCGATTtacta                                                                        <  1:1842585/93‑1 (MQ=255)

                   cccGCTGCCGGTGTTCTGCCTGATGGTTGGTGCGGCAATTGGTCTGGGGAGCATTCTTACACCATACGCCGCCGGCCCCAGTCCGATTtacta                                                                        <  2:389689/93‑1 (MQ=255)

                    ccGCTGCCGGTGTTCTGCCTGATGGTTGGTGCGGCAATTGGTCTGGGGAGCATTCTTACACCATACGCCGCCGGCCCCAGTCCGATTtactac                                                                       >  1:4917060/1‑93 (MQ=255)

                    ccGCTGCCGGTGTTCTGCCTGATGGTTGGTGCGGCAATTGGTCTGGGGAGCATTCTTACACCATACGCCACCGGCCCCAGTCCGATTtactac                                                                       <  2:6591419/93‑1 (MQ=255)

                      gCTGCCGGTGTTCTGCCTGATGGTTGGTGCGGCAATTGGTCTGGGGAGCATTCTTACACCATACGCCGCCGGCCCCAGTCCGATTTACTACgg                                                                     >  1:5245392/1‑93 (MQ=255)

                       cTGCCGGTGTTCTGCCTGATGGTTGGTGCGGCAATTGGTCTGGGGAGCATTCTTACACCATACGCCGCCGGCCCCAGTCCGATTTACTACGGc                                                                    <  2:4409270/93‑1 (MQ=255)

                        tGCCGGTGTTCTGCCTGATGGTTGGTGCGGCAATTGGTCTGGGGAGCATTCTTACACCATACGCCACCGGCCCCAGTCCGATTTACTACGGCa                                                                   >  2:2613220/1‑93 (MQ=255)

                         gCCGGTGTTCTGCCTGATGGTTGGTGCGGCAATTGGTCTGGGGAGCATTCTTACCCCATACGCCGCCGGCCCCAGTCCGATTTACTACGGCAg                                                                  <  1:2728477/93‑1 (MQ=255)

                         gCCGGTGTTCTGCCTGATGGTTGGTGCGGCAATTGGTCTGGGGAGCATTCTTACACCATACGCCGCCGGCCCCAGTCCGATTTACTACGGCAg                                                                  >  1:2093771/1‑93 (MQ=255)

                         gCCGGTGTTCTGCCTGATGGTTGGTGCGGCAATTGGTCTGGGGAGCATTCTTACACCATACGCCGCCGGCCCCAGTCCGATTTACTACGGCAg                                                                  <  1:1641292/93‑1 (MQ=255)

                          ccGGTGTTCTGCCTGATGGTTGGTGCGGCAATTGGTCTGGGGAGCATTCTTACACCATACGCCGCCGGCCCCAGTCCGATTTACTACGGCAGt                                                                 <  2:2229722/93‑1 (MQ=255)

                          ccGGTGTTCTGCCTGATGGTTGGTGCGGCAATTGGTCTGGGGAGCATTCTTACACCATACGCCGCCGGCCCCAGTCCGATTTACTACGGCAGt                                                                 <  2:5065946/93‑1 (MQ=255)

                           cGGTGTTCTGCCTGATGGTTGGTGCGGCAATTGGTCTGGGGAGCATTCTTACACCATACGCCGCCGGCCCCAGTCCGATTTACTACGGCAGTg                                                                <  2:129571/93‑1 (MQ=255)

                           cGGTGTTCTGCCTGATGGTTGGTGCGGCAATTGGTCTGGGGAGCATTCTTACACCATACGCCACCGGCCCCAGTCCGATTTACTACGGCAGTg                                                                >  2:6284595/1‑93 (MQ=255)

                            ggTGTTCTGCCTGATGGTTGGTGCGGCAATTGGTCTGGGGAGCATTCTTACACCATACGCCGCCGGCCCCAGTCCGATTTACTACGGCAGTgg                                                               >  1:1686404/1‑93 (MQ=255)

                            ggTGTTCTGCCTGATGGTTGGTGCGGCAATTGGTCTGGGGAGCATTCTTACACCATACGCCGCCGGCCCCAGTCCGATTTACTACGGCAGTgg                                                               >  1:1384119/1‑93 (MQ=255)

                            ggTGTTCTGCCTGATGGTTGGTGCGGCAATTGGTCTGGGGAGCATTCTTACACCATACGCCGCCGGCCCCAGTCCGATTTACTACGGCAGTgg                                                               >  1:1086713/1‑93 (MQ=255)

                             gtgtTCTGCCTGATGGTTGGTGGGGGAATTGGTCTGGGGAGCATTCTTACACCATACGCCGCCGGCCCCAGTCCGATTTACTACGGCAGTGGt                                                              <  2:3397547/93‑1 (MQ=255)

                              tgtTCTGCCTGATGGTTGGTGCGGCAATTGGTCTGGGGAGCATTCATACACCATACGCCACCGGCCCCAGTCCGATTTACTACGGCAGTGGtt                                                             <  2:5437821/93‑1 (MQ=255)

                              gtttCTGCCTGATGGTTGGTGCGGCAATTGGTCTGGGGAGCATCCTTCCACCATCCGCCGCCGGCCCCAGTCCGA‑TTACTACGGCAGTGGTTa                                                            <  2:803593/91‑1 (MQ=255)

                               gtTCTGCCTGATGGTTGGTGCGGCAATTGGTCTGGGGAGCATTCTTACACCATACGCCGCCGGCCCCAGTCCGATTTACTACGGCAGTGGTTa                                                            <  1:517992/93‑1 (MQ=255)

                                ttCTGCCTGATGGTTGGTGCGGCAATTGGTCTGGGGAGCATTCTTACACCATACGCCGCCGGCCCCAGTCCGATTTACTACGGCAGTGGTTAt                                                           <  1:6221280/93‑1 (MQ=255)

                                ttCTGCCTGATGGTTGGTGCGGCAATTGGTCTGGGGAGCATTCTTACACCATACGCCGCCGGCCCCAGTCCGATTTACTACGGCAGTGGTTAt                                                           <  2:5365512/93‑1 (MQ=255)

                                 tCTGCCTGATGGTTGGTGCGGCAATTGGTCTGGGGAGCATTCTTACACCATACGCCGCCGGCCCCAGTCCGATTTACTACGGCAGTGGTTATc                                                          <  1:2204802/93‑1 (MQ=255)

                                  cTGCCTGATGGTTGGTGCGGCAATTGGTCTGGGGAGCATTCTTACACCATACGCCGCCGGCCCCAGTCCGATTTACTACGGCAGTGGTTATCt                                                         >  1:1924272/1‑93 (MQ=255)

                                  cTGCCTGATGGTTGGTGCGGCAATTGGTCTGGGGAGCATTCTTACACCATACGCCGCCGGCCCCAGTCCGATTTACTACGGCAGTGGTTATCt                                                         <  2:4359224/93‑1 (MQ=255)

                                  cTGCCTGATGGTTGGTGCGGCAATTGGTCTGGGGAGCATTCTTACACCATACGCCGCCGGCCCCAGTCCGATTTACTACGGCAGTGGTTATCt                                                         <  1:975221/93‑1 (MQ=255)

                                   tGCCTGATGGTTGGTGCGGCAATTGGTCTGGGGAGCATTCTTACACCATACGCCGCCGGCCCCATTCCGATTTACTACGGCAGTGGTTATCTg                                                        <  2:7345367/93‑1 (MQ=255)

                                   tGCCTGATGGTTGGTGCGGCAATTGGTCTGGGGAGCATTCTTACACCATACGCCACCGGCCCCAGTCCGATTTACTACGGCAGTGGTTATCTg                                                        <  2:1876478/93‑1 (MQ=255)

                                    gCCTGATGGTTGGTGCGGCAATTGGTCTGGGGAGCATTCTTACACCATACGCCGCCGGCCCCAGTCCGATTTACTACGGCAGTGGTTATCTGc                                                       <  1:4666155/93‑1 (MQ=255)

                                     ccTGATGGTTGGTGCGGCAATTGGTCTGGGGAGCATTCTTACACCATACGCCGCCGGCCCCAGTCCGATTTACTACGGCAGTGGTTATCTGcc                                                      <  1:3797567/93‑1 (MQ=255)

                                     ccTGATGGTTGGTGCGGCAATTGGTCTGGGGAGCATTCTTACACCATACGCCACCGGCCCCAGTCCGATTTACTACGGCAGTGGTTATCTGcc                                                      >  2:3467226/1‑93 (MQ=255)

                                     ccTGATGGTTGGTGCGGCAATTGGTCTGGGGAGCATTCTTACACCATACGCCACCGGCCCCAGTCCGATTTACTACGGCAGTGGTTATCTGcc                                                      <  1:3383842/93‑1 (MQ=255)

                                     ccTGATGGTTGGTGCGGCAATTGGTCTGGGGAGCATTCTTACACCATACGCCACCGGCCCCAGTCCGATTTACTACGGCAGTGGTTATCTGcc                                                      <  1:3166614/93‑1 (MQ=255)

                                       tGATGGTTGGTGCGGCAATTGGTCTGGGGAGCATTCTTACACCATACGCCGCCGGCCCCAGTCCGATTTACTACGGCAGTGGTTATCTGCCaa                                                    >  2:770607/1‑93 (MQ=255)

                                       tGATGGTTGGTGCGGCAATTGGTCTGGGGAGCATTCTTACACCATACGCCGCCGGCCCCAGTCCGATTTACTACGGCAGTGGTTATCTGCCaa                                                    >  2:4755514/1‑93 (MQ=255)

                                        gATGGTTGGTGCGGCAATTGGTCTGGGGAGCATTCTTACACCATACGCCGCCGGCCCCAGTCCGATTTACTACGGCAGTGGTTATCTGCCAAc                                                   <  2:6523204/93‑1 (MQ=255)

                                        gATGGTTGGTGCGGCAATTGGTCTGGGGAGCATTCTTACACCATACGCCACCGGCCCCAGTCCGATTTACTACGGCAGTGGTTATCTGCCAAc                                                   <  2:5016097/93‑1 (MQ=255)

                                          tggttggtGCGGCAATTGGTCTGGGGAGCATTCTTACACCATACGCCGCCGGCCCCAGTCCGATTTACTACGGCAGTGGTTATCTGCCAAcgg                                                 <  2:7053286/93‑1 (MQ=255)

                                          tggttggtGCGGCAATTGGTCTGGGGAGCATTCTTACACCATACGCCACCGGCCCCAGTCCGATTTACTACGGCAGTGGTTATCTGCCAAcgg                                                 >  1:7429584/1‑93 (MQ=255)

                                           ggttggtGCGGCAATTGGTCTGGGGAGCATTCTTACACCATACGCCGCCGGCCCCAGTCCGATTTACTACGGCAGTGGTTATCTGCCAAcggc                                                <  1:2038247/93‑1 (MQ=255)

                                           ggttggtGCGGCAATTGGTCTGGGGAGCATTCTTACACCATACGCCGCCGGCCCCAGTCCGATTTACTACGGCAGTGGTTATCTGCCAAcggc                                                <  1:1918082/93‑1 (MQ=255)

                                             ttggtGCGGCAATTGGTCTGGGGAGCATTCTTACACCATACGCCGCCGGCCCCAGTCCGATTTACTACGGCAGTGGTTATCTGCCAAcggcgg                                              >  1:4287251/1‑93 (MQ=255)

                                              tggtGCGGCAATTGGTCTGGGGAGCATTCTTACACCATACGCCGCCGGCCCCAGTCCGATTTACTACGGCAGTGGTTATCTGCCAACGGCGGa                                             <  1:4220338/93‑1 (MQ=255)

                                               ggtGCGGCAATTGGTCTGGGGAGCATTCTTACACCATACGCCGCCGGCCCCAGTCCGATTTACTACGGCAGTGGTTATCTGCCAACGGCGGAt                                            >  2:3224543/1‑93 (MQ=255)

                                                     gCAATTGGTCTGGGGAGCATTCTTACACCATACGCCACCGGCCCCAGTCCGATTTACTACGGCAGTGGTTTTCTGCCAACGGCGGATTACTgg                                      <  2:5134831/93‑1 (MQ=255)

                                                     gCAATTGGTCTGGGGAGCATTCTTACACCATACGCCACCGGCCCCAGTCCGATTTACTACGGCAGTGGTTATCTGCCAACGGCGGATTACTgg                                      <  1:5403251/93‑1 (MQ=255)

                                                       aaTTGGTCTGGGGAGCATTCTTACACCATACGCCGCCGGCCCCAGTCCGATTTACTACGGCAGTGGTTATCTGCCAACGGCGGATTACTGGCg                                    <  2:3973403/93‑1 (MQ=255)

                                                        aTTGGTCTGGGGAGCATTCTTACACCATACGCCACCGGCCCCAGTCCGATTTACTACGGCAGTGGTTATCTGCCAACGGCGGATTACTGGCGa                                   <  2:4880445/93‑1 (MQ=255)

                                                         ttGGTCTGGGGAGCATTCTTACACCATACGCCGCCGGCCCCAGTCCGATTTACTACGGCAGTGGTTATCTGCCAACGGCGGATTACTGGCGAc                                  <  2:2268943/93‑1 (MQ=255)

                                                         ttGGTCTGGGGAGCATTCTTACACCATACGCCGCCGGCCCCAGTCCGATTTACTACGGCAGTGGTTATCTGCCAACGGCGGATTACTGGCGAc                                  <  2:115481/93‑1 (MQ=255)

                                                         ttGGTCTGGGGAGCATTCTTACACCATACGCCGCCGGCCCCAGTCCGATTTACTACGGCAGTGGTTATCTGCCAACGGCGGATTACTGGCGAc                                  <  2:6158359/93‑1 (MQ=255)

                                                          tGGTCTGGGGAGCATTCTTACACCATACGCCGCCGGCCCCAGTCCGATTTACTACGGCAGTGGTTATCTGCCAACGGCGGATTACTGGCGACt                                 <  1:6047828/93‑1 (MQ=255)

                                                           ggTCTGGGGAGCATTCTTACACCATACGCCGCCGGCCCCAGTCCGATTTACTACGGCAGTGGTTATCTGCCAACGGCGGATTACTGGCGACTg                                >  2:4557752/1‑93 (MQ=255)

                                                            gTCTGGGGAGCATTCTTACACCATACGCCGCCGGCCCCAGTCCGATTTACTACGGCAGTGGTTATCTGCCAACGGCGGATTACTGGCGACTgg                               <  1:7365037/93‑1 (MQ=255)

                                                            gTCTGGGGAGCATTCTTACACCATACGCCGCCGGCCCCAGTCCGATTTACTACGGCAGTGGTTATCTGCCAACGGCGGATTACTGGCGACTgg                               <  1:5137890/93‑1 (MQ=255)

                                                             tCTGGGGAGCATTCTTACACCATACGCCGCCGGCCCCAGTCCGATTTACTACGGCAGTGGTTATCTGCCAACGGCGGATTACTGGCGACTggg                              <  2:2069181/93‑1 (MQ=255)

                                                               tGGGGAGCATTCTTACACCATACGCCGCCGGCCCCAGTCCGATTTACTACGGCAGTGGTTATCTGCCAACGGCGGATTACTGGCGACTggggg                            <  1:3477521/93‑1 (MQ=255)

                                                                ggggAGCATTCTTACACCATACGCCGCCGGCCCCAGTCCGATTTACTACGGCAGTGGTTATCTGCCAACGGCGGATTACTGGCGACTGGGGGc                           <  1:2910408/93‑1 (MQ=255)

                                                                ggggAGCATTCTTACACCATACGCCGCCGGCCCCAGTCCGATTTACTACGGCAGTGGTTATCTGCCAACGGCGGATTACTGGCGACTGGGGGc                           <  1:5202/93‑1 (MQ=255)

                                                                ggggAGCATTCTTACACCATACGCCGCCGGCCCCAGTCCGATTTACTACGGCAGTGGTTATCTGCCAACGGCGGATTACTGGCGACTGGGGGc                           <  1:502836/93‑1 (MQ=255)

                                                                ggggAGCATTCTTACACCATACGCCACCGGCCCCAGTCCGATTTACTACGGCAGTGGTTATCTGCCAACGGCGGATTACTGGCGACTGGGGGc                           <  2:2455612/93‑1 (MQ=255)

                                                                ggggAGCATTCTTACACCATACGCCACCGGCCCCAGTCCGATTTACTACGGCAGTGGTTATCTGCCAACGGCGGATTACTGGCGACTGGGGGc                           <  2:6702348/93‑1 (MQ=255)

                                                                  ggAGCATTCTTACACCATACGCCGCCGGCCCCAGTCCGATTTACTACGGCAGTGGTTATCTGCCAACGGCGGATTACTGGCGACTGGGGGCGa                         >  2:6466237/1‑93 (MQ=255)

                                                                  ggAGCATTCTTACACCATACGCCACCGGCCCCAGTCCGATTTACTACGGCAGTGGTTATCTGCCAACGGCGGATTACTGGCGACTGGGGGCGa                         >  2:3649753/1‑93 (MQ=255)

                                                                    aggaTTCTTACCCCATACGCCGCCGGGCCCAGTTCGATTTTCTTCGGGAGTGGTTTTCTGCCAACGGGGGATTACTGGCGACTGGGGGCGAtt                       <  1:2327336/90‑1 (MQ=255)

                                                                     gCATTCTTACACCATACGCCGCCGGCCCCAGTCCGATTTACTACGGCAGTGGTTATCTGCCAACGGCGGATTACTGGCGACTGGGGGCGAttt                      >  2:1932962/1‑93 (MQ=255)

                                                                     gCATTCTTACACCATACGCCGCCGGCCCCAGTCCGATTTACTACGGCAGTGGTTATCTGCCAACGGCGGATTACTGGCGACTGGGGGCGAttt                      >  2:1668223/1‑93 (MQ=255)

                                                                     gCATTCTTACACCATACGCCGCCGGCCCCAGTCCGATTTACTACGGCAGTGGTTATCTGCCAACGGCGGATTACTGGCGACTGGGGGCGAttt                      <  2:3819122/93‑1 (MQ=255)

                                                                        ttCTTACACCATACGCCGCCGGCCCCAGTCCGATTTACTACGGCAGTGGTTATCTGCCAACGGCGGATTACTGGCGACTGGGGGCGATTTTTg                   >  1:4349194/1‑93 (MQ=255)

                                                                              cacCATACGCCGCCGGCCCCAGTCCGATTTACTACGGCAGTGGTTATCTGCCAACGGCGGATTACTGGCGACTGGGGGCGATTTTTGGGCTGa             >  1:6828925/1‑93 (MQ=255)

                                                                              cacCATACGCCGCCGGCCCCAGTCCGATTTACTACGGCAGTGGTTATCTGCCAACGGCGGATTACTGGCGACTAGGGGCGATTTTTggggtca             >  1:2958704/1‑89 (MQ=255)

                                                                               acCATACGCCGCCGGCCCCAGTCCGATTTACTACGGCAGTGGTTATCTGCCAACGGCGGATTACTGGCGACTGGGGGCGATTTTTGGGCTGat            >  1:4207813/1‑93 (MQ=255)

                                                                                ccATACGCCGCCGGCCCCAGTCCGATTTACTACGGCAGTGGTTATCTGCCAACGGCGGATTACTGGCGACTGGGGGCGATTTTTGGGCTGata           <  1:7102729/93‑1 (MQ=255)

                                                                                ccATACGCCGCCGGCCCCAGTCCGATTTACTACGGCAGTGGTTATCTGCCAACGGCGGATTACTGGCGACTGGGGGCGATTTTTGGGCTGata           >  2:42345/1‑93 (MQ=255)

                                                                                 cATACGCCGCCGGCCCCAGTCCGATTTACTACGGCAGTGGTTATCTGCCAACGGCGGATTACTGGCGACTGGGGGCGATTTTTGGGCTGatat          <  2:6091057/93‑1 (MQ=255)

                                                                                  aTACGCCACCGGCCCCAGTCCGATTTACTACGGCAGTGGTTATCTGCCAACGGCGGATTACTGGCGACTGGGGGCGATTTTTGGGCTGATAtt         >  2:1683604/1‑93 (MQ=255)

                                                                                   tACGCCGCCGGCCCCAGTCCGATTTACTACGGCAGTGGTTATCTGCCAACGGCGGATTACTGGCGACTGGGGGCGATTTTTGGGCTGATATTc        <  2:4077615/93‑1 (MQ=255)

                                                                                   tACGCCGCCGGCCCCAGTCCGATTTACTACGGCAGTGGTTATCTGCCAACGGCGGATTACTGGCGACTGGGGGCGATTTTTGGGCTGATATTc        <  2:1833965/93‑1 (MQ=255)

                                                                                    aCGCCGCCGGCCCCAGTCCGATTTACTACGGCAGTGGTTATCTGCCAACGGCGGATTACTGGCGACTGGGGGCGATTTTTGGGCTGATATTcc       <  2:6793428/93‑1 (MQ=255)

                                                                                     cgccgcCGGCCCCAGTCCGATTTACTACGGCAGTGGTTATCTGCCAACGGCGGATTACTGGCGACTGGGGGCGATTTTTGGGCTGATATTCCt      <  1:6482556/93‑1 (MQ=255)

                                                                                      gccgccGGCCCCAGTCCGATTTACTACGGCAGTGGTTATCTGCCAACGGCGGATTACTGGCGACTGGGGGCGATTTTTGGGCTGATATTCCTc     >  1:5300091/1‑93 (MQ=255)

                                                                                       ccACCGGCCCCAGTCCGATTTACTACGGCAGTGGTTATCTGCCAACGGCGGATTACTGGCGACTGGGGGCGATTTTTGGGCTGATATTCCTCg    >  2:1217366/1‑93 (MQ=255)

CGCTGGCGATGCCGGAAATCCCGCTGCCGGTGTTCTGCCTGATGGTTGGTGCGGCAATTGGTCTGGGGAGCATTCTTACACCATACGCCACCGGCCCCAGTCCGATTTACTACGGCAGTGGTTATCTGCCAACGGCGGATTACTGGCGACTGGGGGCGATTTTTGGGCTGATATTCCTCGTA  >  NC_000913/3206899‑3207080

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