BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	AE015451_phaM1	4,586,057	+C	intergenic (+140/+75)	PP_4061 CDS (conserved protein of unknown function) → / ← PP_4063 CDS (putative Long‑chain‑fatty‑acid‑CoA ligase)	–/–

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AE015451_phaM1	4,586,056	1	.	C	99.4%	618.3 / ‑2.9	160	intergenic (+139/+76)	PP_4061 CDS (conserved protein of unknown function)/PP_4063 CDS (putative Long‑chain‑fatty‑acid‑CoA ligase)	–/–
Reads supporting (aligned to +/- strand): ref base . (1/0); new base C (27/132); total (28/132)
Fisher's exact test for biased strand distribution p-value = 1.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 179 total aligned reads displayed.

TAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTCT‑‑GC‑GGTAAACCCGCTCCTACC‑AGGGA‑CGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGG  >  AE015451_phaM1/4585921‑4586183

tAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑AGGGa                                                                                                                                  <  1:15435616/137‑1 (MQ=39)

  gACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑AGGGA‑cc                                                                                                                               <  2:4263709/137‑2 (MQ=39)

  gACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑AGGGA‑cc                                                                                                                               <  2:18355378/137‑2 (MQ=39)

  gACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑AGGGA‑cc                                                                                                                               <  1:12393250/137‑2 (MQ=39)

  gACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑AGGGA‑cc                                                                                                                               <  2:2370907/137‑2 (MQ=39)

  gACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑AGGGA‑cc                                                                                                                               <  1:5930025/137‑2 (MQ=39)

  gACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑AGGGA‑cc                                                                                                                               <  1:150845/137‑2 (MQ=39)

   aCCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑AGGGA‑ccg                                                                                                                              <  1:367100/137‑3 (MQ=38)

   aCCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑AGGGA‑ccg                                                                                                                              <  2:15524731/137‑3 (MQ=38)

     cACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑AGGGA‑ccgtg                                                                                                                            <  1:2439389/137‑5 (MQ=38)

      aCGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTCTTCGCGGGTAAACCCGCCACTAAC‑AAGGA‑ccat                                                                                                                             >  2:12426252/1‑134 (MQ=255)

       cGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑AGGGACCGTGcc                                                                                                                           <  2:2178149/137‑1 (MQ=38)

          cACCTGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGTTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGc                                                                                                                        <  2:16319719/137‑1 (MQ=38)

          cACCGGCGTGAACAGCTGGTCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGc                                                                                                                        <  2:1987213/137‑1 (MQ=38)

          cACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCAGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCGCCTACC‑AGGGACCGTGCCGGc                                                                                                                        <  2:14624250/137‑1 (MQ=38)

           aCCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGcc                                                                                                                       >  2:12758416/1‑137 (MQ=38)

               gCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGa                                                                                                                   <  1:3676242/137‑1 (MQ=38)

               gCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGa                                                                                                                   <  1:7280551/137‑1 (MQ=38)

               gCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGa                                                                                                                   <  2:16029919/137‑1 (MQ=38)

               gCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGa                                                                                                                   <  1:1885163/137‑1 (MQ=38)

               gCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGa                                                                                                                   <  2:16029880/137‑1 (MQ=38)

               gCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGa                                                                                                                   <  1:11251049/137‑1 (MQ=38)

                cGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGaa                                                                                                                  <  2:4223198/137‑1 (MQ=38)

                   gAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATgc                                                                                                               <  1:4651443/137‑1 (MQ=38)

                    aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑ATGGACCGTGCCGGCCTTGAATgcg                                                                                                              <  1:14061102/137‑1 (MQ=38)

                    aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATgcg                                                                                                              <  2:2753710/137‑1 (MQ=38)

                    aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATgcg                                                                                                              <  1:13650204/137‑1 (MQ=38)

                    aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATgcg                                                                                                              <  2:16958767/137‑1 (MQ=38)

                    aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATgcg                                                                                                              <  1:5550407/137‑1 (MQ=38)

                    aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATgcg                                                                                                              <  1:5572717/137‑1 (MQ=38)

                    aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATgcg                                                                                                              <  1:13527451/137‑1 (MQ=38)

                    aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATgcg                                                                                                              <  2:246654/137‑1 (MQ=38)

                    aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATgcg                                                                                                              <  1:13527412/137‑1 (MQ=38)

                    aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATgcg                                                                                                              <  1:8440283/137‑1 (MQ=38)

                    aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATgcg                                                                                                              <  1:12758416/137‑1 (MQ=38)

                    aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATgcg                                                                                                              <  1:11691004/137‑1 (MQ=38)

                    aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATgcg                                                                                                              <  2:10697267/137‑1 (MQ=38)

                    aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATgcg                                                                                                              <  2:8221450/137‑1 (MQ=38)

                    aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATgcg                                                                                                              <  1:3480853/137‑1 (MQ=38)

                    aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATgcg                                                                                                              <  1:15912504/137‑1 (MQ=38)

                    aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATgcg                                                                                                              <  2:9089059/137‑1 (MQ=38)

                    aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATgcg                                                                                                              <  2:13852447/137‑1 (MQ=38)

                    aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATgcg                                                                                                              <  1:18052642/137‑1 (MQ=38)

                    aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATgcg                                                                                                              <  1:17939383/137‑1 (MQ=38)

                    aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATgcg                                                                                                              <  2:17480120/137‑1 (MQ=38)

                    aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGACGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATgtg                                                                                                              <  1:16201017/137‑3 (MQ=38)

                    aaCAGCTGATCTGCGCTGCTCCATCGAAAGTTTGAGGAGGACCGGGTGCTTGGCAAGTTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑ATGGACCGTACCGGCCTTGAATgcg                                                                                                              <  2:11782704/137‑1 (MQ=37)

                     aCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATgcgg                                                                                                             <  1:16001660/137‑1 (MQ=38)

                     aCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATgcgg                                                                                                             <  1:13361205/137‑1 (MQ=38)

                      cAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATgcggc                                                                                                            <  1:19345172/137‑1 (MQ=38)

                      cAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATgcggc                                                                                                            <  1:5768725/137‑1 (MQ=38)

                       aGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑ATGGACCGTGCCGGCCTTGAATgcggcg                                                                                                           <  2:18576558/137‑1 (MQ=38)

                       aGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATgcggcg                                                                                                           <  1:17060919/137‑1 (MQ=38)

                       aGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATgcggcg                                                                                                           <  2:11691423/137‑1 (MQ=38)

                       aGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATgcggcg                                                                                                           <  1:15376053/137‑1 (MQ=38)

                       aGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATgcggcg                                                                                                           <  2:6590620/137‑1 (MQ=38)

                         cTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATGCGGCGCt                                                                                                         <  1:18404560/137‑1 (MQ=38)

                         cTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATGCGGCGCt                                                                                                         <  2:14742715/137‑1 (MQ=38)

                         cTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATGCGGCGCt                                                                                                         <  2:1361752/137‑1 (MQ=38)

                         cTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATGCGGCGCt                                                                                                         <  1:18125040/137‑1 (MQ=38)

                         cTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATGCGGCGCt                                                                                                         <  2:5853843/137‑1 (MQ=38)

                         cTGATCTGCACTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATGCGGCGCt                                                                                                         <  2:3941837/137‑1 (MQ=38)

                          ttaTCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTGGAAGGCGgcgcgg                                                                                                        <  2:4597407/135‑3 (MQ=37)

                          tGATCTGCGCTGCTCCATTGAAAGGTTTAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑ATGGACCGTGCCGGCCTTGAATGCGGCGCTg                                                                                                        <  2:9671274/137‑1 (MQ=38)

                          tGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATGCGGCGCTg                                                                                                        <  1:16634923/137‑1 (MQ=38)

                               tGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑ATGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCt                                                                                                   <  2:18278688/137‑1 (MQ=38)

                               tGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCt                                                                                                   <  2:2664637/137‑1 (MQ=38)

                               tGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCt                                                                                                   <  2:9750241/137‑1 (MQ=38)

                               tGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCt                                                                                                   <  2:15645242/137‑1 (MQ=38)

                               tGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCt                                                                                                   <  1:3504712/137‑1 (MQ=38)

                               tGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCt                                                                                                   <  1:8677588/137‑1 (MQ=38)

                                   ctgctCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGTTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTa                                                                                               <  2:11441035/137‑1 (MQ=38)

                                   ctgctCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTACACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTa                                                                                               <  1:19128751/137‑1 (MQ=38)

                                   ctgctCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCTGCCTTGAATGCGGCGCTGCCCCTGGTa                                                                                               <  2:17184057/137‑1 (MQ=38)

                                   ctgctCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTa                                                                                               <  1:405755/137‑1 (MQ=38)

                                   ctgctCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTa                                                                                               <  2:6281913/137‑1 (MQ=38)

                                   ctgctCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTa                                                                                               <  1:17400254/137‑1 (MQ=38)

                                   ctgctCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTa                                                                                               <  2:2960959/137‑1 (MQ=38)

                                   ctgctCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTa                                                                                               <  1:9605382/137‑1 (MQ=38)

                                   ctgctCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTa                                                                                               <  2:1225217/137‑1 (MQ=38)

                                   ctgctCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTa                                                                                               <  1:10608784/137‑1 (MQ=38)

                                   ctgctCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTa                                                                                               <  2:16631609/137‑1 (MQ=38)

                                    tgctCCATTGAAAGGTTGGGGGTGGCAGGGTGCTTGGCGAGGGGAAGAGCTGGAATGGATCGGGTCGGCCTCTGC‑‑GCGGGTAAACCCGCTCCTACC‑GTGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGGAg                                                                                              <  1:2291774/137‑1 (MQ=255)

                                    tgctCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAg                                                                                              <  1:18965246/137‑1 (MQ=38)

                                         cATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATgcggcg                                                                                                           >  1:5014335/1‑119 (MQ=37)

                                         cATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATgcggcg                                                                                                           <  2:5014335/119‑1 (MQ=37)

                                         cATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAAACCGCTCCCCCCAAGGGA‑ACTGCCGGACCTGGACGCTGTGCTGTCCCCGGTAGGAgag                                                                                         >  2:5259717/1‑135 (MQ=255)

                                          agtaaaaGGTGGGGGAGGACCGGGTGCTGGACAAGGTGAAGATCTTGCATTGATTGGGCTGTCCTCTTC‑‑GCGGGGAAACCCGCCCCTACC‑AGGGACCGTGCCGTCCTGGAATGCGACGCTGCCCCGGGTAGGAGCgg                                                                                        <  2:14357948/133‑1 (MQ=255)

                                                   ttGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAg                                                                                              <  1:11569844/122‑1 (MQ=37)

                                                                ggTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTC‑‑GCGGGTACAACCACCCCTACA‑AGGGACCGTGCCGGCCTTGAATGCGGCGGTGTCCGTGGGCGGAGCGGGTTTCCCCGAGAAGGggccaga                                                                  >  2:7702213/1‑134 (MQ=255)

                                                                                                                                ccTACC‑AGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCAtac    >  2:13482998/1‑135 (MQ=255)

                                                                                                                                 cTACC‑AGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACAATGATCTCCCGCATGCg   >  2:4914235/1‑137 (MQ=255)

                                                                                                                                  tACC‑GGGGACCGTGCCGGCCTTGAATGTGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCACAGATGCAGCCGCAATCTCCGCCACACTAATCCCACGCATGCgg  >  2:5211795/1‑137 (MQ=38)

                                                                                                                                  tACC‑CGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCAGCAATCTCCGCCACACTGATCTCACGTATGCgg  >  2:3943020/1‑137 (MQ=39)

                                                                                                                                  tACC‑AGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCgg  >  1:16380305/1‑137 (MQ=255)

                                                                                                                                  tACC‑AGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCgg  >  1:6320407/1‑137 (MQ=255)

                                                                                                                                  tACC‑AGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCgg  >  1:1168912/1‑137 (MQ=255)

                                                                                                                                  tACC‑AGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCgg  >  1:8480878/1‑137 (MQ=255)

                                                                                                                                  tACC‑AGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCgg  >  2:14874341/1‑137 (MQ=255)

                                                                                                                                  tACC‑AGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCgg  >  1:10591154/1‑137 (MQ=255)

TAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTCT‑‑GC‑GGTAAACCCGCTCCTACC‑AGGGA‑CGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGG  >  AE015451_phaM1/4585921‑4586183

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