BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	AE015451_phaM1	1,499,480	2 bp→AC	intergenic (+57/‑106)	trpS CDS (Tryptophan‑‑tRNA ligase) → / → zapE CDS (putative enzyme with nucleoside triphosphate hydrolase domain)	–/–

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AE015451_phaM1	1,499,480	0	C	A	92.9%	407.7 / ‑4.9	197	intergenic (+57/‑107)	trpS CDS (Tryptophan‑‑tRNA ligase)/zapE CDS (putative enzyme with nucleoside triphosphate hydrolase domain)	–/–
Reads supporting (aligned to +/- strand): ref base C (10/4); new base A (108/74); total (118/79)
Fisher's exact test for biased strand distribution p-value = 4.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*	AE015451_phaM1	1,499,481	0	A	C	98.0%	468.5 / NA	197	intergenic (+58/‑106)	trpS CDS (Tryptophan‑‑tRNA ligase)/zapE CDS (putative enzyme with nucleoside triphosphate hydrolase domain)	–/–
Reads supporting (aligned to +/- strand): ref base A (3/1); new base C (115/77); total (118/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 198 total aligned reads displayed.

CGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGA‑CTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTG‑T‑T‑T‑CGGGC‑A‑C‑G‑C‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCG‑TTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATC  >  AE015451_phaM1/1499347‑1499608

||

cGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACCGCAGGA‑CTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGGAACAGGCGCCGCGGTGCCTGGCCCCg                                                                                                                                          >  2:15389267/1‑137 (MQ=255)

cGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGA‑CTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGccccg                                                                                                                                          >  2:12152572/1‑134 (MQ=255)

    gTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGA‑CTGAGGCTGCCAAACGGTGGGAGCGGGCGTNCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCtt                                                                                                                                      <  1:17033930/137‑1 (MQ=255)

    gTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGA‑CTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCtt                                                                                                                                      <  1:7910927/137‑1 (MQ=255)

    gTTCGACAACGCGGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGA‑CTGAGGCTGCCAAACGGTGGGAGTGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCtt                                                                                                                                      <  1:55827/137‑1 (MQ=38)

      tCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGA‑CTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTg                                                                                                                                    <  1:4279472/137‑1 (MQ=255)

      tCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGA‑CTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTg                                                                                                                                    <  1:14437767/137‑1 (MQ=255)

          cAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGA‑CTGAGGCTGCCAAACGGTGGGAGCGGGTGTGTCCGCGAAACAAGCGCCGCGGTGTATGGCACCGGCTTTGCCgg                                                                                                                                >  2:7476835/1‑137 (MQ=255)

          cAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGA‑CTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCATGAAACAGGCGCCGCGGTGCATGGCACCAGCTTTGCCgg                                                                                                                                >  2:5418683/1‑137 (MQ=255)

            aCGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGA‑CTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGgtg                                                                                                                              <  2:6474840/137‑1 (MQ=255)

             cGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGA‑CTGAGGCTGCCAAACGGTGGGAGCGGGCGTGTCCGCGAAAAAGGCGCCGCGGTGTACGGCACCGGCCTTGCCGgttt                                                                                                                             >  1:10062162/1‑135 (MQ=255)

                tGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGA‑CTGAGGCTGCCAAACGGTGGGAGCGGGTGTGCCCGCGAAACAGGCGCCGCGGTGTATGGCACCGGCTTTGCCGGTG‑t‑tcg                                                                                                                        >  2:2245373/1‑135 (MQ=38)

                tGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGA‑CTGAGGCTGCCAAACGGTGGGAGCGGGCGTTCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑t‑tcg                                                                                                                        >  1:16631698/1‑135 (MQ=255)

                tGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGA‑CTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAAAAGCCACGGTTCCTGGGCCCCGATTTTCCGGTG‑t‑tcg                                                                                                                        >  1:4120464/1‑135 (MQ=255)

                tGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGA‑CTGAGGCTGCCAAACGGTGGGAACGGGCGTGCCCGCGAAACAGGCGCCACGGTGCATGGAACCGGCTTTGACGGTG‑t‑tcg                                                                                                                        >  2:10360552/1‑135 (MQ=38)

                  aaGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGA‑CTGAGGCTGCCAAACGGTGGGAGCGGGGGTGTCCCCAAAACAGGCGCCGCGGTGCACTGCCCCCGCTTTCCCCGTG‑t‑t‑ttcg                                                                                                                     >  1:2612005/1‑134 (MQ=255)

                  aaGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGA‑CTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAAAAGGCGCCACGGTGCATGGCACCGGATTTGCCGGTG‑t‑tcgcg                                                                                                                      >  2:1916892/1‑133 (MQ=37)

                  aaGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGA‑CTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGAGAAACAGGCGCAGCGGTGCACGACACCGGCTTTGCCGATG‑t‑tcgcg                                                                                                                      >  2:6329780/1‑133 (MQ=37)

                       cGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGA‑CTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAAAAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑T‑TCG‑CGGGC‑A‑c                                                                                                              >  2:15118833/1‑137 (MQ=37)

                          gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGA‑CTGAGGCTGCCAAACGGTGGGAGCGGGGGTGCCCGCGAAACAAGCGCCACGGTGCACGGCACCAGCTTTGCCGGTG‑T‑TCG‑CGGGC‑A‑c‑a‑c‑c                                                                                                        >  1:1785998/1‑137 (MQ=25)

                          gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGA‑CTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCTCGAAACAAGCGCCGCGGTGCATGGCACCAGCTTTGCCGGTGTT‑T‑T‑CGGGC‑A‑c‑g‑c‑c                                                                                                        >  2:13274671/1‑137 (MQ=25)

                          gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGA‑CTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGTGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑T‑T‑CGCGGGC‑A‑c‑g‑c‑c                                                                                                        >  2:8526317/1‑137 (MQ=25)

                          gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGA‑CTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑T‑T‑TGCGGGC‑A‑c‑g‑c‑c                                                                                                        >  2:4158715/1‑137 (MQ=255)

                          gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGA‑CTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAAAAGGCGCCACGGTGCCCGGCCCCCGCTTTTCCGGTG‑T‑TCG‑CGGGC‑A‑c‑t‑c‑c                                                                                                        >  2:9334172/1‑137 (MQ=37)

                          gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGA‑CTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCACGAAAACGGCGCCAAGGTTCACTGCCCCAGATTTGCCGGTG‑T‑TCT‑CAGGC‑A‑a‑g‑c‑c                                                                                                        >  1:16820284/1‑137 (MQ=255)

                            tGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGA‑CTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑T‑TCG‑CGGGC‑A‑c‑g‑c‑cca                                                                                                      >  2:17479742/1‑135 (MQ=25)

                                 ccATTCAGAGCCTGAACGAAGCTTTGGCACCACAGAAGGA‑TGGAGGCTGCCAAACGGTGGGGGTGGGCTTGCCCGAGAAAAGGGCGCCGGGGGGCATGGCACTGGCTTTGCCGGTG‑T‑T‑CGCGGGC‑A‑C‑GCC‑C‑GCTccc                                                                                                 <  2:10680248/137‑1 (MQ=16)

                                   aTTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGA‑CTGAGGCTGCCAAACGGTGGGAGCGGGGGTGCCCACGAAACAAGCGCCACGGTGTATGGCACCAGCTTTGCCGGTG‑T‑TCG‑CGGGC‑A‑C‑GCC‑C‑GCTCCcac                                                                                               >  2:13253548/1‑137 (MQ=18)

                                   aTTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGA‑CTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCACGAAACAAGCGCCACGGTGCATGGCAACGGCTTTGCCGGTG‑T‑TCG‑CGGGC‑A‑C‑GCC‑C‑GCTCCcac                                                                                               >  2:2335353/1‑137 (MQ=21)

                                        gagCCTGAACGAATCCTTGGAACCACAGCGGGA‑CTGAGGCTGCCAAATGGTGGGAGTGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGACACCGGCTTTGCCGGTG‑T‑T‑CGCGGGC‑A‑C‑GCC‑C‑GCTCCCACAGACt                                                                                          <  2:2156795/137‑1 (MQ=18)

                                            cTGAACGAAGCCTTGGCACCACAGCAGGA‑CTGAGGCTGCCAAACGGTGGGAGCGGGCGTTCCCGCGAAACAAGCGCCGCGGTGTATGGCACCGGATTTTCCGGTG‑T‑TCT‑CGGGC‑A‑C‑GCC‑C‑GCTCCCACAGAATGAtt                                                                                      >  1:6997782/1‑137 (MQ=21)

                                             tGAACGAATTCTTGGCACCAAAGCAGGA‑CTTAGGCTGCCAAACGGGGGGAGCGGGCGTGCCCGCGAAACAGGCACCGTGGGGCATGGCACCGGCTTTGCCGGTG‑T‑TCG‑CGGGC‑A‑C‑GCC‑C‑GCTCCCACAGACTGATTg                                                                                     <  2:12827813/137‑1 (MQ=17)

                                              gAACGAAGCCTTGGCACCACAGCAGGA‑CTGAGGCTGCCAAACGGTGGGAGCGGGCGTGTCCACAAAACAAAAGGCAAGGTGCACGGCACCCGCTTTGTCTGTGTT‑T‑T‑CGGGCAA‑A‑C‑C‑C‑ACTCCCCCAAAATGATTGc                                                                                    >  1:8055392/1‑137 (MQ=255)

                                              gAACGAAGCCTTGGCACCACAGCAGGA‑CTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTTCCGGTG‑T‑TCG‑CGGGC‑A‑C‑GCC‑C‑GCTCCCACAGACTGATTGc                                                                                    >  2:9123087/1‑137 (MQ=21)

                                              gAACGAAGCCTTGGCACCACAGCAGGA‑CTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑T‑TCG‑CGGGC‑A‑C‑GCC‑C‑GCTCCCACAAACTGATTGc                                                                                    >  2:7144334/1‑137 (MQ=21)

                                              gAACGAAGCCTTGGCACCACAGCAGGA‑CTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCAGTG‑T‑TCG‑CGGGC‑A‑C‑GCC‑C‑GCTCCCCCAGACTGCTTGc                                                                                    >  2:14304056/1‑137 (MQ=18)

                                                 cGAAGCCTTGGCACCACAGCAGGA‑CTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑‑‑T‑T‑CGGGC‑A‑C‑GCC‑C‑GCTCCCACAGACTGATTGCCATTa                                                                               >  1:16329588/1‑137 (MQ=25)

                                                 cGAAGCCTTGGCACCACAGCAGGA‑CTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCAGATTTGCCGGTG‑T‑TCG‑CGGGC‑A‑C‑GCC‑C‑GCTCCCAAAGAATGAGTGCCag                                                                                 >  2:6745367/1‑136 (MQ=18)

                                                 cGAAGCCTTGGCACCACAGCAGGA‑CTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCACCGCGGTGCATAGCACCGGCTTTGCCGGTG‑T‑TCG‑CGGGC‑A‑A‑A‑C‑CAGTTCCCAAAGAATGATTGCCAt                                                                                 >  2:14535565/1‑137 (MQ=18)

                                                  gAAGCCTTGGCACCACAGCAGGA‑CTGAGGCTGCCAAACGGTGGGAGCGGGTGTGCCCGCGAAAAAGGCGCCGCGGTGCATGGCACCGGCTTTGACGGTG‑T‑TCG‑CGGGC‑A‑C‑GCC‑C‑GCCCCCACAGACTGATTGCCCtt                                                                                >  2:9871581/1‑137 (MQ=18)

                                                  gAAGCCTTGGCACCACAGCAGGA‑CTGAGGCTGCCAAACGGTGGGAGCGGGNGGGTCCGCGAAACAAGAGCCACGGTGTATGGCACCCGATTTTCCTGTG‑T‑T‑TGCGGGA‑A‑CAC‑C‑C‑CCCCCCCCAAACCGATTGCCact                                                                                >  1:14856323/1‑135 (MQ=255)

                                                  gAAGCCTTGGCACCACAGCAGGA‑CTGAGGCTGCCAAACGGTGGGAGCGGGGGTGTCCGCGAAAAAAGCGCCGCAGTGCATGGCCCCAGCTTTTCCGGTG‑T‑TCT‑CGGGC‑A‑CAC‑C‑C‑ACTCCCACAAAATGATTTCCAtt                                                                                >  1:13446285/1‑137 (MQ=255)

                                                  gAAGCCTTGGCACCACAGCAGGA‑CTGAGGCTGCCAAACGGTGGGAGCGGGGGTGTCCGCGAAAAAAACGCCCCGGTGCATGGCACCAACTTTTCCGGTGTT‑T‑T‑TGGGC‑A‑C‑GCC‑C‑ACCCCCCCAGAATGATTGTCAtt                                                                                >  1:11041004/1‑137 (MQ=255)

                                                  gAAGCCTTGGCACCACAGCAGGA‑CTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑T‑TCG‑CGGGC‑A‑C‑GCC‑C‑GCTCCCACAGACCGATTGCCAtt                                                                                >  2:14625092/1‑137 (MQ=21)

                                                  gAAGCCTTGGCACCACAGCAGGA‑CTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCAGCTTTGCCGGTG‑T‑T‑TGCGGGC‑A‑C‑GCC‑C‑GCTCCCACAGAATGATTGCCAtt                                                                                >  2:10076694/1‑137 (MQ=21)

                                                  gAAGCCTTGGCACCACAGCAGGA‑CTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAAGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑T‑TCG‑CGGGC‑A‑C‑GCC‑C‑GCTCCCACAGACTGATTGCCAtt                                                                                >  2:4615842/1‑137 (MQ=18)

                                                  gAAGCCTTGGCACCACAGCAGGA‑CTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCAAAAAAAGAGCCACGGGGCATGGCACCGGCTTTTCCAGTG‑T‑TCG‑CGGGA‑A‑CAC‑C‑C‑GCTCCCCCAAACTTATTGCCact                                                                                >  2:15674160/1‑135 (MQ=12)

                                                      ccTTGTCACCACAACAGAA‑CGGAGGCTTACAAACGGTGGAAGCGGGCGTGCCCGCGAAACAGGCACCACGGTGCATGGCACCGACTTTTCCGTTG‑T‑TCC‑CGGGC‑A‑C‑GCC‑C‑GCTCCCACAGACTGATTGCCATTAAGc                                                                            <  2:12043863/137‑1 (MQ=12)

                                                        ttGGCACCACAGCAGGA‑CTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGGGCCGAGGTGTATTGCACCGGCTTTGCCGGTG‑T‑TCG‑CGGGC‑A‑C‑GCC‑C‑GCTCCCACAGAATGAGTGACATTCAGCgg                                                                          >  1:4206868/1‑137 (MQ=17)

                                                           gCACCACCGAAGGA‑CTGGGGTTGCCAAATTGTGGGAGCGGGCGTGCCTTTGAAGCGGTCGCTGCGGTGCATGGCCCCGTCTTTGCCGTTT‑T‑TCG‑CGTGC‑A‑C‑GCC‑C‑GCGCCCACAGACTGATTGCCATTAAGCGGGcc                                                                       <  2:16840681/137‑1 (MQ=12)

                                                           gCACCACAGCAGGA‑CTGAGTTTGCCAAACGTTGGGAGCGGGCGTGCCCGCGAAAAGGGCGCCGTGGTGCATGGCACCGGTTTTGCCGTTT‑T‑TCG‑CGGGC‑A‑C‑TCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                       <  2:11945443/137‑1 (MQ=18)

                                                           ccaccacAG‑AGGACCTGCGTCTGCAAAACGGTGGGAGCGGGTGTGCCCGTAAAAAAGGTCCCATGGTCCCTGACCCTGGCTTTGCCGGTG‑T‑T‑CGTGGGC‑A‑C‑ACC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                       <  2:5791/136‑1 (MQ=255)

                                                           agacagcagcggGA‑TTGGGGCTGCAAAACGGTGGGGTTGGGCGTGCCTGCGAAACAGGCGCCGTGGGGCGTGGCACTGGCTTTGCTGGTT‑T‑TCG‑CGGGC‑A‑C‑CCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                       <  1:13122719/131‑1 (MQ=255)

                                                             accacAGCAGGA‑CTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑T‑T‑CGCGGGC‑A‑C‑GCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTa                                                                     >  2:5045813/1‑137 (MQ=21)

                                                                acaGCGGGA‑TTGGGGTTGCAAAGTGGTGGGAGTGGGCGTACCTGCGAAAGGGGCGCCGTGGTGCGTGGCACTGGTTTTGCGGGTG‑T‑TCG‑TGGAC‑A‑C‑CCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCg                                                                  <  1:12044453/137‑1 (MQ=255)

                                                                acaGCAGGA‑CTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑T‑TCG‑CGGGC‑A‑C‑G‑C‑CTGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCg                                                                  <  1:10076399/137‑1 (MQ=21)

                                                                   gcagGA‑CTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCCTGGCACCGGCCTTGCCGGTG‑T‑TCG‑CGGGC‑A‑C‑GCC‑C‑GCTTCCACAGAATGATTGCCATTCAGCGGGCCTATCGCTa                                                               >  2:15912533/1‑137 (MQ=18)

                                                                   gcagGA‑CTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑T‑T‑CGCGGGC‑A‑C‑GCC‑C‑GCTCCCACAGACTGATTGCCATTCAGCGGGCCTATCGCTa                                                               >  2:16519715/1‑137 (MQ=21)

                                                                   gcagGA‑CTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAAGAGCCGCGGTGCATGGAACCGGATTTGCCAGTG‑T‑TCG‑CGGGC‑A‑C‑GCC‑C‑AACCCCACAGAATGATTGCCAGTGAGCGGGCCCATCGCTa                                                               >  2:2567899/1‑137 (MQ=12)

                                                                    cagGA‑CTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGGAACAGGCGCCGCGGTGCATGGCACCAGATTTGCCAGTG‑T‑TCG‑CGGGC‑A‑C‑GCC‑C‑GCTCCCCCAGACTGATTGCCATTAAGCGGGTCTATCGCtcc                                                              >  2:6118829/1‑135 (MQ=18)

                                                                    cagGA‑CTGAGGCTGCCAAACGGTGGGAGCGGGAGTGCCCGCGAAAAAAGCGCCACAGGGTATGGCACCAGCTTTGTCGGTG‑T‑T‑TGTGGGT‑A‑CAC‑C‑C‑GCTCCCCCAGAATGATTGCCCTTCAGCGGGCCTATctctca                                                              >  1:14526343/1‑135 (MQ=255)

                                                                     agGA‑CTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑T‑TCG‑CGGGC‑A‑C‑GCC‑C‑GCTCCCACAGACTGATTGACATTCAGCGGGCCTATCGCTACa                                                             >  2:3709845/1‑137 (MQ=18)

                                                                       gA‑CTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑T‑T‑CGCGGGC‑A‑C‑GCC‑C‑GCTCCAAAAGACTGAGTGCAATTAAGCGGGCCTATCGCTACAGt                                                           >  2:6311251/1‑137 (MQ=18)

                                                                        gattGAGGCTACCAAATGGTGGGAGCGGTCGTGCCCGCGAAACAGACGCCGCGGTGCATGGCACCGGCTTTACCGGGG‑T‑T‑CGCGGGC‑A‑C‑GCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGt                                                           <  2:8101926/134‑1 (MQ=18)

                                                                          cTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑T‑TCG‑CGGGC‑A‑A‑GCC‑C‑GCTCCCCCAGAATGATTGCCATTAAACCGGCCCCTCGCTGCAGTgg                                                         >  2:13436845/1‑136 (MQ=14)

                                                                           tGAGTCTACAAAATGGTGGGAGTGGGCGTGCCCGCGAAACGGGCGCCGTGGTGCATGGCACCGGTTTTGCTGGTG‑T‑TCG‑CGGGC‑A‑C‑CCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGAc                                                        <  2:868166/137‑1 (MQ=255)

                                                                           tGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAAAAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑T‑TCG‑CGGGC‑A‑C‑GCC‑C‑GCTCCCAAAGACTGATTGTCATTAAGCGGGCCTATCGCTACAGTGaa                                                        >  2:16407324/1‑136 (MQ=18)

                                                                             gggCTGCCAAATGGGGGGAGTGGGTTTGCCTGGGAAAAGGGTGCCGCGGTGCGGGGCACTGGTTTTGCCGGTG‑T‑T‑TGTGGGC‑A‑C‑GCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACgg                                                      <  1:10344492/136‑1 (MQ=255)

                                                                             aGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACCGGGGCCGCGGTGCCTGGCACTGGTTTTGTTGGTG‑T‑TCG‑TGGGC‑C‑C‑G‑CTC‑GCTTCCACAGAATGATTTTCATTAAGCGGGACCCTCGCTACCGTGAAgg                                                      >  1:12943394/1‑137 (MQ=12)

                                                                                cTGCCACACGGTGGGAGCGGGCGTGCCCGCAAAACCAGCACCGCGGTTCATGGCACCGGCTTTGCCGGTG‑T‑T‑CGCGGGC‑A‑C‑G‑C‑CTACTCCCACAGACTGATTGCCACTAAGCGGGCCTATCGCTACAGTGACGGccc                                                   >  2:4386153/1‑137 (MQ=18)

                                                                                cTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑T‑TCG‑CGGGC‑A‑C‑GCC‑C‑GCTCCCAAAGAATGATTGCCATTAAGCGGGCATATCGCTAAAGTGAAGGacc                                                   >  2:2908079/1‑137 (MQ=17)

                                                                                cTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCAGCTTTGCCGGTG‑T‑TCG‑CGGGC‑C‑C‑GCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGccc                                                   >  2:8418754/1‑137 (MQ=21)

                                                                                cTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCAGCTTTGCCGGTG‑T‑TCG‑CGGGC‑A‑A‑GCC‑C‑GCTCCCACAGAATGATTGCCATTAAGCGGGCCTATCGATACAATGGCGGccc                                                   >  2:1747756/1‑137 (MQ=17)

                                                                                     aaaCGGTGGGAGCGGGCGTGCCCGCGAAACCGGCGCCGCGGTGCATTGCACCGGCTTTGCCGGTG‑T‑T‑CGCGGGC‑A‑C‑GCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCG‑gttt                                             >  2:12886312/1‑137 (MQ=21)

                                                                                     aaaCGGTGGGAGCGGGCGTGCCCGCGAAACAGGTGCCGCGGTGCATGGCCCCGGATTTGCCGGTG‑T‑TCG‑CGGGC‑A‑C‑GCC‑C‑GCTCCCCCAAACTGATTGCCATTCAGCGGGCCTATCGCTACATTGACGGCCCT‑tttt                                             >  2:11132556/1‑137 (MQ=18)

                                                                                     aaaCGGTGGGAGCGGGCGTGCCCGCGAAACAGGAGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑T‑TCG‑CGGGA‑A‑C‑GAC‑C‑GCTCCCAAAGACTGATTGCCATTCAGCGGGCCTATCGCTACAGTGAAGGTCCG‑tttt                                             >  2:11088710/1‑137 (MQ=18)

                                                                                     aaaCGGTGGGAGCGGGCGTGCCCGCGAAACAAGCGCCACGGTGCATGGCACCGGCTTTGCCGGTG‑T‑TCG‑CGGGC‑A‑C‑GCC‑C‑GCTCCCACAAACTGATTGCCATTAAGCGGGACTATCGCTACAGTGACGGCCCG‑tttt                                             >  2:6236668/1‑137 (MQ=17)

                                                                                       agggggggAGCGGGCGTTCCCGAGAAACGGGCGCTGGGGTGCGTGGCATTGGTTTTTCTGGTG‑T‑T‑TGGGGGC‑A‑C‑GCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCG‑tttttt                                           <  1:15747898/132‑1 (MQ=255)

                                                                                       acttttggAGCGGGTTTGCCCGCGAAACAGTCGCCCCTTCCCATGGCACCGGCTTTGCCTGTG‑T‑T‑CGCGGGC‑A‑C‑GCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCG‑tttttt                                           <  2:2605184/131‑1 (MQ=16)

                                                                                           tGGGAGCGGGCGTTCCCGAGAAACAGGCGCCGCGGTGCGTGGCACCGGCTTTGCCGGTG‑T‑TCG‑CGGGC‑A‑C‑GCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCG‑TTTTttgttg                                       <  2:13546489/137‑1 (MQ=21)

                                                                                                  gggCGTGCCTGCGAAAGAGGCCCCGCGGTACATGGCACCGGTTTTGTTGGTT‑T‑T‑CGCGGGC‑A‑C‑GCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCG‑TTTTTTGTTGCCTCGCt                                <  2:10287539/137‑1 (MQ=25)

                                                                                                     cGTGCCCGCGAAACAGGCGCCACGGTGCATTGCACCGGCTTTTCCCATG‑T‑TCT‑CGGGC‑A‑CAC‑C‑C‑ACTCCCACAGACTGATTGTCATTCAGCGGGCCTCCTGCTACAATGACAGCCCG‑TTTTTTGTTTCCCCACTAAc                             >  2:12646643/1‑137 (MQ=255)

                                                                                                        gaccGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑T‑T‑TGTGGGA‑A‑C‑GCC‑C‑GCTTCCACAGAATGATTGCCATTCAGCGGGACTATCACTCCAGTGAAAGCCCG‑TTTTTTGTTTCCCCACCAACGaa                          >  2:16391926/3‑137 (MQ=12)

                                                                                                         cccGCGAAACAGGCGCCGCGGTGCATGGCACCGTCTTTGCAGGTT‑T‑TCG‑CGGGC‑A‑C‑GCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCG‑TTTTTTGTTGCCTCGCTAACGAAt                         <  2:17178478/137‑1 (MQ=25)

                                                                                                           ttcgAAACAGGCACCTCGGGGCGCGGCACCTGCTTTTCCGGTG‑TCC‑G‑TGGGC‑A‑CAC‑C‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCG‑TTTTTTGTTGCCTCGCTAACGAAGCa                       <  2:3393770/135‑1 (MQ=255)

                                                                                                           tgcgAAACAGGCACCGGGGTAAGTGCCCCCGGCTTTTCCGTTG‑T‑TCC‑CGGCC‑ATC‑C‑C‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCG‑TTTTTTGTTGCCTCGCTAACGAATCa                       <  1:6864924/136‑1 (MQ=255)

                                                                                                           cgcgAAACGGGCGCCGTGGTGCGTGGCACCGTTTTTGCTGGTG‑T‑T‑CGCGGGC‑A‑C‑GCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCG‑TTTTTTGTTGCCTCGCTAACGAATCa                       <  2:484541/137‑1 (MQ=25)

                                                                                                           cgcgAAACGGGCGCCGCGGTGAGTGGCACCGGCTTTGCCGGTG‑T‑TCG‑GGGGC‑A‑C‑GCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCG‑TTTTTTGTTGCCTCGCTAACGAATCa                       <  2:2627045/137‑1 (MQ=25)

                                                                                                           cgcgAAACATTCTCTACGGGGCGTGGCGCCGTCTTTGCCGTTT‑T‑T‑TGCGGAT‑A‑C‑GCC‑C‑GCTTCCACAGACTGGTTGCCATTAAGCGGGCCTGTCGCTACAGTGACGGCCCG‑TTTTTTGTTGCCTCGCTAACGAATCa                       <  2:16982801/137‑1 (MQ=255)

                                                                                                           cgcgAAACAGGCGCCGCGTTGCATGGCACCGGTTTTGCCGGTG‑T‑T‑CGCGGGC‑A‑C‑GCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCG‑TTTTTTGTTGCCTCGCTAACGAATCa                       <  2:2016868/137‑1 (MQ=37)

                                                                                                           cgcgAAACAGGCGCCGCGGTGAATGGCACTGGTTTTGCCGGTG‑T‑TCG‑CGGAC‑A‑C‑CCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCG‑TTTTTTGTTGCCTCGCTAACGAATCa                       <  1:3080486/137‑1 (MQ=255)

                                                                                                           cgcgAAACAGGCACCGCGGGGCATGGCACCGGTTTTTCAGGGG‑T‑TCG‑CGGGC‑A‑C‑ACC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCG‑TTTTTTGTTGCCTCGCTAACGAATCa                       <  2:5743031/137‑1 (MQ=255)

                                                                                                           cacgAAAAAGGCGCCGTGGGGCGTGGCACTGTCTTTTCCGGTG‑T‑T‑TGCGGGC‑A‑C‑GCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCG‑TTTTTTGTTGCCTCGCTAACGAATCa                       <  1:6906198/135‑1 (MQ=25)

                                                                                                           caaaaaaCAGGCACCGCGGGGCATGGCACCGTCTTTTCCGGTG‑T‑TCG‑CGGAC‑A‑C‑CCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCG‑TTTTTTGTTGCCTCGCTAACGAATCa                       <  1:6097346/133‑1 (MQ=255)

                                                                                                                aaCAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑T‑TCG‑CGGGC‑A‑C‑GCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCG‑TTTTTTGTTGCCTCGCTAACGAATCATGACt                  >  2:6879348/1‑137 (MQ=37)

                                                                                                                 aCAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑T‑TCG‑CGGGC‑A‑C‑GCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCG‑TTTTTTGTTGCCTCGCTAACGAATCATGACTc                 >  1:12535005/1‑137 (MQ=37)

                                                                                                                      cgccgtGGTGCATGGCACCGGTTTTGTTGGTG‑T‑TCG‑CGGAC‑A‑C‑GCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCG‑TTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa            <  1:9522607/137‑1 (MQ=255)

                                                                                                                      cgccgcGGGGCATGGCACCGGTTTTTATGGTT‑T‑T‑CGCGGGC‑A‑C‑GCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCG‑TTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa            <  1:6711586/137‑1 (MQ=255)

                                                                                                                                cATGGCACCGGTTTTTCTGGTG‑T‑T‑AGCGGGC‑A‑C‑TCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCG‑TTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATc  <  2:1843489/137‑1 (MQ=255)

                                                                                                                                 aTGGCACCGGCTTTGCCGGTG‑T‑TCG‑CGGGC‑A‑C‑CCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATc  <  1:15011177/137‑1 (MQ=37)

||

CGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGA‑CTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTG‑T‑T‑T‑CGGGC‑A‑C‑G‑C‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCG‑TTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATC  >  AE015451_phaM1/1499347‑1499608

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 9 ≤` `ATCG/ATCG` `< 10 ≤` `ATCG/ATCG` `< 24 ≤` `ATCG/ATCG` `< 40 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`