BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	AE015451_phaM1	5,555,608	+GCC	intergenic (‑84/+87)	PP_4887 CDS (conserved protein of unknown function) ← / ← PP_4888 CDS (putative methyl‑accepting chemotaxis transducer)	–/–

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AE015451_phaM1	5,555,608	1	.	G	92.3%	769.6 / 77.2	281	intergenic (‑84/+87)	PP_4887 CDS (conserved protein of unknown function)/PP_4888 CDS (putative methyl‑accepting chemotaxis transducer)	–/–
Reads supporting (aligned to +/- strand): ref base . (1/20); new base G (17/234); total (18/263)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.89e-01
*	AE015451_phaM1	5,555,608	2	.	C	92.5%	887.6 / 72.5	281	intergenic (‑84/+87)	PP_4887 CDS (conserved protein of unknown function)/PP_4888 CDS (putative methyl‑accepting chemotaxis transducer)	–/–
Reads supporting (aligned to +/- strand): ref base . (1/20); new base C (17/242); total (18/263)
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,555,608	3	.	C	81.8%	823.4 / 201.9	281	intergenic (‑84/+87)	PP_4887 CDS (conserved protein of unknown function)/PP_4888 CDS (putative methyl‑accepting chemotaxis transducer)	–/–
Reads supporting (aligned to +/- strand): ref base . (2/49); new base C (16/214); total (18/263)
Fisher's exact test for biased strand distribution p-value = 5.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 294 total aligned reads displayed.

TGCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCAT‑‑A‑‑‑GCG‑‑‑CC‑GGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTG  >  AE015451_phaM1/5555472‑5555730

|||

tgCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCAT‑‑A‑‑‑gcg                                                                                                                                <  1:1142727/137‑1 (MQ=255)

    ggAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCAT‑‑A‑‑‑gcggccc                                                                                                                            <  1:9218097/137‑5 (MQ=255)

                   ccTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCGT‑‑A‑‑‑GCGGCCCC‑GGCAATTTCCGCgt                                                                                                            <  2:16130281/137‑1 (MQ=38)

                   ccTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCAT‑‑A‑‑‑GCGGCCCC‑GGCAATTTCCGCgt                                                                                                            <  2:15513528/137‑1 (MQ=38)

                   ccTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCAT‑‑A‑‑‑GCGGCCCC‑GGCAATTTCCGCgt                                                                                                            <  2:11632324/137‑1 (MQ=38)

                   ccTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCAT‑‑A‑‑‑GCGGCCCC‑GGCAATTTCCGCgt                                                                                                            <  2:10482511/137‑1 (MQ=38)

                   ccTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCAT‑‑A‑‑‑GCGGCCCC‑GGCAATTTCCGCgt                                                                                                            <  2:3346844/137‑1 (MQ=38)

                   ccTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCAT‑‑A‑‑‑GCGGCCCC‑GGCAATTTCCGCgt                                                                                                            <  1:272055/137‑1 (MQ=38)

                     tGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCAT‑‑A‑‑‑GCGGCCCC‑GGCAATTTCCGCgtgt                                                                                                          <  2:17676335/137‑1 (MQ=38)

                        gCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCAT‑‑A‑‑‑GCGGCCCC‑GGCAATTTCCGCGTGTATg                                                                                                       <  1:9314504/137‑1 (MQ=38)

                        gCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCAT‑‑A‑‑‑GCGGCCCC‑GGCAATTTCCGCGTGTATg                                                                                                       <  1:714437/137‑1 (MQ=38)

                        gCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCAT‑‑A‑‑‑GCGGCCCC‑GGCAATTTCCGCGTGTATg                                                                                                       <  2:8205806/137‑1 (MQ=38)

                        gCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCAT‑‑A‑‑‑GCGGCCCC‑GGCAATTTCCGCGTGTATg                                                                                                       <  2:760059/137‑1 (MQ=38)

                         cAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCAT‑‑A‑‑‑GCGGCCCC‑GGCAATTTCCGCGTGTATGc                                                                                                      <  2:17437308/137‑1 (MQ=38)

                         cAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCAT‑‑A‑‑‑GCGGCCCC‑GGCAATTTCCGCGTGTATGc                                                                                                      <  1:7532568/137‑1 (MQ=38)

                         cAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCAT‑‑A‑‑‑GCGGCCCC‑GGCAATTTCCGCGTGTATGc                                                                                                      <  1:13604921/137‑1 (MQ=38)

                         cAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCAT‑‑A‑‑‑GCGGCCCC‑GGCAATTTCCGCGTGTATGc                                                                                                      <  1:17654368/137‑1 (MQ=38)

                         cAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCAT‑‑A‑‑‑GCGGCCCC‑GGCAATTTCCGCGTGTATGc                                                                                                      <  1:13914906/137‑1 (MQ=38)

                         cAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCAT‑‑A‑‑‑GCGGCCCC‑GGCAATTTCCGCGTGTATGc                                                                                                      <  1:12331298/137‑1 (MQ=38)

                         cAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCAT‑‑A‑‑‑GCGGCCCC‑GGCAATTTCCGCGTGTATGc                                                                                                      <  2:18823720/137‑1 (MQ=38)

                         cAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAACAACACCTCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCAT‑‑A‑‑‑GCGGCCCC‑GGCAATTTCCGCGTGTATGc                                                                                                      <  2:11080305/137‑1 (MQ=37)

                              aTTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCAT‑‑A‑‑‑GCGGCCCC‑GGCAATTTCCGCGTGTATGcaacaa                                                                                                 <  1:5576424/137‑1 (MQ=38)

                                 aGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCAT‑‑A‑‑‑GCGGCCCC‑GGCAATTTCCGCGTGTATGCAACAAAcc                                                                                              <  2:11344568/137‑1 (MQ=38)

                                 aGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCAT‑‑A‑‑‑GCGGCCCC‑GGCAATTTCCGCGTGTATGCAACAAAcc                                                                                              <  2:1521995/137‑1 (MQ=38)

                                 aGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCAT‑‑A‑‑‑GCGGCCCC‑GGCAATTTCCGCGTGTATGCAACAAAcc                                                                                              <  2:13567359/137‑1 (MQ=38)

                                 aGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCAT‑‑A‑‑‑GCGGCCCC‑GGCAATTTCCGCGTGTATGCAACAAAcc                                                                                              <  1:11037429/137‑1 (MQ=38)

                                   tCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCAT‑‑A‑‑‑GCGGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTg                                                                                            <  1:16659340/137‑1 (MQ=38)

                                   tCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCAT‑‑A‑‑‑GCGGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTg                                                                                            <  1:14681348/137‑1 (MQ=38)

                                   tCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCAT‑‑A‑‑‑GCGGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTg                                                                                            <  1:2241576/137‑1 (MQ=38)

                                   tCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCAT‑‑A‑‑‑GCGGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTg                                                                                            <  2:8445392/137‑1 (MQ=38)

                                   tCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCAT‑‑A‑‑‑GCGGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTg                                                                                            <  2:14082571/137‑1 (MQ=38)

                                      ggCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTTTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCAT‑‑A‑‑‑GCGGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTgggg                                                                                         <  2:84163/137‑1 (MQ=38)

                                      ggCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCAT‑‑A‑‑‑GCGGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTgggg                                                                                         <  2:14004011/137‑1 (MQ=38)

                                      ggCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCAT‑‑A‑‑‑GCGGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTgggg                                                                                         <  2:12349441/137‑1 (MQ=38)

                                      ggCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCAT‑‑A‑‑‑GCGGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTgggg                                                                                         <  1:9820685/137‑1 (MQ=38)

                                      ggCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCAT‑‑A‑‑‑GCGGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTgggg                                                                                         <  2:19113459/137‑1 (MQ=38)

                                      ggCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCAT‑‑A‑‑‑GCGGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTgggg                                                                                         <  1:8289240/137‑1 (MQ=38)

                                      ggCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCAT‑‑A‑‑‑GCGGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTgggg                                                                                         <  2:15809570/137‑1 (MQ=38)

                                      ggCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCAT‑‑A‑‑‑GCGGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTgggg                                                                                         <  1:6199274/137‑1 (MQ=38)

                                      ggCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCAT‑‑A‑‑‑GCGGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTgggg                                                                                         <  2:3735452/137‑1 (MQ=38)

                                      ggCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCAT‑‑A‑‑‑GCGGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTgggg                                                                                         <  1:3790172/137‑1 (MQ=38)

                                      ggCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCACTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCAT‑‑A‑‑‑GCGGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTgggg                                                                                         <  2:15588912/137‑1 (MQ=38)

                                      ggCGGACATGAGCAACGTCGTGTGGCATGTGTCGCTGATAGCATAGTAGCATCGCTCGGGGGTGGGAGCGGACTTGGGGCGCGAAAGGGCCGTAT‑‑A‑‑‑GCGGT‑CCTGGCAATTTCCGCGTGTATGCAACAAACCTtggg                                                                                         <  1:11774604/137‑1 (MQ=255)

                                           aCATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCAT‑‑A‑‑‑GCGGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTGGGgctgct                                                                                    <  2:5506874/137‑1 (MQ=38)

                                                     cTCCTTTTGTCATGTGTCGCCGACAGCATAACAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGAGAAAGGGCCGCAT‑‑A‑‑‑GCGGCCCC‑GGCAATTTCCGCGTGGAGGAAACAAACCTGGGGCTGCTATGCAGCCCt                                                                          <  2:16399190/137‑1 (MQ=37)

                                                             ttCATGTGTCGCCGACAGCATAGCAACCCCTCTCGGGTGGGGGAGTGGCTTTGTGTTGCGAAAGGGCCGCAG‑‑A‑‑‑GCGTCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTGTGCATCCCTTTCGCGac                                                                  <  2:3504900/136‑1 (MQ=255)

                                                             ttCAGGTGTAGCCGAAAACAGAGCAACCCCGCTCAGGTGTGGGAGAGGCCTTGTGTCGCGAAAGGGCCAAAG‑‑A‑‑‑GCGGCCCC‑GGAAATTTCTTCGGGTATGCAACAAACCTGGGGCTGCTGTGCAGCCCTTTCGCGac                                                                  <  2:9086603/136‑1 (MQ=16)

                                                             gTCATGTGTCGCCGACAGCAGAGAAACACCGCTCAGGGGTGGGAGCGGCTTGGTGTCGAGAAGGGGCCGCAT‑‑A‑‑‑GCGGCCCC‑GGCACTTTCCGCGTGTGTGCAACAAACCTGGGGCTGCGGTGCAGCCCTTTCGCGac                                                                  <  2:18158229/137‑1 (MQ=255)

                                                               cgtttgtCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCAT‑‑A‑‑‑GCGGCCCC‑GGCAATTTCCGCGTGTGTGCAACAAACCTGGGGCTGCTATGCATCCCTTTCGCGacac                                                                <  2:13453772/133‑1 (MQ=37)

                                                               cATGTTTCGCCGACAGCATAGCAACACCGCTCAGGCGTGGGAGCGGCCTTGTGTAGAAAAAGGGCCGCAT‑‑A‑‑‑GCGGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTGGGACGGTGAGGAAGCCCTTTCGAGacac                                                                <  2:7366806/137‑1 (MQ=255)

                                                               cATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCAT‑‑A‑‑‑GCGGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGAGacac                                                                <  2:12498342/137‑1 (MQ=38)

                                                               cATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGACCTTGTTTCGCGAAAGGTCCGCAT‑‑A‑‑‑GCGGC‑CCTGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGacac                                                                <  2:13701849/137‑1 (MQ=255)

                                                                 tttttCGCCGACAGCATAGCAACACCTCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCAT‑‑A‑‑‑GCGGCCCC‑GGCACTTTCCGCGTGTATGCAACAAACCTGGGCCTGCTATGCAGCCCTTTCGCGACACaa                                                              <  2:2526211/133‑1 (MQ=37)

                                                                 tgttccGCCGACAGCATAGCAACACCTCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGACCACAT‑‑A‑‑‑GCGGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACaa                                                              <  2:11800236/132‑1 (MQ=37)

                                                                 tgttTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCAG‑‑A‑‑‑GCGGCCCC‑GGCAATTTCCGCGTGTATACAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACaa                                                              <  2:14140136/137‑1 (MQ=37)

                                                                 tgtgTCGCTGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCAT‑‑A‑‑‑GCGGCCCC‑GGAAATTTCCGTGTGTAGGAAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACaa                                                              <  2:14360328/137‑1 (MQ=25)

                                                                 tgtgTCGCCGAAAGCATAGCAACACCGCTCAGGTGGGGGAGCGGCTTTGTTTCGCGAAAGGGCCACAT‑‑A‑‑‑GCGGCCCC‑GGAAGTTTCCGCGGGGATGCAACAAACTTGGGGCTGCTATGCAGCCCTTTCGCGACACaa                                                              <  2:9737461/137‑1 (MQ=25)

                                                                 tgtgTCGCCAAAAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCAT‑‑A‑‑‑GCGGCCCC‑GGCAATTTCCGCGTGTATACAACAAACCTGGGGCTGCTGTGCAGCCCTTTCGCGACACaa                                                              <  2:13221872/137‑1 (MQ=37)

                                                                 tgtgTCACCGACAGCATAGCAACATCGCTCGGGTGTGGGAGCGGCTTTGTGTCGCGAAAGGGCCAAGG‑‑A‑‑‑GCGGC‑CCTGGCAATTTCCCCGTGTATGCAACAAACTTGGGGCTACTAGGCAGCCCTTTCGCGACACaa                                                              <  2:4952693/137‑1 (MQ=255)

                                                                 ggtgtCGCCGACAACATAGCAACACCGCTCAGGTGTGGGAGTGGCTTTGTGTCGAGAAAGGGCCGCGT‑‑A‑‑‑GCGGCCCC‑GGAAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTGTGCATCCTTTTCTCGACACaa                                                              <  2:2259456/136‑1 (MQ=255)

                                                                      cGCTGACAGCATAGCAACACCGCTCAGGTTTGGGAGCGGCTTTGTGTCGCGAAAGGGCTGCGT‑‑G‑‑‑GCGGCCCC‑GTCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCCATACAGCCCTTTTGCGACACAAGGccc                                                         <  2:4603033/137‑1 (MQ=255)

                                                                                agagaaaCACCGTTAAGGTGGGGGAGAGGCCTTGTGTCGCGAAAGGGCAGCAG‑‑A‑‑‑GCGGC‑CCTGGAACTTTCCGCGTGGAGACAAAAAACCTGGGGTTGCTATGAAGCCCTTTCGCGAAACAAGGCCCACAAAGGAcc                                               <  2:3149130/132‑1 (MQ=255)

                                                                                aTAGCAACACCGCTCGGGTGTGGGAGGGGCTTTGTGTCGCGAAAGGGCCGCAT‑‑A‑‑‑GCGGC‑CCGGGCAGTTTCCGCGTGTGTGAAAAAAACTTGGGGCTGCTATGCAGCCTTTTCGCGACACAAGGCCCACAAAGGAcc                                               <  2:598780/137‑1 (MQ=25)

                                                                                aTAGCAACACCGCTCAGGTGTGGGAGCGACCTTGTGTCGCGAAAGGGCCGCAT‑‑A‑‑‑GCGGC‑CCTGGCACTTTCCGCTTGGCTAAAACAAACCTGGGGCTGCTGTGAATCCCTTTCGCGACACAAGGCCCACAAAGGAcc                                               <  2:2069796/137‑1 (MQ=25)

                                                                                aTAACAACACCCCTAAGGTGTGGGGGCGGCCTTGTGTCGCGAAAGGGCCGAAT‑‑AGTGGCC‑‑‑CC‑GGCAATTTCCGCGTGTAGACAACAAACTTGGGTCTGCTGTGCAGCCTTTTCGAGACACAAGGCCCACAAAGGAcc                                               <  2:12783138/137‑1 (MQ=19)

                                                                                 taacAACACCGCTCGGGTGGGGGAGCGGCTTTGTGTCGTGAAAGGGCCGCAT‑‑A‑‑‑GCGGC‑CCTGGCAATTTCCGCGTGTATACAACAAACCTGGGGCGGCGATGCAGCCTTTTCGCGACACAAGGCCCACAAAGGACcg                                              <  2:9469447/134‑1 (MQ=255)

                                                                                  aGCAACACCTCTAGGGTTGGGGAGCGGCCTTGTGTCGAGAAAGGGCCGCAG‑‑A‑‑‑GCGGC‑CCAGGCAATTTCTGCGTGTATGAAAAAAACTTGGGGCTGCTGTGCAGCCCTTTCGAGACACAAGGCCCACAAAGGACcgc                                             <  1:4812863/137‑1 (MQ=19)

                                                                                  aGCAACACCGCTCAGGTGTGGGAGCGGCCTTTTTCCGCGAAAGGGCCGCAT‑‑AGTGGCC‑‑‑CC‑GGCAATTTCCGCGTTTATGCAACAAACCTGGGGCTGCTATGCATCCCTTTCGCGACACAAGGCCCACAAAGGACcgc                                             <  2:17885359/137‑1 (MQ=37)

                                                                                  aGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCAT‑‑A‑‑‑GCGGCCCC‑GGCAATTTCCGCGTGTATGcaacaa                                                                                                 <  2:5835287/85‑1 (MQ=21)

                                                                                       cacCGCTCGGGTGTGGGAGCGGCTTTGTTTCGCGAAAGGGCTGCAT‑‑A‑‑‑GCGACCCC‑GGCAATTTCCGCGTGTATAAAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAgg                                        <  1:6756754/137‑1 (MQ=37)

                                                                                       cacCGCTCAGGTGTGGGAGCGGCCTTGTTTCGCGAAGGGGCCGCAT‑‑AGTGGCC‑‑‑CC‑GGCAGTTTCCGCGTGTATGCAACAAACCTGGGGTTGCGGTGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAgg                                        <  2:2714501/137‑1 (MQ=37)

                                                                                          cGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCAT‑‑A‑‑‑GCGGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTg                                                                                            >  2:8703284/1‑82 (MQ=18)

                                                                                          cGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCAT‑‑A‑‑‑GCGGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTg                                                                                            <  1:8703284/82‑1 (MQ=18)

                                                                                             taaGGTGGGGGAGCGGCTTTGTGTCGCGAAAGGGCCGAAG‑‑A‑‑‑GCGGCCCC‑GGCAGTTTCCGCGGGTATGCAACAAACTTGGGGCTGCTGTGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAg                                  <  2:9590145/135‑1 (MQ=25)

                                                                                             tCAGTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCAT‑‑A‑‑‑GCGGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTGTGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAg                                  <  2:10183645/137‑1 (MQ=37)

                                                                                             tCAGGTGTGGGAGCGGCTTTGTGTCGCGAAAGGGCCGTAT‑‑A‑‑‑GCGGCCCC‑GGAAATTTCCGCGTGTATGCAACAAACCTGGGGCGGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAg                                  <  2:17283166/137‑1 (MQ=37)

                                                                                             tCAGGTGGGGGAGCGGCCTTGTGTCGAGAAAGGGCCGCAT‑‑A‑‑‑GCGGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACTTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAg                                  <  2:13105603/137‑1 (MQ=37)

                                                                                             ccAGGTGTGGGAGCGGCCTGGGTTCGCGAAAGGGCCGCAT‑‑A‑‑‑GCGGCCCC‑GACAATTTCCGCGGGTAGGAAACAAACCTGGGGCTGCGGTGCAGCCCTTTCGCGAAACAAGGCCCACAAAGGACCGCGCAGGCTTCAg                                  <  2:5637401/136‑1 (MQ=255)

                                                                                                  tgtgGGGGTGGCTTTGTGTCGCGAAAGGGCCGCGT‑‑AGTGGCC‑‑‑CC‑GGCAGTTTCCGTGTGTGTGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCg                             <  2:19418395/137‑1 (MQ=25)

                                                                                                  tgtgGGATCGGCTTTGTTTCGCGAAAGGGCCGAAA‑‑A‑‑‑GAGGC‑CCAGGCAGTTCCCGTGTGCGTACAACAAACCTGGGGGTGCGAGGCAGCCCTTTCGCGANACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCg                             <  1:18816175/137‑1 (MQ=255)

                                                                                                  tgtgGGAGCGGCTTTGTGTCGAGAAGGGGCCACAG‑‑A‑‑‑GCGTCCCC‑GGCATTTTCTGTGTGCAGGAAAAAAACTTGGGGCTGCGAGGAATCCCTTTCGAGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCg                             <  1:10133660/137‑1 (MQ=255)

                                                                                                  tgtgGGAGCGGCCTTGTGTTGAGAAAGTGCCGCAT‑‑A‑‑‑GCGGC‑CCTGGAAGTTTCCGCGTGTATGCAACAAACTTNGGGCGGCCATGCAACCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCg                             <  1:17962740/137‑1 (MQ=255)

                                                                                                  tgtgGGAGCGGCCTTGTGTCGCGAAAGGGCCGCAT‑‑A‑‑‑GCGGCCCC‑GGCACTTTCCGCGTGTATGCAACAACCCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCg                             <  2:6351118/137‑1 (MQ=37)

                                                                                                  tgtgGGAGCGGCCTTGTGTCGCGAAAGGGCCGCAG‑‑A‑‑‑GCGGCCCC‑GGAAATTTCTGTGTGTGTGCAACAAACTTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCg                             <  2:12097926/137‑1 (MQ=25)

                                                                                                  tgtgGGAGCGGCCTTGTGTCACGAAAGGTCCGCAC‑‑A‑‑‑GCGGCCCC‑GAAAATTTCCGCGTGTATGCAACAAACCTGGGGCTACTATGCATCCTTTTCGCAACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCg                             <  1:10151500/137‑1 (MQ=25)

                                                                                                  tgtgGGAGCGGCCTTGTGTAGCGAAAGGTCCGCAT‑‑A‑‑‑GCGGCCCC‑GGCAGTTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCg                             <  2:18401989/137‑1 (MQ=37)

                                                                                                  tgggggAGTGGCCTTGTTTCGAGAAAGGGCCACAT‑‑AGTGGCC‑‑‑CC‑GGCAATTTCCGTGTGTATGCAAAAAACTTGGGGTGGCTATAAATCCTTTTCGCAACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCg                             <  1:13027985/134‑1 (MQ=255)

                                                                                                   gtgGGAGTGGCTTTGTGTCGCGAAAGGTCCATAT‑‑AGTGGCC‑‑‑CC‑GGAAATTTCCGCGTGTAGGCAAAAAACTTGGGGCTGCTAGGCAGCCTTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCgg                            <  1:15285035/137‑1 (MQ=255)

                                                                                                   gtgGGAGCGGCTTTGTGTCGAGAAAGTGCCGCAT‑‑A‑‑‑GCGGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCgg                            <  1:4062155/137‑1 (MQ=38)

                                                                                                   gggggAGCGGCTTTGTGTCGCGAAAGGGCCGCAT‑‑AGTGGCC‑‑‑CC‑GGCAATTTCCGTGGGTAGGAAAAAAACTGGGGGTTGTTGTGCATCCTTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCgg                            <  1:17195252/135‑1 (MQ=16)

                                                                                                   gggggAGAGGCCTTGTGTCGAGAAAGGGCCGAATGGAG‑‑GCC‑‑‑CC‑GGAAGTTTCCGCGTGGAGGAAACAAACTTGGGGCGGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCgg                            <  1:6101186/135‑1 (MQ=255)

                                                                                                                  tgtCGCGAAAGGGCCGCAT‑‑A‑‑‑GCGGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACTTGGGGCTGCGATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTg             <  2:5391115/137‑1 (MQ=255)

                                                                                                                     cgcgAAGGGGCCGCGT‑‑A‑‑‑GCGCCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTTGGGCTGCGAGGCAGCTCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGgcg          <  1:14841956/137‑1 (MQ=255)

                                                                                                                         aaaGGGCCGCAG‑‑A‑‑‑GCGGCCCC‑GGCAATTTCCGCGTGTATGAAACAAACCTGGGGCTGCTATGCATCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCa      <  2:10714789/137‑1 (MQ=255)

                                                                                                                         aaaGGGCCACAT‑‑A‑‑‑GCGGCCCC‑GGCAATTTCCGCGTGCAGAAAACAAACCGGGGGCTGCTAGGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCa      <  2:10923445/137‑1 (MQ=255)

                                                                                                                          aaaggCCGCAT‑‑A‑‑‑GCGGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTCTGCAGCCCTTTTGCGGCACAAGGCCCACAAAGGAACGCGCAGGCTTCAGGCCCGGAACTGGCCCAAGCTGGCGCGTAg     >  2:7133521/4‑137 (MQ=255)

                                                                                                                          aaGGGCCGCAT‑‑A‑‑‑GCGGCCCC‑GGCCATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGTAGCCCCTTCGCGAAACAAGGACCACAAAGGACCGCGCCGGCTTCCGGCCCGGAACCTGCCCAGGCTGGCGCGCAg     >  2:4234362/1‑137 (MQ=255)

                                                                                                                          aaGGGCCGCAT‑‑A‑‑‑ACGGCCCC‑GGCAATTTCTGCGTGTATGCAACAAACCTGGGGTTGCTATGTAGTCCTTTCGCGAAACAAGGCCAACAAAGGACCGCGCCGGCTTCAGGCCCGGAACTGGCCCCGGGCGGCACGCAg     >  2:15683997/1‑137 (MQ=255)

                                                                                                                          aaGGGACGCAT‑‑A‑‑‑GCGGCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCCGGGGCTGCTATGCAGCCCCTTCGCGAAACAAGGCCCACAAAGGACCGCGTAGGATTCAAGCCCGGAACTGGCCCCGGCTGGTGCGCaa     >  2:7428915/1‑136 (MQ=255)

                                                                                                                             ggCCGCAT‑‑A‑‑‑GCGTCCCC‑GGCAATTTCCGCGTGTATGCAACAAACCTGGGTTTGCTGTGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTg  <  2:16602909/137‑1 (MQ=255)

|||

TGCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCAT‑‑A‑‑‑GCG‑‑‑CC‑GGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTG  >  AE015451_phaM1/5555472‑5555730

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