BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	AE015451_phaM1	1,499,508	+C	intergenic (+85/‑79)	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,506	1	.	C	95.5%	840.2 / 37.4	269	intergenic (+83/‑81)	trpS CDS (Tryptophan‑‑tRNA ligase)/zapE CDS (putative enzyme with nucleoside triphosphate hydrolase domain)	–/–
Reads supporting (aligned to +/- strand): ref base . (11/1); new base C (185/68); total (200/69)
Fisher's exact test for biased strand distribution p-value = 1.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 283 total aligned reads displayed.

GATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGC‑CACGGCTTT‑GCCGGT‑G‑TT‑T‑CG‑G‑G‑‑CAC‑G‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTA‑‑‑‑GAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTC  >  AE015451_phaM1/1499373‑1499637

gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGGGTGCCCGCGAAACAGGCGCCACGGTGCATGGA‑ACCGGCTTT‑TTCGGT‑GTTT‑T‑CG‑G‑G‑‑CAc‑a‑cc                                                                                                                                       >  1:10467994/1‑137 (MQ=255)

gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCCTGGC‑ACCGGATTT‑GCCGGT‑G‑TTCG‑CG‑G‑G‑‑CAc‑g‑cc                                                                                                                                       >  2:16218334/1‑137 (MQ=25)

gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGC‑ACCAGATTT‑GCCGGT‑G‑TTCG‑CA‑G‑G‑‑CAc‑g‑cc                                                                                                                                       >  2:310580/1‑137 (MQ=25)

gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAAGCGCCGCGGTGCACGGC‑ACCAGCTTT‑GCCGGT‑G‑TTCG‑CG‑G‑G‑‑CAc‑g‑cc                                                                                                                                       >  2:10816246/1‑137 (MQ=25)

gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGAGTGCCCGCGAAAAAAGAGCCGCGGTGCATGGG‑ACCAGATTT‑TCCGGT‑G‑TTCT‑TG‑G‑G‑‑CAc‑a‑cc                                                                                                                                       >  1:10175655/1‑137 (MQ=255)

 aTGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTTCCCGCGAAACAGGCGCCGCGGGGCATGGC‑ACCGGCTTT‑GCCAGT‑G‑TTCG‑CG‑G‑G‑‑CAc‑g‑ccc                                                                                                                                      >  2:9600394/1‑136 (MQ=37)

  tGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGGGTGGCCGCGAAACAAGAGCCACGGTGCCTGGC‑CCCGGGTTT‑TCCGGT‑GTTT‑T‑CG‑G‑G‑‑CAc‑gaccg                                                                                                                                      >  1:4777954/1‑133 (MQ=255)

  tGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAAGCGCCGCGGTGCATGGC‑AACAGCTTT‑TCCCGT‑G‑TTCG‑CG‑G‑G‑‑CAc‑g‑cccg                                                                                                                                     >  1:10938027/1‑135 (MQ=25)

       ccATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAAGCGCCGCGGTGCATGGC‑ACCGGGTTT‑GCCGGT‑G‑TTCG‑CG‑G‑G‑‑CAC‑GCCCGCTccc                                                                                                                                 >  2:15297598/1‑137 (MQ=18)

              gatccTGAACGAAGCTTTGGCACCACAGCAGGACGGAGGCTGCCAAGCGGTGGGAGCGGGCGTGCCCGTGAAACATGCGCCGCGGTGCATGGC‑ACCGGCTTT‑GCCGGT‑G‑TT‑CGCG‑G‑A‑‑CAC‑GCCCGCTCCCACAGACt                                                                                                                          <  2:15394393/134‑1 (MQ=17)

                 ccTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGTGTTTCCACGAAAAAAGCGCCTCAGTGTATCGC‑ACCAGCTTT‑TTCGGT‑G‑TTCT‑CT‑G‑G‑‑CACAC‑CCCCTCCCCCAGACtact                                                                                                                       >  1:954359/1‑134 (MQ=255)

                    gAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGC‑ACCAGCTTT‑GCCGGT‑G‑TTCG‑CG‑G‑G‑‑CAC‑GCCCGCTCCCACAGACTGATTGc                                                                                                                    >  2:4323706/1‑137 (MQ=21)

                        gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGGGTGTCCGCAAAACAGGCGCCGCGGGGCAGGGC‑ACCGGCTTT‑GTCGGT‑G‑TTCG‑CG‑G‑G‑‑TAC‑GCCCACCCCCACAAACCAATTGCCact                                                                                                                >  1:15265554/1‑135 (MQ=255)

                        gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGTGCCGCGGTGCATGGC‑ACCAGATTT‑GCCGGT‑G‑TTCG‑CG‑G‑G‑‑CAC‑GCCCACTCCCACAGAATGACTGCCact                                                                                                                >  2:5352466/1‑135 (MQ=18)

                        gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGC‑ACCGGCTTT‑GCCGGT‑G‑TTCG‑CG‑G‑G‑‑CAC‑GCCCGCTCCCACAGACTGATTGCCact                                                                                                                >  2:3671859/1‑135 (MQ=21)

                        gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAAGCGCCGCGGTGCATGGC‑ACCGGCTTT‑GCCGGT‑G‑TTCG‑CG‑G‑G‑‑CAC‑GCCCGCTCCCACAGACTGATTGCCAtt                                                                                                                >  2:7218445/1‑137 (MQ=18)

                        gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAAGCACCACGGTGCCTGGC‑ACCCGCCTT‑TCCGGT‑G‑TTCT‑CG‑G‑G‑‑CACAC‑CCGCCCCCACAAACCGATTTcccct                                                                                                                >  2:5264328/1‑134 (MQ=16)

                        gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAAAAGGCGCCGCAGGGTATGGC‑AACGGATTT‑GCCGG‑‑G‑AT‑T‑CGTGAG‑‑CGC‑GCCCGATCCCACAGACCGATTGCCAtt                                                                                                                >  2:279590/1‑137 (MQ=19)

                        gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAAAAGGCGCCGAAGTGCATGGC‑ACCGGGTTT‑GACGGT‑G‑TTCT‑CG‑G‑G‑‑ACC‑GTCCGCTCCCAAAGAACGATTGCCAtt                                                                                                                >  2:7117617/1‑137 (MQ=21)

                        gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAAAAAGCGCCGCAGTGCATGGC‑ACCAGCTTT‑GCCGGT‑G‑TTCG‑CG‑G‑G‑‑CAC‑GCCCACCCCCACCGAATGATTGCCact                                                                                                                >  1:12233814/1‑135 (MQ=17)

                        gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGCGGGAGCGGGCGTGCCCGCGAAACCGGCGCCGCGGTGCACGGT‑ACCGGCTTT‑GCCGGT‑G‑TTCG‑GG‑G‑G‑‑CAC‑GCCCGCTCCCACAGGATGGTTGACAtt                                                                                                                >  2:11254792/1‑137 (MQ=17)

                        gAAGCCTTGGCACCAAAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGGGCCGCGGTGCCTGGC‑CCCGGCTTT‑GCCGGT‑G‑TT‑TGCG‑G‑G‑‑CCCTC‑CCCCTCCCACAGACTGATTGCCact                                                                                                                >  2:3046908/1‑135 (MQ=21)

                         aaGCCTTGGCACCACAGCGGGACGGAGTCCGCCAAACGGGGGGGGTGGGCTTGCCCGCGAAACAGGCGCCGCGGTGCGTGGC‑ACCGGCTTT‑GCCGGT‑G‑TTCG‑TG‑G‑G‑‑CAC‑CCCCGCTCCCACAGACTGATTGCCATTa                                                                                                               <  2:10482490/137‑1 (MQ=12)

                            ctttGGCACCACAGCAGGGTTGAGGCTGCCAAACGGTGGGAGTGGGTGTTCCCTCGAAACGGGCGCCGTGGTGCATGGC‑ACCGGCTTT‑GCCGGT‑G‑TTCG‑CG‑G‑G‑‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGc                                                                                                            <  2:4130727/135‑1 (MQ=18)

                            ccTTGACACAAAAACAGGACTGAGGCGGCAAAATGGGGGGAGCGGGTGTGCCAAAGAAACAGGCACCGCGGTGCATGAC‑ACCGGCTTG‑GCCGGT‑T‑TTCG‑CG‑G‑G‑‑CAC‑CCCCGCTCCCACAGACTGATTGCCATTAAGc                                                                                                            <  1:14084230/137‑1 (MQ=255)

                              ttGGCACCCAACCAGGACTGATCCGGCCAAATGGTGGGAGCGGGCGTACCCGCGACACATGCACCGCAGTCCATAAC‑ACCGTCTTT‑TCCGGT‑G‑TTAA‑AG‑G‑A‑‑CAT‑ACCCGCTCCCACAGACTGATTGCCATTAAGCgg                                                                                                          <  2:2015233/137‑1 (MQ=255)

                                 gCACCACAGAGGGAGTGAGGCTGCGAAGTGGGGGGAGCGGGTGTTCCTGAGAAAAGGGCGCCGTGTGACATGGC‑ACCGGTTTT‑TCCGGT‑G‑TT‑C‑GG‑G‑GCCCCC‑C‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                                                       <  1:13439395/137‑1 (MQ=255)

                                  caccacAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGTCCGCGAAACAGGCGCCGCGGTGCATGGC‑ACCGGCTTT‑GTCGGT‑G‑TTCG‑CG‑G‑G‑‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCt                                                                                                      >  2:6352587/1‑137 (MQ=21)

                                      acaGCAGGACTGAGTCTTCCAAACGGTGGGGGCGGGCGTGCCCGCGAAACAGTCGCCGCGGTGCAGGGC‑ACTGGCTTT‑GCCGGT‑G‑TT‑CGCG‑G‑G‑‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCg                                                                                                  <  2:11701272/137‑1 (MQ=18)

                                      aaagaagGACGGAGTATGCAAAACGGGGGAAGCGGACTTGCCCGCGAACCAGGCACTGCGGTACATGGC‑GCCGGCTTT‑TCCGGT‑G‑TTCG‑CG‑G‑G‑‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCg                                                                                                  <  2:7002967/132‑1 (MQ=9)

                                         gcagGACTGAGGCTGCCAAACGGTGGGAGCGGGGGTGTCCACGAAAAAAGGGCCGCGGTGCATGGCACCCGATTTT‑CCTGTT‑T‑TT‑T‑CG‑G‑G‑‑TAC‑GCCCCCTCCCCCAAACTGATTGTCCCTCCGCGGGCCTATctctc                                                                                               >  1:7370503/1‑136 (MQ=255)

                                              aCTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCCCGAAACAAGCGCCACGGTGCCTGGC‑ACCAGCCTT‑TTCAGT‑G‑TTCG‑CA‑G‑G‑‑CAC‑GCCCGCTCCCACAAAATGATTGCCACTCAGCGGGACCATCGCTCCAATg                                                                                          >  2:3006359/1‑137 (MQ=12)

                                                tGAGGCTACCAACCGGGGGGAGTGGGCGTACCCGCGAAACAGGCGCCGCGGTGCATGAC‑ACCGGCTTT‑GCTGGT‑G‑TT‑CGCG‑G‑G‑‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGAc                                                                                        <  2:13798245/137‑1 (MQ=21)

                                                    gttGCCAAACGGGGGGGGCGGGCGTTCCTGTGAAACGTGCGCCGCGGTGCGTGGC‑ACCGGTTTT‑GCCGGT‑G‑TTCG‑CG‑G‑G‑‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGcc                                                                                    <  2:11912032/135‑1 (MQ=18)

                                                     cTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGC‑ACCGGCTTT‑GCCGGT‑G‑TTCG‑CG‑G‑G‑‑AAC‑GCCCGCTCCCACAGAATGATTGCCCTTAAGCGGGCCTATCGTTAAAGTGACGGccc                                                                                   >  2:6537367/1‑137 (MQ=17)

                                                     cTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCAGTGCATGGC‑ACCGGCTTT‑GCCGGT‑G‑TTCG‑CG‑G‑G‑‑CAC‑GCCCGCTCCCACAGACTAATTGCCATTAAGCGGGCCTATCGATACAGTGACGGccc                                                                                   >  2:7879257/1‑137 (MQ=18)

                                                     cTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCACGGTGCACGGC‑ACCAGCTTT‑GCCGGT‑G‑TTCG‑CG‑G‑G‑‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGccc                                                                                   >  2:15911107/1‑137 (MQ=21)

                                                          aaaCGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGC‑ACCGGCTTT‑GCCGGT‑G‑TTCG‑AG‑G‑G‑‑CAC‑GCCCGCTCCCCCAGACTGATTGCCATTAAGCGGGCCTATCGATACAGTGACGGCCCGtttt                                                                              >  2:15119378/1‑137 (MQ=18)

                                                          aaaCGGTGGGAGCGGGCGTGCCCGCGAAACAAGAGCCACGGTGCACGGC‑ACCAGCTTT‑TCCGGT‑G‑TTCT‑CG‑G‑G‑‑CAC‑GCCCACTCCCACAAAACGATTGTCACTAAGCGGGCCTCTCTCTCCAATGACGGCCCGtttt                                                                              >  2:10640746/1‑137 (MQ=9)

                                                          aaaCGGTGGGAGCGGGCGTGCCCGCAAAACAAGAGCCGCGGTGCACGGC‑ACCGGCTTT‑GCCGGT‑G‑TTCG‑CG‑G‑G‑‑CAC‑GCCCGCTCCCCCCGACTGATTGCCATTCAGGGGGACTATCGCCCCCGTGACGGCCCGtttt                                                                              >  2:5872676/1‑137 (MQ=9)

                                                           aaCGGTGGGAGCGGGCGTGTCCGCGAAACAGGCGCCGCGGTGCATGGC‑ACCGGCTTT‑GCCGGT‑G‑TT‑TGCG‑G‑G‑‑CAC‑GACCGCTCCCAAAAAATGAGTGTCATTCAGCGGGCCCATCACTCCAGTGAAGGCCCGttttt                                                                             >  2:666725/1‑137 (MQ=17)

                                                           aaCGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGC‑ACCGGCTTT‑GCCGGT‑G‑TTCG‑CG‑G‑G‑‑CAC‑GCCCGCTCCCAAAGACTGATTGCCACTCAGCGGGCCTATCGCTCCAGTGACGGCCCGttttt                                                                             >  2:12246450/1‑137 (MQ=18)

                                                           aaCGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGC‑ACCGGCTTT‑GCCGGT‑G‑TT‑CGCG‑G‑G‑‑AAC‑GCCCACTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGttttt                                                                             >  2:8667077/1‑137 (MQ=21)

                                                                     gCGGGTGTGCCCGCGAAACAGGCGCCGCGGTGCATGGC‑ACCGGATTT‑GCCGGT‑G‑TTCG‑CG‑G‑G‑‑CAC‑GCCCGCTCCCACAGAATGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTTTTGCCTCg                                                                   >  2:5525290/1‑137 (MQ=21)

                                                                     gCGGGCGTGCCCGCGAAACAGGGGCCACGGTGCCTGGC‑ACCAGCTTT‑GCCCGT‑G‑TT‑TGCG‑G‑G‑‑AAC‑GCCCCCTCCCCCAAACTGATTGCCATTAAGCAGGCCTACCTCTCCAACGAAAGCCCATTTTTTTTTTCCTcc                                                                   >  2:9746565/1‑136 (MQ=9)

                                                                      tgggTGTGCCCGCGAAACGGTCGCCGCTGTGCATGGC‑ACCGGCTTG‑GCCGGT‑G‑TTCG‑CG‑G‑G‑‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGc                                                                  <  2:14577690/136‑1 (MQ=25)

                                                                      cGGGCTTGCCCGCGAAACGGGCGCCGCGGTACATGGC‑ACCGGTTTT‑TCTGGT‑G‑TT‑TGCG‑G‑G‑‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGc                                                                  <  2:2506023/137‑1 (MQ=25)

                                                                              cccGAGAAACAGTCGCCGCGGTGCATGGC‑ACTGGCTTT‑GCTGGT‑G‑TTCG‑CG‑G‑G‑‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAAt                                                          <  2:4914145/137‑1 (MQ=25)

                                                                                ggcgAAAAAGGCTCCGCGGGGCATGGC‑ACCGTTTTT‑GCTGGT‑G‑TT‑CGCG‑G‑G‑‑AAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                        <  2:4527606/136‑1 (MQ=255)

                                                                                ctcgAAAAAGGCGCCGTGTTGCTTTGC‑ATCGGCTTT‑GCCGGT‑G‑TT‑CGCG‑G‑G‑‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                        <  2:9073355/135‑1 (MQ=25)

                                                                                cgcgAAACAGGCGCAGCGGTGCATGGC‑ACTGGTTTT‑GCCGGT‑G‑TTCG‑CG‑G‑G‑‑CAC‑GCCCGCTCCAACAGACTGATTGCCATTAAGCGGGCCTATCGTTATAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                        <  2:13179536/137‑1 (MQ=25)

                                                                                cgcaAAACATTCGCCGCGGTGCGTGCC‑CCCGTCTTT‑TCCGTT‑G‑TT‑CGCG‑G‑A‑‑CAC‑CCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                        <  1:4407140/137‑1 (MQ=255)

                                                                                cgaaaaaCAGGCGCCGTGGTGCGGGGC‑ACTGGTTTT‑GCCGGT‑G‑TT‑TGCG‑G‑G‑‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                        <  1:200421/133‑1 (MQ=37)

                                                                                cgaaaaaCAGGCGCCGTGGTGCATGGC‑ACCGTCTTT‑GCTGGT‑G‑TTCG‑CG‑G‑A‑‑CAC‑CCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGATAACGAATCa                                                        <  1:14170624/133‑1 (MQ=255)

                                                                                 gagAAACGGGCGCCGTGGTGCATGGC‑ACTGGTTTT‑GCTGGT‑G‑TTCG‑TG‑G‑G‑‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCAt                                                       <  2:10164180/135‑1 (MQ=25)

                                                                                           cgcttcgGTGCATGGC‑ACAAGCTTT‑GCTGGT‑G‑GTCG‑TG‑G‑G‑‑CAC‑GCCCGCTCCCACAGACTGATCGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCGCGCTAACGAATCATGACTCCCTTa                                             <  2:9145510/132‑1 (MQ=25)

                                                                                           cgccgcGGGGCATGGC‑ACCGGTTTT‑GCTGGT‑G‑TTCG‑TG‑G‑G‑‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa                                             <  2:14255641/137‑1 (MQ=37)

                                                                                            gccgcGGTGCATGGC‑ACCGGCTTT‑GCCGGT‑G‑TTCG‑CG‑G‑G‑‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTA‑‑‑‑g                                        >  1:6665899/1‑137 (MQ=37)

                                                                                                cgGTGCATGGC‑ACCGGCTTT‑GCCGGT‑G‑TTCG‑CG‑G‑G‑‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTA‑‑‑‑GAACg                                    >  1:3694153/1‑137 (MQ=37)

                                                                                                cgGTGCATGGC‑ACCGGCTTT‑GCCGGT‑G‑TTCG‑CG‑G‑G‑‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTA‑‑‑‑GAACg                                    >  2:8862138/1‑137 (MQ=37)

                                                                                                cgGCGCATGCC‑ACCGGCTTT‑GCTGGTGG‑TC‑G‑CG‑G‑G‑‑CAC‑CCCCACTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTA‑‑‑‑GAACg                                    <  2:11401572/137‑1 (MQ=255)

                                                                                                     cATGGC‑ACCGGTTTT‑GCCGGT‑G‑TT‑CGCG‑G‑G‑‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTA‑‑‑‑GAACGCTATc                               <  2:6068448/137‑1 (MQ=37)

                                                                                                     cATGGC‑ACCGGCTTT‑TCTGGG‑G‑TTCT‑CG‑G‑G‑‑CAC‑TCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTA‑‑‑‑GAACGCTATc                               <  2:2669011/137‑1 (MQ=255)

                                                                                                     cATGGC‑ACCGGCTTT‑GCTGGT‑G‑TTCG‑CG‑G‑G‑‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTA‑‑‑‑GAACGCTATc                               <  2:4006194/137‑1 (MQ=37)

                                                                                                          tacACCGCTTTGGCCGGT‑GTTT‑T‑CG‑G‑G‑‑CAC‑ACCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTA‑‑‑‑GAACGCTATCAAg                            <  1:7093405/135‑1 (MQ=255)

                                                                                                            accGGCTTT‑GCCGGT‑G‑TTCG‑CG‑G‑G‑‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTA‑‑‑‑GAACGCTATCAAGCAg                         >  1:11577323/3‑137 (MQ=38)

                                                                                                            accGGCTTT‑GCCGGT‑G‑TTCG‑CG‑G‑G‑‑CAC‑GCCCGCTCCCAAAGACTGATTGCAATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAAAGAATCATGACTACCTTA‑‑‑‑GAACGCTATCAAGCAg                         >  2:15548631/3‑137 (MQ=37)

                                                                                                            accGGCTTT‑GCCGGT‑G‑TT‑CGCG‑G‑G‑‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCAGTTTTTTGTTGCATCGCTAACGAATCATGACTCCCTTA‑‑‑‑GAACGCTATCAAGCAg                         >  2:2458510/3‑137 (MQ=38)

                                                                                                               ggCTTT‑GCCGGT‑G‑TTCG‑CG‑G‑G‑‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTA‑‑‑‑GAACGCTATCAAGCAGATc                      <  2:13010545/137‑1 (MQ=39)

                                                                                                                 cTTT‑GCCGGT‑G‑TTCG‑CG‑G‑G‑‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGCATCATGACTCCCTTA‑‑‑‑GAACGCTATCAAGCAGATCTg                    >  1:6305201/1‑137 (MQ=38)

                                                                                                                 cTTT‑GCCGGT‑G‑TTCG‑CG‑G‑G‑‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTA‑‑‑‑GAACGCTATCAAGCAGATCTg                    >  1:258492/1‑137 (MQ=39)

                                                                                                                 cTTT‑GCCGGT‑G‑TTCG‑CG‑G‑G‑‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTA‑‑‑‑GAACGCTATCAAGCAGATCTg                    >  2:5548426/1‑137 (MQ=39)

                                                                                                                 cTTT‑GCCGGT‑G‑TTCG‑CG‑G‑G‑‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTA‑‑‑‑GAACGCTATCAAGCAGATCTg                    >  1:322314/1‑137 (MQ=39)

                                                                                                                 cTTT‑GCCGGT‑G‑TTCG‑CG‑G‑G‑‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTA‑‑‑‑GAACGCTATCAAGCAGATCTg                    >  1:12761366/1‑137 (MQ=39)

                                                                                                                 cTTT‑GCCGGT‑G‑TTCG‑CG‑G‑G‑‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTA‑‑‑‑GAACGCTATCAAGCAGATCTg                    >  2:7612788/1‑137 (MQ=39)

                                                                                                                 cTTT‑GCCGGT‑G‑TTCG‑CG‑G‑G‑‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTA‑‑‑‑GAACGCTATCAAGCAGATCTg                    >  2:14893039/1‑137 (MQ=39)

                                                                                                                 cTTT‑GCCGGT‑G‑TTCG‑CG‑G‑G‑‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTA‑‑‑‑GAACGCTATCAAGCAGATCTg                    >  1:13673630/1‑137 (MQ=39)

                                                                                                                 cTTT‑GCCGGT‑G‑TTCG‑CG‑G‑G‑‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTA‑‑‑‑GAACGCTATCAAGCAGATCTg                    >  1:12269330/1‑137 (MQ=39)

                                                                                                                 cTTT‑GCCGGT‑G‑TTCG‑CG‑G‑G‑‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTA‑‑‑‑GAACGCTATCAAGCAGATCTg                    >  1:13019444/1‑137 (MQ=39)

                                                                                                                 cTTT‑GCCGGT‑G‑TTCG‑CG‑G‑G‑‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTA‑‑‑‑GAACGCTATCAAGCAGATCTg                    >  2:10288156/1‑137 (MQ=39)

                                                                                                                 cTTT‑GCCGGT‑G‑TTCG‑CG‑G‑G‑‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTA‑‑‑‑GAACGCTATCAAGCAGATCTg                    >  2:13143096/1‑137 (MQ=39)

                                                                                                                 cTTT‑GCCGGT‑G‑TTCG‑CG‑G‑G‑‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTA‑‑‑‑GAACGCTATCAAGCAGATCTg                    >  2:12722608/1‑137 (MQ=39)

                                                                                                                 cTTT‑GCCGGT‑G‑TTCG‑CG‑G‑G‑‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTA‑‑‑‑GAACGCTATCAAGCAGATCTg                    >  2:10393720/1‑137 (MQ=39)

                                                                                                                 cTTT‑GCCGGT‑G‑TTCG‑CG‑G‑G‑‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTA‑‑‑‑GAACGCTATCAAGCAGATCTg                    >  2:4435138/1‑137 (MQ=39)

                                                                                                                 cTTT‑GCCGGT‑G‑TTCG‑CG‑G‑G‑‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTA‑‑‑‑GAACGCTATCAAGCAGATCTg                    >  2:3497560/1‑137 (MQ=39)

                                                                                                                 cTTT‑GCCGGT‑G‑TTCG‑CG‑G‑G‑‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTA‑‑‑‑GAACGCTATCAAGCAGATCTg                    >  1:9799626/1‑137 (MQ=39)

                                                                                                                 cTTT‑GCCGGT‑G‑TTCG‑CG‑G‑G‑‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCCACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTA‑‑‑‑GAACGCTATCAAGCAGATCTg                    >  2:16350021/1‑137 (MQ=38)

                                                                                                                 cTTT‑GCCGGT‑G‑TT‑CGCG‑G‑G‑‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTCTCGCTACAGTAACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTA‑‑‑‑GAACGCTATCAAGCAGATCTg                    >  1:5551501/1‑137 (MQ=38)

                                                                                                                   tt‑GCCGGT‑G‑TTCG‑CG‑G‑G‑‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTA‑‑‑‑GAACGCTATCAAGCAGATCTGaa                  >  1:8845510/1‑137 (MQ=39)

                                                                                                                   tt‑GCCGGT‑G‑TTCG‑CG‑G‑G‑‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTA‑‑‑‑GAACGCTATCAAGCAGATCTGaa                  >  1:4015640/1‑137 (MQ=39)

                                                                                                                   tt‑GCCGGT‑G‑TTCG‑CG‑G‑G‑‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTA‑‑‑‑GAACGCTATCAAGCAGATCTGaa                  >  1:7762074/1‑137 (MQ=39)

                                                                                                                   tt‑GCCGGT‑G‑TTCG‑CG‑G‑G‑‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTA‑‑‑‑GAACGCTATCAAGCAGATCTGaa                  >  2:1676835/1‑137 (MQ=39)

                                                                                                                   tt‑GCCGGT‑G‑TTCG‑CG‑G‑G‑‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTA‑‑‑‑GAACGCTATCAAGCAGATCTGaa                  >  1:6024710/1‑137 (MQ=39)

                                                                                                                   tt‑GCCGGT‑G‑TTCG‑CG‑G‑G‑‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACAAATCATGACTCCCTTAGAACGAACGCTATCAAGCAGATc                      >  1:15612728/1‑137 (MQ=37)

                                                                                                                              ttttcgcG‑G‑G‑‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTA‑‑‑‑GAACGCTATCAAGCAGATCTGAAACGTccc           <  1:2533061/131‑1 (MQ=255)

                                                                                                                                 tcgcG‑G‑G‑‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTA‑‑‑‑GAACGCTATCAAGCAGATCTGAAACGTCCCGAc        >  2:9817597/4‑137 (MQ=255)

                                                                                                                                 tcgcG‑G‑G‑‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTA‑‑‑‑GAACGCTATCAAGCAGATCTGAAACGTCCCGAc        >  1:16062923/4‑137 (MQ=255)

                                                                                                                                  cgcG‑G‑G‑‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTA‑‑‑‑GAACGCTATCAAGCAGATCTGAAACGTCCCGAct       >  1:15515814/3‑137 (MQ=255)

                                                                                                                                     g‑g‑g‑‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTA‑‑‑‑GAACGCTATCAAGCAGATCTGAAACGTCCCGActtct    >  2:2424270/1‑137 (MQ=255)

                                                                                                                                     g‑g‑g‑‑CAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTA‑‑‑‑GAACGCTATCAAGCAGATCTGAAACGTCCCGActtct    >  2:16399179/1‑137 (MQ=255)

GATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGC‑CACGGCTTT‑GCCGGT‑G‑TT‑T‑CG‑G‑G‑‑CAC‑G‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTA‑‑‑‑GAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTC  >  AE015451_phaM1/1499373‑1499637

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