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	36.1%	189.3 / 25.8	135	intergenic (+69/‑84)	hyi CDS (hydroxypyruvate isomerase)/glxR CDS (tartronate semialdehyde reductase)	–/–
Reads supporting (aligned to +/- strand): ref base A (46/37); new base G (0/47); total (47/88)
Fisher's exact test for biased strand distribution p-value = 1.59e-12
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 136 total aligned reads displayed.

GGAGTACAAGCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGGAACAACTCGGGACCCTGTGGGAGCGGG‑CTCGCCCGCG‑AACACCGGCAAAGCCGGTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGCCCATGGCT  >  AE015451_phaM1/4889488‑4889757

ggAGTACAAGCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGGAACAACTCGGGACCCTGTGGGGGCGGG‑ATCGCCCCCG‑AAACACGTCAAAGACAGTGCCCTCCACCTAg                                                                                                                                       >  2:5227009/1‑137 (MQ=255)

ggAGTACAAGCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGGAACAACTCGGGACCCTGTGGGAGCGGG‑GTTGGCCCCG‑AACAACAGAAAAACCAGTGCCATCCCCCGaa                                                                                                                                       >  2:17662350/1‑136 (MQ=255)

  aGTACAAGCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGATGAACAACTCGGGACCCTGTGGGAGCGGG‑CTCGCCCGCG‑AACACCGGCAAAGCCGGTGCCATCCACCGAGtt                                                                                                                                     <  2:8490908/137‑1 (MQ=255)

    tACAAGCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGGAACAACTCGGGACCCTGTGGGAGCGGG‑CTCGCCCGCG‑AACACCGGCAAAGCCGGTGCCATCCACCGAGTTGt                                                                                                                                   <  2:15558611/137‑1 (MQ=255)

      cAAGCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGGAACAACTCTGGACCCTGTGGGAGCGGG‑CTCGCCCGCG‑AACACCGGCAAAGCCGGTGCCATCCACCGAGTTGTCt                                                                                                                                 <  2:12462667/137‑1 (MQ=255)

        aGCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGGAACAACTCGGGACCCTGTGGGAGCGGG‑CTCGCCCGCG‑AACACCGGCAAAGCCGGTGCCATCCACCGAGTTGTCTGc                                                                                                                               <  2:5510583/137‑1 (MQ=255)

         gCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGGAACAACTCGGGACCCTGTGGGAGCGGG‑CTCGCCCGCG‑AACACCGGCAAAGCCGGTGCCATCCACCGAGTTGTCTGCt                                                                                                                              <  1:11590431/137‑1 (MQ=255)

         gCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGGAACAACTCGGGACCCTGTGGGAGCGGG‑CTCGCCCGCG‑AACACCGGCAAAGCCGGTGCCATCCACCGAGTTGTCTGCt                                                                                                                              <  2:16555453/137‑1 (MQ=255)

         gCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGGAACAACTCGGGACCCTGTGGGAGCGGG‑CTCGCCCGCG‑AACACCGGCAAAGCCGGTGCCATCCACCGAGTTGTCTGCt                                                                                                                              <  2:10252973/137‑1 (MQ=255)

         gCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGGAACAACTCGGGACCCTGTGGGAGCGGG‑CTCGCCCGCG‑AACACCGGCAAAGCCGGTGCCATCCACCGAGTTGTCTGCt                                                                                                                              <  1:16599980/137‑1 (MQ=255)

         gCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGGAACAACTCGGGACCCTGTGGGAGCGGG‑CTCGCCCGCG‑AACACCGGCAAAGCCGGTGCCATCCACCGAGTTGTCTGCt                                                                                                                              <  1:11590474/137‑1 (MQ=255)

         gCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGGAACAACTCGGGACCCTGTGGGAGCGGG‑CTCGCCCGCG‑AACACCGGCAAAGCCGGTGCCATCCACCGAGTTGTCTGCt                                                                                                                              <  1:15388896/137‑1 (MQ=255)

         gCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGGAACAACTCGGGACCCTGTGGGAGCGGG‑CTCGCCCGCG‑AACACCGGCAAAGCCGGTGCCATCCACCGAGTTGTCTGCt                                                                                                                              <  1:15585508/137‑1 (MQ=255)

                 ccaccaccGAAGTGGGCTTGGGTTGGCGGAAAACCCACAACGCAATTTGAGAAACAACTCGGGACCCTGTGGGAGCGGG‑CTCGCCCGCG‑AACACCGGCAAAGCCGGTGCCATCCACCGATTTGTCTGTTTTGCGGGc                                                                                                                      <  2:17719164/137‑1 (MQ=25)

                 ccaccaccGAAGCGGGCCTGGGCTGGTGGAAAACCCACAACGCAATCTGAGGAACAACTTGGGACCTTGTGGGAGCGGT‑CTCGCCCGCG‑AACACCGGCAAAGCCGGTGCCATCCACTGAGTTGTCTGTTTCGTGGGc                                                                                                                      <  2:13510135/137‑1 (MQ=37)

                 ccaccaccGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGGAACAACTCGGGACCCTGTGGGAACGGG‑ATCGCCCACGAAACACAGACAAA‑ACGGTGCCCCCCCCCGAGGTGTTTGCTTCTCGGGc                                                                                                                      >  2:6314091/1‑137 (MQ=255)

                 ccaccacaGAAGTGGGCTTGGGTGGGCTGAAAACCCACAACGCAATCTGAGGAACAACTTGGGACCCTTTGGGAGCGGG‑CTCGCCCGCG‑AACACCGGCAAAGCCGGTGCCATCCACCGGTTTGTTTGCTTCGCGGGc                                                                                                                      <  2:16135901/137‑1 (MQ=25)

                                  cTGGGCTGGCTGAAAACCCACAACGCAATCTGAGGAACAACTCGGGACCCTGTGGGAGCGGG‑GTCTCCCACA‑CACACCAACAAAACCAGTTTCACCCCCCCAGGTTTCTTCTTCTCTGGTGAAACCCCTCCcccaag                                                                                                     >  2:3902439/1‑135 (MQ=255)

                                  cTGGGCTGGCTGAAAACCCACAACGCAATCTGAGGAACAACTCGGGACCCTGTGGGAGCGGG‑GCCACCCCCA‑AACCACGGGAAAAAAGGGGGCACCCACCAAATTGTCCGCTTCGCGGGGAAGGCCGCTCCCAAAga                                                                                                     >  2:14037204/1‑136 (MQ=255)

                                  cTGGGCTGGCTGAAAACCAACAACGCAATCTGAGGAACAACTCGGGACCCTGTGGGAGCGGG‑CTCGCCCGCG‑AACACCAGCAAAGCAGGTGCCATCCACCAAGGTGTCTGCTTCGTGGGCGAGCCCGCTCCCACCgg                                                                                                     >  2:8797266/1‑137 (MQ=37)

                                        tggcgaaaaaCCCACAACACAATCTGAGAAACACCTCGGGACCTTGGGGGGGCGGG‑CTCGCCCACG‑AACACCGAAAAAGCTGGTGCCACCCACCGGGTTGTTTGTTTTGCGGGCGACCCCGCTCCCACAGGGGGCCa                                                                                               <  2:3584423/131‑1 (MQ=255)

                                        gggtggAAAACCCACAAAGAAGTTTGAGAAACAACTTGGGACCTTGTGGGAGTGGT‑TTCTCCCGAA‑AAACCCGCCAAACCCGGGGCCATCCACTGGTTTGTTTGCTTCTCGGTCAACCCCGCTCCCACAGGGGGCCa                                                                                               <  2:8739231/132‑1 (MQ=255)

                                         ggctgAAACCCCCCAACGCAGTTTGAGAACCACCTTGGGCCTTTGTGGGAGTGGT‑TTCGCCTGAG‑AACACGGGCAAGTTTGGTGCCCTCCACTGGTTTGTTTGCTTCTCGGGCGATCCCGCTCCCACAGGGGGCCaa                                                                                              <  2:18742565/137‑1 (MQ=255)

                                         ggctaAAAACCCACAACGCAATCTGAGAAACAACTCGGGACCCTGTGGGAGTGGG‑CTCCCCCGCA‑AACACCGACAAAGCCGGTGCCACCCACCGGTTTTTTTTTTTCGCGGGCGATCCCGCTCCCACAGGGGGCCaa                                                                                              <  2:16330432/137‑1 (MQ=255)

                                         ggcggAAACCCCCCAACGCAATCGGAGAAACAACTTGGGACCCTTTGGGGGTGGT‑CTCTCCTGCG‑AACACTGGCAAAGCTGGTGCCATCCACCGAGTTTTCTGCTTCGGGGTCGAGCCCGCTCCCACAGGGGGCCaa                                                                                              <  2:12163937/137‑1 (MQ=255)

                                             gAAAACCCACAACGCAATCTGAGGAACAACTCGGGACCCTGTGGGAGCGGG‑GTCGGCCCCGCAAACCCGGAAAAGCCGGGTCACTCCCCCAAGTGTTTTTTTTTTCGGGTGATCCCCCTCCCACAGGGGGCaaaaat                                                                                           >  2:14151453/1‑135 (MQ=255)

                                              aaacCCCAAAACAAAATCTGAGAAAAACCTGGGAACCTTGGGGGAGTGGC‑CTCACCAAAA‑AACACCGAAAAAACCGGTGCCACCCCCCGATTTTTTTGTTCCGTGGCCGACCCCGCTCCCACAGGGGGCCAAAgctt                                                                                         <  2:3122381/137‑1 (MQ=255)

                                              aaacCCCAAAAAGAAATCGGAGAAAAACCTCGGGACCCTTTGGGAGTGGG‑CTCACCCCAA‑AACACCAGCAAAACCGGTGCCCCCCCCCGCTTGGTCTTCTCCTTGGTCGACCCCGCTCCCACAGGGGGCCAAAgctt                                                                                         <  2:12534867/137‑1 (MQ=255)

                                              aaaaCCCACACCACACTCTGAGACACACCCCGGGACCCTGTGGGATTGGT‑CTCACCCGCG‑AACACCGGCAAAGCCGGTTCCTTCCCTCCATTTTTTTTCTTCATGGGCAACCCCGCTCCCACAGGGGGCCAAAgctt                                                                                         <  2:7501666/137‑1 (MQ=255)

                                              aaaaCCCAAAAAGAAATCGGAGAAACACCGTGGAACCTGGGGGGGGGGGG‑CTCGCCAAAA‑AACACCGGAAAAGCTGTTACAATCCACGGGGTTTTTTGTTTTGGGGTCGATCCCGCTCCCACAGGGGGCCAAAgctt                                                                                         <  2:14865181/137‑1 (MQ=255)

                                                aaCCCACAACGAACTCGGGGGAACAATCTGGGACTTTGTGGGAGCGGT‑CCCGCCCGCG‑AACACTGGCAAATCCGGTGCCACCCACCGGGTTTTTGGTTTCGTGGGCGAGCCCGCTCCCACAGGGGGCCAAAgcttgc                                                                                       <  2:4605514/137‑1 (MQ=255)

                                                  cccACAATACACTCGGAGGACCACCTCTGGCCCTTGTGGGACCGGT‑TTCGCCCACA‑AAAACCGAAAAAACCGGTGCAACCCCCCGAGTTGTCGGCTTAGGGGTAGATCCCGCTCCCACAGGGGGCCAAAgcttgctt                                                                                     <  2:2095357/137‑1 (MQ=255)

                                                      aaacaaaaTCGGAGAAACAATTGGGGACTTGGGGGGAGTGGG‑CTCGCCAGAA‑AACACAGGAAAATCTGGTGCCACCCCCTGGTTTTTTTGTTTTGTGGGCGACCCCGCTCCCACAGGGGGCCAAAGCTTGCTTgaga                                                                                 <  2:440174/131‑1 (MQ=255)

                                                            actCTGAGGAACACCTCGGGACCCTGGGGGAGCGGG‑CTCGCCCGCG‑AAAACCGGCAAAGCTGGTGCCAGCCACCGAGTTTTTTGCTTCGCGGGCGGTCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATcac                                                                           <  2:8590530/135‑1 (MQ=255)

                                                                tggggAACAACTCGGGACCTTGTGGGAGCGGT‑CTCGCCCGCG‑AACACCGGCAAACCCGTTACCATCCACCGAGTTGTCTGTTTCGTGGGCGACCCCGCTCCCATAGGGGGCCAAAGCTTGCTTGAGATATCACAAtt                                                                       <  2:17889016/134‑1 (MQ=37)

                                                                tgggaaacgacGCGGGACCTGGGGGGGGTGGG‑TTCGCCCAAG‑AACACTAGCAACCCTGTTGCCCCCCACCGAGTTTTTTGTTTCGGGGGCGATCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGAGATCACAAtt                                                                       <  2:3201015/132‑1 (MQ=255)

                                                                  gggAACAACTCGGAACTCTTTGGGAGCGGT‑CTCGCCCGCG‑AACACCGGAAAAGCCGGTGCCATCCACCGAGTTGTTTGTTCCGTGGGCAGCCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTAc                                                                     <  2:16822620/136‑1 (MQ=25)

                                                                  aGGCACAACTTTGAACCTTGTGGAATGGGG‑CTAGCTGGCG‑AACCCCGGCACACTCGGTGCAATCAAAGTGTTTGTCTGCTTCGTGGGTGATCCCGCGCCCAAAGGGGGCCAAAGCTTGCTTGAGATATCACAATTAc                                                                     <  2:11137908/137‑1 (MQ=255)

                                                                               tactttgtgGGGGTGGG‑CTCGCCTGAA‑AACATTGGCAAAGTGGGTTCCATCCACTTGTTTGTCTTCTCCTCGGCCGCCCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAAttt                                                        <  1:1849595/132‑1 (MQ=255)

                                                                               aaCCTTGTGGGGGGGGG‑CTCACCTGAA‑AAAACCGGAAAGTCTGGTGCCATCCACTGATTTGTCGGCTTTGTGGGCAACCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAAttt                                                        <  2:4639790/136‑1 (MQ=255)

                                                                                  ccctggggAGGGGGTCTCGCCCGAA‑AAAACGGACAAATCTGTTGCAGTCAACCGATTTGTTTGTTCCATGGGCAACCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTtct                                                      <  1:17802670/132‑1 (MQ=255)

                                                                                               tcctttcccGAG‑AACACTGCCAAAGCTGGTCCCCTCCCCTGGTTTTTTTTCTTCTTGGTCGACCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATc                                         <  1:7468112/131‑1 (MQ=255)

                                                                                               gtccctcccACG‑AACACCGGCAAACCCGTTCCCACCCACCGGTTTGTCTGCTCCGCGGGCGAACCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATc                                         <  2:13714575/131‑1 (MQ=255)

                                                                                                     cccGCG‑ACCCCCGTCAGGGCGGGTGCCATCCCTCGGGTTTTTTGTTTCTGGGGCGCCCCCGCTCCCACAGTGGGCCAAAGCTTGCTTGAGATATCACAATTACAACGAGTTAATTTCTCATGGCTAAAATCGGTTTc                                   <  2:509928/137‑1 (MQ=255)

                                                                                                       cgag‑AACACCGGCAAAGCCGTTGCCACCCACTGAGTTGTCTGCTTCGCGGGCGATCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCAt                                 <  2:13450549/134‑1 (MQ=255)

                                                                                                       cctaaaaCACGGGCAAAGTTGGTCCCATCCACCGGTTTTTCTTTTTTGTGGCCAACCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCa                                  <  1:10430851/132‑1 (MQ=255)

                                                                                                            aaCACCGGCAAATCTGTGGCCCCCCACCGGTTTGTTTGTTTCGTGGGCGACCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGc                             <  2:16889154/137‑1 (MQ=255)

                                                                                                             acataggCAAATCTGGTGCCATCCACTGGTTTGTTTGCTTTGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCa                            <  2:7371098/132‑1 (MQ=255)

                                                                                                             acacTGGCAAGGCCGGTGCCGTCCCCCGGGTGGTCGGTTTCGTGGGCGACCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCa                            <  2:11626984/137‑1 (MQ=255)

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

                                                                                                              cacCGGCAAAGCCGGTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCAc                           >  1:9411358/1‑137 (MQ=255)

                                                                                                                  ggaaaaGCCGGTGCCACCCCCCGGTTTGTCTGCTTCGTGGGCGACCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGc                       <  2:5687838/134‑1 (MQ=255)

                                                                                                                  ggCAAAGCGGGGGCCATCCATCGGTTTTTTTGTTTCTCGGGCGACCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGc                       <  1:9206391/137‑1 (MQ=255)

                                                                                                                  gaaaaaTCCGTTGCCCCCCCCCGATTTTTCTTTTTCTCGGTCAACCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGc                       <  1:8681752/134‑1 (MQ=255)

                                                                                                                  aaaaaaGCTGGTGCCACCCACTGGTTTGTCTTCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGTCACCGGc                       <  2:15586636/134‑1 (MQ=255)

                                                                                                                   aaaaaTCCGGGCCCACCCACCGGGTTGCCTGCTTCGTGGGCCATCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCTGTTTCATCGGCACCGGca                      <  1:9433204/135‑1 (MQ=255)

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

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

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

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

                                                                                                                      aaGCCGGTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAACTACAAAGAGGTAATTTCTCATGGCTAACATCGGTTTCATCGGCACCGGcatca                   >  1:12705411/1‑137 (MQ=255)

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

                                                                                                                       gtcaggTACCATCCCCCGTGTTGTCGGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGcatcat                  <  2:9128111/133‑1 (MQ=255)

                                                                                                                        tcctgTGCCATCCACCGATTTGTTTTTTTTGCGGGCGCTCCCGCTCCCACAGTGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATg                 <  2:15694095/133‑1 (MQ=255)

                                                                                                                        tccGTTGCCATCCCCTGGGTTGTTTGTTTAGCGGGAAACCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATg                 <  1:18272970/136‑1 (MQ=255)

                                                                                                                        gCCGTTTCCATCCCCCGGGTGGTCCGCTGCGCGGGCGGTCCCGCTCCCACAGGGGGCCAAAGCTTGCGTGAGATATCACAATTACAAAGAGGTGATTGCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATg                 <  2:6623323/137‑1 (MQ=255)

                                                                                                                        gCCGGTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGTCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATg                 <  1:8605256/137‑1 (MQ=255)

                                                                                                                        cccGGTACCCCCCCCCGCGTTTTCTGCTCCGCGGTCGACCCCGCTCCCACAGGGGGCCAAAGCTTGCCTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATg                 <  1:16298201/136‑1 (MQ=255)

                                                                                                                        accGGGGCAATCCACCGGTTTGTCTGCTTTACGGGCAATCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAAGTACAAAGAGGTAAGTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATg                 <  2:6636080/136‑1 (MQ=255)

                                                                                                                          tggTGCCGTCCACCGGGTTGTCTGCTCCGAGGGCAACCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATggg               <  2:16334666/136‑1 (MQ=255)

                                                                                                                          tggTGCCATCCACCGGGTTTTCTGTTTCGTGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATggg               <  2:2788300/136‑1 (MQ=255)

                                                                                                                           ggTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGc              >  2:14145812/1‑137 (MQ=255)

                                                                                                                              gCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAg           >  2:10566301/1‑137 (MQ=255)

                                                                                                                              gCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAg           >  1:16745719/1‑137 (MQ=255)

                                                                                                                               ccgtCCACTGGGTTGTCTGCTTCGCGGGCGACCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGc          <  2:7818107/134‑1 (MQ=255)

                                                                                                                               ccacccacTGGGTTTTCTGCTTCGCGGGCGATCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGc          <  2:1792686/137‑1 (MQ=255)

                                                                                                                               ccATCCACTGAGTGGTCTGTTTCTCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGc          <  1:2147208/137‑1 (MQ=255)

                                                                                                                               ccATCCACCGGTTTTTTGGTTTCGCGGGAAACCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGc          <  1:16391677/137‑1 (MQ=255)

                                                                                                                               ccATCCACCGGTTTTTCTGCTTCGCGGCCGACCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGc          <  2:17502668/137‑1 (MQ=255)

                                                                                                                               ccATCCACCGGTTTGCTTGCTTCGTGGGCGACCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGc          <  2:861110/137‑1 (MQ=255)

                                                                                                                               ccATCCACCGGGTTGTCTGCTTCGCGGGCGACCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGc          <  2:8806523/137‑1 (MQ=255)

                                                                                                                               ccACCCATCGATTTGTCTTCTTCGCGGGCGACCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGc          <  2:4560442/137‑1 (MQ=255)

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

                                                                                                                                 aTCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGccc        >  1:4109055/1‑137 (MQ=255)

                                                                                                                                 aTCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGccc        >  1:4109083/1‑137 (MQ=255)

                                                                                                                                    cACTGAGTTGTTTACTTCGTGGGCGACCCCGCTCCCACAGGGGGCCAAAGCTTACTTGAGATATCACAATTACAAAGAGATAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGCCCATg     <  2:8335421/137‑1 (MQ=255)

                                                                                                                                    cACCGAGTTGTCTGCTTCGCGGGCGAGTCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGCCCATg     <  2:3705942/137‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

                                                                                                                                      ccGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGCCCATGGc   >  2:12580356/1‑137 (MQ=255)

                                                                                                                                      ccGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGCCCATGGc   >  1:16399314/1‑137 (MQ=255)

GGAGTACAAGCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGGAACAACTCGGGACCCTGTGGGAGCGGG‑CTCGCCCGCG‑AACACCGGCAAAGCCGGTGCCATCCACCGAGTTGTCTGCTTCGCGGGCGAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGCCCATGGCT  >  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: `‑`