BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,354,595	A→G	L534P (CTG→CCG)	cadA ←	lysine decarboxylase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,354,595	0	A	G	85.7%	627.8 / 97.5	266	L534P (CTG→CCG)	cadA	lysine decarboxylase 1
Reads supporting (aligned to +/- strand): ref base A (18/20); new base G (110/118); total (128/138)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.01e-01

Only 100 of 323 total aligned reads displayed.

GCAGGTTCAGGTCGAACGCACGTTTAAAGTCAGTCAGAGCACGCAGCAGGCTCAGTGCTTTGGTCTTATCGATACCGATGCTGAACAGGAACAGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCACGATGCTGGCCGGAATACCAAAGTCGCTCATG  >  NC_000913/4354503‑4354687

  aGGTTCAGGTCGAACGCACGTTTAAAGTCAGTCAGAGCACGCGGCAGGCTCAGTGCTTTGGTCTTATCGATACCGATGCTGAACAGGAAcggc                                                                                            >  1:119488/1‑93 (MQ=255)

  aGGTTCAGGTCGAACGCACGTTTAAAGTCAGTCAGAGCACGCAGCAGGCTCAGTGCTTTGGTCTTATCGATACCGATGCTGAACAGGAAcggc                                                                                            <  2:6554518/93‑1 (MQ=255)

  aGGTTCAGGTCGAACGCACGTTTAAAGTCAGTCAGAGCACGCAGCAGGCTCAGTGCTTTGGTCTTATCGATACCGATGCTGAACAGGAAcggc                                                                                            >  2:4933804/1‑93 (MQ=255)

   ggTTCAGGTCGAACGCACGTTTAAAGTCAGTCAGAGCACGCAGCAGGCTCAGTGCTTTGGTCTTATCGATACCGATGCTGAACAGGAAcggca                                                                                           >  2:3088649/1‑93 (MQ=255)

   ggTTCAGGTCGAACGCACGTTTAAAGTCAGTCAGAGCACGCAGCAGGCTCAGTGCTTTGGTCTTATCGATACCGATGCTGAACAGGAAcggca                                                                                           >  2:3499346/1‑93 (MQ=255)

   ggTTCAGGTCGAACGCACGTTTAAAGTCAGTCAGAGCACGCAGCAGGCTCAGTGCTTTGGTCTTATCGATACCGATGCTGAACAGGAAcggca                                                                                           <  1:5939738/93‑1 (MQ=255)

   ggTTCAGGTCGAACGCACGTTTAAAGTCAGTCAGAGCACGCAGCAGGCTCAGTGCTTTGGTCTTATCGATACCGATGCTGAACAGGAAcggca                                                                                           <  1:3325725/93‑1 (MQ=255)

   ggTTCAGGTCGAACGCACGTTTAAAGTCAGTCAGAGCACGCAGCAGGCTCAGTGCTTTGGTCTTATCGATACCGATGCTGAACAGGAAcggca                                                                                           >  1:6421252/1‑93 (MQ=255)

    gTTCAGGTCGAACGCACGTTTAAAGTCAGTCAGAGCACGCAGCAGGCTCAGTGCTTTGGTCTTATCGATACCGATGCTGAACAGGAAcggcag                                                                                          >  1:6112948/1‑93 (MQ=255)

    gTTCAGGTCGAACGCACGTTTAAAGTCAGTCAGAGCACGCAGCAGGCTCAGTGCTTTGGTCTTATCGATACCGATGCTGAACAGGAAcggcag                                                                                          <  2:5066711/93‑1 (MQ=255)

     ttCAGGTCGAACGCACGTTTAAAGTCAGTCAGAGCACGCAGCAGGCTCAGTGCTTTGGTCTTATCGATACCGATGCTGAACAGGAACGGCAgg                                                                                         <  1:6828877/93‑1 (MQ=255)

         ggTCGAACGCACGTTTAAAGTCAGTCAGAGCACGCAGCAGGCTCAGTGCTTTGGTCTTATCGATACCGATGCTGAACAGGAACGGCAGGTtat                                                                                     >  2:7197434/1‑93 (MQ=255)

         ggTCGAACGCACGTTTAAAGTCAGTCAGAGCACGCAGCAGGCTCAGTGCTTTGGTCTTATCGATACCGATGCTGAACAGGAACAGCAGGTtat                                                                                     <  1:2765939/93‑1 (MQ=255)

          gTCGAACGCACGTTTAAAGTCAGTCAGAGCACGCAGCAGGCTCAGTGCTTTGGTCTTATCGATACCGATGCTGAACAGGAACGGCAGGTtata                                                                                    <  2:6320287/93‑1 (MQ=255)

          gTCGAACGCACGTTTAAAGTCAGTCAGAGCACGCAGCAGGCTCAGTGCTTTGGTCTTATCGATACCGATGCTGAACAGGAACGGCAGGTtata                                                                                    <  2:5919094/93‑1 (MQ=255)

          gTCGAACGCACGTTTAAAGTCAGTCAGAGCACGCAGCAGGCTCAGTGCTTTGGTCTTATCGATACCGATGCTGAACAGGAACAGCAGGTtata                                                                                    <  2:1863785/93‑1 (MQ=255)

          gTCGAACGCACGTTTAAAGTCAGTCAGAGCACGCAGCAGGCTCAGTGCTTTGGTCTTATCGATACCGATGCTGAACAGGAACAGCAGGTtata                                                                                    <  2:552802/93‑1 (MQ=255)

             gAACGCACGTTTAAAGTCAGTCAGAGCACGCAGCAGGCTCAGTGCTTTGGTCTTATCGATACCGATGCTGAACAGGAACAGCAGGTTATACgg                                                                                 <  1:4962936/93‑1 (MQ=255)

              aaCGCACGTTTAAAGTCAGTCAGAGCACGCAGCAGGCTCAGTGCTTTGGTCTTATCGATACCGATGCTGAACAGGAACGGCAGGTTATACGGa                                                                                <  2:5223352/93‑1 (MQ=255)

              aaCGCACGTTTAAAGTCAGTCAGAGCACGCAGCAGGCTCAGTGCTTTGGTCTTATCGATACCGATGCTGAACAGGAACGGCAGGTTATACGGa                                                                                >  2:1949371/1‑93 (MQ=255)

              aaCGCACGTTTAAAGTCAGTCAGAGCACGCAGCAGGCTCAGTGCTTTGGTCTTATCGATACCGATGCTGAACAGGAACGGCAGGTTATACGGa                                                                                <  1:2832737/93‑1 (MQ=255)

               aCGCACGTTTAAAGTCAGTCAGAGCACGCAGCAGGCTCAGTGCTTTGGTCTTATCGATACCGATGCTGAACAGGAACGGCAGGTTATACGGAc                                                                               >  1:5864471/1‑93 (MQ=255)

               aCGCACGTTTAAAGTCAGTCAGAGCACGCAGCAGGCTCAGTGCTTTGGTCTTATCGATACCGATGCTGAACAGGAACGGCAGGTTATACGGAc                                                                               >  2:1526883/1‑93 (MQ=255)

                   aCGTTTAAAGTCAGTCAGAGCACGCAGCAGGCTCAGTGCTTTGGTCTTATCGATACCGATGCTGAACAGGAACGGCAGGTTATACGGACCGGt                                                                           <  2:5679438/93‑1 (MQ=255)

                     gTTTAAAGTCAGTCAGAGCACGCAGCAGGCTCAGTGCTTTGGTCTTATCGATACCGATGCTGAACAGGAACGGCAGGTTATACGGACCGGttt                                                                         >  1:1287371/1‑93 (MQ=255)

                       ttAAAGTCAGTCAGAGCACGCAGCAGGCTCAGTGCTTTGGTCTTATCGATACCGATGCTGAACAGGAACGGCAGGTTATACGGACCGGTTTtc                                                                       <  2:2722087/93‑1 (MQ=255)

                       ttAAAGTCAGTCAGAGCACGCAGCAGGCTCAGTGCTTTGGTCTTATCGATACCGATGCTGAACAGGAACGGCAGGTTATACGGACCGGTTTtc                                                                       <  1:533719/93‑1 (MQ=255)

                       ttAAAGTCAGTCAGAGCACGCAGCAGGCTCAGTGCTTTGGTCTTATCGATACCGATGCTGAACAGGAACGGCAGGTTATACGGACCGGTTTtc                                                                       <  1:5603436/93‑1 (MQ=255)

                         aaaGTCAGTCAGAGCACGCAGCAGGCTCAGTGCTTTGGTCTTATCGATACCGATGCTGAACAGGAACGGCAGGTTATACGGACCGGTTTtctc                                                                     >  1:2327900/1‑93 (MQ=255)

                         aaaGTCAGTCAGAGCACGCAGCAGGCTCAGTGCTTTGGTCTTATCGATACCGATGCTGAACAGGAACAGCAGGTTATACGGACCGGTTTtctc                                                                     <  1:4411701/93‑1 (MQ=255)

                           agtcagtcAGAGCACGCAGCAGGCTCAGTGCTTTGGTCTTATCGATACCGATGCTGAACAGGAACGGCAGGTTATACGGACCGGTTTTCTcaa                                                                   >  1:985044/1‑93 (MQ=255)

                           agtcagtcAGAGCACGCAGCAGGCTCAGTGCTTTGGTCTTATCGATACCGATGCTGAACAGGAACAGCAGGTTATACGGACCGGTTTTCTcaa                                                                   >  1:7469667/1‑93 (MQ=255)

                           agtcagtcAGAGCACGCAGCAGGCTCAGTGCTTTGGTCTTATCGATACCGATGCTGAACAGGAACAGCAGGTTATACGGACCGGTTTTCTcaa                                                                   <  2:3214105/93‑1 (MQ=255)

                            gtcagtcaGAGCACGCAGCAGGCTCAGTGCTTTGGTCTTATCGATACCGATGCTGAACAGGAACGGCAGGTTATACGGACCGGTTTTCTcaac                                                                  >  1:470945/1‑93 (MQ=255)

                            gtcagtcaGAGCACGCAGCAGGCTCAGTGCTTTGGTCTTATCGATACCGATGCTGAACAGGAACGGCAGGTTATACGGACCGGTTTTCTcaac                                                                  >  2:5120784/1‑93 (MQ=255)

                             tcagtcagAGCACGCAGCAGGCTCAGTGCTTTGGTCTTATCGATACCGATGCTGAACAGGAACGGCAGGTTATACGGACCGGTTTTCTcaaca                                                                 >  1:2663730/1‑93 (MQ=255)

                              cagtcagAGCACGCAGCAGGCTCAGTGCTTTGGTCTTATCGATACCGATGCTGAACAGGAACGGCAGGTTATACGGACCGGTTTTCTcaacaa                                                                <  1:1673212/93‑1 (MQ=255)

                                  cagAGCACGCAGCAGGCTCAGTGCTTTGGTCTTATCGATACCGATGCTGAACAGGAACAGCAGGTTATACGGACCGGTTTTCTCAACAACGAt                                                            <  1:6303360/93‑1 (MQ=255)

                                   agagCACGCAGCAGGCTCAGTGCTTTGGTCTTATCGATACCGATGCTGAACAGGAACGGCAGGTTATACGGACCGGTTTTCTCAACAACGATg                                                           <  1:4533480/93‑1 (MQ=255)

                                    gagCACGCAGCAGGCTCAGTGCTTTGGTCTTATCGATACCGATGCTGAACAGGAACGGCAGGTTATACGGACCGGTTTTCTCAACAACGATGc                                                          <  2:1281150/93‑1 (MQ=255)

                                        aCGCAGCAGGCTCAGTGCTTTGGTCTTATCGATACCGATGCTGAACAGGAACGGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATg                                                      >  2:2621431/1‑93 (MQ=255)

                                        aCGCAGCAGGCTCAGTGCTTTGGTCTTATCGATACCGATGCTGAACAGGAACGGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATg                                                      >  2:4071151/1‑93 (MQ=255)

                                        aCGCAGCAGGCTCAGTGCTTTGGTCTTATCGATACCGATGCTGAACAGGAACGGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATg                                                      <  2:5421014/93‑1 (MQ=255)

                                          gcagcaGGCTCAGTGCTTTGGTCTTATCGATACCGATGCTGAACAGGAACGGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGtt                                                    >  2:332768/1‑93 (MQ=255)

                                          gcagcaGGCTCAGTGCTTTGGTCTTATCGATACCGATGCTGAACAGGAACGGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGtt                                                    >  1:824587/1‑93 (MQ=255)

                                           cagcagGCTCAGTGCTTTGGTCTTATCGATACCGATGCTGAACCTGAACGGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTtc                                                   <  1:7305332/93‑1 (MQ=255)

                                              cagGCTCAGTGCTTTGGTCTTATCGATACCGATGCTGAACAGGAACAGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTtcgtc                                                <  1:5484176/93‑1 (MQ=255)

                                                ggCTCAGTGCTTTGGTCTTATCGATACCGATGCTGAACAGGAACGGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGa                                              <  2:5190623/93‑1 (MQ=255)

                                                ggCTCAGTGCTTTGGTCTTATCGATACCGATGCTGAACAGGAACGGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGa                                              >  2:6172439/1‑93 (MQ=255)

                                                  cTCAGTGCTTTGGTCTTATCGATACCGATGCTGAACAGGAACGGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAgg                                            <  1:2547951/93‑1 (MQ=255)

                                                  cTCAGTGCTTTGGTCTTATCGATACCGATGCTGAACAGGAACGGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAgg                                            >  1:70973/1‑93 (MQ=255)

                                                    cAGTGCTTTGGTCTTATCGATACCGATGCTGAACAGGAACAGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTa                                          <  1:773272/93‑1 (MQ=255)

                                                     aGTGCTTTGGTCTTATCGATACCGATGCTGAACAGGAACGGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTAt                                         >  2:3587282/1‑93 (MQ=255)

                                                       tGCTTTGGTCTTATCGATACCGATGCTGAACAGGAACGGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTAttt                                       <  2:6905140/93‑1 (MQ=255)

                                                       tGCTTTGGTCTTATCGATACCGATGCTGAACAGGAACAGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTAttt                                       >  2:1139058/1‑93 (MQ=255)

                                                       tGCTTTGGTCTTATCGATACCGATGCTGAACAGGAACAGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTAttt                                       <  2:4190217/93‑1 (MQ=255)

                                                        gCTTTGGTCTTATCGATACCGATGCTGAACAGGAACGGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTc                                      <  2:7041559/93‑1 (MQ=255)

                                                         cTTTGGTCTTATCGATACCGATGCTGAACAGGAACGGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCg                                     >  1:4832438/1‑93 (MQ=255)

                                                          tttGGTCTTATCGATACCGATGCTGAACAGGAACGGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGc                                    >  1:4131130/1‑93 (MQ=255)

                                                           ttGGTCTTATCGATACCGATGCTGAACAGGAACGGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGcc                                   >  2:7519014/1‑93 (MQ=255)

                                                           ttGGTCTTATCGATACCGATGCTGAACAGGAACGGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGcc                                   >  2:3773480/1‑93 (MQ=255)

                                                           ttGGTCTTATCGATACCGATGCTGAACAGGAACGGCAGGTTATACGGACCGGTTTTCTCAACAACGACGCCATGTTCGTCGAGGTATTTCGcc                                   >  1:3900880/1‑93 (MQ=255)

                                                             ggTCTTATCGATACCGATGCTGAACAGGAACGGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCAc                                 <  2:2478935/93‑1 (MQ=255)

                                                             ggTCTTATCGATACCGATGCTGAACAGGAACGGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGGATTTCGCCAc                                 >  1:1841485/1‑93 (MQ=255)

                                                              gTCTTATCGATACCGATGCTGAACAGGAACGGCAGGTTATACGGACCGGTTTTTTCAACAACGATGCCATGTTCGTCGAGGGATTTCGCCACg                                >  1:3689151/1‑93 (MQ=255)

                                                               tCTTATCGATACCGATGCTGAACAGGAACGGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCACGa                               <  1:3780859/93‑1 (MQ=255)

                                                               tCTTATCGATACCGATGCTGAACAGGAACGGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCACGa                               <  1:6949561/93‑1 (MQ=255)

                                                               tCTTATCGATACCGATGCTGAACAGGAACAGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCACGa                               <  1:3052201/93‑1 (MQ=255)

                                                                cTTATCGATACCGATGCTGAACAGGAACGGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCACGAt                              <  1:6720256/93‑1 (MQ=255)

                                                                cTTATCGATACCGATGCTGAACAGGAACGGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCACGAt                              <  1:3237805/93‑1 (MQ=255)

                                                                  tATCGATACCGATGCTGAACAGGAACGGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCACGATGc                            >  1:5796389/1‑93 (MQ=255)

                                                                   aTCGATACCGATGCTGAACAGGAACGGCAGGTTATACGGACCGGTTTTTTCAAAAACGATGCCAGGTTCGCCAAGGAATTTCCCCACGATCCt                           >  2:2865174/1‑93 (MQ=255)

                                                                   aTCGATACCGATGCTGAACAGGAACGGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCACGATGCt                           <  2:6391827/93‑1 (MQ=255)

                                                                    tCGATACCGATGCTGAACAGGAACGGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCACGATGCTg                          >  1:4209414/1‑93 (MQ=255)

                                                                    tCGATACCGATGCTGAACAGGAACGGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCACGATGCTg                          >  1:1014319/1‑93 (MQ=255)

                                                                    tCGATACCGATGCTGAACAGGAACGGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCACGATGCTg                          >  2:4719596/1‑93 (MQ=255)

                                                                     cGATACCGATGCTGAACAGGAACGGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCACGATGCTgg                         <  2:4425818/93‑1 (MQ=255)

                                                                      gATACCGATGCCGAACAGGAACGGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCACGATGCTGGc                        <  2:4588016/93‑1 (MQ=255)

                                                                       aTACCGATGCTGAACAGGAACGGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCACGATGCTGGcc                       <  1:1605252/93‑1 (MQ=255)

                                                                       aTACCGATGCTGAACAGGAACGGCAGGTTATACGGACCGGTTTTCCCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCACGATGCTGGcc                       >  1:1367271/1‑93 (MQ=255)

                                                                          ccGATGCTGAACAGGAACGGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCACGATGCTGGCCGGa                    >  1:394350/1‑93 (MQ=255)

                                                                          ccGATGCTGAACAGGAACGGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCACGATGCTGGCCGGa                    >  1:1473018/1‑93 (MQ=255)

                                                                            gATGCTGAACAGGAACGGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCACGATGCTGGCCGGAAt                  <  1:5268864/93‑1 (MQ=255)

                                                                            gATGCTGAACAGGAACGGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCACGATGCTGGCCGGAAt                  >  1:5546473/1‑93 (MQ=255)

                                                                             aTGCTGAACAGGAACGGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCACGATGCTGGCCGGAATa                 >  1:5764262/1‑93 (MQ=255)

                                                                             aTGCTGAACAGGAACGGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCACGATGCTGGCCGGAATa                 >  2:2181614/1‑93 (MQ=255)

                                                                             aTGCTGAACAGGAACGGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCACGATGCTGGCCGGAATa                 >  2:2267657/1‑93 (MQ=255)

                                                                                cTGAACAGGAACGGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCACGATGCTGGCCGGAATACCa              <  1:431438/93‑1 (MQ=255)

                                                                                  gaacaggaacggCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCACGATGCTGGCCGGAATACCaaa            <  2:3867592/93‑1 (MQ=255)

                                                                                   aacaggaacagCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCACGATGCTGGCCGGAATACCAAAg           <  2:1768490/93‑1 (MQ=255)

                                                                                    acaggaacggCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCACGATGCTGGCCGGAATACCAAAGt          <  2:7416737/93‑1 (MQ=255)

                                                                                    acaggaacggCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCACGATGCTGGCCGGAATACCAAAGt          <  2:91082/93‑1 (MQ=255)

                                                                                    acaggaacggCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCACGATGCTGGCCGGAATACCAAAGt          >  2:4755148/1‑93 (MQ=255)

                                                                                     caggaacggCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCACGATGCTGGCCGGAATACCAAAGTc         <  2:683198/93‑1 (MQ=255)

                                                                                     caggaacggCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCACGATGCTGGCCGGAATACCAAAGTc         >  2:2821/1‑93 (MQ=255)

                                                                                      aggaacagCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCTTCGAGGTATTTCGCCACGATGCTGGCCGGAATACCAAAGTCg        >  1:7373928/1‑93 (MQ=255)

                                                                                       ggaacggCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCACGATGCTGGCCGGAATACCAAAGTCGc       >  2:6927914/1‑93 (MQ=255)

                                                                                       ggaacggCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCACGATGCTGGCCGGAATACCAAAGTCGc       <  1:1875237/93‑1 (MQ=255)

                                                                                         aacggCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCACGATGCTGGCCGGAATACCAAAGTCGCTc     <  1:287856/93‑1 (MQ=255)

                                                                                            agcagGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCACGATGCTGGCCGGAATACCAAAGTCGCTCATg  >  2:6489051/1‑93 (MQ=255)

GCAGGTTCAGGTCGAACGCACGTTTAAAGTCAGTCAGAGCACGCAGCAGGCTCAGTGCTTTGGTCTTATCGATACCGATGCTGAACAGGAACAGCAGGTTATACGGACCGGTTTTCTCAACAACGATGCCATGTTCGTCGAGGTATTTCGCCACGATGCTGGCCGGAATACCAAAGTCGCTCATG  >  NC_000913/4354503‑4354687

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