BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AE015451_phaM1	4,889,623	0	A	G	40.7%	157.8 / 39.5	141	intergenic (+69/‑84)	hyi CDS (hydroxypyruvate isomerase)/glxR CDS (tartronate semialdehyde reductase)	–/–
Reads supporting (aligned to +/- strand): ref base A (47/33); new base G (0/55); total (48/93)
Fisher's exact test for biased strand distribution p-value = 5.92e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 142 total aligned reads displayed.

GGAGTACAAGCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGGAACAAC‑TCGGGACCCTGTGGGAGC‑GGGCTCGCCCGCGAACACCGGCAAAGCCGGTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGCCCATGGCT  >  AE015451_phaM1/4889488‑4889757

ggAGTACAAGCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGGAACAAC‑TCGGGACCCTGTGGGAGC‑GGGGTCTCCCCCGAACACCAGCAAAACCAGTGTCACCCCCCAAg                                                                                                                                       >  2:13595310/1‑137 (MQ=255)

  aGTACAAGCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGGAACAAC‑TCGGGACCCTGTGGGAGC‑GGGCTCGCCCGCGAACACCGGCAAAGCCGGTGCCATCCACCGAGtt                                                                                                                                     <  1:10632684/137‑1 (MQ=255)

    tACAAGCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGGAACAAC‑TCGGGACCCTGTGGGAGC‑GGGCTCGACCGAGAACACCGGCAAAGCCGGTGTCATCCACCGAGTTGt                                                                                                                                   >  1:8302018/1‑137 (MQ=255)

         gCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGGAACAAC‑TCGGGACCCTGTGGGAGC‑GGGCTCGCCCGCGAACACCGGCAAAGCCGGTGCCATCCACCGAGTTGTCTGCt                                                                                                                              <  1:10725/137‑1 (MQ=255)

         gCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGGAACAAC‑TCGGGACCCTGTGGGAGC‑GGGCTCGCCCGCGAACACCGGCAAAGCCGGTGCCATCCACCGAGTTGTCTGCt                                                                                                                              <  1:8862287/137‑1 (MQ=255)

            gCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGGAACAAC‑TCGGGACCCTGTGGGAGCGGGGGGTGCCCGCGAAAACCGGAAAAAACAAGGCCCCCCACCAAGTTGTCtttttt                                                                                                                            >  1:14191367/1‑132 (MQ=255)

                 ccaccaccaGAGGGGGCCTGGGCGGGTGGAAACCCCAAAACGCAATCTGAGGAACAAC‑TCGGGACCCGGGGGGGGC‑GGGCTTGCCCGCAAACACGGGCAAGGCTGGGGCCGTCCACCGGGTTGTCGGTTTCGTGGGc                                                                                                                      <  2:14534082/137‑1 (MQ=255)

                 ccaccaccaAAGTGGGTCGGGGCTGGCTGAAACCCCACAACACAATCTGAGGACCACC‑TCGGGACCCTGTGGGAAC‑GGGCGCGCCCGCGAACACCGGCAAATCTGTTGCCATCCCCCGATTTTTCTGCTCCGCGGGc                                                                                                                      <  2:17621113/137‑1 (MQ=21)

                 ccaccaccGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGGAACAAC‑TCGGGACCCTGTGGGAGC‑GGGGTCGCCCACAAACAACAGCAAAAACAGTGGCCTCCCCCCAGGTGTTttcttctcaggt                                                                                                                      >  2:13410134/1‑131 (MQ=255)

                 ccaccaccGAAGCGGGCCGGGGCTGGCTGAAAACCCACAACGCAATCTGAGGAACAAC‑TCGGGACCCTGGGGGAGC‑GGGCTCGCCCGCAAACACCGGCAAAGCCGGTGCCATCCACCGAGTTGTTTGCTTCGCGGGc                                                                                                                      <  2:15503130/137‑1 (MQ=255)

                 ccacaacgGAAGTGGGCTTGGGCTGGCGGAAAACCCACAACGCACTCTGAGGAACAAC‑TTGGGACCTTGTGGGAGT‑GGGCTCTCCCGCGAACACCGGCAAATCCGGTGCCACCCACCGAGTTGTCTGCTTCGCGGGc                                                                                                                      <  2:15578870/137‑1 (MQ=25)

                          aaGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGGAACAAC‑TCGGGACCCTGTGGGAGC‑GGGGTCGACCACGCAAACAGGGAAAAACAGTTTCCCCCCCCAAATTGTCTGCTTCCCGGGCGAgaaccct                                                                                                             >  2:7935357/1‑131 (MQ=255)

                                  cTGGGCTGGCTGAAAACCCACAACGCAATCTGAGGAACAAC‑TCGGGACCCTGTGGGAGC‑GGGGTTTTCCGCGAACACCCGCAAAACACGTGTCATCCCCCAAAGCTTCTGCTTCACAGGCGAGCCCGCTccccccga                                                                                                     >  2:15919390/1‑132 (MQ=255)

                                       ctgtctgAAAACCCATAACCCAACCGGAGAAACAAC‑TCGGGACCTTTGGGGAGC‑GGGCTCGCCCGCAAACACCGGCAAAGCTGGTGCCATCCACCGGGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGcc                                                                                                <  2:15290731/137‑1 (MQ=255)

                                        ggttgaaaaaCCCACAACCCAACTTGAGAAACAAC‑TGGGGACCTTTTGGGGGC‑GGGCTCTCCCGCAAACCCCGGCAAATCTGTTGCCATCCACCGGTTTTTTTGTTTCGGGGGCGATCCCGCTCCCACAGGGGGCCa                                                                                               <  2:4462574/131‑1 (MQ=255)

                                         ggcgaaaaaCCCACAACGAAGTTGGAGGAAAAAC‑TGGGGACCTTGGGGGGGC‑GGTTTCTCCCGAGAAAACGGGAAAAGCGGGTGCCACCCACTGGGTTGTTTGTTTCGGGGGCGACCCCGCTCCCACAGGGGGCCaa                                                                                              <  2:2613391/132‑1 (MQ=255)

                                         gacgaaaaaCCCCCAACTCAACTTGAAAAAAAAT‑TCGGGACCTTGTGGGATT‑GGTCTCGCCTGTGAACACTTACAAGACCGTTTCCATCCCCCGGTTTTTTTTTTTCGGGGGCGATCCCGCTCCCACAGGGGGCCaa                                                                                              <  2:12442989/132‑1 (MQ=255)

                                             gAAAACCCACAACGCAATCTGAGGAACAAC‑TCGGGCCCCTGTGGGAGC‑GGGGGCGCCCACCACAACCAACAAAAACAATGCCATCCCCCAAAGTGTCTGTTTCCTCGGGGAGCCCGGTCCCACAGGGGGCCaaaaac                                                                                          >  2:3357611/1‑134 (MQ=255)

                                              aaaaCCAACAACGAAATCGGAGAAACAAC‑TTGGAACCCTGTGGGAGC‑GGGCTCTCCCGAGAACACCGGAAAACCGGGGGCCACCAACCGGGTTGTCGTCTTCGTGGGCGATCCCGCTCCCACAGGGGGCCAAAgctt                                                                                         <  2:16983852/137‑1 (MQ=255)

                                              aaaaCCAAAAACAAAATCTGAGAAAAAAC‑TTGGACCCTTGGGGGAGC‑GGCTTCTCCTGAGAACACCGCAAACTCTGTTCCCGTCCCCTGGTTTTTCTGCTTCGGGGGAAACCCCGCTCCCACAGGGGGCCAAAgctt                                                                                         <  2:4485579/137‑1 (MQ=255)

                                                accccACACCACACTGTGAGAAACACC‑TTGGGACCTTGTGGGAGT‑GGTCCCCCCCGCAAACAACGAAAAGCCCGGTCCAACCCACTGGTGTTTTTGCTTCGCGGGAAATCCCGCTCCCACAGGGGGCCAAAgcttgc                                                                                       <  1:12076437/135‑1 (MQ=255)

                                                 aCCCCCACCACAGTCTGAGGAACAAC‑TTGGGCCCTTGGGGGAGT‑GGGCTCGCCCGAGAACATCGGCAAAGCTGGTGCCGTCCACGGGGTTGTCTGTTTCGTGGGCGACCCCGCTCCCACAGGGGGCCAAAgcttgct                                                                                      <  2:11194805/137‑1 (MQ=255)

                                                  cccACAACGAACTCGGA‑GAACAACTTCGGTCCCTTGGGGGATC‑GGTCCCTCCCGCAAACACCGACAACTCTGTTCCCCTCCCCCGGTTTTTCTTTTTGGAGGGAGACCCCGCTCCCACAGGGGGCCAAAgcttgctt                                                                                     <  2:10705160/137‑1 (MQ=255)

                                                      caACGCACTCGGAGGAACAAC‑TGGGGACCCTGTGGGAGC‑GGTTTCTCCCGGAAACACCGAAAAATCGGGGGCCGCCACCCGGTTTTTTTGTTTCGTGGGCGATCCCGCTCCCACAGGGGGCCAAAGCTTGCTTgaga                                                                                 <  2:1543188/137‑1 (MQ=255)

                                                      caACGCAACTTGAGAACCACT‑TTGGGACCTTTTGGGATC‑GGGCTCTCCCAAGAACACCGGAAAACCTGGTGCCACCCACGGGTTTGTTTTTTTTGGGGGAGATCCCGCTCCCACAGGGGGCCAAAGCTTGCTTgaga                                                                                 <  2:8461990/137‑1 (MQ=255)

                                                        ccGCAATCGGAGAAACAAC‑TTGGGACCTTGTGGGAGC‑GGTCTCTCCCGTAAACACCGGAAAATCTGGTGCCATCCACAGGGTTTTTTGTTTCGGGGGAGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGata                                                                               <  2:13435179/136‑1 (MQ=255)

                                                        acaaagtCTGGGAAACAAC‑CCGAGACCTTGTGGGACA‑GTACCCACCCACAAACACCGGAAAAGCCGGTGACATCCCCCGGGTGGTCTTTTCTGTGTGCGATCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGata                                                                               <  2:11176645/131‑1 (MQ=255)

                                                        aCGCAATTGGAGAAAAAAC‑TCGGGACCCTGTGGGAGT‑GGGCTCGCCCGGAAACACGGGCAAGGCGGGTGCCATCCACTGGGTTTTTTGTTTCGTGGGCGACCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGata                                                                               <  2:13435152/137‑1 (MQ=255)

                                                               gggagaAACAAC‑TTGGGACCTGGGGGGAGG‑GGGTTCGCCGGGGAACACGGAAAAAGCGGGTGCCATCCACCGGTTTGTCTGCTTCGGGGGGGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAAt                                                                        <  2:1185268/135‑1 (MQ=255)

                                                               cggAGGAACACT‑TTGGAACCTTGTGGGGGC‑GGGTTCGCCTACGACCACCGGCAAAGCAGTTGCCACCCCCCGAGTGGCTTTCTTTGCGGGCGATCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAAt                                                                        <  2:10667527/135‑1 (MQ=255)

                                                                tGAGAAACACC‑TTGGGACCCTGGGGGAGT‑GGTCTCGCCAGCGAAAACGGGAAAATCCGGTGCCACCCCCCGGGTTGTTTGTTTCGCGGGCGACCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAAtt                                                                       <  2:3560626/137‑1 (MQ=255)

                                                                 gagaAACAAC‑TTGGGACCCTGTGGGAGG‑TGGCTCGCCTGCGAACACCGGAAAAGCTGGTGCCCCCCCCCGGTTTGTCTTTTTCGCGGGCGATCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTa                                                                      <  2:16906965/137‑1 (MQ=255)

                                                                              gggACCCTGTGGGAGC‑GGGCTCGACCGCGAACACCGGCAAAGCCAGTGCCACCCACCGAGTTGTCCGATTCGAAGGCGAGCCCGCTCCCAAAGGGGGCCAAAACTTTATTGAGACATCAAACTTAAAAAGAGGGAAt                                                          >  2:8921697/1‑137 (MQ=21)

                                                                                gcccttgtgGGAGT‑GGTCTCTCCCGAAAACACCGACAAATCTGGTGCCCCCCCCCGTGTTTTCTGTTTTGCGGCCGACCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAAttt                                                        <  1:4181323/132‑1 (MQ=255)

                                                                                gacctggtgGGAGG‑GGGTTCGCCGGAAAACACCGGAAAAGTGGGGTCCATCCACTGGTTTTTTTGTTTTGTGGGCGATCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAAttt                                                        <  1:4181357/131‑1 (MQ=255)

                                                                                gACCCTGGGGGAGT‑GGGCTAGCCCGCGAACACAGACAAAGCTGGTGCCATCCACCGGTGTGTCTGCTCCGCGGGCGACCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAAttt                                                        <  2:4169273/137‑1 (MQ=255)

                                                                                aaCCCTGGGGGAGC‑GGCCTCTCCCGAAAAACCTGCCAGATCGGGTGCCATCCCTCGGGTTGTTTTTTTCGGGGTCGACCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAAttt                                                        <  1:17589284/136‑1 (MQ=255)

                                                                                 acccgttgGGAAA‑GGGCTCGCTCGCGAACACCGAAAAAGCCGTTGCCATCCACCGCGTTGTTGGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTtc                                                       <  2:1809565/131‑1 (MQ=255)

                                                                                      ttgGGGGT‑GGGCTCTCCCGAGAACACCGGAAAATCCGGTGCCACCCACCGGGTTGTCTGCTCCGCGGGTAATCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATg                                                  <  2:4337427/136‑1 (MQ=255)

                                                                                            atgggtcTCGCCGGAGAAAACTGAAAAACCTGGTTCCACCCACTGGGTTTTCTGCTTCGGGGTAGACCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTaa                                             <  2:12083319/131‑1 (MQ=255)

                                                                                                          agAAAACCGGCAAACCCGGTGCCCCTCACCGATTCGTCTGTTTGGTGGGCGACCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTGCTCATGGCTAAAATCGGTTTCATCg                               <  2:4986619/136‑1 (MQ=255)

                                                                                                           gAACACCGACAAAGCTGGTGCCATCCACCGATTTGTCTGTTTCTCGGGCAACCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAATAGGTAATTTCTCATGGCTAAAATCGGTTTCATCgg                              <  2:8770859/137‑1 (MQ=255)

                                                                                                             acacCGGCAAAGCCGGTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCa                            >  2:1989421/1‑137 (MQ=255)

                                                                                                             acacCGGCAAAGCCGGTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCa                            >  2:14182842/1‑137 (MQ=255)

                                                                                                              cacggacAAGGCTGGTGCAATCGACAGGGTGGTCTGTTTCCTGGGGGATCCCGCCCCCACAGGGGGCCAAAGCTTGCTGGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCAc                           <  2:15157826/131‑1 (MQ=255)

                                                                                                              cacTGGCAAATCCGGGGCCATCCACCGGTTTGTCTGCTTTGTGGGCGATCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCAc                           <  2:8409166/137‑1 (MQ=255)

                                                                                                              cacGGGAAAGGCGGGGGCCATCCACCGGGTTGTTGGCTTCGCGGGCGACCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCAc                           <  2:13175053/137‑1 (MQ=255)

                                                                                                                  ggCAAATCTGGTGCCACCAACTGGGTTGTTTGCTTTGTGGGTGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGc                       <  2:16351411/137‑1 (MQ=255)

                                                                                                                  ggCAAAGCCGGTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGc                       >  2:11372732/1‑137 (MQ=255)

                                                                                                                  ggCAAAGCCGGTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGc                       >  1:5106342/1‑137 (MQ=255)

                                                                                                                  gacAAATCTGGTCCCCTCCCCCGATTTTTCTTTTTCGCGGTCGACCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGc                       <  1:13321080/135‑1 (MQ=255)

                                                                                                                  gaaaaaGCTGGTACCATCCACCGGTTTGTCTGCTTCGCGGGCAACCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGc                       <  2:7557571/134‑1 (MQ=255)

                                                                                                                   gaaaaGCCGGTGCCATCCACCGGGTTTTTTGTTTTGGGGGCGACCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGca                      <  2:6585136/135‑1 (MQ=255)

                                                                                                                    cAAAGCCGGTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGcat                     >  1:1848270/1‑137 (MQ=255)

                                                                                                                    cAAAGCCGGTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGcat                     >  2:18263294/1‑137 (MQ=255)

                                                                                                                    cAAAGCCGGTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGcat                     >  2:4036973/1‑137 (MQ=255)

                                                                                                                    cAAAGCCGGTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCCAAAATCGGTTTCATCGGCACCGGcat                     >  2:5947355/1‑137 (MQ=255)

                                                                                                                      aaGCCGGTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGcatca                   >  2:16163449/1‑137 (MQ=255)

                                                                                                                      aaGCCGGTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCAACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGcatca                   >  2:16163448/1‑137 (MQ=255)

                                                                                                                       agatggTGCCGTCCACCGGGTTGTTTGTTTCGTGGGCGAGCCCGCTCCCACAGGGGGCCAAAGTTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGcatcat                  <  2:11625431/133‑1 (MQ=255)

                                                                                                                        tctggTGCCATCCACGGGGTTTTTTTTTTGGTGGGCAACCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATg                 <  1:7682245/134‑1 (MQ=255)

                                                                                                                        tctggTGCCACCCACCGAGTTGTCTGTTTCGGGGGAGACCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATg                 <  2:5551108/134‑1 (MQ=255)

                                                                                                                        cccGGTGCCCTCCACCGGGTTTTCTGTTTCGCGGGAGACCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATg                 <  2:15535603/136‑1 (MQ=255)

                                                                                                                        ccaggTCCCATCCACCGGGTTGTTTGCTCCCCGGCCGATCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATg                 <  2:15845520/134‑1 (MQ=255)

                                                                                                                         ccGGTCCCCTCCCCCGCTTTGTTTTCTTCTCGAGCGACCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATgg                <  1:13480129/137‑1 (MQ=255)

                                                                                                                           ggTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGc              <  2:14777315/137‑1 (MQ=255)

                                                                                                                            gTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCa             >  1:1877416/1‑137 (MQ=255)

                                                                                                                            gTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCa             >  1:9579487/1‑137 (MQ=255)

                                                                                                                            gTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCa             >  2:5278404/1‑137 (MQ=255)

                                                                                                                            gTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCa             >  1:6727253/1‑137 (MQ=255)

                                                                                                                            gTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCa             >  1:1676028/1‑137 (MQ=255)

                                                                                                                             ttccATCCACCGGTTTTTTTGTTTCGCGGGCAACCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCaa            <  2:5074357/135‑1 (MQ=255)

                                                                                                                               ccctCCACTGGGTTTTTTGCTTCGGGGGCGATCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGc          <  2:1156009/134‑1 (MQ=255)

                                                                                                                               ccctCCACCGAGTTTTCTGTTTCGTGGGCGATCCCGCTCCCACAGTGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGc          <  2:8066242/134‑1 (MQ=255)

                                                                                                                               ccccccACCGGGTTTTCTGATTCGCGGACAACCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGc          <  1:12415694/133‑1 (MQ=255)

                                                                                                                               ccacccacCGGTTTGTCTGTTTCGTGGGAGACCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGc          <  1:8919950/137‑1 (MQ=255)

                                                                                                                               ccacccacCGGGTTTTTTGCTGCGCGGGAGAGCCAGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAATAGGTAATTTTTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGc          <  2:13039655/137‑1 (MQ=255)

                                                                                                                               ccATCCCCCGATTTGTCTTCTTCGTGGGAGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGGGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGc          <  2:12827568/137‑1 (MQ=255)

                                                                                                                               ccATCCACCGATTTGTCTGCTTCGCGGGCGATCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGc          <  2:12827567/137‑1 (MQ=255)

                                                                                                                               ccATCCACCGAGTTGTCTGCTTCAAGGGCGACCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGc          <  2:13597437/137‑1 (MQ=255)

                                                                                                                                cgtCCACCGAGTTTTTTCTTTTTTGGACGAACCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGcc         <  2:2457323/135‑1 (MQ=255)

                                                                                                                                cacccACTGATTTGTCTTTTTCTTGGTCGACCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGcc         <  2:18236653/134‑1 (MQ=255)

                                                                                                                                cATCCACCGAGTTGTCTGCTTCGCGGGTGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGcc         <  2:7976594/137‑1 (MQ=255)

                                                                                                                                cATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGcc         >  1:4573119/1‑137 (MQ=255)

                                                                                                                                cATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGcc         >  1:17111708/1‑137 (MQ=255)

                                                                                                                                cATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGcc         >  1:3676217/1‑137 (MQ=255)

                                                                                                                                cATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGcc         >  1:4909238/1‑137 (MQ=255)

                                                                                                                                cATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGcc         >  2:11743681/1‑137 (MQ=255)

                                                                                                                                cATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGcc         >  2:11743705/1‑137 (MQ=255)

                                                                                                                                  tCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGCCCa       >  1:3876261/1‑137 (MQ=255)

                                                                                                                                    cACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGCCCATg     >  1:67885/1‑137 (MQ=255)

                                                                                                                                     aCCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCCTGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGCCCATgg    >  2:14098705/1‑137 (MQ=255)

                                                                                                                                     aCCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGCCCATgg    >  1:12922994/1‑137 (MQ=255)

                                                                                                                                     aCCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGCCCATgg    >  1:4560257/1‑137 (MQ=255)

                                                                                                                                     aCCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGCCCATgg    >  1:4396183/1‑137 (MQ=255)

                                                                                                                                     aCCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGCCCATgg    >  1:14514701/1‑137 (MQ=255)

                                                                                                                                     aCCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGCCCATgg    >  1:3187630/1‑137 (MQ=255)

                                                                                                                                     aCCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGCCCATgg    >  1:11067638/1‑137 (MQ=255)

                                                                                                                                     aCCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGCCCATgg    >  1:1609953/1‑137 (MQ=255)

                                                                                                                                     aCCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGCCCATgg    >  1:1609926/1‑137 (MQ=255)

GGAGTACAAGCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGGAACAAC‑TCGGGACCCTGTGGGAGC‑GGGCTCGCCCGCGAACACCGGCAAAGCCGGTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGCCCATGGCT  >  AE015451_phaM1/4889488‑4889757

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