BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,613,076	T→C	T276A (ACT→GCT)	perM ←	putative permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,613,076	0	T	C	99.1%	344.4 / ‑4.6	107	T276A (ACT→GCT)	perM	putative permease
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (42/64); total (43/64)
Fisher's exact test for biased strand distribution p-value = 4.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.75e-01

Only 100 of 109 total aligned reads displayed.

ACCGCTTCGGAGAACAACACCGGCACTAACAGGTTGCCGTCCAGCGCCTGAATAATCAGATACACCGCGAAGCAGCTCCAGAATTCAGTACCTGCGCCAAACTGGAATAGCGCCACGCCAACCACCGGAATGGTCACCACAAATGCGCCAATGTACGGAATAAGAACCGAGAAGCCGA  >  NC_000913/2612988‑2613165

aCCGCTTCGGAGAACAACACCGGCACTAACAGGTTGCCGTCCAGCGCCTGAATAATCAGATACACCGCGAAGCAGCTCCAGAATTCGGCACCt                                                                                       >  1:1110700/1‑93 (MQ=255)

aCCGCTTCGGAGAACAACACCGGCACTAACAGGTTGCCGTCCAGCGCCTGAATAATCAGATACACCGCGAAGCAGCTCCAGAATTCAGCACCt                                                                                       >  1:510939/1‑93 (MQ=255)

  cGCTTCGGAGAACAACACCGGCACTAACAGGTTGCCGTCCAGCGCCTGAATAATCAGATACACCGCGAAGCAGCTCCAGAATTCAGCACCTgc                                                                                     <  1:4581148/93‑1 (MQ=255)

  cGCTTCGGAGAACAACACCGGCACTAACAGGTTGCCGTCCAGCGCCTGAATAATCAGATACACCGCGAAGCAGCTCCAGAATTCAGCACCTgc                                                                                     <  2:3861132/93‑1 (MQ=255)

   gCTTCGGAGAACAACACCGGCACTAACAGGTTGCCGTCCAGCGCCTGAATAATCAGATACACCGCGAAGCAGCTCCAGAATTCAGCACCTgcg                                                                                    >  2:1626852/1‑93 (MQ=255)

   gCTTCGGAGAACAACACCGGCACTAACAGGTTGCCGTCCAGCGCCTGAATAATCAGATACACCGCGAAGCAGCTCCAGAATTCAGCACCTgcg                                                                                    >  2:2878384/1‑93 (MQ=255)

        ggAGAACAACACCGGCACTAACAGGTTGCCGTCCAGCGCCTGAATAATCAGATACACCGCGAAGCAGCTCCAGAATTCAGCACCTGCGCCaaa                                                                               <  1:735583/93‑1 (MQ=255)

         gagaACAACACCGGCACTAACAGGTTGCCGTCCAGCGCCTGAATAATCAGATACACCGCGAAGCAGCTCCAGAATTCAGCACCTGCGCCAAAc                                                                              >  2:1665263/1‑93 (MQ=255)

         gagaACAACACCGGCACTAACAGGTTGCCGTCCAGCGCCTGAATAATCAGATACACCGCGAAGCAGCTCCAGAATTCAGCACCTGCGCCAAAc                                                                              <  2:1068426/93‑1 (MQ=255)

         gagaACAACACCCGCACTAACAGGTTGCCGTCCACCGCCTGAATAATCAGATACACCGCGAAGCAGCTCCAGAATTCAGCACCTGCGCCAAAc                                                                              >  2:1665262/1‑93 (MQ=255)

          agaACAACACCGGCACTAACAGGTTGCCGTCCAGCGCCTGAATAATCAGATACACCGCGAAGCAGCTCCAGAATTCAGCACCTGCGCCAAACt                                                                             <  2:884061/93‑1 (MQ=255)

          agaACAACACCGGCACTAACAGGTTGCCGTCCAGCGCCTGAATAATCAGATACACCGCGAAGCAGCTCCAGAATTCAGCACCTGCGCCAAACt                                                                             <  1:1068200/93‑1 (MQ=255)

           gaACAACACCGGCACTAACAGGATGCCGTCCAGCGCCTGAATAATCAGATACACCGCGAAGCAGCTCCAGAATTCAGCACCTGCGCCAAACTg                                                                            <  1:1604580/93‑1 (MQ=255)

              caacaCCGGCACTAACAGGTTGCCGTCCAGCGCCTGAATAATCAGATACACCGCGAAGCAGCTCCAGAATTCAGCACCTGCGCCAAACTGGaa                                                                         <  2:3165363/93‑1 (MQ=255)

               aacaCCGGCACTAACAGGTTGCCGTCCAGCGCCTGAATAATCAGATACACCGCGAAGCAGCTCCAGAATTCAGCACCTGCGCCAAACTGGAAt                                                                        <  1:887145/93‑1 (MQ=255)

               aacaCCGGCACTAACAGGTTGCCGTCCAGCGCCTGAATAATCAGATACACCGCGAAGCAGCTCCAGAATTCAGCACCTGCGCCAAACTGGAAt                                                                        <  2:601632/93‑1 (MQ=255)

               aacaCCGGCACTAACAGGTTGCCGTCCAGCGCCTGAATAATCAGATACACCGCGAAGCAGCTCCAGAATTCAGCACCTGCGCCAAACTGGAAt                                                                        <  2:4142320/93‑1 (MQ=255)

               aacaCCGGCACTAACAGGTTGCCGTCCAGCGCCTGAATAATCAGATACACCGCGAAGCAGCTCCAGAATTCAGCACCTGCGCCAAACTGGAAt                                                                        >  1:457347/1‑93 (MQ=255)

                accccGGCACTAACGGGTTGCCGTCCAGCGCCTGAATACTCAGATACACCGCGAAGCAGCGCCAGAATTCAGCCCCTGCGCCAAACTGGAATa                                                                       <  2:3618078/90‑1 (MQ=255)

                acacCGGCACTAACAGGTTGCCGTCCAGCGCCTGAATAATCAGATACACCGCGAAGCAGCTCCAGAATTCAGCACCTGCGCCAAACTGGAATa                                                                       <  1:4106591/93‑1 (MQ=255)

                acacCGGCACTAACAGGTTGCCGTCCAGCGCCTGAATAATCAGATACACCGCGAAGCAGCTCCAGAATTCAGCACCTGCGCCAAACTGGAATa                                                                       <  1:2825034/93‑1 (MQ=255)

                acacCGGCACTAACAGGTTGCCGTCCAGCGCCTGAATAATCAGATACACCGCGAAGCAGCTCCAGAATTCAGCACCTGCGCCAAACTGGAATa                                                                       <  1:2450031/93‑1 (MQ=255)

                 cacCGGCACTAACAGGTTGCCGTCCAGCGCCTGAATAATCAGATACACCGCGAAGCAGCTCCAGAATTCAGCACCTGCGCCAAACTGGAATAg                                                                      <  2:2335694/93‑1 (MQ=255)

                  acCGGCACTAACAGGTTGCCGTCCAGCGCCTGAATAATCAGATACACCGCGAAGCAGCTCCAGAATTCAGCACCTGCGCCAAACTGGAATAgc                                                                     >  1:3661995/1‑93 (MQ=255)

                     ggCACTAACAGGTTGCCGTCCAGCGCCTGAATAATCAGATACACCGCGAAGCAGCTCCAGAATTCAGCACCTGCGCCAAACTGGAATAGcgcc                                                                  <  1:1956096/93‑1 (MQ=255)

                      gCACTAACAGGTTGCCGTCCAGCGCCTGAATAATCAGATACACCGCGAAGCAGCTCCAGAATTCAGCACCTGCGCCAAACTGGAATAGcgcca                                                                 <  1:4233992/93‑1 (MQ=255)

                       cACTAACAGGTTGCCGTCCAGCGCCTGAATAATCAGATACACCGCGAAGCAGCTCCAGAATTCAGCACCTGCGCCAAACTGGAATAGcgccac                                                                <  1:3206993/93‑1 (MQ=255)

                       cACTAACAGGTTGCCGTCCAGCGCCTGAATAATCAGATACACCGCGAAGCAGCTCCAGAATTCAGCACCTGCGCCAAACTGGAATAGcgccac                                                                <  2:3471955/93‑1 (MQ=255)

                          tAACAGGTTGCCGTCCAGCGCCTGAATAATCAGATACACCGCGAAGCAGCTCCAGAATTCAGCACCTGCGCCAAACTGGAATAGcgccacgcc                                                             >  2:2908269/1‑93 (MQ=255)

                           aaCAGGTTGCCGTCCAGCGCCTGAATAATCAGATACACCGCGAAGCAGCTCCAGAATTCAGCACCTGCGCCAAACTGGAATAGcgccacgcca                                                            >  2:3094944/1‑93 (MQ=255)

                            aCAGGTTGCCGTCCAGCGCCTGAATAATCAGATACACCGCGAAGCAGCTCCAGAATTCAGCACCTGCGCCAAACTGGAATAGagccacgccac                                                           >  2:538199/1‑92 (MQ=255)

                            aCAGGTTGCCGTACAGCGCCTGAATAAACAGATACTCCGCGAAGCAGCACCAGAATTCAGCAACTGCGCCAAACTGGAATAGCGCCACGCCaa                                                           <  2:4358835/93‑1 (MQ=255)

                             cAGGTTGCCGTCCAGCGCCTGAATAATCAGATACACCGCGAAGCAGCTCCAGAATTCAGCACCTGCGCCAAACTGGAATAGCGCCACGCCAac                                                          <  2:3261196/93‑1 (MQ=255)

                             aaGGTTGCCGTCCAGCGCCTGAATAATCAGATACACCGCGAAGCAGCTCCAGAATTCAGCACCTGCGCCAAACTGGAATAGCGCCACGCCAac                                                          <  2:1323153/92‑1 (MQ=255)

                              aGGTTGCCGTCCAGCGCCTGAATAATCAGATACACCGCGAAGCAGCTCCAGAATTCAGCACCTGCGCCAAACTGGAATAGCGCCACGCCAacc                                                         <  1:712865/93‑1 (MQ=255)

                              aGGTTGCCGTCCAGCGCCTGAATAATCAGATACACCGCGAAGCAGCTCCAGAATTCAGCACCTGCGCCAAACTGGAATAGCGCCACGCCAacc                                                         >  2:975395/1‑93 (MQ=255)

                               ggTTGCCGTCCAGCGCCTGAATAATCAGATACACCGCGAAGCAGCTCCAGAATTCAGCACCTGCGCCAAACTGGAATAGCGCCACGCCAacca                                                        >  1:4247251/1‑93 (MQ=255)

                                gTTGCCGTCCAGCGCCTGAATAATCAGATACACCGCGAAGCAGCTCCAGAATTCAGCACCTGCGCCAAACTGGAATAGCGCCACGCCAaccac                                                       <  1:2268661/93‑1 (MQ=255)

                                 ttGCCGTCCAGCGCCTGAATAATCAGATACACCGCGAAGCAGCTCCAGAATTCAGCACCTGCGCCAAACTGGAATAGCGCCACGCCAaccacc                                                      <  1:1138346/93‑1 (MQ=255)

                                 ttGCCGTCCAGCGCCTGAATAATCAGATACACCGCGAAGCAGCTCCAGAATTCAGCACCTGCGCCAAACTGGAATAGCGCCACGCCAaccacc                                                      <  1:1314232/93‑1 (MQ=255)

                                 ttGCCGTCCAGCGCCTGAATAATCAGATACACCGCGAAGCAGCTCCAGAATTCAGCACCTGCGCCAAACTGGAATAGCGCCACGCCAaccacc                                                      <  2:3403625/93‑1 (MQ=255)

                                  tcccGTCCAGCGCCTGAATAATCAGATACACCGCGAAGCAGCTCCAGAATTCAGCACCTGCGCCAAACTGGAATAGCGCCACGCCAACCACCg                                                     <  2:3891426/91‑1 (MQ=255)

                                  tGCCGTCCAGCGCCTGAATAATCAGATACACCGCGAAGCAGCTCCAGAATTCAGCACCTGCGCCAAACTGGAATAGCGCCACGCCAACCACCg                                                     <  1:883537/93‑1 (MQ=255)

                                    ccGTCCAGCGCCTGAATAATCAGATACACCGCGAAGCAGCTCCAGAATTCAGCACCTGCGCCAAACTGGAATAGCGCCACGCCAACCACCGGa                                                   <  2:154619/93‑1 (MQ=255)

                                      gTCCAGCGCCTGAATAATCAGATACACCGCGAAGCAGCTCCAGAATTCAGCACCTGCGCCACACTGGAATAGCGCCACGCCAACCACCGGAAt                                                 >  1:1462891/1‑93 (MQ=255)

                                      gTCCAGCGCCTGAATAATCAGATACACCGCGAAGCAGCTCCAGAATTCAGCACCTGCGCCAAACTGGAATAGCGCCACGCCAACCACCGGAAt                                                 <  2:1089909/93‑1 (MQ=255)

                                      gTCCAGCGCCTGAATAATCAGATACACCGCGAAGCAGCTCCAGAATTCAGCACCTGCGCCAAACTGGAATAGCGCCACGCCAACCACCGGAAt                                                 >  1:982035/1‑93 (MQ=255)

                                       gcgaGCGCCTGAATAAGCAGATACAGCGCGAAGCAGCTCCAGAATTCAGCACCTGCGCCAAACTGGAATAGCGCCACGCCAACCACCGGAATg                                                <  2:4061092/90‑1 (MQ=255)

                                        ccAGCGCCTGAATAATCAGATACACCGCGAAGCAGCTCCAGAATTCAGCACCTGCGCCAAACTGGAATAGCGCCACGCCAACCACCGGAATgg                                               >  2:4337263/1‑93 (MQ=255)

                                         cAGCGCCTGAATAATCAGATACACCGCGAAGCAGCTCCAGAATTCAGCACCTGCGCCAAACTGGAATAGCGCCACGCCAACCACCGGAATGGt                                              <  1:791009/93‑1 (MQ=255)

                                          aGCGCCTGAATAATCAGATACACCGCGAAGCAGCTCCAGAATTCAGCACCTGCGCCAAACTGGAATAGCGCCACGCCAACCACCGGAATGGTc                                             >  2:1498585/1‑93 (MQ=255)

                                          aGCGCCTGAATAATCAGATACACCGCGAAGCAGCTCCAGAATTCAGCACCTGCGCCAAACTGGAATAGCGCCACGCCAACCACCGGAATGGTc                                             >  1:3190243/1‑93 (MQ=255)

                                           gcgcCTGAATAATCAGATACACCGCGAAGCAGCTCCAGAATTCAGCACCTGCGCCAAACTGGAATAGCGCCACGCCAACCACCGGAATGGTca                                            >  1:3998967/1‑93 (MQ=255)

                                            cgcCTGAATAATCAGATACACCGCGAAGCAGCTCCAGAATTCAGCACCTGCGCCAAACTGGAATAGCGCCACGCCAACCACCGGAATGGTcac                                           <  1:3483747/93‑1 (MQ=255)

                                             gcCTGAATAATCAGATACACCGCGAAGCAGCTCCAGAATTCAGCACCTGCGCCAAACTGGAATAGCGCCACGCCAACCACCGGAATGGTcacc                                          >  1:1167006/1‑93 (MQ=255)

                                             gcCTGAATAATCAGATACACCGCGAAGCAGCTCCAGAATTCAGCACCTGCGCCAAACTGGAATAGCGCCACGCCAACCACCGGAATGGTcacc                                          >  1:1141248/1‑93 (MQ=255)

                                              ccTGAATAATCAGATACACCGCGAAGCAGCTCCAGAATTCAGCACCTGCGCCAAACTGGAATAGCGCCACGCCAACCACCGGAATGGTcacca                                         >  2:3968983/1‑93 (MQ=255)

                                              ccTGAATAATCAGATACACCGCGAAGCAGCTCCAGAATTCAGCACCTGCGCCAAACTGGAATAGCGCCACGCCAACCACCGGAATGGTcacca                                         >  2:3968984/1‑93 (MQ=255)

                                                  aataatCAGATACACCGCGAAGCAGCTCCATAATTCAGCACCTGCGCCAAACTGGAATAGCGCCACGCCAACCACCGGAATGGTCACCACaaa                                     <  1:3251198/93‑1 (MQ=255)

                                                  aataatCAGATACACCGCGAAGCAGCTCCAGAATTCAGCACCTGCGCCAAACTGGAATAGCGCCACGCCAACCACCGGAATGGTCACCACaaa                                     <  2:3600530/93‑1 (MQ=255)

                                                   ataatCAGATACACCGCGAAGCAGCTCCAGAATTCAGCACCTGCGCCAAACTGGAATAGCGCCACGCCAACCACCGGAATGGTCACCACAAAt                                    >  2:2413051/1‑93 (MQ=255)

                                                     aatCAGATACACCGCGAAGCAGCTCCAGAATTCAGCACCTGCGCCAAACTGGAATAGCGCCACGCCAACCACCGGAATGGTCACCACAAATgc                                  >  2:1292406/1‑93 (MQ=255)

                                                        cAGATACACCGCGAAGCAGCTCCAGAATTCAGCACCTGCGCCAAACTGTAATAGCGCCACGCCAACCACCGGAATGGTCACCACAAATGCGcc                               <  2:484440/93‑1 (MQ=255)

                                                        cAGATACACCGCGAAGCAGCTCCAGAATTCAGCACCTGCGCCAAACTGGAATAGCGCCACGCCAACCACCGGAATGGTCACCACAAATGCGcc                               <  2:694873/93‑1 (MQ=255)

                                                        cAGATACACCGCGAAGCAGCTCCAGAATTCAGCACCTGCGCCAAACTGGAATAGCGCCACGCCAACCACCGGAATGGTCACCACAAATGCGcc                               <  1:1747082/93‑1 (MQ=255)

                                                          gATACACCGCGAAGCAGCTCCAGAATTCAGCACCTGCGCCAAACTGGAATAGCGCCACGCCAACCACCGGAATGGTCACCACAAATGCGCCaa                             <  2:1111455/93‑1 (MQ=255)

                                                           aTACACCGCGAAGCAGCTACAGAATTCAGCACCTGCGCCAAACTGGAATAGCGCCACGCCAACCACCGGAATGGTCACCACAAATGCGCCAAt                            <  1:2774573/93‑1 (MQ=255)

                                                            tACACCGCGAAGCAGCTCCAGAATTCAGCACCTGCGCCAAACTGGAATAGCGCCACGCCAACCACCGGAATGGTCACCACAAATGCGCCAATg                           >  1:1334609/1‑93 (MQ=255)

                                                             acacCGCTAAGCAGCTCCAGAATTCAGCACCTGCGCCAAACTGGAATAGCGCCACGCCAACCACCGGAATGGTCACCACAAATGCGCCAATGt                          >  2:4747165/1‑93 (MQ=255)

                                                             acacCGCGAAGCAGCTCCAGAATTCAGCACCTGCGCCAAACTGGAATAGCGCCACGCCAACCACCGGAATGGTCACCACAAATGCGCCAATGt                          >  2:939734/1‑93 (MQ=255)

                                                             acacCGCGAAGCAGCTCCAGAATTCAGCACCTGCGCCAAACTGGAATAGCGCCACGCCAACCACCGGAATGGTCACCACAAATGCGCCAATGt                          <  1:3733803/93‑1 (MQ=255)

                                                               acCGCGAAGCAGCTCCAGAATTCAGCACCTGCGCCAAACTGGAATAGCGCCACGCCAACCACCGGAATGGTCACCACAAATGCGCCAATGTAc                        <  1:2119191/93‑1 (MQ=255)

                                                               acCGCGAAGCAGCTCCAGAATTCAGCACCTGCGCCAAACTGGAATAGCGCCACGCCAACCACCGGAATGGTCACCACAAATGCGCCAATGTAc                        >  1:4559035/1‑93 (MQ=255)

                                                               acCGCGAAGCAGCTCCAGAATTCAGCACCTGCGCCAAACTGGAATAGCGCCACGCCAACCACCGGAATGGTCACCACAAATGCGCCAATGTAc                        >  1:1387177/1‑93 (MQ=255)

                                                                 cgcgAAGCAGCTCCAGAATTCAGCACCTGCGCCAAACTGGAATAGCTCCGCGCCAACCACCGGAATGGTCACCACAAATGCGCCAATGTACgg                      <  1:3233329/93‑1 (MQ=255)

                                                                  gcgAAGCAGCTCCAGAATTCAGCACCTGCGCCAAACTGGAATAGCGCCACGCCAACCACCGGAATGGTCACCACAAATGCGCCAATGTACGGa                     >  1:4297139/1‑93 (MQ=255)

                                                                  gcgAAGCAGCTCCAGAATTCAGCACCTGCGCCAAACTGGAATAGCGCCACGCCAACCACCGGAATGGTCACCACAAATGCGCCAATGTACGGa                     >  1:3323155/1‑93 (MQ=255)

                                                                   cgAAGCAGCTCCAGAATTCAGCACCTGCGCCAAACTGGAATAGCTCCACGCCAACCACCGGAATGGTCACCACAAATGCGCCAATGTACGGaa                    <  1:1981442/93‑1 (MQ=255)

                                                                    gAAGCAGCTCCAGAATTCAGCACCTGCGCCAAACTGGAATAGCGCCACGCCAACCACCGGAATGGTCACCACAAATGCGCCAATGTACGGAAt                   >  2:1200804/1‑93 (MQ=255)

                                                                       gcagcTCCAGAATTCAGCACCTGCGCCAAACTGGAATAGCGCCACGCCAACCACCGGAATGGTCACCACAAATGCGCCAATGTACGGAATAAg                <  1:3013172/93‑1 (MQ=255)

                                                                       gcagcTCCAGAATTCAGCACCTGCGCCAAACTGGAATAGCGCCACGCCAACCACCGGAATGGTCACCACAAATGCGCCAATGTACGCAATAAg                <  1:225240/93‑1 (MQ=255)

                                                                         agcTCCAGAATTCAGCACCTGCGCCAAACTGGAATAGCGCCACGCCAACCACCGGAATGGTCACCACAAATGCGCCAATGTACGGAATAAGaa              >  2:811728/1‑93 (MQ=255)

                                                                          gcTCCAGAATTCAGCACCTGCGCCAAACTGGAATAGCGCCACGCCAACCACCGGAATGGTCACCACAAATGCGCCAATGTACGGAATAAGAAc             <  2:2043780/93‑1 (MQ=255)

                                                                          gcTCCAGAATTCAGCACCTGCGCCAAACTGGAATAGCGCCACGCCAACCACCGGAATGGTCACCACAAATGCGCCAATGTACGGAATAAGAAc             <  1:4027859/93‑1 (MQ=255)

                                                                           cTCCAGAATTCAGCACCTGCGCCAAACTGGAATAGCGCCACGCCAACCACCGGAATGGTCACCACAAATGCGCCAATGTACGGAATAAGAAcc            >  1:3494946/1‑93 (MQ=255)

                                                                           cTCCAGAATTCAGCACCTGCGCCAAACTGGAATAGCGCCACGCCAACCACCGGAATGGTCACCACAAATGCGCCAATGTACGGAATAAGAAcc            <  1:2133125/93‑1 (MQ=255)

                                                                            tCCAGAATTCAGCACCTGCGCCAAACTGGAATAGCGCCACGCCAACCACCGGAATGGTCACCACAAATGCGCCAATGTACGGAATAAGAACCg           >  1:2790874/1‑93 (MQ=255)

                                                                             ccAGAATTCAGCACCTGCGCCAAACTGGAATAGCGCCACGCCAACCACCGGAATGGTCACCACAAATGCGCCAATGTACGGAATAAGAACCga          >  1:3191252/1‑93 (MQ=255)

                                                                             ccAGAATTCAGCACCTGCGCCAAACTGGAATAGCGCCACGCCAACCACCGGAATGGTCACCACAAATGCGCCAATGTACGGAATAAGAACCga          <  2:2627126/93‑1 (MQ=255)

                                                                             ccAGAATTCAGCACCTGCGCCAAACTGGAATAGCGCCACGCCAACCACCGGAATGGTCACCACAAATGCGCCAATGTACGGAATAAGAACCga          <  2:4426030/93‑1 (MQ=255)

                                                                              cAGAATTCAGCACCTGCGCCAAACTGGAATAGCGCCACGCCAACCACCGGAATGGTCACCACAAATGCGCCAATGTACGGAATAAGAACCgag         >  1:849574/1‑93 (MQ=255)

                                                                               aGAATTCAGCACCTGCGCCAAACTGGAATAGCGCCACGCCAACCACCGGAATGGTCACCACAAATGCGCCAATGTACGGAATAAGAACCgaga        >  1:860093/1‑93 (MQ=255)

                                                                               aGAATTCAGCACCTGCGCCAAACTGGAATAGCGCCACGCCAACCACCGGAATGGTCACCACAAATGCGCCAATGTACGGAATAAGAACCgaga        <  1:3242000/93‑1 (MQ=255)

                                                                               aGAATTCAGCACCTGCGCCAAACTGGAATAGCGCCACGCCAACCACCGGAATGGTCACCACAAATGCGCCAATGTACGGAATAAGAACCgaga        <  1:2269277/93‑1 (MQ=255)

                                                                                gAATTCAGCACCTGCGCCAAACTGGAATAGCGCCACGCCAACCACCGGAATGGTCACCACAAATGCGCCAATGTACGGAATAAGAACCGAGaa       >  1:4572472/1‑93 (MQ=255)

                                                                                  aTTCAGCACCTGCGCCAAACTGGAATAGCGCCACGCCAACCACCGGAATGGTCACCACAAATGCGCCAATGTACGGAATAAGAACCGAGAAGc     <  2:974234/93‑1 (MQ=255)

                                                                                   ttCAGCACCTGCGCCAAACTGGAATAGCGCCACGCCAACCACCGGAATGGTCACCACAAATGCGCCAATGTACGGAATAAGAACCGAGAAGcc    >  1:984472/1‑93 (MQ=255)

                                                                                    tCAGTACCTGCGCCAAACTGGAATAGCGCCACGCCAACCACCGGAATGGTCACCACAAATGCGCCAATGTACGGAATAAGAACCGAGAAGCCg   >  2:4112638/1‑93 (MQ=255)

                                                                                     cAGCACCTGCGCCAAACTGGAATAGCGCCACGCCAACCACCGGAATGGTCACCACAAATGCGCCAATGTACGGAATAAGAACCGAGAAGCCGa  <  2:3813010/93‑1 (MQ=255)

                                                                                     cAGCACCTGCGCCAAACTGGAATAGCGCCACGCCAACCACCGGAATGGTCACCACAAATGCGCCAATGTACGGAATAAGAACCGAGAAGCCGa  <  1:1089152/93‑1 (MQ=255)

ACCGCTTCGGAGAACAACACCGGCACTAACAGGTTGCCGTCCAGCGCCTGAATAATCAGATACACCGCGAAGCAGCTCCAGAATTCAGTACCTGCGCCAAACTGGAATAGCGCCACGCCAACCACCGGAATGGTCACCACAAATGCGCCAATGTACGGAATAAGAACCGAGAAGCCGA  >  NC_000913/2612988‑2613165

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