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.7%	2280.6 / 42.4	598	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/11); new base G (280/303); total (283/315)
Fisher's exact test for biased strand distribution p-value = 5.91e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*	AE015451_phaM1	1,070,246	2	.	A	97.6%	2091.5 / 37.8	598	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/11); new base A (274/304); total (283/315)
Fisher's exact test for biased strand distribution p-value = 6.15e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 604 total aligned reads displayed.

ATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGT‑‑G‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGG  >  AE015451_phaM1/1070111‑1070373

||

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

            gTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGcgcc                                                                                                                      <  2:17442444/137‑1 (MQ=255)

                tGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCTAAGGGCTGATCGGGTTTTTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCcgcgc                                                                                                                  <  2:182379/137‑1 (MQ=255)

                tGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCcgcgc                                                                                                                  <  1:14322213/137‑1 (MQ=255)

                     aTCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCgcgggcg                                                                                                             <  1:3981604/137‑1 (MQ=255)

                       cccGGCGTGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGCCTCACCAGGCCAAGGGCTGATCGGGTTTGTTCGGAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGggcggc                                                                                                           <  1:10002448/137‑1 (MQ=38)

                        ccGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGgcggcg                                                                                                          >  1:3066899/1‑137 (MQ=255)

                        ccGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGgcggcg                                                                                                          >  2:4385739/1‑137 (MQ=255)

                        ccGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGgcggcg                                                                                                          >  2:249971/1‑137 (MQ=255)

                        ccGGCATGGATCGCTGTGCTGTGGGCGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGATTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGggggcg                                                                                                          >  1:17539950/1‑137 (MQ=255)

                        ccGGCATGGATCGCCGTGCTGTGGGTGGCTTTCAAGCTGAAGAGGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGgcggcg                                                                                                          >  2:15327921/1‑137 (MQ=39)

                          ggCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCt                                                                                                        >  2:10870069/1‑137 (MQ=255)

                             aTGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGa                                                                                                     >  2:16511150/1‑137 (MQ=255)

                              tGGATCGCTGTGCTGTGGGTGTCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGAt                                                                                                    >  2:7062127/1‑137 (MQ=38)

                               ggATCGCTGTGCTGTGGTTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTTGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGAtc                                                                                                   <  2:8955668/137‑1 (MQ=38)

                                  tCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTca                                                                                                >  2:8162067/1‑137 (MQ=255)

                                   cGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTcac                                                                                               <  2:11902193/137‑1 (MQ=255)

                                        tgctgtgGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCg                                                                                          >  1:12817242/1‑137 (MQ=255)

                                         gctgtgGGTGGCTTTCAAGCTGAAGACGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGa                                                                                         <  1:13969872/137‑1 (MQ=255)

                                         gctgtgGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAATGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGa                                                                                         <  1:10509072/137‑1 (MQ=255)

                                         gctgtgGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGa                                                                                         <  2:9260267/137‑1 (MQ=255)

                                         gctgtgGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGa                                                                                         <  2:12021248/137‑1 (MQ=255)

                                         gctgtgGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGa                                                                                         <  1:5965839/137‑1 (MQ=255)

                                         gctgtgGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGa                                                                                         <  1:4788614/137‑1 (MQ=255)

                                         gctgtgGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGa                                                                                         <  1:2152141/137‑1 (MQ=255)

                                         gctgtgGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGa                                                                                         <  1:13779468/137‑1 (MQ=255)

                                                gtggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGAtc                                                                                                   <  1:12393617/120‑1 (MQ=255)

                                                gtggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGc                                                                                  <  2:12277449/137‑1 (MQ=255)

                                                gtggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGc                                                                                  >  2:4965301/1‑137 (MQ=255)

                                                gtggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGc                                                                                  <  2:7746665/137‑1 (MQ=255)

                                                gtggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGAGGCGCTCGATCTCACCGGCGATGAAAGc                                                                                  >  1:6190816/1‑137 (MQ=38)

                                                 tggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGcc                                                                                 >  1:14762093/1‑137 (MQ=255)

                                                  ggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCt                                                                                >  1:5531501/1‑137 (MQ=255)

                                                     tttCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTaaa                                                                             >  1:8404007/1‑137 (MQ=255)

                                                     tttCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTaaa                                                                             >  1:5710455/1‑137 (MQ=255)

                                                     tttCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGATCGATCTCACAGGCGATGAAAGCCTaaa                                                                             >  1:11035652/1‑137 (MQ=255)

                                                       tCAAGCTGGAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAgg                                                                           >  1:8151783/1‑137 (MQ=38)

                                                       tCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAgg                                                                           >  1:11273021/1‑137 (MQ=255)

                                                       tCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAgg                                                                           >  1:18312033/1‑137 (MQ=255)

                                                       tCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGACGCGCGGGAGGCGCTCGATCTCACAGGCAATGAAAGCACAAAgg                                                                           >  1:18658839/1‑137 (MQ=37)

                                                        cAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGTGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGc                                                                          >  1:7527652/1‑137 (MQ=255)

                                                               aagaagGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTc                                                                   <  1:17324930/137‑1 (MQ=255)

                                                               aagaagGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTc                                                                   <  1:13297615/137‑1 (MQ=255)

                                                                      cTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTggcggc                                                            <  2:11200077/137‑1 (MQ=255)

                                                                        cGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGccc                                                          >  2:18540570/1‑137 (MQ=255)

                                                                          ccAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTg                                                        <  1:15441361/137‑1 (MQ=255)

                                                                            aGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGAc                                                      >  2:1294114/1‑137 (MQ=255)

                                                                            aGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGAc                                                      >  1:3508861/1‑137 (MQ=255)

                                                                             ggCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGAcg                                                     <  2:18891003/137‑1 (MQ=255)

                                                                                 aaGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGtt                                                 >  2:6498689/1‑137 (MQ=255)

                                                                                  aGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTc                                                <  1:7952733/137‑1 (MQ=255)

                                                                                      cTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGCG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCt                                            >  2:16804703/1‑137 (MQ=255)

                                                                                         aTCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTcc                                         <  1:15073466/137‑1 (MQ=255)

                                                                                         aTCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTcc                                         <  1:10803481/137‑1 (MQ=255)

                                                                                         aTCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTcc                                         <  1:11146321/137‑1 (MQ=255)

                                                                                           cGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCAt                                       >  2:15039195/1‑137 (MQ=255)

                                                                                           cGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCAt                                       >  2:2151081/1‑137 (MQ=255)

                                                                                           cGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCACTGACGCGTTCCCCTTCCAt                                       >  1:1951393/1‑137 (MQ=255)

                                                                                            gggTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCCCCGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATg                                      >  1:9192056/1‑137 (MQ=38)

                                                                                            gggTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATg                                      >  1:8712268/1‑137 (MQ=255)

                                                                                            gggTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATg                                      >  2:10589066/1‑137 (MQ=255)

                                                                                             ggTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCTGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGc                                     <  2:15822677/137‑1 (MQ=255)

                                                                                                   ttCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCa                               <  1:16926500/137‑1 (MQ=255)

                                                                                                      gTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAAtg                            >  2:5751637/1‑137 (MQ=255)

                                                                                                      gTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAAtg                            >  1:5161988/1‑137 (MQ=255)

                                                                                                      gTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGACCTCACCGGCGATGAACACCTACAGGCATGCGTCTGGCGGCCCTAACGCATTCCCCTTCCATGCCCCTCCAtg                            >  1:2736989/1‑137 (MQ=37)

                                                                                                           gggTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGcc                       >  2:3597953/1‑137 (MQ=255)

                                                                                                           gggTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGACCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGcc                       >  1:12044490/1‑137 (MQ=255)

                                                                                                                  ccGTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCa                <  1:15670913/137‑1 (MQ=255)

                                                                                                                    gTTAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAg              >  1:3258514/1‑137 (MQ=255)

                                                                                                                     ttAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgc             <  2:8447229/137‑1 (MQ=255)

                                                                                                                     ttAGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgc             <  1:7092852/137‑1 (MQ=255)

                                                                                                                     ttAGGTCTTCGACGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgc             <  2:17788676/137‑1 (MQ=38)

                                                                                                                       aGGTTTTCGGCGCTTGTGAG‑‑ATTGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcac           <  1:11619777/137‑1 (MQ=38)

                                                                                                                       aGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcac           <  1:6821095/137‑1 (MQ=255)

                                                                                                                       aGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcac           <  1:14498610/137‑1 (MQ=255)

                                                                                                                       aGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcac           <  1:6556091/137‑1 (MQ=255)

                                                                                                                       aGGTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcac           <  1:7328077/137‑1 (MQ=255)

                                                                                                                        ggTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCTCGCGAGCGGCTCCCGATCTCCCAAGCAATGAAAGCCTCAAGGCATGCGCCTCTCGGCCCCAACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCACacacg          >  1:16656565/1‑137 (MQ=255)

                                                                                                                        ggTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACGGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATACCCCTCAATGCTGCCTCCGCCAAgcacg          >  2:6097222/1‑137 (MQ=38)

                                                                                                                        ggTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcacg          >  1:16012370/1‑137 (MQ=255)

                                                                                                                        ggTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcacg          >  2:6845067/1‑137 (MQ=255)

                                                                                                                        ggTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcacg          >  2:3997923/1‑137 (MQ=255)

                                                                                                                        ggTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCATAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcacg          >  2:13532273/1‑137 (MQ=255)

                                                                                                                         gTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGTTGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcacgc         <  2:17071053/137‑1 (MQ=38)

                                                                                                                         gTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcacgc         <  2:11415621/137‑1 (MQ=255)

                                                                                                                         gTTTTCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGACCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcacgc         <  1:13545127/137‑1 (MQ=255)

                                                                                                                          ttttCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGTTCGATCTCACCAGCGATGAAAGCCTAAAAGCATGCATCTGGTGGCCATGACGCGTTCCCCTTCCATGCCCCTAAATGCTGCCTCCGCCAAgcacgca        >  1:18851225/1‑137 (MQ=37)

                                                                                                                          ttttCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcacgca        <  1:4357338/137‑1 (MQ=255)

                                                                                                                          ttttCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcacgca        <  1:10651724/137‑1 (MQ=255)

                                                                                                                          ttttCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcacgca        >  1:3718298/1‑137 (MQ=255)

                                                                                                                          ttttCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcacgca        <  2:3126096/137‑1 (MQ=255)

                                                                                                                          ttttCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcacgca        <  2:14602954/137‑1 (MQ=255)

                                                                                                                          ttttCGGCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcacgca        <  2:14168875/137‑1 (MQ=255)

                                                                                                                          ttttCGGCGCCTGT‑‑GAAAGCGAGCGCCGCGCGGGCGGCGCTCGACCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcacgca        <  2:10208577/137‑1 (MQ=38)

                                                                                                                          tttgCGGCGCCTGTGAG‑‑ATCGAGCACCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTGCCATGCCCCTCAATGCTGCCTCCGCCAAgcacgca        <  2:8060256/137‑1 (MQ=38)

                                                                                                                               ggCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAgtg   >  2:3813100/1‑137 (MQ=255)

                                                                                                                               ggCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAgtg   <  2:16127442/137‑1 (MQ=255)

                                                                                                                               ggCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAgtg   <  1:964922/137‑1 (MQ=255)

                                                                                                                               ggCGCCTGTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAgtg   >  2:12612864/1‑137 (MQ=255)

||

ATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGT‑‑G‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGG  >  AE015451_phaM1/1070111‑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: `‑`