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	91.2%	474.8 / ‑4.3	240	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 (17/4); new base A (115/102); total (132/108)
Fisher's exact test for biased strand distribution p-value = 1.97e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*	AE015451_phaM1	1,499,481	0	A	C	93.6%	527.2 / ‑0.6	240	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 (12/3); new base C (120/99); total (132/108)
Fisher's exact test for biased strand distribution p-value = 6.41e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 240 total aligned reads displayed.

GTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTG‑CCCGCGAAACAGGCGCCGCGGTGCA‑TGGCCACGGCTTTGCCGGT‑GTT‑T‑‑CGGG‑CA‑C‑G‑CC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATC  >  AE015451_phaM1/1499351‑1499608

||

gTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTG‑CCCGCGAAACAGGCGCCGCGGTGCA‑TGGCACCGGCtt                                                                                                                                    <  2:14726836/137‑1 (MQ=255)

gTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTG‑CCCGCGAAACAGACGCCGCGGTGCA‑TGGCACCGGCtt                                                                                                                                    >  2:13677157/1‑137 (MQ=255)

  tCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTG‑CCCGCGAAACAGGCGCCGCGGTGCA‑TGGCACCGGCTTTg                                                                                                                                  <  1:7923167/137‑1 (MQ=255)

  tCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTG‑CCCGCGAAACAGGCGCCGCGGTGCA‑TGGCACCGGCTTTg                                                                                                                                  <  2:9093304/137‑1 (MQ=255)

  tCGACAACGAGGAAGCCCGTGATGCGGCTACTTAGAGCCTGAACGAAGCCTTGGCACCACGGCAGGACGGAGGCTGCTAAACGGTGGGAGCGGGCGTG‑CCCGAGAAACAGGTGTCGCGGTGCA‑TGGCGCCGGcggtg                                                                                                                                  <  2:7328666/137‑5 (MQ=25)

       aaCGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTG‑CCCACGAAACAGGCGCCGCGGTGTA‑TGGCACCGGCTTTGCCGgt                                                                                                                             >  2:9731538/1‑137 (MQ=255)

            tGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGGGTG‑CCCCCGAAAAAAGAGCCACGGTGCA‑CGGAACCAGCTTTTCCGGT‑Gtt‑tg                                                                                                                      >  1:5626501/1‑136 (MQ=255)

            tGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGGGTG‑ACCACGAAAAAGGCGCCACAGTGTA‑TTGCCCCGGCTTTTTCGGT‑ttt‑tt                                                                                                                      >  1:8514989/1‑136 (MQ=255)

            tGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTG‑CCCGTGAAACAGGCGCCGCGGTGCA‑TGGCACCAGCTTTTCCAGT‑Gttca                                                                                                                       >  2:5775966/1‑135 (MQ=255)

              aaGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCTTG‑CCCGCGAAACAAGTGCCGCAGTGCA‑TGGCACCGGATTTTCCGGT‑Gttcgcg                                                                                                                     >  1:12341991/1‑133 (MQ=37)

                      gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGGGTG‑CCCGCGAAAAAAGAGCCACGGTGCC‑CGGCACCGGCTTTGCCAGT‑GTT‑TG‑CGGG‑CA‑c‑g‑cc                                                                                                       >  1:16450621/1‑137 (MQ=255)

                      gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTG‑CCCGCGAAACAGGCGCCACGGTGCA‑TGGCACCGGCGTTGCCGGT‑GTT‑TG‑CGGG‑CA‑c‑g‑cc                                                                                                       >  2:13032706/1‑137 (MQ=37)

                      gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTG‑CCCGCAAAACAAGCGCCACGGTGCA‑TGGCACCGGTTTTTCCGGT‑TTTCC‑‑CGGG‑CA‑c‑a‑cc                                                                                                       >  1:8776118/1‑137 (MQ=25)

                      gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGAGTG‑ACCGCGAAACAGGAGCCACAGTGCA‑CGGCACCAAATTTTCCGGT‑GTTCT‑‑CAGG‑CA‑c‑g‑cc                                                                                                       >  1:13290437/1‑137 (MQ=255)

                       aTGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTG‑CCCGCGAAACAGGCGCCGCGGTGCA‑TGGCACCGGCTTTGCCGGT‑GTTCG‑‑CGGG‑CA‑c‑g‑ccc                                                                                                      >  2:14191680/1‑136 (MQ=25)

                        tGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTG‑CCCTCGAAACAAGCGCCGCGGTGCA‑CGGCACCAGATTTGCCAGT‑GTTCG‑‑CGGG‑CA‑c‑g‑cccg                                                                                                     >  1:14786143/1‑135 (MQ=255)

                              cATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGGGTG‑TCCGCGAAACAAGCGCCACGGTGCA‑CGGCACCAAAGTTGCCTGT‑GTTCT‑‑CGGG‑TA‑C‑GCCC‑GCTcccc                                                                                               >  1:3614345/1‑136 (MQ=255)

                               aTTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGTGTG‑TCCGCGAAAAAGGCGCCGCGGGGCA‑TGGCACCGGATTTGACAGT‑GTTCG‑‑CGGG‑CA‑C‑GCCC‑GCTCCcaa                                                                                              >  2:2557943/1‑136 (MQ=25)

                               aTTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGAGTG‑CCCGCGAAACAGGCGCCGCGGTGCA‑TGGCACCGGCTTTGCCGGT‑GTTCG‑‑CGGG‑CACC‑C‑CC‑GCTCCcaa                                                                                              >  2:1979024/1‑136 (MQ=21)

                                  cAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTG‑CCCGCGAAACAGGCGCCGCGGTGCC‑TGGCACCGGATTTGCCGGT‑GTTCG‑‑AGGG‑CA‑C‑GCCC‑GCTCCCACAGa                                                                                           >  1:13817184/1‑137 (MQ=21)

                                        cTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTG‑CCCGCGAAACAGGCGCCGCGGTGCA‑TGGCACCGGCTTTTCCGGT‑GTTCG‑‑CGGG‑CA‑C‑GCCC‑GCTCCCACAGACTGAtt                                                                                     >  2:13396118/1‑137 (MQ=21)

                                         tGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTG‑ACCGCGAAACAGGCGCCACGGTGCA‑CGGCACCGGATTTGGCGGT‑GTTCG‑‑AGGG‑CA‑T‑GCCC‑GCTCCCACAGAATGATTg                                                                                    >  2:7349217/1‑137 (MQ=21)

                                          gAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTG‑CCCGCGAAACAGGCGCCGCGGTGCA‑TGGCAACGGATTTGCCGGT‑GTT‑CG‑CTGG‑CA‑C‑G‑CCAGCTCCCAAAGAATGATTGc                                                                                   >  2:236899/1‑137 (MQ=21)

                                          gAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTG‑CCCAGGAAACAGGCGCCGCGGTGCA‑TGGCAACGGTCTTGCCGGT‑GTTCG‑‑CGGG‑CA‑C‑GCCC‑GATCCCAAAGACTGATTga                                                                                   >  2:16402286/1‑136 (MQ=25)

                                             cGAAGCTTTGGCACCACAGCAGGACTGAGGCTACCAAATGGTGGGGGCGGGCGTG‑CCCGCGAAACAGGCGCTGCGGTGCA‑TGGCACCGGTTTTGCCGGT‑GTT‑CG‑CGGG‑CA‑C‑GCCC‑GCTCCCACAGACTGATTGCCAt                                                                                <  2:16539234/137‑1 (MQ=18)

                                             cGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTG‑CCCGCGGAACAGGAGCCGAGGTGCC‑TGGCACCGGCTTTGCCGGT‑GTTCC‑‑CGGG‑CA‑C‑GCCC‑GCTCCCACAGAATGATTGCCAt                                                                                >  2:8047876/1‑137 (MQ=21)

                                             cGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTG‑CCCGCGAAACAAGCGCCGCGGTGTA‑TGGCACCGGCTTTTCCAGT‑GTT‑TG‑CGGG‑CA‑C‑GCCC‑GCTCCCACAAACCGATTGCCac                                                                                >  2:667958/1‑136 (MQ=18)

                                             cGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTG‑ACCGCGAAAAAAGCGCCACGGTGTC‑TGGCCACGGCTTTGCCGGT‑GTTCG‑‑CGGG‑CA‑CAC‑CC‑GCTCCCACAGAATGGTTGCCAt                                                                                >  1:4027971/1‑137 (MQ=21)

                                              gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAATGGTGGGAGCGGGCGTG‑TCCGCGAAAAAGGCGCCAAAGTGCA‑TGGCCACGGATTTGACAGT‑GTTCA‑‑CGGG‑GA‑C‑GCCC‑ACTTCCACAGGATGATTGCCAtt                                                                               >  1:18713819/1‑137 (MQ=255)

                                              gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTG‑TCCGCAAAAAAAGCGCCGCAGTGCA‑TGGAACCAGCTTTGCCGGT‑GTTCA‑‑CGGG‑CA‑C‑GCCC‑GCCCCCACAAACCGATTGCCAtt                                                                               >  2:3068787/1‑137 (MQ=255)

                                              gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTG‑TCCACAAAACAGGCGCCACAGTGCC‑CTGCAACAGCCTTTCCGGT‑GTTCT‑‑CGGGACC‑C‑C‑CC‑GCCCCCCCAGAACGATTGcccct                                                                               >  1:3128091/1‑134 (MQ=255)

                                              gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTG‑GCCGCGAAACAAGCACCGCGGTGCA‑TGGCAACAGATTTGCCGGT‑GTTCT‑‑CGGG‑CA‑C‑GACC‑ACTCCCACAGAATGACTGCCagt                                                                               >  1:17589470/1‑135 (MQ=255)

                                              gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTG‑CCCGCGAAACAGGCGCCGCGGTGTA‑TGGCACCGGCTTTGCCGGT‑GTT‑TG‑CGGG‑CA‑C‑GCCC‑GCTCCCACAGACTGATTGCCAtt                                                                               >  2:15096648/1‑137 (MQ=21)

                                              gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTG‑CCCGCGAAACAGGCGCCGCAGTGCA‑TGGCACCGGCTTTGCCGGTGGTC‑G‑‑CGGG‑CA‑C‑GACC‑GCTCACAAAGAATGAGTGACAtt                                                                               >  2:4689714/1‑137 (MQ=18)

                                              gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTG‑CCCGCAAAACAGGCGCCGCAGTGTA‑TGGCCCCGGCTTTGCCGGT‑GTT‑TG‑TGGG‑CA‑C‑GCCC‑GCTCCCACCGACTGATTTCCAtt                                                                               >  1:10128531/1‑137 (MQ=25)

                                              gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTG‑CCCGAAAAACAAGCGCCGCGGTGCA‑TGGCACCAGATTTTCCGGT‑GTTCT‑‑CGGG‑CA‑C‑GCCC‑ACCCCCACAAAATGATTGCCAtt                                                                               >  1:12761755/1‑137 (MQ=18)

                                              aaatccTTGGCACAACAGAAGGACGGAGGCGACCAACAGGTGGGAGCGGGCTTA‑CCCCAGAAACAGGCGCCGCGGGGAA‑TGGCACCGGCTTGGCCGGG‑GTTCA‑‑CGGG‑CA‑C‑GCCC‑GCTCCCACAGACTGATTGCCAtt                                                                               <  2:4332713/133‑1 (MQ=9)

                                                  ccTTTGCACCACATAAGGCCGGAGGCTCCCAAACGGGGGAAGCGGTCGTG‑CCCGCGAAACAGTCTCCGCGGTGCA‑GGGAACCGGCTTTTCC‑GT‑GTT‑TCGCGGG‑CA‑C‑TCCC‑GCTCCCACAGACTGATTGCCATTAAGc                                                                           <  2:4246188/137‑1 (MQ=255)

                                                  ccTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTG‑CCCGCGCAACAGGCGCCGCTGTGTA‑TGGCCCCGGCTTTGCCAGT‑GTTCT‑‑CGGG‑CA‑CAC‑CC‑GCTCCCAAAGCCTGCTTTCCATTaaac                                                                           >  2:1080183/1‑135 (MQ=21)

                                                  ccTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTG‑CCCGCGAAAAAGGAGCCGCAGGGCA‑GGGCAACGGATTTGCCGGT‑GTTCG‑‑CGGG‑CA‑C‑GCCC‑GCTCCCCCAGACTGATTGCCATTAAga                                                                           >  2:3359670/1‑136 (MQ=21)

                                                   cTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTG‑CCCGCGAAACAGGGGCCCCGGTGTA‑TGGCCCCTGATTTGCCAGT‑GTTCG‑‑CTGG‑CA‑C‑GCCC‑GCTCCCCCAGACTGACTGCCCttcaaca                                                                          >  2:12122412/1‑132 (MQ=16)

                                                    ttgtccccacATCAGGATTGAGTTTGCCGAATGGGGGGATTGGGTGTC‑CCTGCGAAACAGGCGCTGCGGTACA‑TGGCACCGGCTTTTCTGGT‑TTTCG‑‑CGGG‑CA‑C‑TCCC‑GCTCCCACAGACTGATTGCCATTAAGCgg                                                                         <  1:17111333/131‑1 (MQ=255)

                                                    ttGGCACCACAGCAGGACTGAGGCTGCCAAATGTTGTGAGCGGGTGTT‑CCCGCGAAGCAGGTGCCGTGGTACA‑TGGCACCGGTTTTGCCGGT‑GTT‑CG‑CGGG‑CA‑C‑TCCC‑GCTCCCACAGACTGATTGCCATTAAGCgg                                                                         <  1:13493050/137‑1 (MQ=16)

                                                    ttGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTG‑CCCGCGAAACAGGCGCCGCGGTGCA‑TGGCACCGGCTTTGCCGGT‑GTT‑CG‑CGGG‑CA‑C‑GCCC‑GCTCCCACAGACTGATTGCCATTCAGCgg                                                                         >  2:3511199/1‑137 (MQ=18)

                                                    ttGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTG‑CCCGCGAAACAAGCGCCGCGGTGCA‑TGGCAACGGCTTTGCCGGT‑GTT‑CG‑AGGG‑CA‑C‑GCCC‑GATTCCAAAGAATGATTGCAATTAAGCgg                                                                         >  2:17827242/1‑137 (MQ=18)

                                                    tgggCACCACAACAGGACGGAGTCTACCAAATGGTGGGAGTGGGCGTG‑CCCGCGAAACATGCGCCGCGGTGCA‑TGGCACTGGTTTTGCTGTT‑GTT‑CG‑TGGG‑CA‑C‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCgg                                                                         <  1:11569460/135‑1 (MQ=12)

                                                       tcacaacAGCGGGACTGGGGCTACAAAATGGGGGGAGTGGGCGTG‑CCCGTGAAACGGTCGCTGTGGGGCA‑TGGCACTGGTTTTGCTGGT‑GTT‑GG‑CGGG‑CA‑C‑CCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                      <  1:15597757/136‑1 (MQ=255)

                                                       gCACCACAGCAGGACTGGGGCTGCAAAATGGGGGGAGCGGGTGTG‑CCCGCGAAACGGGCTCCGCGGTGCA‑GGGCACTGGCTTTGCCGGT‑GTT‑CG‑CGGA‑CA‑C‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                      <  2:1246081/137‑1 (MQ=18)

                                                       gCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGTGTG‑CCCGCGAAACAGGCACCGCGGTGCA‑TGGCACCGGCTTTGCCGGT‑GTT‑CG‑CGGG‑CA‑C‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                      <  2:1703177/137‑1 (MQ=21)

                                                       gCACCAAAGCAGGATGGAGGTTGCAAAATGGTGGGAGCGGGTGTG‑CCCGAGAAACAGGCGCCGCGGTGCA‑TGGCACCGTTTTTGCTGGT‑GTTCG‑‑CGGG‑CA‑C‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                      <  2:16235435/137‑1 (MQ=12)

                                                       ccaccaaAGAAGGATGGGGGTTGCAAACAGGTGGGAGAGGGTGTG‑CCCGAGAAACGGTTGCCGTGGGGCG‑GGGAACCGGTTTTTCCGGG‑GTT‑TG‑TGGG‑CA‑C‑CCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                      <  1:12478977/136‑1 (MQ=255)

                                                        caccacAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTG‑CCCGCGAAACAGGCGCCGCGGTGCA‑TGGCACCGGCTTTGCCGGT‑GTTCG‑‑CGGG‑CA‑C‑GCCC‑GCTCCCAAAGACTGATTGCCATTAAGCGGGCCt                                                                     >  2:18609087/1‑137 (MQ=21)

                                                        caccacAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTG‑CCCGCGAAACAGGCGCCGCGGTGCA‑TGGCAACGGCTTTGCCGGT‑GTT‑CG‑CGGG‑CC‑C‑GCCC‑GCTCCCAAAGACTGATTGACATTAAGCGGGCCt                                                                     >  2:17677708/1‑137 (MQ=21)

                                                        caccacAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTG‑CCCGCGAAAAAAGCGCCGCAGTGCA‑CGGCACCGGCTTTTCCGGT‑GTTCG‑‑CGGG‑CA‑C‑GCCC‑GCTCCCACAGACTGATTGCCATTAAACGGGGCt                                                                     >  2:13256619/1‑137 (MQ=17)

                                                        cacaacAGCGGGATTGGTGCTGAAAAGGGGTGGGAGAGGGTGTG‑CCCGAGAAACGGGCGCCGTGGGGCA‑TGGCACTGGTTTTGCTGGT‑GTTCG‑‑TGGG‑CC‑C‑TCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCt                                                                     <  1:11055775/137‑1 (MQ=255)

                                                               gcagGATGGGGGCTACAAAATGGTGGGAGGGGGTGTG‑CCCGAAAAAGGGGCGCCGTGGTGCG‑GGGCATTGGTTTTTTTGGT‑GTTCG‑‑TGGA‑CA‑C‑CCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTa                                                              <  2:1730777/137‑1 (MQ=255)

                                                               gcagGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTG‑CCCGCGAAACAGGCACCGCGGTGCA‑TGGCACCGGCTTTGCCGGT‑GTTCG‑‑CGGG‑CA‑C‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCtc                                                              >  2:3700121/1‑136 (MQ=21)

                                                               gcagGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTG‑CCCGCGAAAAAGGCGCCGCGGTGCA‑TGGCACCGGCTTTGCCGGT‑GTTCG‑‑CGGG‑CA‑C‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTa                                                              >  2:8330667/1‑137 (MQ=21)

                                                                cagGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTG‑CCCGCGAAACAGGCGCCGCGGTGCA‑TGGCACCGGCTTTGCCGGT‑GTTCG‑‑CGGG‑CA‑C‑GCCC‑GCTCCCACAGAATGATTGCCATTAAGCGGGCCTATCGCTAc                                                             >  2:11235976/1‑137 (MQ=18)

                                                                cagGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTG‑CCCGCGAAACAGGCGCCGCGGTGCA‑TGGCACCAGCTTTGCCGGT‑GTTCG‑‑CGGG‑CA‑C‑GCCC‑GCTCCCACAGACTGAGTGCCATTCAGCGGGCCTATCGCTAc                                                             >  2:15719263/1‑137 (MQ=18)

                                                                     cTGAGGCTGCCAAACGGTGGGAGCGGGCGTG‑CCCACGAAACAAGCGCCACGGTGCA‑TGGCACCGGCTTTGCCGGT‑GTT‑TG‑CGGG‑CA‑C‑GCCC‑GCTCCCACAAACTGATTGCCATTAAGCGGGCCTCTCGCTACAGTGa                                                        >  2:497398/1‑137 (MQ=21)

                                                                        aGGCTGCCAAACGGTGGGAGCGGGCGTG‑CCCGCGAAAAAAGCGCCGCTGTGCA‑TGGCACCGGCTTTGCCGGT‑GTTCG‑‑CGGG‑CC‑C‑GACC‑GATCCCAAAGACCGCTCTCCATTAAACGGGCACATCGCTACAATGACgg                                                     >  2:18264003/1‑137 (MQ=9)

                                                                        aGGCTGCCAAACGGTGGGAGCGGGAGTG‑CCCGCGAAACAGGAGCCACGGTGCA‑TGGCCCCCGCTTTTCCGGT‑TTTCG‑‑CCGG‑GA‑C‑GCCC‑GCTCCCACAGAACGGGTGCCATTCAGCAGGACTAGCGCTACCGTGAcag                                                     >  1:7474550/1‑135 (MQ=255)

                                                                           cTGCCAAACGGTGGGAGCGGGCGTGCCCCGCGAAACAGGCGCCGCGGTGCA‑TGGCACCGGCTTTGCCGGT‑GTTCG‑‑CGGG‑CA‑C‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGcc                                                   >  2:11570827/1‑137 (MQ=21)

                                                                           cTGCCAAACGGTGGGAGCGGGCGTG‑CCCGTGAAAAAAGCGCCACGGCGCA‑TGGCACCGGCTTTGCCTGT‑GTTCG‑‑CGGG‑CA‑C‑GCCC‑GCTCCCCCAGACTGAGTGCCATTCAGCGGGCCTCTCGCTACAATGACGGacc                                                  >  2:2656877/1‑137 (MQ=12)

                                                                           cTGCCAAACGGTGGGAGCGGGCGTG‑CCCGCGAAACCAGCGCCCCGGTGCA‑CGGCACCCGCTTTTCCGGT‑GTTCT‑‑CTGG‑TA‑C‑GCCC‑ACTCCCCCAGACCGATTGCCATTCAGCGGGCCTATCGCTACAGTGACGGacc                                                  >  2:16304470/1‑137 (MQ=16)

                                                                           cTGCCAAACGGTGGGAGCGGGCGTG‑CCCGCGAAACAGGCGCCGCGGTGCA‑TGGCACCGGCTTTGCCAGT‑GTT‑TG‑CGGG‑CA‑A‑GCCC‑GCTCCCACAGAATGATTGCCATTAAGCGGGACTCTCGCTAAAATGaaggacc                                                  >  2:16887610/1‑131 (MQ=17)

                                                                           cTGCCAAACGGTGGGAGCGGGCGTG‑CCCGCGAAACAGGCGCCACGGGGCA‑TGGCACCGGCTTTGACGGT‑GTTCG‑‑CAGG‑CA‑C‑GCCC‑GCTCCCACAGACTGATTGCCATTCAGCGGGACTATCGCTACAGTGGCGGtcc                                                  >  2:12641295/1‑137 (MQ=18)

                                                                           cTGCCAAACGGTGGGAGCGGGCGTG‑CCCGCGAAACAAGCGCCGCAGTGCA‑TGGCACCAGCTTTGCCGGT‑GTTCG‑‑CGGG‑CA‑C‑GCCC‑GCTCCCACAGACTGATTGCACTTCCGCGGGCCTACCGCTACAGTGACGGcac                                                  >  2:18149894/1‑135 (MQ=14)

                                                                           cTGCCAAACGGTGGGAGCGGGCGTG‑CCCGCGAAAAAGGCGCCGTGGTGCA‑TGGCACCGGCTTTGCCGGT‑GTTCG‑‑CGGG‑TA‑A‑GCCC‑GCTCCCACAAGCTGATTGCCATTAAGCGGGCCTATCACCACAGTGACAGccc                                                  >  2:11018482/1‑137 (MQ=17)

                                                                           cTGCCAAACGGTGGGAGCGGGCGTG‑CCCGCAAAACAGGCGCCACGGTGCA‑TGGCACCGGCTTTGCCGGT‑GTTCG‑‑TGGG‑CA‑C‑GCCC‑GCTCCCACAGACTGATTGCCCTTAAGAGGGCATATCGCTAAAGTGACGGccc                                                  >  2:18249873/1‑137 (MQ=18)

                                                                           cTGCCAAACGGTGGGAGCGGGCGTG‑CCCACGAAACAGGCACCACGGTGCA‑TGGCACCGGATTTGCCGGT‑GTTCA‑‑CGGG‑CA‑A‑GACC‑CATCCCACAAAATGATTGCCATTAAGCGGGCCTATCGCTCCAATGACGGccc                                                  >  2:2116738/1‑137 (MQ=255)

                                                                           cTGCCAAACGGTGGGAGCGGGAGTG‑CCCACAAAACAAGCGCCACGGTGCA‑CGGCCCCAGCTTTGCCGGT‑GTTCG‑‑CGGG‑CA‑C‑GCCC‑ACTCCCACAGACTGATTGCCCCTCAGCGGGCCTCCCACTCCCGTGACGGccc                                                  >  2:731398/1‑137 (MQ=12)

                                                                             gCCAAACGGTGGGAGCGGGCGTG‑CCCGAGAAACAGGCGCCGCGGTGCA‑TGGCACCGGCTTTGCCGGT‑GTTCG‑‑CGGG‑CA‑C‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGt                                                >  1:13120974/1‑137 (MQ=21)

                                                                                aaaCGGTGGGAGCGGGCGTG‑CCCGCGAAACAAGTGTCGCGGTGCC‑TGGCACCGGCTTTTCCTGT‑GTTCG‑‑CGGG‑CC‑C‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGACTCTCGCTACCGTGACGGCCCGtttt                                             >  2:7967727/1‑137 (MQ=17)

                                                                                aaaCGGTGGGAGCGGGCGTG‑CCCGCGAAAAAAGCGCCGCGGTGCA‑TGGCACCGGATTTGCCAGT‑GTTCG‑‑CGGG‑CA‑CAC‑CC‑ACTCCCACAGAACGATTGCCATTAAGCGGGCCTATCGCTCCAGTGGCGGCCCGtttt                                             >  1:14160223/1‑137 (MQ=17)

                                                                                aaaCGGTGGGAGCGGGCGTG‑CCCACGAAACAGGCGCCGCGGTGCA‑TGGCACCAGCTTTGTCGGT‑GTTCG‑‑CGGG‑CA‑A‑GCCC‑GCTCCCACAAAACGATTGCCATTCCACGGGCCTATCGCTCCAATGACGGCCCGgttt                                             >  2:12708698/1‑137 (MQ=9)

                                                                                 aaCGGTGGGAGCGGGCGTG‑TCCGCAAAACAGGCGCCGCGGTGCA‑TGGCACCGGATTTGTCGGT‑GTTCG‑‑CGGG‑CC‑C‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGGCTATCGCTACAGTGACGGCACGttttt                                            >  2:14991187/1‑137 (MQ=255)

                                                                                 aaCGGTGGGAGCGGGCGTG‑CCCGCGAAACAAGCGCCACGGTGTC‑TGGCACCGGCTTTTCCGGT‑GTTCT‑‑CGGG‑CA‑CAC‑CC‑GCTCCCCCCAACCGATTGCCCTTCAGCGGGTCTCTCGCTCCAATGACGGCCCCttttt                                            >  2:16812437/1‑137 (MQ=9)

                                                                                           gCGGGCGTG‑CCCGCGAAACAAGCACCGCGGTGCACTGCCCCCGCCTTTCCCGTT‑TTC‑T‑‑CGGG‑CA‑CAC‑CC‑CCTCCCACAGAACGACTGTCCCTCAGCGGGCCTATCGTTCCAGTGACGGCCCTTTTTTTTTTTCCTcc                                  >  2:533098/1‑136 (MQ=12)

                                                                                                   tccccGCGAAAAAGGCGCCGCGGTGCA‑TGGCACCTGCTTTGCCGGT‑GTTCG‑‑TGGG‑CA‑C‑CCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGa                           <  1:10522407/135‑1 (MQ=255)

                                                                                                       cgcgAAACGGGTGCCGTGGTGCA‑TGGCACTGGTTTTTCTGGT‑GTTCG‑‑CGGG‑CA‑C‑TCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                       <  2:14419425/137‑1 (MQ=255)

                                                                                                       cgcgAAACAGGCGCCGCGGTGCA‑TGGCACCGGTTTTTCCGGT‑GTT‑CG‑CGGG‑CA‑C‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                       <  2:13788660/137‑1 (MQ=25)

                                                                                                       cgcgAAACAGGCGCCGCGGTGCA‑TGGCACCGGCTTTGCTGGG‑GTTCG‑‑CGGG‑CA‑C‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                       <  2:13772008/137‑1 (MQ=25)

                                                                                                       cgcgAAACAGGCGCCGCGGTACA‑TGACACTGGCTTTACTGGT‑GTTCG‑‑CGGG‑CA‑C‑CCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                       <  1:18942267/137‑1 (MQ=255)

                                                                                                       cgcgAAACAGGCCCCGCGGTACG‑TGGCACCGGCTTTGCCGGT‑GTT‑CG‑CGGG‑CA‑C‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                       <  2:3914407/137‑1 (MQ=25)

                                                                                                       cgaaaaaCAGGCGCCGCGGTACA‑TGGCACCGGCTTTGCCGGT‑GTTAG‑‑CGGG‑CA‑C‑ACCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                       <  2:326564/133‑1 (MQ=37)

                                                                                                       agaaaaaCGGGCGCCGCGGGAAA‑TGGCACTGGTTTTGCTGGT‑GTTCG‑‑TGGG‑CA‑C‑ACCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                       <  1:6544651/133‑1 (MQ=255)

                                                                                                        gcgAAACAGGCTCCGCGGTGCA‑TGGCACCTGCTTTGCCGGT‑GTT‑CG‑CGGG‑CA‑C‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCAt                      <  2:7213849/137‑1 (MQ=25)

                                                                                                        gcgAAACAGGCGCCGGGGTGCG‑TGGCACCGGTTTTGCTGGG‑GTTCG‑‑TGGG‑CA‑C‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCAt                      <  1:8970433/137‑1 (MQ=25)

                                                                                                          aaaaaAGGCTCCGCGGTGCA‑TGACATCGGTTTTGCCGGT‑TTTCG‑‑CGGC‑CA‑C‑CCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGa                    <  1:16885027/136‑1 (MQ=255)

                                                                                                             aCAGGCGCCGCGGTGCA‑TGGCACCGGCTTTGCCGGT‑GTTCG‑‑CGGG‑CA‑C‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTc                 >  2:14101033/1‑137 (MQ=37)

                                                                                                             aCAGGCGCCGCGGTGCA‑TGGCACCGGCTTTGCCGGT‑GTTCG‑‑CGGG‑CA‑C‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTc                 >  2:8713361/1‑137 (MQ=37)

                                                                                                                  cgctgtggTGCG‑TGGCACCGGTTTTGCTGGT‑GTTCT‑‑CGGA‑CA‑C‑CCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa            <  1:13180200/131‑1 (MQ=255)

                                                                                                                  cgctgtggTGCA‑TGGCACTGTTTTTGCTGGT‑GTTCG‑‑TGGG‑CA‑C‑TCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa            <  1:712035/131‑1 (MQ=255)

                                                                                                                  cgccgcGGGGCA‑GGGCACCGGCTTTGCCGGT‑GTT‑CG‑CGGG‑CA‑C‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa            <  2:10290123/137‑1 (MQ=37)

                                                                                                                  caccgcGGGGCA‑TGGCACCGGCTTTGCTGGT‑GTT‑CG‑CGGG‑CA‑CTC‑CC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa            <  2:5912540/135‑1 (MQ=255)

                                                                                                                    ccgcGGTGCA‑GGGCACCGGCTTGGCCGGT‑GTT‑CG‑CGGG‑CA‑C‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGa          <  2:9732185/137‑1 (MQ=37)

                                                                                                                       cgGTGCA‑TGGCACCGGCTTTTCTGGT‑GTT‑CG‑CGGG‑CA‑C‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACg       <  2:16595538/137‑1 (MQ=37)

                                                                                                                            cA‑TGGCATCGGTTTTGCTGGT‑GTTCG‑‑TGGG‑CA‑C‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATc  <  1:13986620/137‑1 (MQ=38)

||

GTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTG‑CCCGCGAAACAGGCGCCGCGGTGCA‑TGGCCACGGCTTTGCCGGT‑GTT‑T‑‑CGGG‑CA‑C‑G‑CC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATC  >  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: `‑`