BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,696,865	T→C	P53P (CCT→CCC)	yfhH →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,696,865	0	T	C	99.7%	1039.4 / ‑5.3	328	P53P (CCT→CCC)	yfhH	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (159/167); total (161/167)
Fisher's exact test for biased strand distribution p-value = 4.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.33e-01

Only 100 of 346 total aligned reads displayed.

AGATATGAACGGCTTACTTCGTATCCGTCAGCGTTACCAGGGGCTTGCCCAGAGTGATAAAAAACTGGCGGATTATCTGCTGCTACAACCTGATACGGCGCGCCATTTAAGCTCTCAGCAACTGGCCAACGAAGCCGGAGTCAGTCAGTCCAGCGTCGTGAAGTTCGCGCAAAAACTCGG  >  NC_000913/2696775‑2696954

aGATATGAACGGCTTACTTCGTATCCGTCAGCGTTACCAGGGACTTGCCCAGAGTGATAAAAAACTGGCGGATTATCTGCTGCTACAACCCGa                                                                                         <  1:6474161/93‑1 (MQ=255)

 gATATGAACGGCTTACTTCGTATCCGTCAGCGTTACCAGGGGCTTGCCCAGAGTGATAAAAAACTGGCGGATTATCTGCTGCTACAACCCGAt                                                                                        <  1:8475537/93‑1 (MQ=255)

 aatatGAACGGCTTACTTCGTNTCCGTCAGCGTTACCAGGGGCTTGCCCAGAGTGATAAAAAACTGGCGGATTATCTGCTGCTACAACCCGAt                                                                                        <  1:7643081/92‑1 (MQ=255)

   tataAACGGCTTACTTCGTATCCGTCAGCGTTACCAGGGGCTTGCCCAGAGTGATAAAAAACTGGCGGATTATCTGCTGCTACAACCCGATAc                                                                                      >  2:664962/1‑93 (MQ=255)

   tatGAACGGCTTACTTCTTATCCGTCAGCGTTACCAGGGGCTTGCCCAGAGTGATAAAAAACTGGCGGATTATCTGCTGCTACAACCCGATAc                                                                                      <  2:9496972/93‑1 (MQ=255)

   tatGAACGGCTTACTTCGTATCCGTCAGCGTTACCAGGGGCTTGCCCAGAGTGATAAAAAACTGGCGGATTATCTGCTGCTACAACCCGATAc                                                                                      <  1:912970/93‑1 (MQ=255)

     tGAACGGCTTACTTCGTATCCGTCAGCGTTACCAGGGGCTTGCCCAGAGTGATAAAAAACTGGCGGATTATCTGCTGCTACAACCCGATACgg                                                                                    <  2:8474844/93‑1 (MQ=255)

     tGAACGGC‑TACTTCGTATCCGTCAGCGTTACCAGGGGCTTGCCCAGAGTGATAAAAAACTGGCGGATTATCTGCTGCTACAACCCGATACGgc                                                                                   <  1:3293008/93‑1 (MQ=255)

       aaCGGCTTACTTCGTATCCGTCAGCGTTACCAGGGGCTTGCCCAGAGTGATAAAAAACTGGCGGATTATCTGCTGCTACAACCCGATACGgcg                                                                                  >  1:2602708/1‑93 (MQ=255)

       aaCGGCTTACTTCGTATCCGTCAGCGTTACCAGGGACTTGCCCAGAGTGATAAAAAACTGGCGGATTATCTGCTGCTACAACCCGATACGgcg                                                                                  >  2:7729024/1‑93 (MQ=255)

        aCGGCTTACTTCGTATCCGTCAGCGTTACCAGGGGCTTGCCCAGAGTGATAAAAAACTGGCGGATTATCTGCTGCTACAACCCGATACGgcgc                                                                                 <  1:9745486/93‑1 (MQ=255)

         cGGCTTACTTCGTATCCGTCAGCGTTACCAGGGGCTTGCCCAGAGTGATAAAAAACTGGCGGATTATCTGCTGCTACAACCCGATACGgcgcg                                                                                >  2:7692285/1‑93 (MQ=255)

          ggCTTACTTCGTATCCGTCAGCGTTACCAGGGGCTTGCCCAGAGTGATAAAAAACTGGCGGATTATCTGCTGCTACAACCCGATACGgcgcgc                                                                               <  1:3968973/93‑1 (MQ=255)

          ggCTTACTTCGTATCCGTCAGCGTTACCAGGGGCTTGCCCAGAGTGATAAAAAACTGGCGGATTATCTGCTGCTACAACCCGATACGgcgcgc                                                                               >  1:5702612/1‑93 (MQ=255)

            cTTACTTCGTATCCGTCAGCGTTACCAGGGGCTTGCCCAGAGTGATAAAAAACTGGCGGATTATCTGCTGCTACAACCCGATACGGCGCGCCa                                                                             <  1:1031291/93‑1 (MQ=255)

             ttACTTCGTATCCGTCAGCGTTACCAGGGGCTTGCCCAGAGTGATAAAAAACTGGCGGATTATCTGCTGCTACAACCCGATACGGCGCGCCAt                                                                            >  1:6355533/1‑93 (MQ=255)

              tACTTCGTATCCGTCAGCGTTACCAGGGGCTTGCCCAGAGTGATAAAAAACTGGCGGATTATCTGCTGCTACAACCCGATACGGCGCGCCAtt                                                                           >  1:814375/1‑93 (MQ=255)

               aCTTCGTATCCGTCAGCGTTACCAGGGGCTTGCCCAGAGTGATAAAAAACTGGCGGATTATCTGCTGCTACAACCCGATACGGCGCGCCAttt                                                                          <  2:7435329/93‑1 (MQ=255)

               aCTTCGTATCCGTCAGCGTTACCAGGGGCTTGCCCAGAGTGATAAAAAACTGGCGGATTATCTGCTGCTACAACCCGATACGGCGCGCCAttt                                                                          >  2:118688/1‑93 (MQ=255)

               aCTTCGTATCCGTCAGCGTTACCAGGGGCTTGCCCAGAGTGATAAAAAACTGGCGGATTATCTGCTGCTACAACCCGATACGGCGCGCCAttt                                                                          <  2:6174459/93‑1 (MQ=255)

                cTTCGTATCCGTCAGCGTTACCAGGGGCTTGCCCAGAGTGATAAAAAACTGGCGGATTATCTGCTGCTACAACCCGATACGGCGCGCCATTTa                                                                         <  1:8299584/93‑1 (MQ=255)

                cTTCGTATCCGTCAGCGTTACCAGGGGCTTGCCCAGAGTGATAAAAAACTGGCGGATTATCTGCTGCTACAACCCGATACGGCGCGCCATTTa                                                                         <  2:1621506/93‑1 (MQ=255)

                cTTCGTATCCGTCAGCGTTACCAGGGGCTTGCCCAGAGTGATAAAAAACTGGCGGATTATCTGCTGCTACAACCCGATACGGCGCGCCATTTa                                                                         <  1:3204244/93‑1 (MQ=255)

                 ttCGTATCCGTCAGCGTTACCAGGGGCTTGCCCCGAGTGATAAAAAACTGGCGGATTATCTGCTGCTACAACCCGATACGGCGCGCCATTTaa                                                                        >  2:8075045/1‑93 (MQ=255)

                 ttCGTATCCGTCAGCGTTACCAGGGGCTTGCCCAGAGTGATAAAAAACTGGCGGATTATCTGCTGCTACAACCCGATACGGCGCGCCAtttga                                                                        >  1:7296551/1‑91 (MQ=255)

                 ttCGTATCCGTCAGCGTTACCAGGGGCTTGCCCAGAGTGATAAAAAACTGGCGGATTATCTGCTGCTACAACCCGATACGGCGCGCCATTTaa                                                                        >  2:4084633/1‑93 (MQ=255)

                 ttCGTATCCGTCAGCGTTACCAGGGGCTTGCCCAGAGTGATAAAAAACTGGCGGATTATCTGCTGCTACAACCCGATACGGCGCGCCATTTaa                                                                        <  1:1694703/93‑1 (MQ=255)

                   cGTATCCGTCAGCGTTACCAGGGGCTTGCCCAGAGTGATAAAAAACTGGCGGATTATCTGCTGCTACAACCCGATACGGCGCGCCATTTAAGc                                                                      <  2:7980040/93‑1 (MQ=255)

                   cGTATCCGTCAGCGTTACCAGGGGCTTGCCCAGAGTGATAAAAAACTGGCGGATTATCTGCTGCTACAACCCGATACGGCGCGCCATTTAAGc                                                                      <  2:4678354/93‑1 (MQ=255)

                    gTATCCGTCAGCGTTACCAGGGGCTTGCCCAGAGTGATAAAAAACTGGCGGATTATCTGCTGCTACAACCCGATACGGCGCGCCATTTAAGct                                                                     >  1:4372582/1‑93 (MQ=255)

                     tATCCGTCAGCGTTACCAGGGGCTTGCCCAGAGTGATAAAAAACTGGCGGATTATCTGCTGCTACAACCCGATACGGCGCGCCATTTAAGctc                                                                    <  1:7069045/93‑1 (MQ=255)

                       tCCGTCAGCGTTACCAGGGGCTTGCCCAGAGTGATAAAAAACTGGCGGATTATCTGCTGCTACAACCCGATACGGCGCGCCATTTAAGCtctc                                                                  >  2:3598128/1‑93 (MQ=255)

                       tCCGTCAGCGTTACCAGGGGCTTGCCCAGAGTGATAAAAAACTGGCGGATTATCTGCTGCTACAACCCGATACGGCGCGCCATTTAAGCtctc                                                                  >  1:9307001/1‑93 (MQ=255)

                        ccGTCAGCGTTACCAGGGGCTTGCCCAGAGTGATAAAAAACTGGCGGATTATCTGCTGCTACAACCCGATACGGCGCGCCATTTAAGCTCTCa                                                                 <  2:3104048/93‑1 (MQ=255)

                            cAGCGTTACCAGGGGCTTGCCCAGAGTGATAAAAAACTGGCGGATTATCTGCTGCTACAACCCGATACGGCGCGCCATTTAAGCTCTCAGCaa                                                             <  2:2813723/93‑1 (MQ=255)

                             aGCGTTACCAGGGGCTTGCCCAGAGTGATAAAAAACTGGCGGATTATCTGCTGCTACAACCCGATACGGCGCGCCATTTAAGCTCTCAGCAAc                                                            <  1:631964/93‑1 (MQ=255)

                             aGCGTTACCAGGGGCTTGCCCAGAGTGATAAAAAACTGGCGGATTATCTGCTGCTACAACCCGATACGGCGCGCCATTTAAGCTCTCAGCAAc                                                            <  2:9753220/93‑1 (MQ=255)

                              gCGTTACCAGGGGCTTGCCCAGAGTGATAAAAAACTGGCGGATTATCTGCTGCTACAACCCGATACGGCGCGCCATTTAAGCTCTCAGCAACt                                                           <  2:6556649/93‑1 (MQ=255)

                              gCGTTACCAGGGGCTTGCCCAGAGTGATAAAAAACTGGCGGATTATCTGCTGCTACAACCCGATACGGCGCGCCATTTAAGCTCTCAGCAACt                                                           <  1:2504691/93‑1 (MQ=255)

                              gCGTTACCAGGGGCTTGCCCAGAGTGATAAAAAACTGGCGGATTATCTGCTGCTACAACCCGATACGGCGCGCCATTTAAGCTCTCAGCAACt                                                           <  1:9578728/93‑1 (MQ=255)

                               cGTTACCAGGGGCTTGCCCAGAGTGATAAAAAACTGGCGGATTATCTGCTGCTACAACCCGATACGGCGCGCCATTTAAGCTCTCAGCAACTg                                                          <  2:2506093/93‑1 (MQ=255)

                                gTTACCAGGGGCTTGCCCAGAGTGATAAAAAACTGGCGGATTATCTGCTGCTACAACCCGATACGGCGCGCCATTTAAGCTCTCAGCAACTgg                                                         >  2:276101/1‑93 (MQ=255)

                                 ttACCAGGGGCTTGCCCAGAGTGATAAAAAACTGGCGGATTATCTGCTGCTACAACCCGATACGGCGCGCCATTTAAGCTCTCAGCAACTGGc                                                        <  1:8411731/93‑1 (MQ=255)

                                 ttACCAGGGGCTTGCCCAGAGTGATAAAAAACTGGCGGATTATCTGCTGCTACAACCCGATACGGCGCGCCATTTAAGCTCTCAGCAACTGGc                                                        <  1:7895062/93‑1 (MQ=255)

                                 ttACCAGGGGCTTGCCCAGAGTGATAAAAAACTGGCGGATTATCTGCTGCTACAACCCGATACGGCGCGCCATTTAAGCTCTCAGCAACTGGc                                                        <  2:2548127/93‑1 (MQ=255)

                                  ttccAGGGGCTTGCCCAGAGTGATAAAAAACTGGCGGATTATCTGCTGCTACAACCCGATACGGCGCGCCATTTAAGCTCTCAGCAACTGGcc                                                       <  1:2976018/91‑1 (MQ=255)

                                    ccAGGGGCTTGCCCAGAGTGATAAAAAACTGGCGGATTATCTGCTGCTACAACCCGCTACGGCGCGCCATTTAAGCTCTCAGCAACTGGCCaa                                                     >  1:7546424/1‑93 (MQ=255)

                                    ccAGGGGCTTGCCCAGAGTGATAAAAAACTGGCGGATTATCTGCTGCTACAACCCGATACGGCGCGCCATTTAAGCTCTCAGCAACTGGCCaa                                                     >  2:635427/1‑93 (MQ=255)

                                    ccAGGGGCTTGCCCAGAGTGATAAAAAACTGGCGGATTATCTGCTGCTACAACCCGATACGGCGCGCCATTTAAGCTCTCAGCAACTGGCCaa                                                     <  2:7901703/93‑1 (MQ=255)

                                    ccAGGGGCTTGCCCAGAGTGATAAAAAACTAGCGGATTATCTGCTGCTACAACCCGATACGGCGCGCCATTTAAGCTCTCAGCAACTGGCCaa                                                     >  1:4154244/1‑93 (MQ=255)

                                     cAGGGGCTTGCCCAGAGTGATAAAAAACTGGCGGATTATCTGCTGCTACAACCCGATACGGCGCGCCATTTAAGCTCTCAGCAACTGGCCAAc                                                    >  2:7074135/1‑93 (MQ=255)

                                     cAGGGGCTTGCCCAGAGTGATAAAAAACTGGCGGATTATCTGCTGCTACAACCCGATACGGCGCGCCATTTAAGCTCTCAGCAACTGGCCAAc                                                    >  2:8731219/1‑93 (MQ=255)

                                      aGGGGCTTGCCCAGAGTGATAAAAAACTGGCGGATTATCTGCTGCTACAACCCGATACGGCGCGCCATTTAAGCTCTCAGCAACTGGCCAACg                                                   >  2:9065476/1‑93 (MQ=255)

                                       ggggCTTGCCCAGAGTGATAAAAAACTGGCGGATTATCTGCTGCTACAACCCGATACGGCGCGCCATTTAAGCTCTCAGCAACTGGCCAACGa                                                  >  2:4151768/1‑93 (MQ=255)

                                        gggCTTGCCCAGAGTGATAAAAAACTGGCGGATTATCTGCTGCTACAACCCGATACGGCGCGCCATTTAAGCTCTCAGCAACTGGCCAACGaa                                                 >  1:6187629/1‑93 (MQ=255)

                                         ggCTTGCCCAGAGTGATAAAAAACTGGCGGATTATCTGCTGCTACAACCCGATACGGCGCGCCATTTAAGCTCTCAGCAACTGGCCAACGAAg                                                >  1:9386926/1‑93 (MQ=255)

                                            ttGCCCAGAGTGATAAAAAACTGGCGGATTATCTGCTGCTACAACCCGATACGGCGCGCCATTTAAGCTCTCAGCAACTGGCCAACGAAGCCg                                             >  2:8541310/1‑93 (MQ=255)

                                             tGCCCAGAGTGATAAAAAACTGGCGGATTATCTGCTGCTACAACCCGATACGGCGCGCCATTTAAGCTCTCAGCAACTGGCCAACGAAGCCgg                                            >  2:4516903/1‑93 (MQ=255)

                                              gCCCAGAGTGATAAAAAACTGGCGGATTATCTGCTGCTACAACCCGATACGGCGCGCCATTTAAGCTCTCAGCAACTGGCCAACGAAGCCGGa                                           <  1:9002386/93‑1 (MQ=255)

                                              gCCCAGAGTGATAAAAAACTGGCGGATTATCTGCTGCTACAACCCGATACGGCGCGCCATTTAAGCTCTCAGCAACTGGCCAACGAAGCCGGa                                           <  1:1835887/93‑1 (MQ=255)

                                                 cAGAGTGATAAAAAACTGGCGGATTATCTGCTGCTACAACCCGATACGGCGCGCCATTTGAGCTCTCAGCAACTGGCCAACGAAGCCGGagtc                                        <  2:3431562/93‑1 (MQ=255)

                                                   gagTGATAAAAAACTGGCGGATTATCTGCTGCTACAACCCGATACGGCGCGCCATTTAAGCTCTCAGCAACTGGCCAACGAAGCCGGagtcag                                      <  1:1148040/93‑1 (MQ=255)

                                                     gTGATAAAAAACTGGCGGATTATCTGCTGCTACAACCCGATACGGCGCGCCATTTAAGCTCTCAGCAACTGGCCAACGAAGCCGGagtcagtc                                    <  1:5048330/93‑1 (MQ=255)

                                                      tGATAAAAAACTGGCGGATTATCTGCTGCTACAACCCGATACGGCGCGCCATTTAAGCTCTCAGCAACTGGCCAACGAAGCCGGagtcagtca                                   <  1:3764040/93‑1 (MQ=255)

                                                      tGATAAAAAACTGGCGGATTATCTGCTGCTACAACCCGATACGGCGCGCCATTTAAGCTCTCAGCAACTGGCCAACGAAGCCGGagtcagtca                                   >  2:5784695/1‑93 (MQ=255)

                                                      tGATAAAAAACTGGCGGATTATCTGCTGCTACAACCCGATACGGCGCGCCATTTAAGCTCTCAGCAACTGGCCAACGAAGCCGGagtcagtca                                   >  2:1363678/1‑93 (MQ=255)

                                                      tGATAAAAAACTGGCGGATTATCTGCTGCTACAACCCGATACGGCGCGCCATTTAAGCTCTCAGCAACTGGCCAACGAAGCCGGagtcagtca                                   <  1:2330840/93‑1 (MQ=255)

                                                      tGATAAAAAACTGGCGGATTATCTGCTGCTACAACCCGATACGGCGCGCCATTTAAGCTCTCAGCAACTGGCCAACGAAGCCGGagtcagtca                                   >  1:5326469/1‑93 (MQ=255)

                                                      tGATAAAAAACTGGCGGATTATCTGCTGCTACAACCCGATACGGCGCGCCATTTAAGCTCTCAGCAACTGGCCAACGAAGCCGGagtcagtca                                   <  1:4795176/93‑1 (MQ=255)

                                                       gATAAAAAACTGGCGGATTATCTGCTGCTACAACCCGATACGGCGCGCCATTTAAGCTCTCAGCAACTGGCCAACGAAGCCGGagtcagtcag                                  >  1:272169/1‑93 (MQ=255)

                                                       gATAAAAAACTGGCGGATTATCTGCTGCTACAACCCGATACGGCGCGCCATTTAAGCTCTCAGCAACTGGCCAACGAAGCCGGagtcagtcag                                  >  2:1290165/1‑93 (MQ=255)

                                                          aaaaaaCTGGCGGATTATCTGCTGCTACAACCCGATACGGCGCGCCATTTAAGCTCTCAGCAACTGGCCAACGAAGCCGGAGTCAGTCAGTcc                               <  2:7199574/93‑1 (MQ=255)

                                                          aaaaaaCTGGCGGATTATCTGCTGCTACAACCCGATACGGCGCGCCATTTAAGCTCTCAGCAACTGGCCAACGAAGCCGGAGTCAGTCAGTcc                               >  1:3412383/1‑93 (MQ=255)

                                                           aaaaaCTGGCGGATTATCTGCTGCTACAACCCGATACGGCGCGCCATTTAAGCTCTCAGCAACTGGCCAACGAAGCCGGAGTCAGTCAGTCCa                              >  2:4983911/1‑93 (MQ=255)

                                                           aaaaaCTGGCGGATTATCTGCTGCTACAACCCGATACGGCGCGCCATTTAAGCTCTCAGCAACTGGCCAACGAAGCCGGAGTCAGTCAGTCCa                              <  2:5373772/93‑1 (MQ=255)

                                                            aaaaCTGGCGGATTATCTGCTGCTACAACCCGATACGGCGCGCCATTTAAGCTCTCAGCAACTGGCCAACGAAGCCGGAGTCAGTCAGTCCAg                             >  2:5664646/1‑93 (MQ=255)

                                                            aaaaCTGGCGGATTATCTGCTGCTACAACCCGATACGGCGCGCCATTTAAGCTCTCAGCAACTGGCCAACGAAGCCGGAGTCAGTCAGTCCAg                             <  2:3285467/93‑1 (MQ=255)

                                                              aaCTGGCGGATTATCTGCTGCTACAACCCGATACGGCGCGCCATTTAAGCTCTCAGCAACTGGCCAACGAAGCCGGAGTCAGTCAGTCCAGcg                           >  1:2189801/1‑93 (MQ=255)

                                                               aCTGGCGGATTATCTGCTGCTACAACCCGATACGGCGCGCCATTTAAGCTCTCAGCAACTGGCCAACGAAGCCGGAGTCAGTCAGTCCAGcgt                          >  2:3076568/1‑93 (MQ=255)

                                                                cTGGCGGATTATCTGCTGCTACAACCCGATACGGCGCGCCATTTAAGCTCTCAGCAACTGGCCAACGAAGCCGGAGTCAGTCAGTCCAGcgtc                         <  1:9902695/93‑1 (MQ=255)

                                                                      gATTATCTGGTGCTCCAACCCGATACGGCGCGCCATTTAAGCGCTGAGCAACTGGCCAACGAAGCCGGAGTCAGTCAGTCCAGCGTCGTGAAg                   <  1:7774803/93‑1 (MQ=255)

                                                                      gATTATCTGCTGCTACAACCCGATACGGCGCGCCATTTAAGCTCTCAGCACCTGGCCACCGCAGCCGGAGTCAGTCAGTCCACCGCCGTGAAg                   >  1:6060694/1‑93 (MQ=255)

                                                                      gATTATCTGCTGCTACAACCCGATACGGCGCGCCATTTAAGCTCTCAGCAACTGGCCAACGAAGCCGGAGTCAGTCAGTCCAGCGTCGTGAAg                   >  1:7335210/1‑93 (MQ=255)

                                                                        ttATCTGCTGCTACAACCCGATACGGCGCGCCATTTAAGCTCTCAGCAACTGGCCAACGAAGCCGGAGTCAGTCAGTCCAGCGTCGTGAAGtt                 <  1:1295483/93‑1 (MQ=255)

                                                                         tATCTGCTGCTACAACCCGATACGGCGCGCCATTTAAGCTCTCAGCAACTGGCCAACGAAGCCGGAGTCAGTCAGTCCAGCGTCGTGAAGTTc                <  2:4434026/93‑1 (MQ=255)

                                                                         tATCTGCTGCTACAACCCGATACGGCGCGCCATTTAAGCTCTCAGCAACTGGCCAACGAAGCCGGAGTCAGTCAGTCCAGCGTCGTGAAGTTc                >  2:8323792/1‑93 (MQ=255)

                                                                         tATCTGCTGCTACAACCCGATACGGCGCGCCATGTAAGCTCTCAGCAACTGGCCAACTAAGCCGGAGTCAGTCAGTCCAGCGTCGTGAAGTTc                <  2:5870025/93‑1 (MQ=255)

                                                                          aTCTGCTGCTACAACCCGATACGGCGCGCAATTTAAGCTCACAGCAACTGGCCAACGCCGCCGGAGTCAGTCTGTCCAGCGTCGAAAAGTTcg               <  2:9245663/93‑1 (MQ=255)

                                                                           tCTGCTGCTACAACCCGATACGGCGCGCCATTTAAGCTCTCAGCAACTGGCCAACGAAGCCGGAGTCAGTCAGTCCAGCGTCGTGAAGTTcgc              >  2:8593398/1‑93 (MQ=255)

                                                                             tgctgcTACAACCCGATACGGCGCGCCATTTAAGCTCTCAGCAACTGGCCAACGAAGCCGGAGTCAGTCAGTCCAGCGTCGTGAAGTTCgcgc            <  1:8768025/93‑1 (MQ=255)

                                                                              gctgctACAACCCGATACGGCGCGCCATTTAAGCTCTCAGCAACTGGCCAACGAAGCCGGAGTCAGTCAGTCCAGCGTCGTGAAGTTCGCGCa           <  1:6935782/93‑1 (MQ=255)

                                                                              gctgctACAACCCGATACGGCGCGCCATTTAAGCTCTCAGCAACTGGCCAACGAAGCCGGAGTCAGTCAGTCCAGCGTCGTGAAGTTCGCGCa           >  2:5164813/1‑93 (MQ=255)

                                                                                 gctACAACCCGATACGGCGCGCCATTTAAGCTCTCAGCAACTGGCCAACGAAGCCGGAGTCAGTCAGTCCAGCGTCGTGAAGTTCGCGCaaaa        >  2:2075112/1‑93 (MQ=255)

                                                                                    aCAACCCGATACGGCGCGCCATTTAAGCTCTCAGCAACTGGCCAACGAAGCCGGAGTCAGTCAGTCCAGCGTCGTGAAGTTCGCGCAAAAACt     >  2:1555731/1‑93 (MQ=255)

                                                                                     cAACCCGATACGGCGCGCCATTTAAGCTCTCAGCAACTGGCCAACGAAGCCGGAGTCAGTCAGTCCAGCGTCGTGAAGTTCGCGCAAAAACTc    <  1:8682291/93‑1 (MQ=255)

                                                                                      aaCCCGATACGGCGCGCCATTTAAGCTCTCAGCAACTGGCCAACGa                                                  >  1:4505066/1‑46 (MQ=255)

                                                                                      aaCCCGATACGGCGCGCCATTTAAGCTCTCAGCAACTGGCCAACGAAGCCGGAGTCAGTCAGTCCAGCGTCGTGAAGTTCGCGCAAAAACTCg   >  1:424652/1‑93 (MQ=255)

                                                                                       aCCCGATACGGCGCGCCATTTAAGCTCTCAGCAACTGGCCAACGAAGCCGGAGTCAGTCAGTGCAGCGTCGTGAAGTTCGCGCACAAACTCgg  >  2:6062580/1‑93 (MQ=255)

                                                                                       aCCCGATACGGCGCGCCATTTAAGCTCTCAGCAACTGGCCAACGAAGCCGGAGTCAGTCAGTCCAGCGTCGTGAAGTTCGCGCAAAAACTCgg  <  1:77237/93‑1 (MQ=255)

                                                                                       aCCCGATACGGCGCGCCATTTAAGCTCTCAGCAACTGGCCAACGAAGCCGGAGTCAGTCAGTCCAGCGTCGTGAAGTTCGCGCAAAAACTCgg  >  2:9914047/1‑93 (MQ=255)

AGATATGAACGGCTTACTTCGTATCCGTCAGCGTTACCAGGGGCTTGCCCAGAGTGATAAAAAACTGGCGGATTATCTGCTGCTACAACCTGATACGGCGCGCCATTTAAGCTCTCAGCAACTGGCCAACGAAGCCGGAGTCAGTCAGTCCAGCGTCGTGAAGTTCGCGCAAAAACTCGG  >  NC_000913/2696775‑2696954

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