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	95.1%	314.5 / ‑3.9	144	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 (6/1); new base A (67/70); total (73/71)
Fisher's exact test for biased strand distribution p-value = 1.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*	AE015451_phaM1	1,499,481	0	A	C	96.5%	349.5 / ‑5.5	144	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 (4/1); new base C (69/70); total (73/71)
Fisher's exact test for biased strand distribution p-value = 3.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 144 total aligned reads displayed.

GTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGC‑CACGGCTT‑TGCCGGT‑G‑T‑T‑T‑CG‑GGC‑A‑C‑G‑C‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATC  >  AE015451_phaM1/1499351‑1499608

||

gTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACATGCGCCGCGGTGCATGGC‑ACCGGCtt                                                                                                                                       <  1:2130891/137‑1 (MQ=255)

 ttCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCGGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGC‑ACCGGCtt‑t                                                                                                                                     <  1:10442525/137‑1 (MQ=255)

  tCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGC‑ACCGGCTT‑Tg                                                                                                                                    <  1:9536/137‑1 (MQ=255)

  tCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGC‑ACCGGCTT‑Tg                                                                                                                                    <  1:8354360/137‑1 (MQ=255)

  tCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGC‑ACCGGCTT‑Tg                                                                                                                                    <  1:5154432/137‑1 (MQ=255)

  tCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGC‑ACCGGCTT‑Tg                                                                                                                                    <  1:2913308/137‑1 (MQ=255)

  tCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGC‑ACCGGCTT‑Tg                                                                                                                                    <  2:2736040/137‑1 (MQ=255)

         cGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGTGTGCCCGCGAAAAAAGCGCCGCGGTGTACGGC‑ACCTGTTT‑TGCCGgt‑gg                                                                                                                            >  2:2864049/1‑136 (MQ=255)

            tGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGTCCCCGAAACAAGCGCCGCGGTGCCTTGC‑CCCGGCTC‑TTCCGGT‑G‑t‑t‑tg                                                                                                                      >  1:11151487/1‑136 (MQ=255)

            tGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGC‑ACCGGATT‑TGCCGGT‑G‑t‑tcg                                                                                                                       >  2:6187724/1‑135 (MQ=255)

                      gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGGGTGCCCCCGAAACAAGAGCCACTGTGCACGGC‑ACCAGACT‑TTTCGGT‑G‑T‑T‑TGCG‑GGA‑A‑c‑g‑c‑c                                                                                                      >  1:10365934/1‑137 (MQ=255)

                      gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGC‑ACCGGCTT‑TGCCGGT‑G‑T‑TCG‑CG‑GGC‑A‑c‑g‑c‑c                                                                                                      >  2:4007115/1‑137 (MQ=25)

                      gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAAAAGGCGCCGCGGGGCATGGC‑CCCGGATT‑TGTCGGT‑G‑T‑T‑TGCG‑GGC‑A‑c‑g‑c‑c                                                                                                      >  1:8788177/1‑137 (MQ=37)

                      gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCAAAACAGGGGCCGCGGTGCATGGC‑ACCGGCTT‑TGTCGGT‑G‑T‑TCG‑CA‑GGC‑A‑c‑g‑c‑c                                                                                                      >  2:5692553/1‑137 (MQ=25)

                      gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCAAAACAAGCGCCAAGGTGCCTGGC‑ACCGGATT‑TGCCGGT‑G‑T‑TCG‑CG‑GGC‑A‑c‑g‑c‑c                                                                                                      >  1:8788219/1‑137 (MQ=25)

                      gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCACAAAACAAGCGCCGCTGTGCATGGC‑ACCAGCTT‑TTCCCGT‑G‑T‑TCG‑CG‑GGC‑A‑c‑g‑c‑c                                                                                                      >  1:14698481/1‑137 (MQ=25)

                      gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGAGTGTCCGCGAAACAAGCGACACGGTGTATGGA‑ACCTGACT‑TGTCGGT‑G‑T‑TCG‑CG‑GGC‑A‑c‑g‑c‑c                                                                                                      >  1:6920361/1‑137 (MQ=255)

                      gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGAGTGTCCACGAAACAAGCGCCACGGTGCCTTGCACCCGCCTT‑TCCCGTT‑T‑T‑C‑T‑CG‑GGcccacc                                                                                                           >  1:7926179/1‑132 (MQ=255)

                       aTGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGAGTGCCCACGAAAAAAGCGCCACGGTGCCTGGC‑ACCAGCTT‑TTCCGGT‑G‑T‑T‑TGCG‑GGG‑A‑c‑g‑c‑cc                                                                                                     >  1:3409503/1‑136 (MQ=255)

                        tGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGTGTGCCCACAAAACAGGCGTCGCGGTGTATGGC‑ACCGGCTT‑TGCCTTT‑GTT‑T‑T‑CG‑GGT‑A‑c‑g‑c‑cca                                                                                                    >  1:1574440/1‑135 (MQ=255)

                        tGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGAGCCACGGTGCACGGC‑CCCAGACT‑TTTCGGT‑G‑T‑T‑TGCG‑GGG‑A‑c‑g‑c‑ccg                                                                                                    >  1:2876421/1‑135 (MQ=255)

                        tGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGAGTGTCCGCGAAACAGGCGCCGCGGTGCCTGGC‑ACCGGCTT‑TGCCGGT‑G‑T‑TCG‑CG‑GGC‑A‑c‑g‑c‑cca                                                                                                    >  1:18477335/1‑135 (MQ=37)

                               aTTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGC‑ACCGGCTT‑TGCCGGT‑G‑T‑TCG‑CG‑GGC‑A‑A‑GCC‑CGCTCCcac                                                                                              >  2:16161304/1‑137 (MQ=18)

                               aTTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGAGAAAAAGGCGCCGAGGTGCATGGG‑AACGGCTT‑TGCCGGT‑G‑T‑T‑TGCG‑GGG‑A‑C‑GAC‑CGCTCCcac                                                                                              >  1:13416872/1‑137 (MQ=25)

                                   agggCCTGAACGAAGCCTTGGCACCACAACATGATTGAGGCTGCCAAACGGTGGGAGCGGGTGGGCCCGCGAAACGGGCGCCGCGGTGCATGGC‑ACTGGTTT‑TGCCGGT‑G‑T‑TCG‑TG‑GGA‑A‑C‑GCC‑CGCTCCCACAGAc                                                                                          <  2:12209045/134‑1 (MQ=255)

                                    gagCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCACGAAAAAGGCGCCGCGGTGCATGGC‑AACGGCTT‑TGCCGGT‑G‑T‑TCG‑CG‑GGC‑A‑C‑GCC‑CGCTCCCACAGACt                                                                                         >  2:12727240/1‑137 (MQ=25)

                                    aagCCTGAACAAAACCTTGGCACCAAGGAAGGACGGAGGCTGCAAAACGGGGGGAGCGGGCGTGCCCGCGAAACAGACACCGCGGTGCATGGC‑ACCGGCTT‑TGCCGGT‑G‑T‑TCG‑CG‑GAC‑A‑C‑GCC‑CGCTCCCACAGACt                                                                                         <  2:17107008/136‑1 (MQ=9)

                                          gAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCTCGAAACAAGCGCCGCGGTGCATGGC‑ACCGACTT‑TTCCGGT‑G‑T‑TCG‑CG‑GGC‑ACC‑C‑C‑CACTCCCACAAACTCATTGc                                                                                   >  2:16403143/1‑137 (MQ=255)

                                          gAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCACGGTGCATGGC‑ACCGGCTT‑TGCCGGT‑G‑T‑T‑CGCG‑GGC‑A‑C‑GCC‑CGCTCCCACAGAATGATTGc                                                                                   >  2:8638097/1‑137 (MQ=21)

                                          gAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAAGCGCCGCGGTGCATGGC‑ACCGGCTT‑TGCCGGT‑G‑T‑T‑TGCG‑GGC‑A‑C‑GCC‑CGCTCCCACAGACTGATTGc                                                                                   >  2:17673899/1‑137 (MQ=21)

                                             cGAAGCTTTGTCACCACAGCAGGGTTGAGGCTACCAAACGGTGGGGGCGGGCGTGCGCGCGAAACGGGCGCCGTTGTGCGTGGC‑ACCGGCTT‑TGCTGGT‑G‑T‑T‑TGCG‑GGC‑A‑T‑GCC‑CGCTCCCACAGACTGATTGCCAt                                                                                <  2:15808892/137‑1 (MQ=16)

                                             cGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCAAAACAGGTGCCGCGGTGCATGGC‑ACCGGCTT‑TGCCGGT‑G‑T‑TCG‑CG‑GGC‑A‑C‑ACC‑CGCTCCCAAAGACTGATTgtcct                                                                                >  2:15640509/1‑133 (MQ=18)

                                              gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACCGGCGCCACGGTTCATGGC‑ACCAGCTT‑TGCCGGT‑G‑T‑TCT‑CA‑GGT‑A‑C‑ACC‑CCCTTCCACAGACTGATTGGCAtt                                                                               >  1:13167018/1‑137 (MQ=21)

                                              gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCCTGGC‑ACCAGCTT‑TGCCGGT‑G‑T‑TCT‑CG‑GGC‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCAtt                                                                               >  2:16580361/1‑137 (MQ=21)

                                              gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGC‑ACCGGCTT‑TGCCGGT‑G‑T‑TCT‑CG‑GGC‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCAtt                                                                               >  2:15614547/1‑137 (MQ=25)

                                              gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGC‑ACCGGCTT‑TGCCGGT‑G‑T‑TCG‑CG‑GGC‑A‑C‑GCC‑CGCTCCCACAGAATGATTGCCAtt                                                                               >  2:9156560/1‑137 (MQ=21)

                                              gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCACGGC‑ACCGGCTT‑TGCCGGT‑G‑T‑TCG‑CG‑GGC‑A‑C‑GCC‑CGATGCCACAGAACGATTGCCAtt                                                                               >  2:1205746/1‑137 (MQ=18)

                                              gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCACGAAACAAGCGCCGAGGTGTATGGC‑ACCTGATT‑TTCCGGT‑G‑T‑TCG‑CA‑GGCAA‑A‑C‑C‑CGCTCCCACAGAATGATTGCCCtt                                                                               >  1:3845641/1‑137 (MQ=12)

                                              gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCACGAAAAAAGAGCCACCGTTTCTGGC‑ACCCGATT‑TGCCTGT‑G‑T‑TCT‑CG‑GGA‑A‑CAC‑C‑CCCTCCCCCCAAATGATTGCCact                                                                               >  1:5817128/1‑135 (MQ=255)

                                              gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGAGTGCCCGCGAAACAAGCGCCGCGGTGCATGGC‑ACCGGCTT‑TGCCGGT‑G‑T‑TCG‑CG‑GGC‑A‑C‑GCC‑CGCTCCCCCAGAATGATTGCCAtt                                                                               >  1:17133334/1‑137 (MQ=18)

                                              gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGAGTGCCCGCAAAAAAAGAACCTCAGTGCATGGC‑ACCAGCTT‑TTCCAGT‑GTT‑T‑T‑CG‑GGC‑A‑CAC‑C‑CCCTCCCACCAACCTATTGcccct                                                                               >  2:5433856/1‑134 (MQ=255)

                                              aaaGCTTTGGCACCACAGCAGGATGGGGGCTGCAAAATGGTGGGAGGGGTCGGCCCTGCAAAACGGGCGCCGTGGGGAATGGC‑ACTGGTTT‑TGCTGGT‑G‑T‑TCG‑TG‑GGC‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCAtt                                                                               <  1:12792947/136‑1 (MQ=12)

                                                  ccTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCACGGTGCCTGGC‑ACCAGCTT‑TGCCGGT‑G‑T‑TCA‑CG‑GGC‑A‑C‑GCC‑CACTCCCACAGACTGATTGCCATTAAGc                                                                           >  2:12269267/1‑137 (MQ=21)

                                                  ccTTGGCACCACAGAAGGGTTGAGTCGGCAAAATGGTGGGAGCGGGTGGGCCCGCGAAAAAGACACCGCGGGGCAGGGA‑ACCGGCTT‑TGCTGGTGG‑T‑C‑G‑CG‑GGC‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGc                                                                           <  2:11044341/137‑1 (MQ=9)

                                                    ttGGCACCACAGAAGGATTGGGGTTGCAAAACGGGGGGAGCGGGCGTGCCCGCGAAAAGGGCGCCGGGGGGCAGGGC‑ACTGGTTT‑TGCAGGT‑G‑T‑TCG‑CG‑GGC‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCgg                                                                         <  2:6735418/137‑1 (MQ=12)

                                                    ttGACACCACAGAAGAACTGAGGCTGCCAAACGGTGGGAGCGGGCTTTCCTAAGAAAAAGTCGCCGTGGGACAGGAC‑ACCGGTTT‑TCCCGGT‑G‑T‑T‑CGCG‑GAC‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCgg                                                                         <  2:77658/137‑1 (MQ=255)

                                                      ggCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGTCCACGAAACAGGCGCCACGGTGTATGGC‑ACCGGCTT‑TGCCGGT‑G‑T‑TCG‑CG‑GGC‑A‑C‑GCC‑CGCTCCCCCCAACTGATTTCCATTCAACGGGc                                                                       >  2:10027655/1‑137 (MQ=17)

                                                       gCACCACGGCGGGATTGAGGTGACCAAATGGTGGGAGGGGTTGTGCCCGCGAAACAGGAGCCAAGTGGCATGGA‑ACTGGTTT‑TTCCGGT‑T‑T‑TCG‑CG‑GGC‑A‑C‑ACC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                      <  2:14700943/137‑1 (MQ=255)

                                                       gCACCACAGCAGGGCTGAGGCTGCCAAACGTTGGGAGCGGGCGGCCCCGCGAAACAGGCCCCAAGGGGCAGGGC‑ACTGGCTT‑TGCCGGT‑G‑T‑T‑TGCG‑GGC‑A‑C‑ACC‑CGCTCCCGCAGACTGATTGCCATTAAGCGGGcc                                                                      <  2:3931265/137‑1 (MQ=255)

                                                       gCACCACAGCAGGACTGAGTCTCCCAAACGGTGGGAGCGGGTGTGCCCGCGAAACGGGCGCCGCGGTGCAGGAA‑ACCGGCTT‑GGCCGTT‑G‑T‑TCG‑CG‑GGC‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                      <  2:2851360/137‑1 (MQ=18)

                                                       gCACCACAGCAGGACGGAGGCTGCCAAATGGTGGGAGCGGGTGTGCCCGCGAAACAGGCGCCGCGGTGCATGGC‑ACCGGTTT‑TGCCGGT‑G‑T‑TCG‑CG‑GGC‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                      <  2:8857660/137‑1 (MQ=21)

                                                       gCACCACAGAAGGATTGAGGCGGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAAAAGGCGCAGCGGGGCATGGC‑ACCGGCTT‑TGCCGTT‑G‑T‑T‑CGCG‑GGC‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                      <  2:8853627/137‑1 (MQ=18)

                                                       gCACCACAGAAGGATGGGGGTTGCAAAATGGGGGGAGGGGGTGTCCCCGAAAAACAGGCGCCGCGGTGAAGGGC‑ACTGGTTT‑TGCTGGT‑G‑T‑T‑CGCG‑GGC‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                      <  2:1483340/137‑1 (MQ=12)

                                                       acaccacATCAGGACTGAGGCGGCAAAATGGGGGGAGCGGGTGTGCCAGAGAAACAGGCACAGCGGGACATGGC‑ACCGGTTT‑TGCCGGT‑T‑T‑TCG‑CG‑GGA‑A‑C‑CCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                      <  2:9717769/136‑1 (MQ=255)

                                                        caccacAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCCTGGC‑ACCGGCTT‑TGCCGGT‑G‑T‑TCG‑CG‑GGC‑A‑C‑GCC‑CGCTCCCACAGAATGATTGCCATTAAGCGGGCCt                                                                     >  2:15651896/1‑137 (MQ=21)

                                                            acaGCAGGACTGAGGCTGCCAAATGGTGGGAGCGGGTGTGCCCGCGAAACAGGCGCCGTGGGGCATGGC‑ACTGGCTT‑TGCCGGT‑G‑T‑TCG‑CG‑GGC‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCg                                                                 <  2:16224513/137‑1 (MQ=21)

                                                               gcagGACTGAGGCTGCCAAACGGTGGGAGCGGGAGTGCCCGCGAAACAAGCACCACGGTGCCCGGC‑ACCCGCTT‑TTCCGGT‑G‑T‑TCT‑CG‑GGC‑A‑CAC‑C‑CCCCCCCCCAAAATGATTGCCACTCCACGGGACTATCACtc                                                              >  1:11976457/1‑136 (MQ=255)

                                                                   gACTGAGGCTGCCAAACGGTGGGAGCGGGGGTGCCCACGAAACAAGTACCACGGTGCATGGC‑ACCAGGTT‑TGCCGGT‑G‑T‑TCG‑CG‑GGC‑A‑CAC‑C‑CACCCCCCCAGACCGAGTGCCCTTCAGCGGGACTATCGCTaaaat                                                          >  1:1816872/1‑133 (MQ=9)

                                                                       gAGGCTGACAAATGGTGGGAGCGGGTGTGCCCGCGAAAAAGGCGCCGCGGTGCATGGC‑ACCGGCTT‑TGCCGGT‑G‑T‑TCG‑CG‑GGC‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACg                                                      <  2:16166808/137‑1 (MQ=18)

                                                                           cTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCCTGGC‑ACCGGCTT‑TGCCGGT‑G‑T‑TCG‑CG‑GGC‑A‑C‑GCC‑CGCTCCCAAAGAATGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGccc                                                  >  2:17849984/1‑137 (MQ=21)

                                                                           cTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGC‑ACCGGCTT‑TGCCGGT‑G‑T‑T‑TGCG‑GGC‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCCTTCAGCGGGCCTATCGCTACAGTGACGGccc                                                  >  2:7075035/1‑137 (MQ=21)

                                                                           cTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGC‑ACCGGCTT‑TGCCGGT‑G‑T‑T‑CGCG‑GGC‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGtcc                                                  >  2:12217738/1‑137 (MQ=21)

                                                                           cTGCCAAACGGTGGGAGCGGGAGTGCCCGCAAAACAAGAGCCGCAGTTCCTGGC‑CCCAGCTT‑TTCCCGT‑GTT‑T‑T‑TG‑GGC‑A‑CAC‑C‑CGCTCCCCCAAAACGATTGCCATTCAGCGGGTCTATCCCTCCAATGACGGacc                                                  >  2:941820/1‑137 (MQ=255)

                                                                           cTGCCAAACGGTGGGAGCGGGAGTGCCCACGAAACAAGCGCCACGGTGCATGGC‑ACCGGATT‑TGTCGGT‑G‑T‑TCT‑CG‑GGC‑A‑C‑GCC‑CACTCCCACAAAACGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGccc                                                  >  2:2668583/1‑137 (MQ=255)

                                                                                gaaCGGTGGGAGTGGGTGTGGCCGAGGAAAAGGCGCCGCGGTGCATGGA‑ACCGGATT‑TGAAGGT‑G‑T‑TCG‑CG‑GGC‑A‑A‑G‑CACGCTCCCACAAACTTATTGCCATTAAGTGGGCCTATCACTACAGTGACAGCCCGtttt                                             >  2:15439402/2‑137 (MQ=255)

                                                                                aaaCGGTGGGAGCGGGCGTTCCCGCGAAACAGGCGCCACAGTGCACGGC‑ACCGGCTT‑TGCCTGT‑G‑T‑TCA‑CG‑GGC‑A‑C‑GCC‑CGCCCCCCCCGACTGATTGCCCTTCAACAGGCCTATCGCTACAGTGACGGCCCGtttt                                             >  2:8547585/1‑137 (MQ=255)

                                                                                aaaCGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGC‑ACCGGCTT‑TGCCGGT‑G‑T‑T‑CGCG‑GGC‑A‑C‑GCC‑CGATCCCACAGACTGATTGCCATTAAGCGGGCCTCTCGCTACAGTTCCGGCCCGtttt                                             >  2:17127893/1‑137 (MQ=18)

                                                                                aaaCGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCACGGC‑ACCAGCTT‑TGCCGGT‑G‑T‑TCG‑CG‑GGC‑A‑CAC‑C‑CGCTCCCACAGACTGATTGCCACTAAGCGGGAATATCGATACAGTGACGGAACGtttt                                             >  2:8547538/1‑137 (MQ=17)

                                                                                        ggggtgggTGTGCCTGAGAAACGGGCGCTGTGGTGCAGGGC‑ACTGGTTT‑GGTTGGT‑G‑T‑TAG‑TG‑GAC‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGcc                                     <  1:7834872/132‑1 (MQ=255)

                                                                                          aGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGC‑ACCGGCTT‑TGCCGGT‑G‑T‑TCG‑CG‑GGC‑A‑CAC‑C‑CGCTCCCACAGACCGATTGCCATTAAGCGGGCCTATCGCTAAAATGACGGCCCGTTTTTTGTTGCCTc                                   >  2:18046116/1‑137 (MQ=18)

                                                                                           gCGGGCGTGCCCACGAAACAAGCGCCGCGGTGCACGGC‑ACCGGATT‑TGCCGGT‑G‑T‑TCG‑CG‑GGC‑A‑C‑GCC‑CGCTCCCACAAACTGATTGCCATTCAGCGGGCCTACCGCTACAGTGACGGCCCGTTTTTTTTTGCCTCg                                  >  2:4021771/1‑137 (MQ=18)

                                                                                            tggttgTGCCCGCAAAACGGGCGCCGTGTTGCATGGC‑ACTGTTTT‑TTTCGGT‑G‑T‑TCG‑TG‑GAC‑A‑C‑CCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGc                                 <  1:7107358/132‑1 (MQ=255)

                                                                                            cGGGCATTCCCGCAAAACAGGCGCCGCGGCGCATGGC‑ACCTACTT‑TACCGGT‑G‑T‑TCG‑CG‑GGC‑A‑CTC‑C‑CACTCCCACAGACTGATTACCATTAAACGGGCCTATCGCTACACTGACGGCCCGTTTTTTGTTGCCTCGc                                 <  2:15756961/137‑1 (MQ=255)

                                                                                             gggcgggcCCACGAAACAGGCACCGCGGTGCGTGG‑‑CACCGTTTGTACTGTT‑G‑T‑TCG‑TG‑GGC‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCt                                <  1:3975625/137‑1 (MQ=255)

                                                                                                   gCCCGCGAAACAGGCGCCGCGGTGCATGGC‑ACCGGCTT‑TGCCGGT‑G‑T‑TCG‑CG‑GGC‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGaa                          >  1:18377998/1‑137 (MQ=25)

                                                                                                      cgtgAAACGGGCGCCGGGGTGCATGGC‑ACCGGTTT‑TGCCGGT‑G‑T‑TCG‑CG‑GGC‑A‑C‑ACC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                       <  2:6299534/134‑1 (MQ=25)

                                                                                                      cgtgAAACAGGCGCCGTGGGGCGTGGC‑ACCGGCTT‑TGCTGGT‑G‑T‑T‑TGCG‑GGC‑A‑C‑TCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                       <  1:16280599/134‑1 (MQ=255)

                                                                                                      cgtgAAACAGGCGCCGAGGTGCCTTGC‑ACCGGCTG‑TGCTGGT‑G‑T‑T‑C‑GGTGGC‑A‑C‑GCC‑CGCTCTCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                       <  2:4228083/134‑1 (MQ=255)

                                                                                                      cgcgAAATAGGCGCCGTGGTGCGTGTA‑ACTGGCTT‑TGCTGGT‑G‑T‑TCG‑CG‑GGC‑A‑C‑CCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                       <  2:1736802/137‑1 (MQ=255)

                                                                                                      cgcgAAACGGGCTCCACGGTGCATGGC‑ACCGGCTT‑TACTGTT‑T‑T‑T‑CGTG‑GGC‑A‑C‑TCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                       <  1:3391724/137‑1 (MQ=255)

                                                                                                      cgcgAAACGGGCGCCGTGGGGCAGGAC‑ACCGGTTT‑TGTAGGT‑G‑T‑TCG‑CG‑GAC‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                       <  2:1746102/137‑1 (MQ=255)

                                                                                                      cgcgAAACGGGCGCCGCGGTACGTGTC‑ACCGGCTT‑TTCCGGT‑G‑T‑TCG‑CG‑GGC‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                       <  2:11616370/137‑1 (MQ=25)

                                                                                                      cgcgAAACAGGTGCCGCGGTACATGGC‑ACCGGCTT‑TGCCGGT‑G‑T‑TCG‑CG‑GGC‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                       <  2:6945442/137‑1 (MQ=37)

                                                                                                      cgcgAAACAGGCGCCGTGGTGCGTGGC‑ACCGGCTT‑TTCCGGT‑G‑T‑TCG‑CG‑GGC‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                       <  2:11412587/137‑1 (MQ=25)

                                                                                                      cgcgAAACAGGCGCCGCGGTGCATGGC‑ACCGGCTT‑TGCCGGT‑G‑T‑T‑CGCG‑GAC‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                       <  2:2118104/137‑1 (MQ=25)

                                                                                                      cgcgAAACAGCCCCCGTGGGGAATGGC‑ACTGGCTT‑TGCTGGT‑G‑T‑TCA‑GG‑GGC‑A‑T‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                       <  1:8479008/137‑1 (MQ=255)

                                                                                                      cgagAAACGGGCGCGGTGGTGCATGGC‑ACCGGTTT‑GGCGGGT‑T‑T‑TCG‑CG‑GGC‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                       <  2:4806553/134‑1 (MQ=25)

                                                                                                       gcgAAACAGTCGCCGCGGTGCATGGC‑ACCGGCTT‑TGCCGGT‑G‑T‑T‑CGTG‑GGC‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCAt                      <  2:7635300/137‑1 (MQ=25)

                                                                                                       gcgAAACAGGCGCCGTGTTCCGTGTC‑ACTGGCTT‑GGCTGTG‑G‑T‑T‑CGCG‑GGC‑A‑C‑TCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGCGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCAt                      <  2:14250331/137‑1 (MQ=255)

                                                                                                       gcgAAACAGGCGCCACGGTGCATGGC‑ACCGGTTT‑TGCTGGT‑G‑T‑T‑CGCG‑GAC‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCAt                      <  2:8053469/137‑1 (MQ=255)

                                                                                                       acaaaaCAGGCGCCCCGATGCAGGGC‑ACCGCCTT‑TGCAGGT‑G‑T‑CCG‑CG‑GGC‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTTACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCAt                      <  2:3343609/134‑1 (MQ=25)

                                                                                                         gAAAAAGGCGCCGCGGTGCGTGGC‑ACCGGCTT‑TGCGGGG‑G‑T‑T‑CGCG‑GGC‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGa                    <  2:17958142/137‑1 (MQ=37)

                                                                                                           aataGGCGCCGCAGTGCATGGC‑ACCGGCTT‑TGCCGGT‑G‑T‑TCG‑CG‑GGC‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGATACAGTGACGGCCAGTTGTTTGTTGCCTCGCTAACGAATCGTGCCt                  >  2:10389638/4‑137 (MQ=21)

                                                                                                            aCAGGCGTCGCGGTGCATGGC‑ACCGGCTT‑TGCCGGT‑G‑T‑TCG‑CG‑GGC‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTc                 >  2:3378910/1‑137 (MQ=25)

                                                                                                            aCAGGCGCCGCGGTGCATGGC‑ACCGGGTT‑TGCCGGT‑G‑T‑TCG‑CG‑GGC‑A‑C‑GCC‑CGCTCCCACAAACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGAATc                 >  2:9215777/1‑137 (MQ=25)

                                                                                                                 cgctgtggTGCGTGGC‑ACCGGCTT‑TGCCGTT‑G‑T‑TCG‑CG‑GGC‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa            <  1:13832102/131‑1 (MQ=25)

                                                                                                                 cgccgtGGTGCATGGC‑ACCGGTTT‑GGCTGGT‑G‑T‑TCG‑CG‑GGC‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa            <  2:7873477/137‑1 (MQ=37)

                                                                                                                 cgccgcGGTGCCTGGC‑ACCGGTTT‑TGCCGTT‑G‑TCC‑G‑TG‑GAC‑A‑C‑TCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa            <  1:13987396/137‑1 (MQ=255)

                                                                                                                           cATGGC‑ACCGGCTT‑TGTCGGT‑G‑T‑TAG‑CG‑GGC‑A‑C‑TCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATc  <  1:6980392/137‑1 (MQ=255)

                                                                                                                           cATGGC‑ACCGGCTT‑TGCCGGT‑G‑T‑C‑CGCG‑GGC‑A‑C‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGAATCCCTTAGAACGCTATc  <  2:12989992/137‑1 (MQ=37)

||

GTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGC‑CACGGCTT‑TGCCGGT‑G‑T‑T‑T‑CG‑GGC‑A‑C‑G‑C‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATC  >  AE015451_phaM1/1499351‑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: `‑`