BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	AE015451_phaM1	1,499,498	1 bp→CG	intergenic (+75/‑89)	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,497	1	.	C	77.0%	502.8 / 226.5	245	intergenic (+74/‑90)	trpS CDS (Tryptophan‑‑tRNA ligase)/zapE CDS (putative enzyme with nucleoside triphosphate hydrolase domain)	–/–
Reads supporting (aligned to +/- strand): ref base . (25/31); new base C (150/38); total (175/70)
Fisher's exact test for biased strand distribution p-value = 1.06e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*	AE015451_phaM1	1,499,498	0	T	G	83.2%	452.6 / 36.6	245	intergenic (+75/‑89)	trpS CDS (Tryptophan‑‑tRNA ligase)/zapE CDS (putative enzyme with nucleoside triphosphate hydrolase domain)	–/–
Reads supporting (aligned to +/- strand): ref base T (23/12); new base G (137/36); total (175/70)
Fisher's exact test for biased strand distribution p-value = 1.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 263 total aligned reads displayed.

AAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGA‑CTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCAC‑GGCTTTGCCGGTG‑TT‑T‑CGGGC‑A‑C‑G‑C‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTT  >  AE015451_phaM1/1499365‑1499633

|||

aaGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGA‑CTGAGGCTGCCAAACGGTGGGAGCGGGGGTGTCCTCGAAACAAGAGCCGCGGTGTATGGCACC‑AGCCTTGCTAGTG‑tt‑ttcg                                                                                                                                             >  1:7812040/1‑134 (MQ=255)

   cccGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGA‑CTGAGGCTGCCAAACGGTGGGAGCGGGGGTGCCCGCGAAACAAGAGCCGAGGGGCACGGCACC‑AGACTTGCCGGTt‑tt‑tttgggt                                                                                                                                          >  1:13244592/1‑131 (MQ=255)

     cGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGA‑CTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAAAAGGCGCCGCGGTGCATGGCACC‑GGCTTTGCCGGTG‑TTCG‑CGGGC‑A‑c                                                                                                                                      >  2:15118833/1‑137 (MQ=37)

        gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGA‑CTGAGGCTGCCAAACGGTGGGAGCGGGGGTGCCCGCGAAACAAGCGCCACGGTGCACGGCACC‑AGCTTTGCCGGTG‑TTCG‑CGGGC‑A‑c‑a‑c‑c                                                                                                                                >  1:1785998/1‑137 (MQ=25)

        gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGA‑CTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAAAAGGCGCCACGGTGCCCGGCCCC‑CGCTTTTCCGGTG‑TTCG‑CGGGC‑A‑c‑t‑c‑c                                                                                                                                >  2:9334172/1‑137 (MQ=37)

        gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGA‑CTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCACGAAAACGGCGCCAAGGTTCACTGCCCC‑AGATTTGCCGGTG‑TTCT‑CAGGC‑A‑a‑g‑c‑c                                                                                                                                >  1:16820284/1‑137 (MQ=255)

          tGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGA‑CTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCACGAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTGCCGGTG‑TTCG‑CGGGC‑A‑c‑g‑c‑ccg                                                                                                                              >  2:16520937/1‑135 (MQ=25)

                 aTTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGA‑CTGAGGCTGCCAAACGGTGGGAGCGGGGGTGCCCACGAAACAAGCGCCACGGTGTATGGCACC‑AGCTTTGCCGGTG‑TTCG‑CGGGC‑A‑C‑GCC‑C‑GCTCCcac                                                                                                                       >  2:13253548/1‑137 (MQ=18)

                 aTTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGA‑CTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCCTGGCACC‑GGCTTTGCCGGTG‑TTCG‑CGGGC‑A‑C‑GCC‑C‑GCTccccc                                                                                                                       >  2:2943963/1‑135 (MQ=21)

                 aTTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGA‑CTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCACGAAACAAGCGCCACGGTGCATGGCAAC‑GGCTTTGCCGGTG‑TTCG‑CGGGC‑A‑C‑GCC‑C‑GCTCCcac                                                                                                                       >  2:2335353/1‑137 (MQ=21)

                      gagCCTGAACGAATCCTTGGCACCACAGCAGGA‑CGGAGGCTGCCAAATGGTGGAAGCGGGCGTGCCCGCAAAACAGGCGCCGTGGTGCGTGACACC‑GGCTTTTCCGGTG‑TTCG‑CGGGC‑A‑CTC‑C‑C‑GCTCCCACAGACt                                                                                                                  <  2:4344904/137‑1 (MQ=255)

                           tGAACGAATTCTTGGCACCAAAGCAGGA‑CTTAGGCTGCCAAACGGGGGGAGCGGGCGTGCCCGCGAAACAGGCACCGTGGGGCATGGCACC‑GGCTTTGCCGGTG‑TTCG‑CGGGC‑A‑C‑GCC‑C‑GCTCCCACAGACTGATTg                                                                                                             <  2:12827813/137‑1 (MQ=17)

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

                            gAACGAAGCCTTGGCACCACAGCAGGA‑CTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTGCCAGTG‑TTCG‑CGGGC‑A‑C‑GCC‑C‑GCTCCCCCAGACTGCTTGc                                                                                                            >  2:14304056/1‑137 (MQ=18)

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

                                gAAGCCTTGGCACCACAGCAGGA‑CTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTGCCGGTG‑TTCG‑CGGGC‑A‑C‑GCC‑C‑GCTCCCACAGACCGATTGCCAtt                                                                                                        >  2:14625092/1‑137 (MQ=21)

                                gAAGCCTTGGCACCACAGCAGGA‑CTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCAAC‑GGCTTTGCCGGTG‑TTCG‑CGGGC‑A‑C‑GCC‑C‑GCTCCCACAGACTGATTGCCAtt                                                                                                        >  2:11718330/1‑137 (MQ=21)

                                gAAGCCTTGGCACCACAGCAGGA‑CTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCACCGCGGTGCATGGCAAC‑GGATTTGCCAGTG‑TTCG‑CGGGC‑A‑C‑GCC‑C‑GCTCCCAAAAAATGATTGCCagt                                                                                                        >  2:9164625/1‑135 (MQ=18)

                                gAAGCCTTGGCACCACAGCAGGA‑CTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAAGCGCCGCGGTGCATGGCACC‑GGCTTTGCCGGTG‑TTCG‑CGGGC‑A‑C‑GCC‑C‑GCTCCCACAGACTGATTGCCAtt                                                                                                        >  2:4615842/1‑137 (MQ=18)

                                gAAGCCTTGGCACCACAGCAGGA‑CTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCAAAAAAAGAGCCACGGGGCATGGCACC‑GGCTTTTCCAGTG‑TTCG‑CGGGA‑A‑CAC‑C‑C‑GCTCCCCCAAACTTATTGCCact                                                                                                        >  2:15674160/1‑135 (MQ=12)

                                  aGCCTTGTCACCACAGCAGGA‑CTGGGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACGGGCGCCGCGGTGCATGGCACC‑GGCTTTGCCGGTG‑TTCG‑CGGGC‑A‑C‑GCC‑C‑GCTCCCACAGACTGATTGCCATTaa                                                                                                      <  2:16948129/137‑1 (MQ=18)

                                    ccTTGTCACCACAACAGAA‑CGGAGGCTTACAAACGGTGGAAGCGGGCGTGCCCGCGAAACAGGCACCACGGTGCATGGCACC‑GACTTTTCCGTTG‑TTCC‑CGGGC‑A‑C‑GCC‑C‑GCTCCCACAGACTGATTGCCATTAAGc                                                                                                    <  2:12043863/137‑1 (MQ=12)

                                      ttGGCACCACAGCAGGA‑CTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGGGCCGAGGTGTATTGCACC‑GGCTTTGCCGGTG‑TTCG‑CGGGC‑A‑C‑GCC‑C‑GCTCCCACAGAATGAGTGACATTCAGCgg                                                                                                  >  1:4206868/1‑137 (MQ=17)

                                         ccaccacAG‑AGGACCTGCGTCTGCAAAACGGTGGGAGCGGGTGTGCCCGTAAAAAAGGTCCCATGGTCCCTGACCCT‑GGCTTTGCCGGTG‑TT‑CGTGGGC‑A‑C‑ACC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                                               <  2:5791/136‑1 (MQ=255)

                                         agacagcagcggGA‑TTGGGGCTGCAAAACGGTGGGGTTGGGCGTGCCTGCGAAACAGGCGCCGTGGGGCGTGGCACT‑GGCTTTGCTGGTT‑TTCG‑CGGGC‑A‑C‑CCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                                               <  1:13122719/131‑1 (MQ=255)

                                          caccacAGCAGGA‑CTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTGCCGGTG‑TTCG‑CGGGC‑A‑C‑GCC‑C‑GATCCCACAGACTGATTGCCATTAAGCGGGACt                                                                                              >  2:10402044/1‑137 (MQ=18)

                                          caccacAGCAGGA‑CTGAGGCTGCCAAACGGTGGGAGCGGGAGTGTCCACGAAAAAAGAACCACGGTGTCTGGCACC‑CGATTTTCCGGTGTTT‑T‑CGGGG‑A‑CAC‑C‑C‑CCCCCCCCAAAATGATTGTCACTCAGCGGGGCt                                                                                              >  1:14649028/1‑137 (MQ=255)

                                             ccaagcagGA‑TTGAGGCTGCAAAACGGTGGGAGCGGGCGTGCCCGTGAAACAGGCGCCGCTGTACATGGCACC‑TGCTTTGCCGGTG‑TTCG‑CGGGC‑A‑CTC‑C‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATc                                                                                           <  2:16053061/134‑1 (MQ=18)

                                              acaGCGGGA‑TTGGGGTTGCAAAGTGGTGGGAGTGGGCGTACCTGCGAAAGGGGCGCCGTGGTGCGTGGCACT‑GGTTTTGCGGGTG‑TTCG‑TGGAC‑A‑C‑CCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCg                                                                                          <  1:12044453/137‑1 (MQ=255)

                                              acaGCAGGA‑CTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTGCCGGTG‑TTCG‑CGGGC‑A‑C‑G‑C‑CTGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCg                                                                                          <  1:10076399/137‑1 (MQ=21)

                                                 gcagGA‑CTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAAGAGCCGCGGTGCATGGAACC‑GGATTTGCCAGTG‑TTCG‑CGGGC‑A‑C‑GCC‑C‑AACCCCACAGAATGATTGCCAGTGAGCGGGCCCATCGCTa                                                                                       >  2:2567899/1‑137 (MQ=12)

                                                  cagGA‑CTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGGAACAGGCGCCGCGGTGCATGGCACC‑AGATTTGCCAGTG‑TTCG‑CGGGC‑A‑C‑GCC‑C‑GCTCCCCCAGACTGATTGCCATTAAGCGGGTCTATCGCtcc                                                                                      >  2:6118829/1‑135 (MQ=18)

                                                      gattGAGGCTACCAAATGGTGGGAGCGGTCGTGCCCGCGAAACAGACGCCGCGGTGCATGGCACC‑GGCTTTACCGGGG‑TT‑CGCGGGC‑A‑C‑GCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGt                                                                                   <  2:8101926/134‑1 (MQ=18)

                                                          gAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAAGCGCCACAGTGCCTGGCACC‑AGCTTTTCCGGTG‑TTCG‑CGGGC‑A‑C‑GCC‑C‑GCTCCCACAAACCGATTGCCCTTCAGCGGGCCTCTCGCTCCAATGGCg                                                                               >  2:13415938/1‑137 (MQ=9)

                                                           aGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACCGGGGCCGCGGTGCCTGGCACT‑GGTTTTGTTGGTG‑TTCG‑TGGGC‑C‑C‑G‑CTC‑GCTTCCACAGAATGATTTTCATTAAGCGGGACCCTCGCTACCGTGAAgg                                                                              >  1:12943394/1‑137 (MQ=12)

                                                              cTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTGCCGGTG‑TTCG‑CGGGC‑A‑C‑GCC‑C‑GCTCCCAAAGAATGATTGCCATTAAGCGGGCATATCGCTAAAGTGAAGGacc                                                                           >  2:2908079/1‑137 (MQ=17)

                                                              cTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑AGCTTTGCCGGTG‑TTCG‑CGGGC‑C‑C‑GCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGccc                                                                           >  2:8418754/1‑137 (MQ=21)

                                                              cTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACC‑AGCTTTGCCGGTG‑TTCG‑CGGGC‑A‑A‑GCC‑C‑GCTCCCACAGAATGATTGCCATTAAGCGGGCCTATCGATACAATGGCGGccc                                                                           >  2:1747756/1‑137 (MQ=17)

                                                                   aaaCGGTGGGAGCGGGCGTGCCCGCGAAACCGGCGCCGCGGTGCATTGCACC‑GGCTTTGCCGGTG‑TT‑CGCGGGC‑A‑C‑GCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGgttt                                                                      >  2:12886312/1‑137 (MQ=21)

                                                                   aaaCGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGAGGTGCAGGGCACC‑GGCTTTGCCGGTG‑TTCG‑AGGGC‑A‑A‑GCC‑C‑GCCCCCACAGAATGATTGCCATTAAGCGGGACTCTCGCTACAATGACGGCCCGtttt                                                                      >  2:16489255/1‑137 (MQ=17)

                                                                   aaaCGGTGGGAGCGGGCGTGCCCGCGAAACAGGAGCCGCGGTGCATGGCACC‑GGCTTTGCCGGTG‑TTCG‑CGGGA‑A‑C‑GAC‑C‑GCTCCCAAAGACTGATTGCCATTCAGCGGGCCTATCGCTACAGTGAAGGTCCGtttt                                                                      >  2:11088710/1‑137 (MQ=18)

                                                                   aaaCGGTGGGAGCGGGCGTGCCCGCGAAACAAGCGCCACGGTGCATGGCACC‑GGCTTTGCCGGTG‑TTCG‑CGGGC‑A‑C‑GCC‑C‑GCTCCCACAAACTGATTGCCATTAAGCGGGACTATCGCTACAGTGACGGCCCGtttt                                                                      >  2:6236668/1‑137 (MQ=17)

                                                                       ggggggAGTGGGTGTGCCCGAGAAAAGGGCGCTGGGGGGCGT‑GCCACTGGTTTTGCTGGTG‑TTCG‑GGGGC‑A‑C‑CCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTttgt                                                                  <  1:6448523/134‑1 (MQ=255)

                                                                       ggTGGGAGAGGGCGTGCCTGCGAAACGGGCGCTGCGGTGAATGGCACC‑GGCTTTGCGGGTG‑TT‑CGCGGGA‑ACC‑C‑C‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTttgt                                                                  <  2:631284/137‑1 (MQ=255)

                                                                         tGGGAGCGGGCGTTCCCGAGAAACAGGCGCCGCGGTGCGTGGCACC‑GGCTTTGCCGGTG‑TTCG‑CGGGC‑A‑C‑GCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTttgttg                                                                <  2:13546489/137‑1 (MQ=21)

                                                                           ggAGCGGGCGTGCCCGCGAAACAGGCGCCGTGGTGCATGGCACC‑GTCTTTGCCGGTG‑TT‑AGCGGAC‑A‑C‑GCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGcc                                                              <  2:6774309/137‑1 (MQ=21)

                                                                                   cGTGCCCGCGAAACAGGCGCCACGGTGCATTGCACC‑GGCTTTTCCCATG‑TTCT‑CGGGC‑A‑CAC‑C‑C‑ACTCCCACAGACTGATTGTCATTCAGCGGGCCTCCTGCTACAATGACAGCCCGTTTTTTGTTTCCCCACTAAc                                                      >  2:12646643/1‑137 (MQ=255)

                                                                                       cccGCGAAACAGGCGCCGCGGTGCATGGCACC‑GTCTTTGCAGGTT‑TTCG‑CGGGC‑A‑C‑GCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAAt                                                  <  2:17178478/137‑1 (MQ=25)

                                                                                         cgcgAAATAGTCGCCGTGGTGCCTGGCACC‑GTTTTTGCTGTTT‑TT‑TGTGGGC‑A‑C‑TCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                <  2:8637738/137‑1 (MQ=255)

                                                                                         cgcgAAACGGGCGCCGTGGTGCGTGGCACC‑GTTTTTGCTGGTG‑TT‑CGCGGGC‑A‑C‑GCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                <  2:484541/137‑1 (MQ=25)

                                                                                         cgcgAAACAGGCGCCGCGTTGCATGGCACC‑GGTTTTGCCGGTG‑TT‑CGCGGGC‑A‑C‑GCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                <  2:2016868/137‑1 (MQ=37)

                                                                                         cgcgAAACAGGCGCCGCGGTGCATGGCACC‑GGCTTTTCCGGTG‑TTCG‑TGGGC‑A‑C‑TCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                <  1:6864963/137‑1 (MQ=25)

                                                                                         cgcgAAACAGGCGCCGCGGTGAATGGCACT‑GGTTTTGCCGGTG‑TTCG‑CGGAC‑A‑C‑CCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                <  1:3080486/137‑1 (MQ=255)

                                                                                         cgcgAAACAGGCACCGCGGGGCATGGCACC‑GGTTTTTCAGGGG‑TTCG‑CGGGC‑A‑C‑ACC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                <  2:5743031/137‑1 (MQ=255)

                                                                                         cgcgAAAAAGTCGCTGTGGTGAGTGGCACT‑GGTTTTGCGGGTG‑TT‑TGCGGAC‑A‑C‑CCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                <  1:558772/137‑1 (MQ=255)

                                                                                         cgcaAAAAAGACGCCGCGGGACATGACACC‑GGTTTTGCCGGTG‑TTCG‑TGGGA‑A‑C‑CCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                <  2:2627044/137‑1 (MQ=255)

                                                                                               aCAGGCGCCGCGGTGCATGGCACC‑GGCTTTGCCGGTG‑TTCG‑CGGGC‑A‑C‑GCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTc                                          >  1:12535005/1‑137 (MQ=37)

                                                                                               aCAGGCGCCGCGGTGCATGGCACC‑GGCTTTGCCGGTG‑TTCG‑CGGGC‑A‑C‑GCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTc                                          >  1:6988725/1‑137 (MQ=37)

                                                                                                    ctccgtGGTGCATGGCACC‑GTTTTTTCTGGTG‑TTCG‑CGGAC‑A‑C‑GCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa                                     <  1:2750198/135‑1 (MQ=255)

                                                                                                    cgccggGGGGCATGGCACT‑GGTTTTGCTGGTG‑TT‑CGCGGGC‑A‑C‑CCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa                                     <  1:14915996/137‑1 (MQ=255)

                                                                                                    caccgcGGCGCATGGCACC‑GGTTTTGCTGGTG‑TT‑TGGGGGC‑A‑C‑GCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa                                     <  2:10198965/135‑1 (MQ=37)

                                                                                                     gccgcGGTGCATGGCACC‑GGCTTTGCCGGTG‑TTCG‑CGGGC‑A‑C‑GCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAg                                    >  1:1035865/1‑137 (MQ=37)

                                                                                                         cgGTGCATGGCACC‑GGCTTTGCCGGTG‑TTCG‑CGGGC‑A‑C‑GCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACg                                >  2:3950671/1‑137 (MQ=37)

                                                                                                              cATGGCACC‑GGTTTTTCTGGTG‑TT‑AGCGGGC‑A‑C‑TCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATc                           <  2:1843489/137‑1 (MQ=255)

                                                                                                                 ggcacc‑GGCTTTGCCGGTG‑TTCG‑CGGGC‑A‑C‑GCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAg                        >  2:2222709/6‑137 (MQ=37)

                                                                                                                    acc‑GGCTTTGCCGGTG‑TTCG‑CGGGC‑A‑C‑GCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAg                     >  2:8961495/3‑137 (MQ=38)

                                                                                                                    acc‑GGCTTTGCCGGTG‑TTCG‑CGGGC‑A‑C‑GCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAg                     >  1:2375555/3‑137 (MQ=38)

                                                                                                                      c‑GGCTTTGCCGGTG‑TTCG‑CGGGC‑A‑C‑GCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGAt                   >  1:16772680/1‑137 (MQ=38)

                                                                                                                          cTTTGCCGGTG‑TTCG‑CGGGC‑A‑C‑GCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTATAACGCTATCAAGCAGATCTg                >  2:4080456/1‑137 (MQ=38)

                                                                                                                          cTTTGCCGGTG‑TTCG‑CGGGC‑A‑C‑GCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                >  2:9808887/1‑137 (MQ=39)

                                                                                                                          cTTTGCCGGTG‑TTCG‑CGGGC‑A‑C‑GCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                >  2:6969436/1‑137 (MQ=39)

                                                                                                                          cTTTGCCGGTG‑TTCG‑CGGGC‑A‑C‑GCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                >  2:7205293/1‑137 (MQ=39)

                                                                                                                          cTTTGCCGGTG‑TTCG‑CGGGC‑A‑C‑GCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                >  2:5083728/1‑137 (MQ=39)

                                                                                                                          cTTTGCCGGTG‑TTCG‑CGGGC‑A‑C‑GCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                >  1:6843935/1‑137 (MQ=39)

                                                                                                                          cTTTGCCGGTG‑TTCG‑CGGGC‑A‑C‑GCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                >  1:10226596/1‑137 (MQ=39)

                                                                                                                          cTTTGCCGGTG‑TTCG‑CGGGC‑A‑C‑GCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                >  1:10728007/1‑137 (MQ=39)

                                                                                                                          cTTTGCCGGTG‑TTCG‑CGGGC‑A‑C‑GCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                >  1:13705614/1‑137 (MQ=39)

                                                                                                                          cTTTGCCGGTG‑TTCG‑CGGGC‑A‑C‑GCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                >  1:14392878/1‑137 (MQ=39)

                                                                                                                          cTTTGCCGGTG‑TTCG‑CGGGC‑A‑C‑GCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                >  1:15093673/1‑137 (MQ=39)

                                                                                                                          cTTTGCCGGTG‑TTCG‑CGGGC‑A‑C‑GCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                >  1:15601554/1‑137 (MQ=39)

                                                                                                                          cTTTGCCGGTG‑TTCG‑CGGGC‑A‑C‑GCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                >  1:3392203/1‑137 (MQ=39)

                                                                                                                          cTTTGCCGGTG‑TTCG‑CGGGC‑A‑C‑GCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                >  1:5025392/1‑137 (MQ=39)

                                                                                                                          cTTTGCCGGTG‑TTCG‑CGGGC‑A‑C‑GCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                >  1:9667168/1‑137 (MQ=39)

                                                                                                                          cTTTGCCGGTG‑TTCG‑CGGGC‑A‑C‑GCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                >  2:12197051/1‑137 (MQ=39)

                                                                                                                          cTTTGCCGGTG‑TTCG‑CGGGC‑A‑C‑GCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                >  2:14858723/1‑137 (MQ=39)

                                                                                                                          cTTTGCCGGTG‑TTCG‑CGGGC‑A‑C‑GCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCACGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                >  2:17556756/1‑137 (MQ=38)

                                                                                                                          cTTTGCCGGTG‑TTCG‑CGGGC‑A‑C‑GCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGACCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCACTTAGAACGCTATCAAGCAGATCTg                >  2:8762547/1‑137 (MQ=38)

                                                                                                                          cTTTGCCGGTG‑TTCG‑CGGGC‑A‑C‑GCC‑C‑GCTCACACAGACTGATTGCCATTAAGCGGGCCTATCGCTAAAGTGACGGCACGTTTTTTGTTGCCTCGCTAACGAATCATGGCTCCCTTAGAACGCTATCAAGCAGATCTg                >  2:2633704/1‑137 (MQ=38)

                                                                                                                          cTTTGCCGGTG‑TT‑CGTGGGC‑A‑C‑G‑C‑CAGCTCCCACAGACTGATTGCCATTAAGCGGGCATATCGCTACAGTGACGGCCCGTTTTTTATTGCCTCGCTAACGAATCATGACTCACTTAGAACGCTATCGAGCAGATCTg                >  2:10855645/1‑137 (MQ=255)

                                                                                                                          cTTTGCCGGTG‑TT‑CGCGGGC‑A‑C‑GCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGCCTCCCTTAGAACGCTATCAAGCAGATCTg                >  2:12351114/1‑137 (MQ=38)

                                                                                                                          atttGCCGGTG‑TTCG‑CGGGC‑A‑C‑GCC‑C‑GCTCCCACAGAATGATTGCCATTAAGCGGGCCCATCGCTACAGTGACGGCCCGTTTTTTGTTGTCTCGCTAACGAATCATGACTCCCTTAGGACGCTATCAAGCAGGTCTg                >  2:2133051/2‑137 (MQ=37)

                                                                                                                            ttGCCGGTG‑TTCG‑CGGGC‑A‑C‑GCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATc                           <  1:9042898/124‑1 (MQ=38)

                                                                                                                            ttGCCGGTG‑TTCG‑CGGGC‑A‑C‑GCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGaa              >  2:7814361/1‑137 (MQ=39)

                                                                                                                            ttGCCGGTG‑TTCG‑CGGGC‑A‑C‑GCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGaa              >  2:7768815/1‑137 (MQ=39)

                                                                                                                            ttGCCGGTG‑TTCG‑CGGGC‑A‑C‑GCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTAAAGTGACGGCCCGTTTTTTGTTGCCTCGATAACGAATCATGACTCCCTTAGAACGCTATCAAGTAGATCTGaa              >  2:3112548/1‑137 (MQ=38)

                                                                                                                            ttGCCGGTG‑TTCG‑CGGGC‑A‑C‑GCC‑C‑GCTCCCACAAACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGGCTCCCTTAGAACGCTATAAAGCAGGTCTGaa              >  2:7114176/1‑137 (MQ=38)

                                                                                                                                        tcgcGGGC‑A‑C‑GCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGACCTGAAACGTCCCGAc    >  2:1789650/4‑137 (MQ=255)

                                                                                                                                        tcgcGGGC‑A‑C‑GCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCATTAGAACGCTATCAAGCAGATCTGAAACGTCCCGAc    >  1:4625966/4‑137 (MQ=255)

                                                                                                                                         cgcGGGC‑A‑C‑GCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGAct   >  2:1613518/3‑137 (MQ=255)

                                                                                                                                         cgcGGGC‑A‑C‑GCC‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGAct   >  2:3444218/3‑137 (MQ=255)

|||

AAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGA‑CTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCAC‑GGCTTTGCCGGTG‑TT‑T‑CGGGC‑A‑C‑G‑C‑C‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTT  >  AE015451_phaM1/1499365‑1499633

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