BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	451,404	A→G	F456L (TTC→CTC)	ampG ←	AmpG

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	451,404	0	A	G	99.5%	660.6 / ‑5.1	206	F456L (TTC→CTC)	ampG	AmpG
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (118/87); total (119/87)
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.20e-01

Only 100 of 216 total aligned reads displayed.

GTAGCGCCAGATAATCCAGCAAACCACCAAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAAGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCAGTGTCCACATGGCAAAGG  >  NC_000913/451314‑451493

gTAGCGCCAGATAATCCAGCAAACCACCAAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGt                                                                                         <  1:4959516/93‑1 (MQ=255)

gTAGCGCCAGATAATCCAGCAAACCACCAAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGt                                                                                         <  1:6129777/93‑1 (MQ=255)

 tAGCGCCAGATAATCCAGCAAACCACCAAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTg                                                                                        <  2:5012479/93‑1 (MQ=255)

  aGNGCCNAATAATCCAGCAAACNACCAAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGc                                                                                       <  2:5609317/93‑1 (MQ=255)

  aGCGCCAGATAATCCAGCAAACCACCAAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGc                                                                                       >  1:4804694/1‑93 (MQ=255)

  aGCGCCAGATAATCCAGCAAACCACCAAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGc                                                                                       >  1:1491189/1‑93 (MQ=255)

  aGCGCCAGATAATCCAGCAAACCACCAAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGc                                                                                       >  1:3280673/1‑93 (MQ=255)

     gcCAGATAATCCAGCAAACCACCAAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTc                                                                                    <  1:1843753/93‑1 (MQ=255)

     gcCAGATAATCCAGCAAACCACCAAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTc                                                                                    <  2:4494460/93‑1 (MQ=255)

      ccAGATAATCCAGCAAACCACCAAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCa                                                                                   <  1:5712567/93‑1 (MQ=255)

      ccAGATAATCCAGCAAACCACCAAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCa                                                                                   >  1:2438237/1‑93 (MQ=255)

       cAGATAATCCAGCAAACCACCAAGCACGACGCCAGAAAGGGTGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGCAAGAGAGGTGCGGCaa                                                                                  >  1:3743445/1‑93 (MQ=255)

       cAGATAATCCAGCAAACCACCAAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCaa                                                                                  >  2:3139682/1‑93 (MQ=255)

        aGATAATCCAGCAAACCACCAAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCaaa                                                                                 <  2:2270879/93‑1 (MQ=255)

         gATAATCCAGCAAACCACCAAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGGGCGTCAAAt                                                                                >  2:1669659/1‑93 (MQ=255)

           tAATCCAGCAAACCACCAAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATcc                                                                              >  1:299781/1‑93 (MQ=255)

           tAATCCAGCAAACCACCAAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATcc                                                                              <  1:80102/93‑1 (MQ=255)

             aTCCAGCAAACCACCAAGCACGACGCCAGAAAGGGCGACTAAACCCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAg                                                                            <  1:1059353/93‑1 (MQ=255)

                cAGCAAACCACCAAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAgcgc                                                                         <  1:2821497/93‑1 (MQ=255)

                 aGCAAACCACCAAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGcgcg                                                                        <  1:5233570/93‑1 (MQ=255)

                 aGCAAACCACCAAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGcgcg                                                                        <  2:6542411/93‑1 (MQ=255)

                    aaaCCACCAAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGACAGAGAGGTGCGTCCAATCCAGCGCGTcc                                                                     >  1:169802/1‑93 (MQ=255)

                    aaaCCACCAAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTcc                                                                     >  1:5976747/1‑93 (MQ=255)

                     aaCCACCAAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCa                                                                    >  1:3772372/1‑93 (MQ=255)

                     aaCCACCAAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCa                                                                    <  1:4477488/93‑1 (MQ=255)

                     aaCCACCAAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCa                                                                    <  2:5848216/93‑1 (MQ=255)

                      accaccAAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCAt                                                                   <  2:1822991/93‑1 (MQ=255)

                               cACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAg                                                          <  1:3964806/93‑1 (MQ=255)

                                acgacgCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGc                                                         <  2:6627985/93‑1 (MQ=255)

                                 cgacgCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCa                                                        <  2:4346495/93‑1 (MQ=255)

                                  gacgCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAg                                                       >  2:5798765/1‑93 (MQ=255)

                                  gacgCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAg                                                       <  1:2263086/93‑1 (MQ=255)

                                   acgCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAga                                                      >  2:1475731/1‑92 (MQ=255)

                                     gCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTaa                                                    <  2:4965405/93‑1 (MQ=255)

                                      ccAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAAc                                                   >  1:1018074/1‑93 (MQ=255)

                                      ccAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAAc                                                   <  2:4720776/93‑1 (MQ=255)

                                      ccAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAAc                                                   <  2:640353/93‑1 (MQ=255)

                                       cAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAAcc                                                  >  2:5116197/1‑93 (MQ=255)

                                       cAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAAGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAAcc                                                  >  2:1566198/1‑93 (MQ=255)

                                        aGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACca                                                 >  1:1908776/1‑93 (MQ=255)

                                        aGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACca                                                 <  1:3854294/93‑1 (MQ=255)

                                        aGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGAAGTAACca                                                 >  1:5318630/1‑93 (MQ=255)

                                        aGAAAGGACGACTAAAACCCCACCTTCCAGCAGAGCAGTCAGGAAAGAGAGGTGCGTCAACTCCAGCGCGTCCATCGTCAATAGCAGTAACca                                                 <  2:1236364/93‑1 (MQ=255)

                                         gAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACcac                                                >  1:5546852/1‑93 (MQ=255)

                                           aaGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCacac                                              <  2:6327225/93‑1 (MQ=255)

                                           aaGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCacac                                              >  1:2758519/1‑93 (MQ=255)

                                           aaGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCacac                                              <  1:6246897/93‑1 (MQ=255)

                                            aGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACg                                             >  1:894496/1‑93 (MQ=255)

                                             gggCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACgg                                            <  1:636724/93‑1 (MQ=255)

                                              ggCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAAGCACACGGc                                           >  1:1876024/1‑93 (MQ=255)

                                               gCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGcc                                          <  2:628606/93‑1 (MQ=255)

                                                cGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCa                                         >  1:2565027/1‑93 (MQ=255)

                                                  aCTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAAc                                       <  2:3907656/93‑1 (MQ=255)

                                                    tAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCCCACGGCCAACAg                                     >  2:4635365/1‑93 (MQ=255)

                                                    tAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAg                                     <  1:5471344/93‑1 (MQ=255)

                                                    tAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAg                                     >  2:5761752/1‑93 (MQ=255)

                                                    caaaaCCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAg                                     <  2:2472330/92‑1 (MQ=255)

                                                      aaaCCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCGATAGCAGTAACCACACGGCCAACAGGc                                   <  2:1255227/93‑1 (MQ=255)

                                                      aaaCCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGc                                   <  2:6497279/93‑1 (MQ=255)

                                                         ccccGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGa                                <  1:3929783/93‑1 (MQ=255)

                                                         ccccGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGa                                >  1:1532642/1‑93 (MQ=255)

                                                         ccccGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGa                                >  1:1208237/1‑93 (MQ=255)

                                                          cccGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGAc                               >  1:6061466/1‑93 (MQ=255)

                                                           ccGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGAcg                              >  2:3960854/1‑93 (MQ=255)

                                                            cGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGAcgc                             >  1:4904390/1‑93 (MQ=255)

                                                              aCTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGAcgccc                           >  2:2620613/1‑93 (MQ=255)

                                                              aCTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGAcgccc                           >  1:3169895/1‑93 (MQ=255)

                                                              aCTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGAcgacc                           >  2:5163138/1‑93 (MQ=255)

                                                               cTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGAcgcccg                          <  1:1856441/93‑1 (MQ=255)

                                                                ttCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGAcgcccgc                         <  1:1384663/93‑1 (MQ=255)

                                                                 tCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGAcgcccgcc                        <  1:775028/93‑1 (MQ=255)

                                                                  ccAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCccgccg                       >  1:3307756/1‑93 (MQ=255)

                                                                    agcagagcagGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCgccgcc                     <  1:2048106/93‑1 (MQ=255)

                                                                    agcagagcagGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCgccgcc                     <  2:42228/93‑1 (MQ=255)

                                                                     gcagagcagGCAGGAAAGGGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAGTAGCAGTAACCACACAGCCAACAGGCTGACGCCCGCCGCCa                    >  2:77451/1‑93 (MQ=255)

                                                                       agagcagGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCGAACAGGCTGCCGgccgccgcccgg                  >  2:321959/1‑90 (MQ=255)

                                                                        gagcagGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCAgtg                 >  2:2662687/1‑93 (MQ=255)

                                                                        gagcagGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCAgtg                 >  2:6316333/1‑93 (MQ=255)

                                                                        gagcagGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCAgtg                 >  2:1141738/1‑93 (MQ=255)

                                                                        gagcagGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCAgtg                 >  2:733271/1‑93 (MQ=255)

                                                                         agcagGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCAgtgt                >  1:3694356/1‑93 (MQ=255)

                                                                          gcaggcagGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCAGTGTc               >  1:4223832/1‑93 (MQ=255)

                                                                          gcaggcagGAAAGAGAGGAGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCAGTGTc               >  1:6644516/1‑93 (MQ=255)

                                                                           caggcaggAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCAGTGTcc              <  1:5036712/93‑1 (MQ=255)

                                                                             ggcaggAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCCATAGCAGTAACCACACGGCCAACAGGCTGCCGCCCGCCGCCAGTGTCcac            <  1:6726021/93‑1 (MQ=255)

                                                                             ggcaggAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCAGTGTCcac            >  1:4296413/1‑93 (MQ=255)

                                                                              gcaggAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCGGTGTCcaca           >  2:3049230/1‑93 (MQ=255)

                                                                              gcaggAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCAGTGTCcaca           <  1:3538819/93‑1 (MQ=255)

                                                                              gcaggAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCAGTGTCcaca           <  1:3195515/93‑1 (MQ=255)

                                                                               caggAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCAGTGTCCACAt          <  1:2872232/93‑1 (MQ=255)

                                                                                 ggAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCAGTGTCCACATgg        >  1:4623303/1‑93 (MQ=255)

                                                                                  gAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCAGTGTCCACATGGc       >  2:4417166/1‑93 (MQ=255)

                                                                                   aaaGAGAGGTGCGTCAAATCCAGCGCGTCCATCCGCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCAGTGTCCACATGGCa      <  2:2884827/93‑1 (MQ=255)

                                                                                    aaGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCAGTGTCCACATGGCaa     >  1:1601081/1‑93 (MQ=255)

                                                                                    aaGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCAGTGTCCACATGGCaa     <  2:6690681/93‑1 (MQ=255)

                                                                                     agagagGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCAGTGTCCACATGGCaaa    <  2:791841/93‑1 (MQ=255)

                                                                                      gagaGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACTCCCGCCGCCAGTGTCCACATGGCAAAg   <  2:4198487/93‑1 (MQ=255)

                                                                                      gagaGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCAGTGTCCACATGGCAAAg   >  1:538113/1‑93 (MQ=255)

                                                                                      gagaGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCAGTGTCCACATGGCAAAg   >  2:1036782/1‑93 (MQ=255)

                                                                                       agagGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCAGTGTCCACATGGCAAAgg  <  2:6661738/93‑1 (MQ=255)

GTAGCGCCAGATAATCCAGCAAACCACCAAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAAGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCAGTGTCCACATGGCAAAGG  >  NC_000913/451314‑451493

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