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.9%	2022.0 / 32.1	528	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 . (3/8); new base G (250/264); total (253/275)
Fisher's exact test for biased strand distribution p-value = 2.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*	AE015451_phaM1	1,070,246	2	.	A	97.9%	1903.4 / 28.0	528	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 . (3/8); new base A (250/267); total (253/275)
Fisher's exact test for biased strand distribution p-value = 2.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 533 total aligned reads displayed.

GATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGT‑‑G‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGG  >  AE015451_phaM1/1070110‑1070373

||

 aTTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCtgt‑‑g                                                                                                                                    >  2:11893224/1‑137 (MQ=255)

  ttCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGT‑‑G‑‑a                                                                                                                                 <  2:10907282/137‑1 (MQ=255)

             gTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGcgcc                                                                                                                      <  1:10200358/137‑1 (MQ=255)

             gTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGcgcc                                                                                                                      <  1:11028659/137‑1 (MQ=255)

             gTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGcgcc                                                                                                                      <  1:11571164/137‑1 (MQ=255)

                 tGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCcgcgc                                                                                                                  <  1:14988389/137‑1 (MQ=255)

                    tGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGNGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCgcggg                                                                                                               >  1:14276357/1‑137 (MQ=255)

                    tGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCgcggg                                                                                                               >  1:7979655/1‑137 (MQ=255)

                    tGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCgcggg                                                                                                               <  2:5155858/137‑1 (MQ=255)

                    tGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCgcggg                                                                                                               >  2:4865224/1‑137 (MQ=255)

                      aTCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCgcgggcg                                                                                                             <  2:1900619/137‑1 (MQ=255)

                       tCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAGGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCTTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCgcgggcgg                                                                                                            <  2:7715328/137‑1 (MQ=39)

                        cccGGCAGGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGggcggc                                                                                                           <  2:7055562/137‑1 (MQ=255)

                         ccGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTCGGGGGTATGCCGTTAGGTTTTAGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGgcggcg                                                                                                          >  1:13727933/1‑137 (MQ=39)

                         ccGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGgcggcg                                                                                                          >  1:13270988/1‑137 (MQ=255)

                         ccGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGgcggcg                                                                                                          >  1:12337883/1‑137 (MQ=255)

                         ccGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGgcggcg                                                                                                          >  2:2971427/1‑137 (MQ=255)

                         ccGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGgcggcg                                                                                                          >  1:7681896/1‑137 (MQ=255)

                         ccGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGgcggcg                                                                                                          >  1:1723171/1‑137 (MQ=255)

                            gCATGGGTCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTc                                                                                                       <  2:9801600/137‑1 (MQ=255)

                              aTGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGAGTCACGCGGGCGGCGGTCGa                                                                                                     >  1:1188454/1‑137 (MQ=255)

                                ggATCGCTGTGCTGTGGGTGGCTTGCAAGCTGAAGAAGGCGCGCCATGTCAAGGGCTGATCGGGGGTGTTTGGAAGGGGTAGGCCGGTAGGTTTGCGGTGCCTGT‑‑G‑‑AGGTCGAGTGCCGCGTGGGTGGCGCTCGAtc                                                                                                   <  1:2342574/137‑1 (MQ=255)

                                   tCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCCca                                                                                                >  1:3727967/1‑137 (MQ=255)

                                    cGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTcac                                                                                               <  1:668339/137‑1 (MQ=255)

                                     gctgtgctgtGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTcaca                                                                                              >  1:17409398/1‑137 (MQ=255)

                                     gctgtgctgtGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTcaca                                                                                              >  2:2160686/1‑137 (MQ=255)

                                     gctgtgctgtGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTcaca                                                                                              >  2:7285511/1‑137 (MQ=255)

                                         tgctgtgGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTTAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCg                                                                                          <  2:2725841/137‑1 (MQ=255)

                                          gctgtgGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGa                                                                                         <  1:15420117/137‑1 (MQ=255)

                                          gctgtgGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGa                                                                                         <  1:10559168/137‑1 (MQ=255)

                                          gctatgGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCATGTGAG‑‑A‑‑TCGAGCACCGCGCGGGCGGCGCTCGATCTCACAGGCGa                                                                                         <  1:16884325/137‑1 (MQ=38)

                                           cggtgggtggCTTTCAGGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGAGCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGAt                                                                                        <  2:2371113/135‑1 (MQ=37)

                                            tgtgGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATg                                                                                       >  2:6543084/1‑137 (MQ=255)

                                                 gtggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGc                                                                                  >  2:6012482/1‑137 (MQ=255)

                                                 gtggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGc                                                                                  >  1:8367065/1‑137 (MQ=255)

                                                 gtggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGc                                                                                  <  1:1520889/137‑1 (MQ=255)

                                                 gtggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGc                                                                                  >  1:7166881/1‑137 (MQ=255)

                                                 gtggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGc                                                                                  >  1:9173051/1‑137 (MQ=255)

                                                   ggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCt                                                                                >  2:1201825/1‑137 (MQ=255)

                                                      tttCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTaaa                                                                             >  2:16509218/1‑137 (MQ=255)

                                                        tCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAgg                                                                           >  1:6379262/1‑137 (MQ=255)

                                                        tCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAgg                                                                           >  2:10409547/1‑137 (MQ=255)

                                                        tCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGATTTGTTCGTAAGAGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TTGAGCGCCGCGCGGGAGGCGCTCGATCTCAAAGGCGATGAACGCCTAAAgg                                                                           >  2:11609654/1‑137 (MQ=255)

                                                                aagaagGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTc                                                                   <  1:11992747/137‑1 (MQ=255)

                                                                aagaagGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTc                                                                   <  1:3568619/137‑1 (MQ=255)

                                                                aagaagGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTc                                                                   <  1:12241473/137‑1 (MQ=255)

                                                                aagaagGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTc                                                                   <  1:17616316/137‑1 (MQ=255)

                                                                       cTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTggcggc                                                            <  1:3308830/137‑1 (MQ=255)

                                                                           ccAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTg                                                        <  1:7510310/137‑1 (MQ=255)

                                                                           ccAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTg                                                        <  1:959584/137‑1 (MQ=255)

                                                                            cAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGa                                                       <  1:9869619/137‑1 (MQ=255)

                                                                             aGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGAc                                                      <  1:8965731/137‑1 (MQ=255)

                                                                             aGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGAc                                                      >  2:12448190/1‑137 (MQ=255)

                                                                                       cTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCt                                            >  1:9266003/1‑137 (MQ=255)

                                                                                       cTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCt                                            >  2:4525502/1‑137 (MQ=255)

                                                                                        tGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCtt                                           <  1:8214125/137‑1 (MQ=255)

                                                                                          aTCGGGTTTTTTCTTAATGGGTATGCCGTTGGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCTAGCGCCGCGCGTGCGGCGCTCGATTTAACAGGCGATGAAGGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTccc                                         <  2:8502742/137‑1 (MQ=255)

                                                                                          aTCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTcc                                         <  1:16731139/137‑1 (MQ=255)

                                                                                             gggTTTTTTCGTCAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCCCTCGGGCGGCGCTCGGTCTCACATGCAATGAACGCCTAACGGCATGTGTCTGGCGACCCTGACGCGTTCCCCTTACAtt                                      >  2:13159948/1‑136 (MQ=25)

                                                                                             gggTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGCAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATg                                      >  2:4327671/1‑137 (MQ=255)

                                                                                             gggTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGGCGCGTTCCCCTTCCATg                                      >  2:15870422/1‑137 (MQ=255)

                                                                                             gggTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATg                                      >  2:13723072/1‑137 (MQ=255)

                                                                                             gggTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATg                                      >  2:12739282/1‑137 (MQ=255)

                                                                                              ggTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGc                                     >  2:9283054/1‑137 (MQ=255)

                                                                                                       gTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAAtg                            >  1:4895180/1‑137 (MQ=255)

                                                                                                         aaGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAAtgct                          <  1:2551319/137‑1 (MQ=255)

                                                                                                              gTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTc                     <  1:16070367/137‑1 (MQ=255)

                                                                                                                aTGTCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCg                   <  1:461941/137‑1 (MQ=255)

                                                                                                                aTGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCg                   <  2:10079383/137‑1 (MQ=255)

                                                                                                                   ccGTTAGGTTTTCGGTGCCCGTGAG‑‑A‑‑TCGAACGCCGCGCGGGCGGCACTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCa                <  2:13888516/137‑1 (MQ=255)

                                                                                                                   ccGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCa                <  1:8814045/137‑1 (MQ=255)

                                                                                                                    cGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCATCAATGCTGCCTCCGCCaa               >  2:14184958/1‑137 (MQ=255)

                                                                                                                    cGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCCCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCaa               >  1:14114834/1‑137 (MQ=255)

                                                                                                                      ttgggTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGTGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgc             <  1:11738406/134‑1 (MQ=38)

                                                                                                                        aGGTTTTCTGCGCCTGT‑‑GATA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACATGCGATGAAAGCCTAAAGGCATGCGCCCGGCCGCCCCGACGCATTCCCCTTCCATGCCCCTCAGTGCTGCATCCGCCAAgtat           <  2:8949895/137‑4 (MQ=25)

                                                                                                                        aGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcac           <  2:5476582/137‑1 (MQ=255)

                                                                                                                        aGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcac           <  1:8627239/137‑1 (MQ=255)

                                                                                                                        aGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcac           <  1:7026347/137‑1 (MQ=255)

                                                                                                                        aGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcac           <  2:8570785/137‑1 (MQ=255)

                                                                                                                        aGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcac           <  1:593749/137‑1 (MQ=255)

                                                                                                                        aGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcac           <  1:4166255/137‑1 (MQ=255)

                                                                                                                         ggTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcacg          >  2:9550495/1‑137 (MQ=255)

                                                                                                                          gTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAggacgc         >  1:17585744/1‑137 (MQ=255)

                                                                                                                          gTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcacgc         >  2:11044109/1‑137 (MQ=255)

                                                                                                                          gTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcacgc         >  2:6937051/1‑137 (MQ=255)

                                                                                                                          gTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcacgc         >  2:7963503/1‑137 (MQ=255)

                                                                                                                          gTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcacgc         <  1:14755051/137‑1 (MQ=255)

                                                                                                                          gTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcacgc         >  1:5409423/1‑137 (MQ=255)

                                                                                                                          gTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcacgc         <  1:5128456/137‑1 (MQ=255)

                                                                                                                          gTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCACCGCGCGGGCGGCGCTCGATCTCCCAGACGCTGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCCCCGCACAgcacga         >  2:15040218/1‑136 (MQ=37)

                                                                                                                           ttttCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcacgca        <  1:12564905/137‑1 (MQ=255)

                                                                                                                           ttttCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcacgca        <  2:17434165/137‑1 (MQ=255)

                                                                                                                           ttttCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcacgca        <  2:15557211/137‑1 (MQ=255)

                                                                                                                           ttttCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcacgca        <  2:14539705/137‑1 (MQ=255)

                                                                                                                           ttttCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcacgca        <  1:12850207/137‑1 (MQ=255)

                                                                                                                           ttttCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcacgca        <  2:11217119/137‑1 (MQ=255)

                                                                                                                           ttttCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcacgca        <  1:4452351/137‑1 (MQ=255)

                                                                                                                                ggCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAgtg   <  1:6266048/137‑1 (MQ=255)

                                                                                                                                ggCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAgtg   >  1:10028723/1‑137 (MQ=255)

                                                                                                                                 gcgcCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAgtgg  >  1:15699183/1‑137 (MQ=255)

||

GATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGT‑‑G‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGG  >  AE015451_phaM1/1070110‑1070373

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