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.5%	507.4 / ‑2.7	139	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); major base C (95/40); minor base G (1/1); total (98/41)
Fisher's exact test for biased strand distribution p-value = 5.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*	AE015451_phaM1	5,681,415	2	.	G	99.3%	516.9 / ‑3.2	139	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 (97/41); total (98/41)
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.98e-01
*	AE015451_phaM1	5,681,415	3	.	G	99.3%	530.1 / ‑3.2	139	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 (97/41); total (98/41)
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	4	.	G	97.8%	520.0 / 4.1	139	intergenic (‑133/+71)	PP_4986 CDS (Channel protein_ hemolysin III family)/PP_4987 CDS (putative Chemotaxis protein)	–/–
Reads supporting (aligned to +/- strand): ref base . (2/1); new base G (96/40); total (98/41)
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 161 total aligned reads displayed.

ACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑‑ACCGC‑‑‑‑GCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCTTTCCAA  >  AE015451_phaM1/5681279‑5681533

                                                                                                                                             ||||

aCATGCAAACACCTACCTCCTGAGGTGGGCTGGGTTCCTGTGTTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑‑ACcgc                                                                                                                            <  2:13959682/137‑1 (MQ=255)

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

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

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

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

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

   tGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑‑ACcgccgg                                                                                                                         <  1:9726367/137‑4 (MQ=255)

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

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

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

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

                   ctGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑‑ACCGCCGGGGCTGCTGCGCAGccc                                                                                                         <  1:10334751/137‑1 (MQ=37)

                     tagttGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCGTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTa                                                                                                       <  1:8446860/133‑1 (MQ=37)

                     taGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTa                                                                                                       <  2:17447513/136‑1 (MQ=37)

                     taGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTa                                                                                                       <  2:13889847/136‑1 (MQ=37)

                     gAGGTGGGCTTGGTTCCTGTTGTGATACTGATGGCTTCGGTCGTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCATGCTTTTATTGTCATTGTATGA‑‑‑‑ACCGCCGGGGCTGCTGCGCAGCTCTa                                                                                                       <  1:4199079/137‑1 (MQ=21)

                     gAGGTGGGCTGTGTTCCTGGTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTa                                                                                                       <  1:392980/137‑1 (MQ=37)

                     gAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGTCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTa                                                                                                       <  2:6335197/137‑1 (MQ=37)

                     gAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGNAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTa                                                                                                       <  1:16267562/137‑1 (MQ=37)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                     gAGGTGGGCTGGGTTCCGGTTGTGATGCTGATGGCTTCGTTCTGGGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTa                                                                                                       <  1:11768052/137‑1 (MQ=37)

                        aagggCTGGGTTCCTGTTGTGGTGCGGATGGCTTTGTTTTTTGTGGGAGTGGCCTTGTGTCGCGAAGGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCAGGA‑‑‑‑ACCGCCGGGGTTGCTGCGCAGCCCTATcg                                                                                                    <  1:16499798/135‑1 (MQ=25)

                              tGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑‑ACcgccggg                                                                                                                        >  1:9250962/1‑107 (MQ=255)

                              tGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑‑ACcgccggg                                                                                                                        <  2:9250962/111‑5 (MQ=255)

                                        ttgtGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGACCTTTGTCGCGAAAGGGGTGCGAAGCGGCCCCAGGCTTTCCTTGCCAATGCCCGAACCCACCGG‑‑‑‑GCTGCTGCGCAGACCTATCACGACACAAAGCCACTc                                                                                    >  2:8964528/1‑137 (MQ=255)

                                                                                         aaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTgcgc                                                          <  1:13151186/114‑1 (MQ=37)

                                                                                         aaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGTGCCTCACATCCTAACCAAGGTCAgg                                   >  1:16709624/1‑137 (MQ=38)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                                                                         aaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCAACAGGCCCTGCGCCTCACATCCTAACCAAAGTCAgg                                   >  1:4807289/1‑137 (MQ=38)

                                                                                         aaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTATCGCGAAACAAGGCCGCTCCCACTCCAGACCAACAGGCCCTGCGCCTCACATCCTAAACAAGGTCAgg                                   >  2:18254930/1‑137 (MQ=38)

                                                                                         aaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑‑ACCGCCGGGGCTGCTGCGCAGCACTATCGCGACACGAGGCAGCTCCCACCCCAGACCACCAGACCCTGCGCATCACATCCTAACCACGGTCAgg                                   >  2:16692421/1‑137 (MQ=255)

                                                                                         aaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑‑ACCGCCGGGGCTGCTGCGCAGCACTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACCTCCTAACCAAGGTCAgg                                   >  2:13234188/1‑137 (MQ=38)

                                                                                         aaGGGGCGCGAAGCGGCCCAAGGCTTTCATTGCCACTGCACGA‑‑‑‑ACCGCCGGGGCTGCTGCGCAGCACTATCGCGAAAAAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACAAAGGTAAgg                                   >  2:12517405/1‑137 (MQ=255)

                                                                                           ggggCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGtt                                 >  1:7948607/1‑137 (MQ=38)

                                                                                                gcgAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAgg                                        >  2:18515839/1‑125 (MQ=38)

                                                                                                gcgAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAgg                                        <  1:18515839/125‑1 (MQ=38)

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

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

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

                                                                                                   aaGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATcc                         >  2:14724273/1‑137 (MQ=38)

                                                                                                      cGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCcgcg                      >  2:3571037/1‑137 (MQ=38)

                                                                                                      cGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCcgcg                      >  2:3571069/1‑137 (MQ=38)

                                                                                                         ccccAGGCTTTCATTGCCACTGCATGA‑‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCa                   >  1:18595946/1‑137 (MQ=38)

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

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

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

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

                                                                                                            cAGGCTTTCATTGCCACTGCATGA‑‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCg                                                                                       <  1:15919308/66‑1 (MQ=17)

                                                                                                                 tttCATTGCCACTGCATGA‑‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAg           >  1:5916394/1‑137 (MQ=38)

                                                                                                                 tttCATTGCCACTGCATGA‑‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAg           >  2:2270321/1‑137 (MQ=38)

                                                                                                                 tttCATTGCCACTGCATGA‑‑‑‑ACCGCCGGGGCTGCTGCGCAGACCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAg           >  1:16019773/1‑137 (MQ=38)

                                                                                                                      ttGCCACTGCATGA‑‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTATCGCGGCACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCttt      >  1:12692635/1‑137 (MQ=38)

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

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

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

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

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

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

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

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

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

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

                                                                                                                      ttGCCACTGCATGA‑‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCAACAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCttt      >  1:13772854/1‑137 (MQ=38)

                                                                                                                      ttGCCACTGCATGA‑‑‑‑ACCGCAGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGCTCCGCGTCCCGCACCAGCCttt      >  1:3553542/1‑137 (MQ=38)

                                                                                                                       tGCCACTGCATGA‑‑‑‑ACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCTTTc     >  1:14912999/1‑137 (MQ=38)

                                                                                                                                             ||||

ACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGA‑‑‑‑ACCGC‑‑‑‑GCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCTTTCCAA  >  AE015451_phaM1/5681279‑5681533

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