BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,533,160	A→G	L116P (CTG→CCG)	envZ ←	protein histidine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,533,160	0	A	G	99.5%	1145.5 / ‑3.8	369	L116P (CTG→CCG)	envZ	protein histidine kinase
Reads supporting (aligned to +/- strand): ref base A (1/1); new base G (153/214); total (154/215)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.61e-01

Only 100 of 387 total aligned reads displayed.

CGTACCCAGATATTGGGCGACAGCCAGGTTTTCAGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCAGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAG  >  NC_000913/3533070‑3533249

cgTACCCAGATATTGGGCGACAGCCAGGTTTTCAGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGt                                                                                         <  2:1110155/93‑1 (MQ=255)

cgTACCCAGATATTGGGCGACAGCCAGGTTTTCAGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGt                                                                                         <  2:3642022/93‑1 (MQ=255)

 gTACCCAGATATTGGGCGACAGCCAGGTTTTCAGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGtt                                                                                        <  1:32722/93‑1 (MQ=255)

   aCCCAGATATTGGGCGACAGCCAGGTTTTCAGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTtgc                                                                                      >  1:4446140/1‑93 (MQ=255)

   aCCAAGATATTGGGCGACAGCCAGGTTTTCAGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTtgc                                                                                      <  1:5395552/93‑1 (MQ=255)

     ccAGATATTGGGCGACAGCCAGGTTTTCAGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTtgctg                                                                                    <  2:1052445/93‑1 (MQ=255)

     ccAGATATTGGGCGACAGCCAGGTTTTCAGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTtgctg                                                                                    <  2:6810899/93‑1 (MQ=255)

      cAGATATTGGGCGACAGCCAGGTTTTCAGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTtgctgc                                                                                   >  2:3882426/1‑93 (MQ=255)

      cAGATATTGGGCGACAGCCAGGTTTTCAGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTtgctgc                                                                                   >  2:876689/1‑93 (MQ=255)

        gATATTGGGCGACAGCCAGGTTTTCAGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGGCCACTTCCGTCGGGCCGCCCGGTTGCTgcgc                                                                                 >  2:5510158/1‑93 (MQ=255)

           atTGGGCGACAGCCAGGTTTTCAGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCAt                                                                              >  1:7560395/1‑93 (MQ=255)

              gggCGACAGCCAGGTTTTCAGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTg                                                                           <  1:2342620/93‑1 (MQ=255)

               ggCGACAGCCAGGTTTTCAGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGa                                                                          <  2:2661952/93‑1 (MQ=255)

               ggCGACAGCCAGGTTTTCAGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGa                                                                          >  1:7164801/1‑93 (MQ=255)

                gCGACAGCCAGGTTTTCAGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGAt                                                                         >  1:1631940/1‑93 (MQ=255)

                  gACAGCCAGGTTTTCAGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATgg                                                                       >  2:2905009/1‑93 (MQ=255)

                   aCAGCCAGGTTTTCAGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGc                                                                      >  2:3212192/1‑93 (MQ=255)

                    cAGCCAGGTTTTCAGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCt                                                                     <  1:4000054/93‑1 (MQ=255)

                     aGCCAGGTTTTCAGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCtt                                                                    >  2:6345335/1‑93 (MQ=255)

                     aGCCAGGTTTTCAGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCtt                                                                    <  1:1281170/93‑1 (MQ=255)

                      gCCAGGTTTTCAGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTa                                                                   >  2:5962463/1‑93 (MQ=255)

                           gTTTTCAGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAAt                                                              <  1:5509889/93‑1 (MQ=255)

                           gTTTTCAGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAAt                                                              <  2:937495/93‑1 (MQ=255)

                            ttttCAGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAAtt                                                             >  2:2934629/1‑93 (MQ=255)

                            gtttCAGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAAtt                                                             >  2:7880778/2‑93 (MQ=255)

                               tCAGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCAt                                                          <  1:1113816/93‑1 (MQ=255)

                                taGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATa                                                         <  1:6116382/92‑1 (MQ=255)

                                cAGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATa                                                         <  1:183560/93‑1 (MQ=255)

                                cAGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATa                                                         <  2:7417502/93‑1 (MQ=255)

                                cAGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATa                                                         <  1:3345624/93‑1 (MQ=255)

                                 aGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAg                                                        <  1:4988221/93‑1 (MQ=255)

                                 aGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAg                                                        <  2:162687/93‑1 (MQ=255)

                                  gCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAgt                                                       >  1:1501789/1‑93 (MQ=255)

                                  gCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAgt                                                       <  1:2561699/93‑1 (MQ=255)

                                    cAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAgtgt                                                     >  2:7589306/1‑93 (MQ=255)

                                    cAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAgtgt                                                     <  2:6026939/93‑1 (MQ=255)

                                    cAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAgtgt                                                     <  1:7340077/93‑1 (MQ=255)

                                     cgacgacAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGtt                                                    <  1:1944598/92‑1 (MQ=255)

                                     aGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGtt                                                    <  2:1416/93‑1 (MQ=255)

                                      gacgacAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTg                                                   <  1:5670368/93‑1 (MQ=255)

                                      gacgacAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTg                                                   >  2:7139665/1‑93 (MQ=255)

                                       acgacAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTgc                                                  <  2:6123154/93‑1 (MQ=255)

                                       acgacAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTgc                                                  <  2:6459010/93‑1 (MQ=255)

                                          acAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGcc                                               >  1:3702619/1‑93 (MQ=255)

                                           cAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTgcgctc                                              <  1:1711582/93‑3 (MQ=255)

                                            aGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCa                                             >  1:7238382/1‑93 (MQ=255)

                                            aGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTGAAGAATTCATAGTGTTGCGCCCa                                             >  2:4320163/1‑93 (MQ=255)

                                              gCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAAc                                           <  2:2348724/93‑1 (MQ=255)

                                              gCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAAc                                           <  2:7054377/93‑1 (MQ=255)

                                              gCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAAc                                           <  2:3154999/93‑1 (MQ=255)

                                                gAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGc                                         <  2:3387435/93‑1 (MQ=255)

                                                gAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGc                                         >  1:5101938/1‑93 (MQ=255)

                                                 aaCTTTTGTTGACCTCAACGCGCTCTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCa                                        <  1:418418/93‑1 (MQ=255)

                                                 aaCTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCa                                        <  1:2670213/93‑1 (MQ=255)

                                                 aaCTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCa                                        <  2:1713565/93‑1 (MQ=255)

                                                 aaCTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCa                                        >  2:6546557/1‑93 (MQ=255)

                                                  aCTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAg                                       >  1:6789064/1‑93 (MQ=255)

                                                        gttgCCCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGc                                 <  2:404999/93‑1 (MQ=255)

                                                        gttgACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGc                                 >  2:3927630/1‑93 (MQ=255)

                                                          tgACCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCct                               >  2:4659835/1‑93 (MQ=255)

                                                            aCCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCctct                             >  1:4650886/1‑93 (MQ=255)

                                                            aCCTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCctct                             <  1:6420972/93‑1 (MQ=255)

                                                             ccTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCtt                            <  1:3067074/93‑1 (MQ=255)

                                                             ccTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCtt                            >  1:4547208/1‑93 (MQ=255)

                                                             ccTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCtt                            <  1:2138075/93‑1 (MQ=255)

                                                              cTCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTc                           >  2:183168/1‑93 (MQ=255)

                                                               tCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCg                          >  2:2210242/1‑93 (MQ=255)

                                                               tCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCg                          <  2:5179951/93‑1 (MQ=255)

                                                               tCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCg                          >  1:5967090/1‑93 (MQ=255)

                                                               tCAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCg                          >  1:7419489/1‑93 (MQ=255)

                                                                cAACGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGACTCTTCgg                         <  1:6466508/93‑1 (MQ=255)

                                                                 caCGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGc                        <  2:2731700/92‑1 (MQ=255)

                                                                 aaCGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGc                        >  1:5775648/1‑93 (MQ=255)

                                                                 aaCGCGCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGc                        <  2:475424/93‑1 (MQ=255)

                                                                   cgcgCACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAg                      >  1:7709599/1‑93 (MQ=255)

                                                                    gcgcACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGc                     >  2:4140268/1‑93 (MQ=255)

                                                                     cgcACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGcc                    <  2:4490561/93‑1 (MQ=255)

                                                                      gcACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCt                   <  2:5073433/93‑1 (MQ=255)

                                                                      gcACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCt                   >  1:1011004/1‑93 (MQ=255)

                                                                       cACTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTc                  <  2:7328461/93‑1 (MQ=255)

                                                                         cTTCGGTCGGCCCGCCGGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCGGACCTGCCTCTTCGGCAGCCTCGt                <  2:3027103/93‑1 (MQ=255)

                                                                         cTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGt                >  2:4465941/1‑93 (MQ=255)

                                                                         cTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGt                <  2:6908968/93‑1 (MQ=255)

                                                                         cTTCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGt                <  2:2028944/93‑1 (MQ=255)

                                                                          ttCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGtt               >  1:4850672/1‑93 (MQ=255)

                                                                           tCCGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTg              >  1:3857490/1‑93 (MQ=255)

                                                                           cccGTCGGGCTGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGGGTTGCGCCCAACTCAGACCTGCCTCCTTGGCAGCCTCGTTg              >  1:7046702/2‑93 (MQ=255)

                                                                            cagTCGGGCCGCCGGGTTGCTGCACCACCTGATGGCTAAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCCCTCCGGAAGCCTCGTTgg             <  2:3741861/91‑1 (MQ=255)

                                                                             cGTCGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGa            >  1:7838713/1‑93 (MQ=255)

                                                                              gTCGGGCCGCCCGGTTGCTGCGCCCTCTGATGGCTCAAGAATTCATAGTGCTGCGACCAACGCAGACCTGCCTCTTCCGCAGCCTCGGTGGAg           >  2:3525467/1‑93 (MQ=255)

                                                                                cGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTa         >  1:2901844/1‑93 (MQ=255)

                                                                                cGGGCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTa         <  1:968891/93‑1 (MQ=255)

                                                                                 gggCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTag        >  2:1799193/1‑93 (MQ=255)

                                                                                 gggCCGCCCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTag        <  1:2211082/93‑1 (MQ=255)

                                                                                   gccgccCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTagag      <  1:3576306/93‑1 (MQ=255)

                                                                                    ccgccCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTagaga     <  2:1228242/93‑1 (MQ=255)

                                                                                    ccgccCGGTTGCTCCGCCATCTGATGCCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTagaga     <  2:2620297/93‑1 (MQ=255)

                                                                                     cgccCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTagagag    <  2:5339179/93‑1 (MQ=255)

                                                                                      gccCGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTagagaga   <  2:5672390/93‑1 (MQ=255)

                                                                                       cccGGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTagagagag  >  1:5345668/1‑93 (MQ=255)

CGTACCCAGATATTGGGCGACAGCCAGGTTTTCAGCCAGACGACAGGCGAACTTTTGTTGACCTCAACGCGCACTTCCGTCGGGCCGCCCAGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAG  >  NC_000913/3533070‑3533249

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