BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	AE015451_phaM1	5,681,415	+CGGG	intergenic (‑133/+71)	PP_4986 CDS (Channel protein_ hemolysin III family) ← / ← PP_4987 CDS (putative Chemotaxis protein)	–/–

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AE015451_phaM1	5,681,415	1	.	C	98.9%	353.1 / ‑4.3	93	intergenic (‑133/+71)	PP_4986 CDS (Channel protein_ hemolysin III family)/PP_4987 CDS (putative Chemotaxis protein)	–/–
Reads supporting (aligned to +/- strand): ref base . (0/0); major base C (74/18); minor base G (1/0); total (75/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01
*	AE015451_phaM1	5,681,415	2	.	G	98.9%	353.7 / ‑2.9	93	intergenic (‑133/+71)	PP_4986 CDS (Channel protein_ hemolysin III family)/PP_4987 CDS (putative Chemotaxis protein)	–/–
Reads supporting (aligned to +/- strand): ref base . (1/0); new base G (74/18); total (75/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*	AE015451_phaM1	5,681,415	3	.	G	98.9%	348.2 / ‑2.9	93	intergenic (‑133/+71)	PP_4986 CDS (Channel protein_ hemolysin III family)/PP_4987 CDS (putative Chemotaxis protein)	–/–
Reads supporting (aligned to +/- strand): ref base . (1/0); new base G (74/18); total (75/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01
*	AE015451_phaM1	5,681,415	4	.	G	98.9%	354.1 / ‑2.9	93	intergenic (‑133/+71)	PP_4986 CDS (Channel protein_ hemolysin III family)/PP_4987 CDS (putative Chemotaxis protein)	–/–
Reads supporting (aligned to +/- strand): ref base . (1/0); new base G (74/18); total (75/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 111 total aligned reads displayed.

ACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑ACCGC‑‑‑‑GCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCTTT  >  AE015451_phaM1/5681279‑5681529

                                                                                                                                            ||||

aCATGCAAACACATACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGTGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑ACcgc                                                                                                                        <  2:1812551/137‑1 (MQ=255)

aCATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑ACcgc                                                                                                                        <  1:13490590/137‑1 (MQ=255)

aCATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑ACcgc                                                                                                                        <  1:245636/137‑1 (MQ=255)

aCATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑ACcgc                                                                                                                        <  2:6212126/137‑1 (MQ=255)

aCATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑ACcgc                                                                                                                        <  1:9444823/137‑1 (MQ=255)

  aTGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑ACcgccg                                                                                                                      <  1:9861148/137‑3 (MQ=255)

    gCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑ACcgccggg                                                                                                                    <  2:521324/137‑5 (MQ=255)

     cAAACACATACCTCCTGAGTTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑ACcgccgggg                                                                                                                   <  2:6315287/137‑6 (MQ=255)

     cAAACACATACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑ACcgccgggg                                                                                                                   <  2:10289412/137‑6 (MQ=255)

     cAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTTTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCGTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑ACcgccgggg                                                                                                                   <  2:6623300/137‑6 (MQ=255)

     cAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑ACcgccgggg                                                                                                                   <  2:6705204/137‑6 (MQ=255)

     cAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑ACcgccgggg                                                                                                                   <  1:12936926/137‑6 (MQ=255)

     cAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑ACcgccgggg                                                                                                                   <  1:12936927/137‑6 (MQ=255)

     cAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑ACcgccgggg                                                                                                                   <  2:14855230/137‑6 (MQ=255)

     cAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑ACcgccgggg                                                                                                                   <  2:6315332/137‑6 (MQ=255)

     cAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑ACcgccgggg                                                                                                                   <  1:15208321/137‑6 (MQ=255)

       aaCACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑ACcgccgggg                                                                                                                   >  1:8087422/1‑130 (MQ=255)

       aaCACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑ACcgccgggg                                                                                                                   <  2:8087422/135‑6 (MQ=255)

                    tGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCt                                                                                                    <  1:3015250/137‑1 (MQ=37)

                     gAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTa                                                                                                   <  1:6128310/137‑1 (MQ=37)

                     gAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTa                                                                                                   <  2:10651571/137‑1 (MQ=37)

                     gAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTa                                                                                                   <  1:3019790/137‑1 (MQ=37)

                     gAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTa                                                                                                   <  2:10217448/137‑1 (MQ=37)

                     gAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTa                                                                                                   <  1:8302788/137‑1 (MQ=37)

                     gAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTa                                                                                                   <  2:11337162/137‑1 (MQ=37)

                     gAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTa                                                                                                   <  2:4689462/137‑1 (MQ=37)

                     gAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTa                                                                                                   <  1:16198568/137‑1 (MQ=37)

                     gAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTa                                                                                                   <  1:13843194/137‑1 (MQ=37)

                     gAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTa                                                                                                   <  2:6600388/137‑1 (MQ=37)

                     gAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTa                                                                                                   <  1:13026406/137‑1 (MQ=37)

                     gAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTa                                                                                                   <  2:7869929/137‑1 (MQ=37)

                     gAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTa                                                                                                   <  2:7869930/137‑1 (MQ=37)

                     gAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTa                                                                                                   <  1:1017556/137‑1 (MQ=37)

                     gAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGCGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTa                                                                                                   <  2:984910/137‑1 (MQ=37)

                     gAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGCGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTa                                                                                                   <  2:984911/137‑1 (MQ=37)

                                ggTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTTGCGAAAAAGGAACGAAGCGGCCCCCAACCTTCCCTTCCACCGCACGAACCACCGCG‑‑‑GCTGCTGCGCAACCCCACCGCAACACaa                                                                                        >  2:7343777/1‑137 (MQ=255)

                                                                                       gAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCa                                 >  2:8659386/1‑137 (MQ=38)

                                                                                        aaaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACCCAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAg                                >  1:3982213/1‑137 (MQ=38)

                                                                                         aaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑ACCGCCGGGGCTGTTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAgg                               >  1:5214632/1‑137 (MQ=38)

                                                                                         aaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑ACCGCCGGGGCTGCTGCGCAGTCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAgg                               >  2:7243030/1‑137 (MQ=38)

                                                                                         aaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACGCAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAgg                               >  2:2026372/1‑137 (MQ=38)

                                                                                         aaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCCgg                               >  2:12856123/1‑137 (MQ=38)

                                                                                         aaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAgg                               >  2:5323300/1‑137 (MQ=38)

                                                                                         aaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAgg                               >  1:1605792/1‑137 (MQ=38)

                                                                                         aaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAgg                               >  2:1176498/1‑137 (MQ=38)

                                                                                         aaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAgg                               >  1:12737039/1‑137 (MQ=38)

                                                                                         aaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAgg                               >  2:9583946/1‑137 (MQ=38)

                                                                                         aaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAgg                               >  2:3200513/1‑137 (MQ=38)

                                                                                         aaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAgg                               >  2:13216774/1‑137 (MQ=38)

                                                                                         aaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAgg                               >  1:17312991/1‑137 (MQ=38)

                                                                                         aaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAgg                               >  2:16173140/1‑137 (MQ=38)

                                                                                         aaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAgg                               >  1:11450999/1‑137 (MQ=38)

                                                                                         aaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAgg                               >  2:17814916/1‑137 (MQ=38)

                                                                                         aaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAgg                               >  2:17814954/1‑137 (MQ=38)

                                                                                         aaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAgg                               >  2:17918025/1‑137 (MQ=38)

                                                                                         aaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAgg                               >  2:2686540/1‑137 (MQ=38)

                                                                                         aaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAgg                               >  1:8864026/1‑137 (MQ=38)

                                                                                         aaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAgg                               >  1:4546973/1‑137 (MQ=38)

                                                                                         aaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAgg                               >  1:6407537/1‑137 (MQ=38)

                                                                                         aaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAgg                               >  1:6407574/1‑137 (MQ=38)

                                                                                         aaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAgg                               >  1:6693475/1‑137 (MQ=38)

                                                                                         aaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAgg                               >  1:10677463/1‑137 (MQ=38)

                                                                                         aaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAgg                               >  1:3611381/1‑137 (MQ=38)

                                                                                         aaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAgg                               >  1:18476760/1‑137 (MQ=38)

                                                                                         aaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAgg                               >  1:9141785/1‑137 (MQ=38)

                                                                                         aaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAgg                               >  1:14628725/1‑137 (MQ=38)

                                                                                         aaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAgg                               >  1:3258255/1‑137 (MQ=38)

                                                                                         aaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAgg                               >  1:2883638/1‑137 (MQ=38)

                                                                                         aaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAgg                               >  1:15915613/1‑137 (MQ=38)

                                                                                         aaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGACCCTGCGCCTCACATCCTAACCAAGGTCAgg                               >  1:7674085/1‑137 (MQ=38)

                                                                                         aaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTATCACGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAgg                               >  2:4795272/1‑137 (MQ=38)

                                                                                         aaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑ACCGCCGGGGCTGCTGCGCAGCACTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAgg                               >  1:7674060/1‑137 (MQ=38)

                                                                                           ggggCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCCTGA‑‑‑ACCGCCGGGGCCACTGCTCAGCCCCATCGCGACCCAATGCCGCTCCCACTCCAGCCCACCCGGTCCCTCACCTCACATTCTAAACCAGGTCCGGtt                             >  1:330316/1‑137 (MQ=255)

                                                                                              gcgcGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCa                          >  1:6064533/1‑137 (MQ=38)

                                                                                                  gAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATc                      >  2:11883946/1‑137 (MQ=38)

                                                                                                  gAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATc                      >  2:10623531/1‑137 (MQ=38)

                                                                                                  gAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATc                      >  1:325854/1‑137 (MQ=38)

                                                                                                    aGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCCCTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCcg                    >  1:16743814/1‑137 (MQ=38)

                                                                                                          cccAGGCTTTCATTGCCACTGCATGA‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCAc              >  2:10301610/1‑137 (MQ=38)

                                                                                                           ccAGGCTTTCATTGCCACTGCATGA‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACg             >  2:11586067/1‑137 (MQ=38)

                                                                                                           ccAGGCTTTCATTGCCACTGCATGA‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACg             >  2:11586029/1‑137 (MQ=38)

                                                                                                           ccAGGCTTTCATTGCCACTGCATGA‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACg             >  1:1746331/1‑137 (MQ=38)

                                                                                                                      ttGCCACTGCATGA‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGTTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCttt  >  1:3453071/1‑137 (MQ=38)

                                                                                                                      ttGCCACTGCATGA‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCttt  >  1:12468483/1‑137 (MQ=38)

                                                                                                                      ttGCCACTGCATGA‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCttt  >  2:888961/1‑137 (MQ=38)

                                                                                                                      ttGCCACTGCATGA‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCttt  >  1:10432337/1‑137 (MQ=38)

                                                                                                                      ttGCCACTGCATGA‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCttt  >  2:782292/1‑137 (MQ=38)

                                                                                                                      ttGCCACTGCATGA‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCttt  >  1:12802936/1‑137 (MQ=38)

                                                                                                                      ttGCCACTGCATGA‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCttt  >  1:13677459/1‑137 (MQ=38)

                                                                                                                      ttGCCACTGCATGA‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCttt  >  1:13834968/1‑137 (MQ=38)

                                                                                                                      ttGCCACTGCATGA‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCttt  >  1:16624889/1‑137 (MQ=38)

                                                                                                                      ttGCCACTGCATGA‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCttt  >  1:16624890/1‑137 (MQ=38)

                                                                                                                      ttGCCACTGCATGA‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCttt  >  2:3396854/1‑137 (MQ=38)

                                                                                                                      ttGCCACTGCATGA‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCttt  >  2:2307383/1‑137 (MQ=38)

                                                                                                                      ttGCCACTGCATGA‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCttt  >  2:17050754/1‑137 (MQ=38)

                                                                                                                      ttGCCACTGCATGA‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCttt  >  2:11957412/1‑137 (MQ=38)

                                                                                                                      ttGCCACTGCATGA‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCttt  >  1:4746004/1‑137 (MQ=38)

                                                                                                                      ttGCCACTGCATGA‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTATCGCGAAACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCttt  >  2:11935775/1‑137 (MQ=38)

                                                                                                                      ttGCCACTGCATGA‑‑‑ACCGCCGGGGCTGCTGCACAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTAACATCCTACCCAAGGTCAGGTTGCAGATCCGCGTCAAGCACCAGCCttt  >  1:8078429/1‑137 (MQ=37)

                                                                                                                      ttGCCAATGCATGA‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTATCACGACACAAGGCCGCTCACACTCCCGACCACCAGGCCCTGCGCCTCACATCCTAACAAAGGTCAGGTTACAGATCCGCCTCACGCACCAGCCttt  >  2:13191059/1‑137 (MQ=25)

                                                                                                                                            ||||

ACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑ACCGC‑‑‑‑GCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCTTT  >  AE015451_phaM1/5681279‑5681529

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