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	98.1%	325.4 / ‑2.1	106	F456L (TTC→CTC)	ampG	AmpG
Reads supporting (aligned to +/- strand): ref base A (1/1); new base G (45/59); total (46/60)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.25e-01

Only 100 of 114 total aligned reads displayed.

GTAGCGCCAGATAATCCAGCAAACCACCAAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAAGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCAGTGTCCACATGGCAAA  >  NC_000913/451314‑451491

gTAGCGCCAGATAATCCAGCAAACCACCAAGCACGACGCCAGAAAGGGCGACTAAACCCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGt                                                                                       <  1:2631211/93‑1 (MQ=255)

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

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

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

  aGCGCCAGATAATCCAGCAAACCACCAAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGc                                                                                     <  1:3442429/93‑1 (MQ=255)

   gcgcCAGATAATCCAGCAAACCACCAAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCg                                                                                    >  2:2518422/1‑93 (MQ=255)

    cgcCAGATAATCCAGCAAACCACCAAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGt                                                                                   <  2:4038175/93‑1 (MQ=255)

    cgcCAGATAATCCAGCAAACCACCAAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGt                                                                                   <  1:868574/93‑1 (MQ=255)

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

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

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

      ccAGATAATCCAGCAAACCACCAAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAAGTGCGTCa                                                                                 >  1:4527194/1‑93 (MQ=255)

         gATAATCCAGCAAACCACCAAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAAt                                                                              >  1:3004216/1‑93 (MQ=255)

         gATAATCCAGCAAACCACCAAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAAt                                                                              >  1:273596/1‑93 (MQ=255)

          aTAATCCAGCAAACCACCAAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATc                                                                             <  1:107138/93‑1 (MQ=255)

          aTAATCCAGCAAACCACCAAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATc                                                                             >  1:4236316/1‑93 (MQ=255)

             aTCCAGCAAACCACCAAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAg                                                                          <  2:963110/93‑1 (MQ=255)

              tCCAGCAAACCACCAAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAgc                                                                         >  2:247061/1‑93 (MQ=255)

              tCCAGCAAACCACCAAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAgc                                                                         <  1:2066688/93‑1 (MQ=255)

              tCCAGCAAACCACCAAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAgc                                                                         <  1:2066687/93‑1 (MQ=255)

               ccAGCAAACCACCAAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAgcg                                                                        >  1:456113/1‑93 (MQ=255)

                cAGCAAACCACCAAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAgcgc                                                                       >  2:738952/1‑93 (MQ=255)

                  gCAAACCACCAAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGt                                                                     >  1:4162322/1‑93 (MQ=255)

                   cAAACCACCAAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTc                                                                    <  1:3284909/93‑1 (MQ=255)

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

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

                      accaccAAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCAt                                                                 >  2:3156866/1‑93 (MQ=255)

                       ccaccaAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATc                                                                <  2:2013509/93‑1 (MQ=255)

                       ccaccaAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATc                                                                <  1:3057995/93‑1 (MQ=255)

                       ccaccaAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATc                                                                <  2:2965764/93‑1 (MQ=255)

                         accaAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGt                                                              >  2:4370051/1‑93 (MQ=255)

                         accaAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGt                                                              >  1:4696166/1‑93 (MQ=255)

                          ccaAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTc                                                             <  2:965746/93‑1 (MQ=255)

                          ccaAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTc                                                             <  1:526770/93‑1 (MQ=255)

                            aaGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCaa                                                           >  2:3064723/1‑93 (MQ=255)

                            aaGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCaa                                                           >  1:85196/1‑93 (MQ=255)

                             aGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAAt                                                          >  1:3709527/1‑93 (MQ=255)

                              gCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATa                                                         <  1:1190993/93‑1 (MQ=255)

                              gCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATa                                                         <  2:698119/93‑1 (MQ=255)

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

                                acgacgCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGc                                                       >  1:4085869/1‑93 (MQ=255)

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

                                acgacgCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGc                                                       >  2:3260849/1‑93 (MQ=255)

                                acgacgCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGAGCGTCAAATCAAGCGCGTCCATCGTAAACAGc                                                       >  2:1731668/1‑93 (MQ=255)

                                 cgacgCCAGAAAGGGCGACTAAACCCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCa                                                      <  1:2938334/93‑1 (MQ=255)

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

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

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

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

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

                                         gAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACcac                                              <  2:3965573/93‑1 (MQ=255)

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

                                         gAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACcac                                              <  1:1452159/93‑1 (MQ=255)

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

                                            aGGGCGCCTAACACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAAGTGGGTCAAATCAAGCACGTCCACCGTCAATAGCAGTCACCACACg                                           >  2:1863625/1‑93 (MQ=255)

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

                                             gggCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAAACCAGCGCGTCCATCGTCAATAGCCGTAACCAAACgc                                          >  1:4616992/1‑92 (MQ=255)

                                              ggCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGc                                         <  1:149048/93‑1 (MQ=255)

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

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

                                                 gACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCaa                                      >  2:1154455/1‑93 (MQ=255)

                                                 gACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCaa                                      >  2:3173829/1‑93 (MQ=255)

                                                 gACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCaa                                      <  2:3916674/93‑1 (MQ=255)

                                                  aCTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAAc                                     >  2:642166/1‑93 (MQ=255)

                                                  aCTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAAc                                     >  2:2158863/1‑93 (MQ=255)

                                                  aCTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAAc                                     >  1:1117692/1‑93 (MQ=255)

                                                  aCTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAAc                                     >  1:2273125/1‑93 (MQ=255)

                                                   cTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACa                                    <  1:4090231/93‑1 (MQ=255)

                                                   cTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACa                                    >  2:2645428/1‑93 (MQ=255)

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

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

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

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

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

                                                        aCCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTg                               <  1:4463000/93‑1 (MQ=255)

                                                        aCCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTg                               >  1:2510669/1‑93 (MQ=255)

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

                                                         acccGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGa                              <  1:773603/92‑1 (MQ=255)

                                                          cccGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGAc                             <  2:709868/93‑1 (MQ=255)

                                                          cccGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGAc                             <  1:1708143/93‑1 (MQ=255)

                                                          cccGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGAc                             <  1:136642/93‑1 (MQ=255)

                                                          cccGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGAc                             <  2:4550727/93‑1 (MQ=255)

                                                          cacGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTGAAATCCAGCGCGTCCATCGTCTATAGCAGTAACCACCCGGCCAACAGGCTGAc                             >  2:2532394/3‑93 (MQ=255)

                                                             gACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCAAACAGGCTGAcgcc                          >  2:1587929/1‑93 (MQ=255)

                                                                ttCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGAcgcccgc                       >  1:2765614/1‑93 (MQ=255)

                                                                    agcagagcagGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCgccgcc                   >  1:2430390/1‑93 (MQ=255)

                                                                     gcagagcagGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCa                  <  1:225018/93‑1 (MQ=255)

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

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

                                                                          gcaggcagGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCAGTGTc             <  2:479182/93‑1 (MQ=255)

                                                                          gcaggcagAAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACTCCCGCCGCCAGTGTc             <  2:479183/93‑1 (MQ=255)

                                                                           gaggcaggAAAGAGAGGAGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCAGTGTcc            <  1:4310885/92‑1 (MQ=255)

                                                                            aggcaggAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCAGTGTCca           <  1:2644959/93‑1 (MQ=255)

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

                                                                               caggAAAGAGAGGTGCGGCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCAGTGTCCACAt        <  1:4111508/93‑1 (MQ=255)

                                                                               caggAAACAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCAGTGTCCACAt        <  2:890024/93‑1 (MQ=255)

                                                                                 ggAAAGAGCGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTACCCCCACGGCCAACAGGCTGACGCCCGCCGCCAGTGTCCACATgg      <  2:478151/93‑1 (MQ=255)

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

                                                                                    aaGAGAGGGGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCAGTGTCCACATGGCaa   <  1:2676076/93‑1 (MQ=255)

                                                                                     agagagGTGAGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCATCAGGCTGACGCCCGCCGCCAGTGTCCACATGGCaaa  <  2:2736829/93‑1 (MQ=255)

GTAGCGCCAGATAATCCAGCAAACCACCAAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAAGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCAGTGTCCACATGGCAAA  >  NC_000913/451314‑451491

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