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.4%	501.9 / NA	158	F456L (TTC→CTC)	ampG	AmpG
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (79/78); total (79/79)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 174 total aligned reads displayed.

GTAGCGCCAGATAATCCAGCAAACCACCAAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAAGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCAGTGTCCACATGGCAAAGG  >  NC_000913/451314‑451493

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

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

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

 tAGCGCCAGATAATCCAGCAAACCACCAAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTg                                                                                        >  2:2082011/1‑93 (MQ=255)

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

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

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

        aGATAATCCAGCAAACCACCAAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCaaa                                                                                 >  1:1598214/1‑93 (MQ=255)

        aGATAATCCAGCAAACCACCAAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGCGCGTCaaa                                                                                 >  1:3525004/1‑93 (MQ=255)

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

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

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

           tAATCCAGCAAACCACCAAGCCCGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATcc                                                                              >  1:2480727/1‑93 (MQ=255)

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

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

                cAGCAAACCACCAAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAgtgc                                                                         >  1:4822330/1‑93 (MQ=255)

                  gCAAACCACCAAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGTAAGAGAGGTGCGTCAAATGCAGCGCGt                                                                       >  2:2336075/1‑93 (MQ=255)

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

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

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

                    aaaCCACCAAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGAGCGTCAAATCCAGCGCGTcc                                                                     >  2:4738116/1‑93 (MQ=255)

                       ccaccaAGCACGACGCCAGAAAGGGCGACTAAACCCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATc                                                                  <  1:1120735/93‑1 (MQ=255)

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

                       ccaccaAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAAGTGCGTCAAATCCAGCGCGTCCATc                                                                  <  2:3053901/93‑1 (MQ=255)

                         accaAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGt                                                                <  2:3445007/93‑1 (MQ=255)

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

                          ccaAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTc                                                               >  1:1295057/1‑93 (MQ=255)

                          ccaAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTc                                                               >  1:1523141/1‑93 (MQ=255)

                           caAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCa                                                              <  2:5588351/93‑1 (MQ=255)

                           caAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCa                                                              >  1:2260531/1‑93 (MQ=255)

                            aaGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCaa                                                             <  1:2040643/93‑1 (MQ=255)

                            aaGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCaa                                                             <  1:1279070/93‑1 (MQ=255)

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

                              gCACGACGCCAGAAAGGGAGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATa                                                           >  2:5193146/1‑93 (MQ=255)

                               cACGACGCCAGAAAGGGCGACTAAACCCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAg                                                          <  2:591122/93‑1 (MQ=255)

                                 cgacgCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCa                                                        >  1:6236237/1‑93 (MQ=255)

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

                                    cgCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAAACCAGCGCGTCCATCGTCAATAGCAGTa                                                     >  1:2429319/1‑93 (MQ=255)

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

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

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

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

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

                                        aGAAAGCGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACca                                                 <  1:5925541/93‑1 (MQ=255)

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

                                          aaaGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACcaca                                               >  2:1849826/1‑93 (MQ=255)

                                          aaaGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGAGCGTCCATCGTCAATAGCAGTAACcagc                                               >  1:1859130/1‑91 (MQ=255)

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

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

                                              ggCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCTAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGc                                           <  2:5136758/93‑1 (MQ=255)

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

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

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

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

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

                                                   cTAAAACCACGAAATCCAGCAGAACAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCACGTCCATAGACAAAAGCAGAAACAACACGGCCAACa                                      <  1:6300156/93‑1 (MQ=255)

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

                                                     aaaaCCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAgg                                    >  1:2052526/1‑93 (MQ=255)

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

                                                      aaaCCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGc                                   >  2:3439722/1‑93 (MQ=255)

                                                      aaaCCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGc                                   >  2:5046575/1‑93 (MQ=255)

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

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

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

                                                           ccGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCTAACAGACTGAcg                              >  1:4122279/1‑93 (MQ=255)

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

                                                             gACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGAcgcc                            <  2:1568767/93‑1 (MQ=255)

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

                                                               cTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGAcgcccg                          >  2:330964/1‑93 (MQ=255)

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

                                                                 tCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGAcgcccgcc                        >  1:3666016/1‑93 (MQ=255)

                                                                   cagcagAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCcgccgc                      <  2:3979116/93‑1 (MQ=255)

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

                                                                     gcagagcagGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCa                    >  2:4506184/1‑93 (MQ=255)

                                                                       agagcagGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCAgt                  >  1:6052824/1‑93 (MQ=255)

                                                                         agcagGCAGGAAAGAGAGGTGCGTCAAATCGAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCAgtgg                >  2:1849666/1‑92 (MQ=255)

                                                                         agcagGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCAgtgt                <  1:1943726/93‑1 (MQ=255)

                                                                         agcagGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCAgtgt                <  1:4506655/93‑1 (MQ=255)

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

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

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

                                                                                aggAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCAGTGTCCACATg         <  2:5791408/93‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

                                                                                      gagaGGTGCGTCAAATCCAGCGCGTCCGTTGTCAATAGCAGGATCCTCACGGCCAACAGGCTGAGGCCCTCCGCCAGTGTCCACATGGCAAAg   <  2:3173887/93‑1 (MQ=255)

                                                                                      gagaGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCAGTGTCCACATGGCAAAg   <  1:4904354/93‑1 (MQ=255)

                                                                                      gagaGGCGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCAGTGTCCACATGGCAAAg   <  1:5945729/93‑1 (MQ=255)

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

                                                                                       agagGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCAGTGTCCACATGGCAAAgg  >  2:674112/1‑93 (MQ=255)

GTAGCGCCAGATAATCCAGCAAACCACCAAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAAGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCAGTGTCCACATGGCAAAGG  >  NC_000913/451314‑451493

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