BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	AE015451_phaM1	1,070,248	+GA	intergenic (+48/+94)	aroP‑I CDS (aromatic amino acid transport protein) → / ← PP_0928 CDS (K+‑dependent Na+_Ca+ exchanger related‑protein)	–/–

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AE015451_phaM1	1,070,246	1	.	G	97.4%	1857.2 / 37.5	492	intergenic (+46/+96)	aroP‑I CDS (aromatic amino acid transport protein)/PP_0928 CDS (K+‑dependent Na+_Ca+ exchanger related‑protein)	–/–
Reads supporting (aligned to +/- strand): ref base . (4/9); new base G (238/239); total (242/250)
Fisher's exact test for biased strand distribution p-value = 2.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*	AE015451_phaM1	1,070,246	2	.	A	97.3%	1782.3 / 34.1	492	intergenic (+46/+96)	aroP‑I CDS (aromatic amino acid transport protein)/PP_0928 CDS (K+‑dependent Na+_Ca+ exchanger related‑protein)	–/–
Reads supporting (aligned to +/- strand): ref base . (4/9); new base A (236/240); total (242/250)
Fisher's exact test for biased strand distribution p-value = 2.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 504 total aligned reads displayed.

GATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GT‑‑G‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTG  >  AE015451_phaM1/1070110‑1070372

||

gATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCt‑‑gt                                                                                                                                    >  2:11171580/1‑137 (MQ=255)

 aTTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCt‑‑gt‑‑g                                                                                                                                 >  2:8456310/1‑137 (MQ=255)

 aTTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCt‑‑gt‑‑g                                                                                                                                 >  2:8093752/1‑137 (MQ=255)

 aTTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCt‑‑gt‑‑g                                                                                                                                 >  2:14733883/1‑137 (MQ=255)

             gTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAATGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGcgcc                                                                                                                     <  1:11228203/137‑1 (MQ=255)

             gTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGcgcc                                                                                                                     <  1:6883313/137‑1 (MQ=255)

             gTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTGCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGcgcc                                                                                                                     <  2:4822333/137‑1 (MQ=255)

              tAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGcgccg                                                                                                                    >  1:944084/1‑137 (MQ=255)

                 tGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCcgcgc                                                                                                                 <  1:6982416/137‑1 (MQ=255)

                 tGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCcgcgc                                                                                                                 <  1:2430802/137‑1 (MQ=255)

                  tcTGATCCCTGCATGGATCGCTGTGCTGTGGGGGGCTTTCAAGCTGAAGAAGGCTCGCAAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTGGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCcgcgcg                                                                                                                <  2:7556647/136‑1 (MQ=255)

                    tGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCgcggg                                                                                                              >  2:1283112/1‑137 (MQ=255)

                      aTCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCGGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCgcgggcg                                                                                                            <  1:4419078/137‑1 (MQ=255)

                         ccGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCGCGGgcggcg                                                                                                         >  2:6044938/1‑137 (MQ=255)

                         ccGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCGCGGgcggcg                                                                                                         >  1:10740482/1‑137 (MQ=255)

                         ccGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCGCGGgcggcg                                                                                                         >  1:7008487/1‑137 (MQ=255)

                         ccGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCAGCGCGGgcggcg                                                                                                         >  1:5132679/1‑137 (MQ=255)

                         ccGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAACGCCGCGCGGgcggcg                                                                                                         <  1:6090179/137‑1 (MQ=255)

                         ccGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGACTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCGCGGgcggcg                                                                                                         >  2:11607878/1‑137 (MQ=255)

                              aTGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGa                                                                                                    >  1:14830943/1‑137 (MQ=255)

                                ggATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGAtc                                                                                                  >  2:16517514/1‑137 (MQ=255)

                                ggATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGGNAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGAtc                                                                                                  <  1:15547370/137‑1 (MQ=38)

                                ggATCGCTGCGCTGTGGTTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGAtc                                                                                                  <  1:3111711/137‑1 (MQ=38)

                                    cGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTcac                                                                                              <  1:11322207/137‑1 (MQ=255)

                                     gctgtgctgtGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTcaca                                                                                             >  1:10421875/1‑137 (MQ=255)

                                         tgctgtgGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCg                                                                                         >  2:5986217/1‑137 (MQ=255)

                                          gctgtgGGTGGCTTTCAAGCTGAAGANGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGa                                                                                        <  1:16661781/137‑1 (MQ=255)

                                          gctgtgGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGa                                                                                        <  2:14831953/137‑1 (MQ=255)

                                          gctgtgGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGa                                                                                        <  2:1066485/137‑1 (MQ=255)

                                          gctgtgGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGa                                                                                        <  1:4588607/137‑1 (MQ=255)

                                          gctgtgGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGa                                                                                        <  1:8538812/137‑1 (MQ=255)

                                           ctgtgGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCGGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGAt                                                                                       <  1:3790232/137‑1 (MQ=255)

                                                ggttgCTTTCAAGTTGAAGAAGGCTCACCATGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGTTTTTCGGCGCCT‑‑GT‑‑GATATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAg                                                                                  <  1:5887201/137‑1 (MQ=37)

                                                ggtggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAg                                                                                  <  2:15292215/137‑1 (MQ=255)

                                                   ggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCt                                                                               >  2:14230640/1‑137 (MQ=255)

                                                    gCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTa                                                                              >  1:16366112/1‑137 (MQ=255)

                                                        tCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCGCGTGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAgg                                                                          <  2:7229145/137‑1 (MQ=255)

                                                        tCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGTGGCGCTCGATCTCACAGGCGATGAAAGACTCAAgg                                                                          >  1:4205300/1‑137 (MQ=38)

                                                        tCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCGCAGGCGATGAAAGCCTAAAgg                                                                          >  1:9064037/1‑137 (MQ=255)

                                                        tCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCCCAGGCGATGAAAGCCTAAAgg                                                                          >  1:13374691/1‑137 (MQ=255)

                                                        tCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAgg                                                                          >  1:6663966/1‑137 (MQ=255)

                                                        tCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAgg                                                                          >  2:13671636/1‑137 (MQ=255)

                                                          aaGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAATGGGTATGCCGTTAGGTTTTCTGCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCGCGTGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCa                                                                        <  2:12438297/137‑1 (MQ=38)

                                                           aGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCAt                                                                       >  1:10528915/1‑137 (MQ=255)

                                                                aagaacGCTCGCCATGCCAAGGGCTGCTCGGGTTTGTTCGTAGGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTc                                                                  <  2:2667551/137‑1 (MQ=37)

                                                                  gaagGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTg                                                                <  2:5672579/137‑1 (MQ=255)

                                                                   aagGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTgg                                                               >  2:3361514/1‑137 (MQ=255)

                                                                      gCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTTAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTggcgg                                                            <  1:7232626/137‑1 (MQ=255)

                                                                       cTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTggcggc                                                           <  1:16572947/137‑1 (MQ=255)

                                                                          gCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCt                                                        >  2:13926904/1‑137 (MQ=255)

                                                                           ccAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTg                                                       <  1:12512753/137‑1 (MQ=255)

                                                                           ccAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTg                                                       <  2:5888870/137‑1 (MQ=255)

                                                                           ccAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTg                                                       <  1:2097671/137‑1 (MQ=255)

                                                                           ccAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTg                                                       <  1:8166786/137‑1 (MQ=255)

                                                                            cAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGa                                                      <  2:5348557/137‑1 (MQ=255)

                                                                             aGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCACTCGATCCCACAGGCGATGAAAACCTAAAAGCATGCGTCTGGCGGCCCTGAc                                                     >  1:181706/1‑137 (MQ=38)

                                                                              ggCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGAcg                                                    <  2:7857965/137‑1 (MQ=255)

                                                                               gCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGAcgc                                                   <  1:1000089/137‑1 (MQ=255)

                                                                                  aaGGGCTGATAGGGGTTGTTCGTAAGGGGTATGCTGTTGGGTTTTCGGCGCCTGCGA‑‑G‑‑ATCGAGCGCTGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGtt                                                <  1:11902114/137‑1 (MQ=37)

                                                                                   aGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCAT‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCAGCCCTGACGCGTTc                                               >  2:2015138/1‑137 (MQ=38)

                                                                                       cTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCt                                           >  2:1324379/1‑137 (MQ=255)

                                                                                        tGATTGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCtt                                          <  2:6485535/137‑1 (MQ=255)

                                                                                          aTCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTcc                                        <  2:16199558/137‑1 (MQ=255)

                                                                                          aTCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTcc                                        <  1:3539213/137‑1 (MQ=255)

                                                                                          aTCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTcc                                        <  2:5137232/137‑1 (MQ=255)

                                                                                          aTCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTcc                                        <  1:5293852/137‑1 (MQ=255)

                                                                                            cGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTcac                                                                                              >  1:3266853/1‑81 (MQ=25)

                                                                                            cGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCAt                                      >  2:10097260/1‑137 (MQ=255)

                                                                                            cGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCAt                                      >  1:14978456/1‑137 (MQ=255)

                                                                                            cGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCAt                                      >  2:2757096/1‑137 (MQ=255)

                                                                                             gggTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATg                                     >  1:5582542/1‑137 (MQ=255)

                                                                                             gggTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATg                                     >  1:15309083/1‑137 (MQ=255)

                                                                                             gggTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATg                                     >  2:10232548/1‑137 (MQ=255)

                                                                                             gggTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATg                                     >  2:2307706/1‑137 (MQ=255)

                                                                                             gggTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATg                                     >  1:14387344/1‑137 (MQ=255)

                                                                                                  tGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGccccc                                >  2:12225534/1‑136 (MQ=255)

                                                                                                    ttCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCa                              <  2:15649248/137‑1 (MQ=255)

                                                                                                     tCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCaa                             >  1:12798707/1‑137 (MQ=255)

                                                                                                     tCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCATAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCaa                             >  2:13512512/1‑137 (MQ=255)

                                                                                                        tAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAAtgc                          <  2:2876934/137‑1 (MQ=255)

                                                                                                            gggTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGcc                      >  2:9066607/1‑137 (MQ=255)

                                                                                                            gggTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGcc                      >  2:6898113/1‑137 (MQ=255)

                                                                                                              gTATGCCGTTAGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTc                    <  1:9782897/137‑1 (MQ=255)

                                                                                                                        aGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcac          <  2:11834474/137‑1 (MQ=255)

                                                                                                                        aGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcac          <  2:14370887/137‑1 (MQ=255)

                                                                                                                        aGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcac          <  1:13991285/137‑1 (MQ=255)

                                                                                                                        aGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcac          <  1:9245445/137‑1 (MQ=255)

                                                                                                                        aGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCATGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcac          <  2:3587898/137‑1 (MQ=255)

                                                                                                                        aGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCACTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcac          <  2:9662108/137‑1 (MQ=255)

                                                                                                                        aGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCGCATGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcac          <  2:16037353/137‑1 (MQ=38)

                                                                                                                        aGGTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCCACGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcac          <  1:15837697/137‑1 (MQ=255)

                                                                                                                         ggTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcacg         <  2:6325307/137‑1 (MQ=255)

                                                                                                                         ggTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcacg         >  1:7627205/1‑137 (MQ=255)

                                                                                                                         ggTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcacg         >  1:484064/1‑137 (MQ=255)

                                                                                                                          gTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcacgc        >  1:13060838/1‑137 (MQ=255)

                                                                                                                          gTTTTCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCCCCGCCAAgcacgc        <  2:4267790/137‑1 (MQ=255)

                                                                                                                           ttttCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcacgca       <  2:12734165/137‑1 (MQ=255)

                                                                                                                           ttttCGGCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcacgca       <  1:2873459/137‑1 (MQ=255)

                                                                                                                                ggCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAgtg  <  1:8855220/137‑1 (MQ=255)

                                                                                                                                ggCGCCT‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAgtg  <  1:5839603/137‑1 (MQ=255)

||

GATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCT‑‑GT‑‑G‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTG  >  AE015451_phaM1/1070110‑1070372

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