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	29.1%	221.7 / 23.5	141	intergenic (+69/‑84)	hyi CDS (hydroxypyruvate isomerase)/glxR CDS (tartronate semialdehyde reductase)	–/–
Reads supporting (aligned to +/- strand): ref base A (57/38); new base G (0/39); total (59/82)
Fisher's exact test for biased strand distribution p-value = 1.95e-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 141 total aligned reads displayed.

AGTACAAGCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCT‑GAGGAA‑CAACTCGGGACCCTGTGGGAGCGGG‑CTCGCCCGCGAACACCGGCAAAGCCGGT‑GCCATCCACCGAGTTGTCTG‑CTTCGCGGGC‑GAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGCCCATGGC  >  AE015451_phaM1/4889490‑4889756

aGTACAAGCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCT‑GAGGAA‑CAACTCGGGACCCTGTGGGAGCGGG‑CTCGCCCGCGAACACCGGCAAAGCCGGT‑GCCATCCACCGAGtt                                                                                                                                      <  1:8651199/137‑1 (MQ=255)

  tACAAGCCGCTTACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACACACAACGCAATCT‑GAGGAA‑CAACGCGGGACCCTGTGGGAGCGAG‑ATCGCCCGCGAACACCGGCAAAGCCGGG‑GCCAGACACCGAGGTgg                                                                                                                                    <  2:12434650/137‑2 (MQ=37)

  tACAAGCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCT‑GAGGAA‑CAACTCGGGACCCTGTGGGAGCGGG‑GTCGGCCGCGAAAACCAGCCAAAACCGT‑GTCATCCCCCCAGTTGt                                                                                                                                    >  1:14913171/1‑137 (MQ=255)

  tACAAGCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCT‑GAGGAA‑CAACTCGGGACCCTGTGGGAGCGGG‑GGTGACCACGAACAACCAAAAAGACAGT‑GCCACCCCCCAAGTTGt                                                                                                                                    >  1:15957287/1‑137 (MQ=255)

  tACAAGCCGCCGACCACCACCGAAGCGGGCCTGGGCTTGCTGAAAACCCACAACGCAATCT‑GAGGAA‑CAACTCGGGACCCTTTGGGAGCGGG‑CTCGCCCTCGAACACCGGCAAAGCCGGT‑GCCATCCACCGAGTTGt                                                                                                                                    <  2:12434611/137‑1 (MQ=255)

    cAAGCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCT‑GAGGAA‑CAACTCGGGACCCTGTGGGAGCGGG‑CTCGCCCGCGAACACCGGCAAAGCCGGT‑GCCATCCACCGAGTTGTCt                                                                                                                                  <  2:11313480/137‑1 (MQ=255)

       gCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCT‑GAGGAA‑CAACTCGGGACCCTGTGGGAGCGGG‑CTCGCCCGCGAACACCGGCAAAGCCGGT‑GCCATCCACCGAGTTGTCTG‑Ct                                                                                                                              <  1:10049690/137‑1 (MQ=255)

       gCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCT‑GAGGAA‑CAACTCGGGACCCTGTGGGAGCGGG‑CTCGCCCGCGAACACCGGCAAAGCCGGT‑GCCATCCACCGAGTTGTCTG‑Ct                                                                                                                              <  2:3844573/137‑1 (MQ=255)

        ccGCTGACCACCACCGAAGCGGGTCTGGGCTGGCTGAAAACCCACAACGCAATCT‑GAGGAA‑CAACTCGTGACCCTTTGTGAGCGGG‑CTCGCCCGCGAACACTGGCAAAGCCGGT‑GCCATCCACCTAGTTGTCTG‑Ctt                                                                                                                             <  2:1422800/137‑1 (MQ=39)

            tGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCT‑GAGAAA‑CAACTCGGGACCCTGTGGGAGCGGG‑CTCGCCCGCGAACACCGGCAAAGCCGGT‑GCCATCCACCGAGTTGTCTG‑CTTcgcg                                                                                                                         <  1:14036116/137‑1 (MQ=255)

               ccaccaccGAATCGGGCCTGGGTTGGCTGAAAACCCACAACACAATTT‑GAGGAA‑CAATTCGGGACCTTGTGGGAGCGGG‑CTCGCCCGAGAAAACCGGAAAAGCCGGG‑GCCATCCACCGGGTGGTTGG‑CTTCGTggga                                                                                                                      <  2:2813510/137‑2 (MQ=255)

               ccaccaccGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCT‑GAGGAA‑CAACTTGGGACCCTTTGGGAGCGGG‑CTCTCCCGCGAACACAGGAAAAGCCGGG‑GCCATCCACCGGGTTGTTTG‑TTTCGTGGGc                                                                                                                      <  2:8353695/137‑1 (MQ=255)

               ccaccaccGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCT‑GAGGAA‑CAACTCGGGACCCTGTGGGAGCGGG‑ATCGCCCGCGAAAACCGGAAAAGAAGGT‑GACACCCACCAAAGTGTCTT‑CTTCACgggt                                                                                                                      >  2:13116985/1‑136 (MQ=255)

               ccaccaccGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCT‑GAGGAA‑CAACTCGGGACCCTGTGGGAGCGGG‑ATCGCCCACGAACACCTGAAAAGCCAGT‑GTCCTCCCCCCAGTTTTCTT‑CTTTTCgggt                                                                                                                      >  2:16190994/1‑136 (MQ=255)

                      cGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCT‑GAGGAA‑CAACTCGGGACCCTGTGGGAGCGGG‑GTCGCCCGCGAAAACCAGCAAAAGAGGT‑GTCATCCCCCAAGTTGTTTT‑TTTCTCGGGG‑GAGCCCg                                                                                                              >  2:13815745/1‑137 (MQ=255)

                          gCGGGCCTGGGCTGGCGGAAAACCCACAACGCAACCT‑GAGGAA‑CAACTCGGGACCCTGTGGGAGCGGG‑CTCGCCCGCGAACACCGGCAAAGCCGGT‑GCCATCCACCGAGTTGTTTG‑CTTCGCGGGC‑GAGCCCGcgaa                                                                                                          <  1:13851165/137‑4 (MQ=255)

                           cGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCT‑GAGGAA‑CAACTCGGGACCCTGTGGGAGCGGG‑GTCACCCACAAACCACGGCAAAGAAAGTGGCAACCCACCAAGTT‑TCTT‑CTTCCCGGGC‑AAAGCCCCTccc                                                                                                         >  2:3676080/1‑137 (MQ=255)

                           cGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCT‑GAGGAA‑CAACTCGGGACCCTGTGGGAGCGGG‑ATTCCCCCCGAACACCGGCAAAGACGGA‑GCCCTCCCACGCGGTTTTGG‑CTTCCCTGGC‑GAGAACGCTccc                                                                                                         >  2:4285627/1‑137 (MQ=255)

                                cTGGGCTGGCTGAAAACCCACAACGCAATCT‑GAGGAA‑CAACTCGGGACCCTGTGGGAGCGGG‑AGCGCCCACGAAAAACAAAAAAACAAAG‑GCCACCCACCAAGTTGTCTT‑CTTCTCGGGC‑GAGACCGATCCcccaaa                                                                                                    >  2:3367349/1‑135 (MQ=255)

                                      tggctgAAAAACCACAACGCAATCT‑GAGGAA‑CAACTCGGGACCCTGTGGGAGCGGG‑CTCGCTAACGAACACCGGCAAAGACGGT‑GTCACCCCCCGAGTTGTTTT‑CTTCACGGGCGGGGCCCGCCCCCCAAggggcaa                                                                                               >  2:16630461/1‑134 (MQ=255)

                                      tggcgaaaacCCCACAACACACTTT‑GAGGAA‑AACCTGGGGACCTTGGGGGAGTGGT‑CTCCCCCAAAAAAACCGACAAACCTGGT‑TCCGTCCCCCGGTTTTTCTT‑CTTGGTGGGC‑GACCCCGCTCCCACAGGGGGCCa                                                                                              <  2:15065677/131‑1 (MQ=255)

                                      gggcggAAAACCAAAAAC‑CAAACTGGAGAAA‑CAACTTGGGACCTTGGGGGGGAGGG‑TTCGCCCCCGAACACCGAAAAATCTGGT‑CCCACCCACCG‑GTTTTTTGTCTTCGCGGGC‑AATCCCGCTCCCACAGGGGGCCa                                                                                              <  2:10607256/136‑1 (MQ=255)

                                       tgctaAAAACCTACAACACACTCT‑GGGAAA‑CAACTCGGGACCCTGGGGGAGCGGT‑CTCCCCCGTGAACACCGGCAAATCCGGT‑ACCATCCCCCGATTTGTCTC‑CTTCGCGGAC‑GATCCCGCTCCCACAGGGGGCCaa                                                                                             <  2:15451201/136‑1 (MQ=255)

                                       ggcggAAAACCCACACCAC‑GGCT‑TAGGAACCACCTTGGGCCCTTGTGGGAGCGGG‑CTCGCCCGAAAACACTGGCAAGGCGGTT‑GCCCCCCCCTGGGTTGTTTG‑CTTCGCGGGC‑AACCCCGCTCCCACAGGGGGCCaa                                                                                             <  2:8714634/137‑1 (MQ=255)

                                       gacaaaaaaCCCACACCACACCCT‑GGGGAA‑AAACTCGGAACCCGTTGGGAGAGGT‑TTCGCCCGCGAACACCGGAAAACCTGTT‑GCCATCCCCCGAGTTGTTTG‑CTTTGTGGAC‑GAACCCGCTCCCACAGGGGGCCaa                                                                                             <  2:12816687/132‑1 (MQ=255)

                                        gcggAAACCCCACACCGCAATCG‑GAGAAA‑C‑ACTTTGGAACCTGGGGGGGGGGGTCCCCCCCGAGAACACTGGCAACGTTGGT‑TCCACCCCCTGGTTTGTCTT‑CTTCGTGGGC‑CCGCCCGCTCCCACAGGGGGCCaaa                                                                                            <  2:128888/134‑1 (MQ=255)

                                           gAAAATCCACAACGCAATCT‑GAGGAA‑CAACTCGGGACCCTGTGGGAGCGGG‑CTCGCCCTCAAACACCGGCAAAAACAGT‑GCCCCCCCCCGCATTTTTTT‑TTTTGCGGGT‑GAGCCCGCTCCCACAGGGGGccccaact                                                                                         >  2:11775076/1‑131 (MQ=255)

                                            aaacCCCACAACGCAATCT‑GAGGAA‑CACCTCGGGACCTTGTGGGAGAGGG‑CTCGCCCGCGAACACCGGCAAAGCCGGT‑GCCATCCCCCGGGTTGTCTG‑CTTCGCGGGC‑GAGCCCGCTCCCACAGGGGGCCAAAgctt                                                                                        <  2:12045260/137‑1 (MQ=25)

                                            aaaaCCCACACCCCACTCT‑GAGAAC‑CACCTCGGGACCCTGTGGGAGCGGG‑CTCCCCTGCAAACACCGGAAAACTCGGT‑GCCATCCACTGGTTTTTTTG‑TTTCGTGGGC‑AACCCCGCTCCCACAGGGGGCCAAAgctt                                                                                        <  2:7176505/137‑1 (MQ=255)

                                               accaaaaaCGCACTCT‑GAGGAA‑CAACTTGGGACCTTGGGGGAGCGGG‑CTCGCCCAAGAACACAGAAAAAGCGGGT‑GCCACCCACTGGGTTGTTTG‑TTTCGAGGGC‑GACCCCGCTCCCACAGGGGGCCAAAgcttgct                                                                                     <  2:5995637/131‑1 (MQ=255)

                                                    caACTCACTTT‑GGGGAA‑CAACTCGGAACCCTGGGGGAACGTT‑CTCCCCTGAAAACACCTACAAATCTGGT‑GCCCTCCACTGATTTTTCTC‑TTCCGTGGGC‑AAACCCGCTCCCACAGGGGGCCAAAGCTTGCTTgaga                                                                                <  2:2938140/137‑1 (MQ=255)

                                                      aCGCACCCT‑GAGAAA‑CACCTCGGAACCTTGTGGGGGTGGT‑CTCCCCCGAGAACACCGGCAAACCTGGT‑TCCCCCCACCGGTTTGTTGG‑TTTTGGGGGA‑GACCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGata                                                                              <  2:6326920/137‑1 (MQ=255)

                                                               agAGGAA‑CAACTTGGGACCCTGTGGGAGCGGG‑CTCACCCACAAAAACCGACAAAGCCGGT‑GCCATCCACTGATTTGTCTG‑CTTCGCGGGC‑AAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAAtt                                                                      <  2:15075146/136‑1 (MQ=25)

                                                                    aa‑caacTCGGGACCCTGTGGGAGCGGG‑CTCGCCCGGAAACACTGGCAAGTCGGGT‑GCCACCCACCGGGTTTTTTG‑CTTCGGGGGC‑AACCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACaaa                                                                 <  2:709974/137‑1 (MQ=25)

                                                                               gACCCTTTGGGAGCGGG‑CTTGCCCGGGAACACCGGCAAATCCGGT‑GCCATCCACCGAGTTGTTGG‑CTTCGTGGGC‑GAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAAttt                                                       <  2:16498287/137‑1 (MQ=255)

                                                                               gACCCTGTGGGAGCGGA‑CTCGCCCGCGAACACCGGCAAAGCCGGT‑GCCATCCACCGAGTTGTATG‑CTTCGCGGGC‑GAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAAttt                                                       >  2:7798797/1‑137 (MQ=255)

                                                                                aCCCTGGGGGAGTGGG‑CTCTCCCGCAAACACCGGAAAAGCCGGT‑GCAATCCACAGAGTTGTCGT‑TTTCGCGGGC‑GACCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTATAAAGAGGTAATTtc                                                      <  2:3416538/137‑1 (MQ=255)

                                                                                         gAGTGGG‑CTCGCCCGAAAACACGGGCAAAGCCGGT‑GCCACCCACCGGGTTGTTTG‑CTTCGCGGGC‑AACCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTa                                             <  2:3545853/137‑1 (MQ=255)

                                                                                                       cgagAACACAGGCAAGGCGGGT‑GCCATCCACTGGGTTGTCGG‑TTTCGTGGGA‑GGCCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCAt                                <  1:1841652/134‑1 (MQ=255)

                                                                                                       cgaaaaCACCGACAAAGCCGGT‑ACCACCCACCGGTTTGTCTG‑CTTCGCGGGC‑AACCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCAt                                <  1:15605315/133‑1 (MQ=255)

                                                                                                         ggAACACCGGCAAGGCCGGG‑GCCATCCACCAGGTTGTCTG‑TTTCGCGGGC‑AATCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCg                              <  2:14610008/136‑1 (MQ=255)

                                                                                                         cgAACACCGGCAAAGCCGGT‑GCCATCCACCGAGTTGTCTG‑CTTCGCGGGC‑GAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCg                              >  1:10261182/1‑137 (MQ=255)

                                                                                                         agAACAGGGGAAAAGAGGGT‑GCCACCCCCTGGGTTGTTTG‑TTTCGCGGGA‑GAACCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCg                              <  1:5560576/136‑1 (MQ=255)

                                                                                                           aaCACCGACAAGTCCGGT‑CCCCCCCACCGGTTTGTTGT‑CTCCGCGGGT‑GACCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGc                            <  1:2452505/137‑1 (MQ=255)

                                                                                                            gcacCGGCAAAGCCGGT‑GCCATCCACCGAGTTGTCTG‑CTTCGCGGGC‑GAGCCCGCTCCCACAGGGGGCCCAAGCTTGCTTTAGATCTCCCCATTCCCCCGAGGTCCTTTCTCCTGGCTCCACTCTGTTTCCTCGGcc                           >  1:14770484/2‑136 (MQ=19)

                                                                                                            acacCGGCAAAGCCGGT‑GCCCTCCACTGCTTTGCTTG‑TTTCGCGGTC‑GATCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCa                           <  1:4029420/137‑1 (MQ=37)

                                                                                                            acacCGGCAAAGCCGGT‑GCCATCCACCGAGTTGTCTG‑CTTCGCGGGC‑GAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCa                           >  2:1111708/1‑137 (MQ=255)

                                                                                                            aaacCGGAAAAGCTGGT‑GCACCCCCCCGGGTTGTCTG‑CTTCGTGGTC‑GGGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCa                           <  1:4029376/135‑1 (MQ=255)

                                                                                                             cacCGGCAAAGCCGGT‑GCCATCCACCGAGTTGTCTG‑CTTCGCGGGC‑GAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTAAAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCAc                          >  2:7843001/1‑137 (MQ=255)

                                                                                                                 ggCAAAGCTGGT‑GCCATCCACCGGGTTTCCGG‑CTTAGCGGGC‑AACCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGc                      <  2:11409638/137‑1 (MQ=255)

                                                                                                                 ggCAAAGCCGGG‑GCCATCCACCTCGTTTCCTT‑CTCCACGGGC‑GCCCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAACTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGc                      <  2:9808276/137‑1 (MQ=255)

                                                                                                                 gacAAACCCGTT‑GCCATCCACCGAGTTTTCTG‑CTTCGCGGGA‑GAACCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGc                      <  2:5445762/135‑1 (MQ=255)

                                                                                                                  gCAAGGCTGGT‑TCCCTCCACCGGGTTTTCGG‑TTTCGGGTGT‑AAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTATAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGca                     <  2:13585312/137‑1 (MQ=255)

                                                                                                                   cAAAGCCGGT‑GCCATCCACCGAGTTGTCTG‑CTTCGCGGGC‑GAGCCCGCTCCCACAGGGGGCCAAAGCTTNCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGcat                    >  1:14561821/1‑137 (MQ=255)

                                                                                                                   cAAAGCCGGT‑GCCATCCACCGAGTTGTCTG‑CTTCGCGGGC‑GAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGcat                    >  1:8168064/1‑137 (MQ=255)

                                                                                                                   cAAAGCCGGT‑GCCATCCACCGAGTTGTCTG‑CTTCGCGGGC‑GAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGcat                    >  2:10309974/1‑137 (MQ=255)

                                                                                                                   cAAAGCCGGT‑GCCATCCACCGAGTTGTCTG‑CTTCGCGGGC‑GAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGcat                    >  2:6097205/1‑137 (MQ=255)

                                                                                                                   cAAAGCCGGT‑GCCATCCACCGAGTTGTCTG‑CTTCGCGGGC‑GAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGcat                    >  2:6553789/1‑137 (MQ=255)

                                                                                                                     aaGCCGGT‑GCCATCCACCGAGTTGTCTG‑CTTCGCGGGC‑GAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGcatca                  >  2:4071442/1‑137 (MQ=255)

                                                                                                                     aaGCCGGT‑GCCATCCACCGAGTTGTCTG‑CTTCGCGGGC‑GAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGcatca                  >  2:4619227/1‑137 (MQ=255)

                                                                                                                     aaGCCGGT‑GCCATCCACCGAGTTGTCTG‑CTTCGCGGGC‑GAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGcatca                  >  1:14128210/1‑137 (MQ=255)

                                                                                                                     aaGCCGGT‑GCCATCCACCGAGTTGTCTG‑CTTCGCGGGC‑GAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGcatca                  >  2:7172255/1‑137 (MQ=255)

                                                                                                                      atccGGT‑GCCACCCCCCGGTTTGTCTG‑TTGCGCGGGC‑GAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGcatcat                 <  2:1047090/135‑1 (MQ=255)

                                                                                                                       tctggT‑GCCACCCACTGTTTTTTCTT‑CTTCGCGAAC‑GACCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATg                <  1:5398689/134‑1 (MQ=255)

                                                                                                                       gctggT‑GCCATCCACCGGTTTTTCTG‑TTTCGCGGGC‑GACCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATg                <  2:8972339/134‑1 (MQ=255)

                                                                                                                       cctggT‑GCCATCCACCGAGTTGTCTG‑CTTCGCGGGC‑GAACCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATg                <  2:8144905/134‑1 (MQ=255)

                                                                                                                        ctgtt‑GCCCTCCACCGGGTTGTCTG‑TTTCGCGGGC‑AACCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATgg               <  2:2357560/133‑1 (MQ=255)

                                                                                                                         cGGT‑GCCATCCACCGAGTTGTCTG‑CTTCGCGGGC‑GAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATggg              >  2:12567334/1‑137 (MQ=255)

                                                                                                                           gT‑GCCATCCACCGAGTTGTCTG‑CTTCGCGGGC‑GAGCCCGCTCCCACAGGGGGCCACCGCTTGCCTGAGATATCACACTTACACAGAGGTAATTTCTCATGGCTAAAATCAGTTTCATCGGCACCGGCATCATGGGCa            >  1:16039177/1‑137 (MQ=255)

                                                                                                                           gT‑GCCATCCACCGAGTTGTCTG‑CTTCGCGGGC‑GAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCa            >  1:10007016/1‑137 (MQ=255)

                                                                                                                           gT‑GCCATCCACCGAGTTGTCTG‑CTTCGCGGGC‑GAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCa            >  1:7849860/1‑137 (MQ=255)

                                                                                                                            t‑GCCATCCCCCGGGTTGTTGG‑TTTTGCGGGA‑GAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCaa           <  2:16420463/137‑1 (MQ=255)

                                                                                                                              gCCATCCACCGAGTTGTCTG‑CTTCGCGGGC‑GAGCCCGCTCCCACCGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAg          >  1:12180927/1‑137 (MQ=255)

                                                                                                                               ccctCCACGGAGTTGTTTG‑CTTCATGGAT‑AACCCCACTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGc         <  1:8823423/134‑1 (MQ=255)

                                                                                                                               ccctCCACCGGGTTGTTTT‑TTTCGGGTGC‑AACCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGc         <  2:1379039/134‑1 (MQ=255)

                                                                                                                               ccccccACCGGTTTGTTTG‑CTTTGTGGGT‑AATCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGc         <  1:13786101/133‑1 (MQ=255)

                                                                                                                               ccccccACAGAGTTTTTTT‑CTTCACGGAA‑AACCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGc         <  1:8216124/133‑1 (MQ=255)

                                                                                                                               ccacccacCGGGTTGTCTG‑TTTCGGGGGC‑GAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAga         <  2:3151922/137‑2 (MQ=255)

                                                                                                                               ccacccacCAAGTTGCCGG‑CTTCGCGGCC‑AATCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGATGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGc         <  2:1959259/137‑1 (MQ=255)

                                                                                                                               ccATCCCCCGATTTGTCTG‑CTTCGCGGGC‑GATCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGc         <  2:16042840/137‑1 (MQ=255)

                                                                                                                               ccATCCACTGGTTTTTTTG‑TTTCGCGGGC‑AACCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGc         <  1:15255394/137‑1 (MQ=255)

                                                                                                                               ccATCCACGGGGTTGTTTG‑TTTTGTGGGC‑GATCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGTCACCGGCATCATGGGCAAGc         <  2:4704453/137‑1 (MQ=255)

                                                                                                                               ccATCCACCGGGTTGTCTG‑TTTCGGGGGA‑GATCCCGCTCCCACAGGGGGCTAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCTGCATCATGGGCAAGc         <  2:1882868/137‑1 (MQ=255)

                                                                                                                               ccATCCACCGGGTTGTCTG‑CTTCGCGTGC‑GAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGc         <  1:4483833/137‑1 (MQ=255)

                                                                                                                               ccATCCACCGAGTTGTCTG‑CTTCGCGGGC‑GAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACACAGAGGTAATTTCTCCTGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGc         >  1:5686048/1‑137 (MQ=255)

                                                                                                                               ccATCCACCGAGTTGTCTG‑CTTCGCGGGC‑GAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGc         >  1:5686124/1‑137 (MQ=255)

                                                                                                                               ccATCCACCGAGTTGTCTG‑CTTCGCGGGC‑GAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGc         >  1:5780461/1‑137 (MQ=255)

                                                                                                                               ccATCCACCGAGTTGTCTG‑CTTCGCGGGC‑GAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGc         <  1:15749139/137‑1 (MQ=255)

                                                                                                                                cATCCACCGAGTTGTCTG‑TTTCGCGGGC‑GAACCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGcc        <  2:9778504/137‑1 (MQ=255)

                                                                                                                                cATCCACCGAGTTGTCTG‑CTTCGCGGGC‑GAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTCAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGcc        >  1:16602687/1‑137 (MQ=255)

                                                                                                                                cATCCACCGAGTTGTCTG‑CTTCGCGGGC‑GAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGcc        >  2:1698414/1‑137 (MQ=255)

                                                                                                                                cATCCACCGAGTTGTCTG‑CTTCGCGGGC‑GAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGcc        >  1:16602711/1‑137 (MQ=255)

                                                                                                                                cATCCACCGAGTTGTCTG‑CTTCGCGGGC‑GAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGcc        >  2:15939158/1‑137 (MQ=255)

                                                                                                                                cATCCACCGAGTTGTCTG‑CTTCGCGGGC‑GAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGcc        >  1:12780639/1‑137 (MQ=255)

                                                                                                                                cATCCACCGAGTTGTCTG‑CTTCGCGGGC‑GAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGcc        >  1:4624798/1‑137 (MQ=255)

                                                                                                                                cATCCACCGAGTTGTCTG‑CTTCGCGGGC‑GAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGcc        >  1:4624763/1‑137 (MQ=255)

                                                                                                                                     aCCGAGTTGTCTG‑CTTCGCGGGC‑GAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGCCCATgg   >  2:7773884/1‑137 (MQ=255)

                                                                                                                                     aCCGAGTTGTCTG‑CTTCGCGGGC‑GAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGCCCATgg   >  1:14097478/1‑137 (MQ=255)

                                                                                                                                     aCCGAGTTGTCTG‑CTTCGCGGGC‑GAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGCCCATgg   >  2:1424949/1‑137 (MQ=255)

                                                                                                                                      ccGAGTTGTCTG‑CTTCGCGGGC‑GAGCCCGCTCCCACAGGGGGCCCAAGCTTGCTTGCGATATCACAATTAAAAAGAGGTAATTTCTCATGGATAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGCCCATGGc  >  2:2469799/1‑137 (MQ=255)

AGTACAAGCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCT‑GAGGAA‑CAACTCGGGACCCTGTGGGAGCGGG‑CTCGCCCGCGAACACCGGCAAAGCCGGT‑GCCATCCACCGAGTTGTCTG‑CTTCGCGGGC‑GAGCCCGCTCCCACAGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGCCCATGGC  >  AE015451_phaM1/4889490‑4889756

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: `‑`