Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451_phaM1445,1120CT31.8% 167.1 / 11.9 89intergenic (+85/‑203)PP_0367 CDS (conserved protein of unknown function)/PP_0368 CDS (putative Acyl‑CoA dehydrogenase)–/–
Reads supporting (aligned to +/- strand):  ref base C (31/29);  new base T (4/24);  total (35/54)
Fisher's exact test for biased strand distribution p-value = 9.38e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAG‑ACAACAATTTTGTCTGGT‑CGGG‑CCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAA  >  AE015451_phaM1/444976‑445248
                                                                                                                                         |                                                                                                                                          
cgtcAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAG‑ACAACAATTTTGTc                                                                                                                                            <  1:8177158/137‑1 (MQ=255)
cgtcAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAG‑ACAACAATTTTGTc                                                                                                                                            <  1:8177205/137‑1 (MQ=255)
 gtcAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAG‑ACAACAATTTTGTCt                                                                                                                                           <  1:3970228/137‑1 (MQ=255)
 gtcAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAG‑ACAACAATTTTGTCt                                                                                                                                           <  1:3970267/137‑1 (MQ=255)
 gtcAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAG‑ACAACAATTTTGTCt                                                                                                                                           <  2:2195801/137‑1 (MQ=255)
  tcAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAG‑ACAACAATTTTGTCTg                                                                                                                                          <  2:6067758/137‑1 (MQ=255)
  tcAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAG‑ACAACAATTTTGTCTg                                                                                                                                          <  2:6067804/137‑1 (MQ=255)
    aaGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGGTTGGCCGCGAAGGCGTACGATCCG‑AAAAAAACTTTGTTTGGt                                                                                                                                        >  2:13189842/1‑137 (MQ=255)
      ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAG‑ACAACAATTTTGTCTGGT‑Cg                                                                                                                                     <  2:12288370/137‑1 (MQ=255)
       gCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAG‑ACAACAATTTTGTCTGGT‑Cgg                                                                                                                                    <  2:16126604/137‑1 (MQ=255)
        cATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAG‑ACAACAATTTTGTCTGGT‑Cggg                                                                                                                                   <  2:17618769/137‑1 (MQ=255)
           ccGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGGTTGTCCCTGAACGCGTCCAATCCA‑ACAACAAATTTTTTTGGT‑Ctggacc                                                                                                                                >  2:3954878/1‑134 (MQ=255)
            ctacgaAATGCTTTGCCGGATCACCGACACCTCTTCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAG‑ACAACAATTTTGTCTGGT‑CGGG‑Cctc                                                                                                                              <  2:12452641/135‑1 (MQ=255)
            cgacgaAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAG‑ACAACAATTTTGTCTGGT‑CGGG‑Cctc                                                                                                                              <  1:5154403/137‑1 (MQ=255)
            cgacgaAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAG‑ACAACAATTTTGTCTGGT‑CGGG‑Cctc                                                                                                                              <  2:10385170/137‑1 (MQ=255)
             gacgaAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAG‑ACAACAATTTTGTCTGGT‑CGGG‑Cctct                                                                                                                             <  2:129612/137‑1 (MQ=255)
             gacgaAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAG‑ACAACAATTTTGTCTGGT‑CGGG‑Cctct                                                                                                                             <  1:17668359/137‑1 (MQ=255)
             gacgaAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAG‑ACAACAATTTTGTCTGGT‑CGGG‑Cctct                                                                                                                             <  1:2063131/137‑1 (MQ=255)
             gacgaAATGCTTGGCCGGATCAGCGACACCTGCGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTATGAGTGGGCTTGCCCGCGAATGCGTCCGATCAG‑ACAACAATTTTGTCTGGT‑CGGG‑Cctct                                                                                                                             <  2:16326420/137‑1 (MQ=255)
              acgaAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGATTTTCCTCTAACGCGTCCGCCCAG‑AAAAAAAATTTGTTCGGT‑Ctggtacctt                                                                                                                             >  2:4816065/1‑131 (MQ=255)
                  agggtttGGCTGGATCATCGACACCTACGCCTGAATATTTGCGTCATCCCCTGTTTCAGCGCCTGCTTTGGAGGAGCGGGCTTGCCCGTGAGTGCTCCCGGTAAG‑ACAACAGTTTTGTTTGGT‑GGGG‑CCTTTTcgcg                                                                                                                        <  2:18837935/132‑1 (MQ=255)
                     gCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAG‑ACAACAATTTTGTCt                                                                                                                                           >  1:8182349/1‑117 (MQ=255)
                     gCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAG‑ACAACAATTTTGTCt                                                                                                                                           <  2:8182349/117‑1 (MQ=255)
                     gCTGGGCCGGACCACAAACACCTACCCCTGAATACTTGCTCCCCCCCCTATTTCAACTCCTGCTTTGTAGGAATGGTCTTGCCTGCGAATTCGTCCGATCAA‑ACAACACTTTTGTCTGGT‑CGGG‑CCTCTTCGCggga                                                                                                                     <  2:14093294/137‑2 (MQ=255)
                        tGGCCGGACCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGAGTGCCCACGAACGTGTCTAACCCG‑AAAAAAACTTTGTTTGTT‑TGGG‑CCTCCTCTCGGACAAg                                                                                                                  >  2:9790552/1‑137 (MQ=255)
                          gCCGGATCATCGACACCTACGCCTGAATACTTTCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAT‑ACAACAATTTTGCCTGGT‑CGGG‑CCTCTTCGCGGGCAAGcc                                                                                                                <  2:8504750/137‑1 (MQ=255)
                          gCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTTCCCGCGAATGCGTCCGATCAG‑ACAACAATTTTGTCTGGT‑CGGG‑CCTCTTCGCGGGCAAGcc                                                                                                                <  2:7076484/137‑1 (MQ=255)
                            cGGATCATCGACCCCGACGCCGAAATACTTGCTTCCTCCCCTCTTCCATCTCTTGCTTTGTAGGAGTGGGCTTGCCCGCGAATGCTCCCGCCCAG‑ACAACACTTTTTTTTGGT‑CGGG‑CCTCTTCGCGGACAACCCCg                                                                                                              <  2:12932661/137‑1 (MQ=255)
                                                       cTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCACGCGAATGCGTCCGATCAG‑ACAACAATTTTGTCTGGT‑CGGG‑CCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTggcatggcat                                                                                   >  1:15003563/1‑137 (MQ=255)
                                                             tCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAA‑ACAACAATTTTGTCTGGT‑CGGG‑CCTCCTCGCGGGCAAAACCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGa                                                                             >  1:10512709/1‑137 (MQ=255)
                                                                  cccTATTTCAGCGCCTGCCTTGTAGGAGCGGGGTTTTCCCCGAATGCGTCCTACCCA‑ACAACAAATCTTTCTTGT‑TGGG‑TCTCTTCTCTAGCAAACCCCCCCCCCCCCAAGGTTGTGGCATGGAATGACTGAtcttc                                                                        >  2:8424201/1‑133 (MQ=255)
                                                                        ttCATCACCTGCTTTGGGGGGGTGGGCTTCCCCGCGAGTGCTTCTGCTCAG‑ACAACGTTTTTTTTTGTT‑CTTT‑CTTCTTCGCGGGCAACCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCAtttt                                                                  <  1:17633657/137‑1 (MQ=255)
                                                                               ctttctttGGAGGGGCGGTCTTACCCAAGAGTGTTTCCAACCAA‑AAAACACTTTT‑TTTGGT‑TGGGTCCTTTTCGCGAACAACCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGc                                                           <  2:6380874/131‑1 (MQ=255)
                                                                               cctgctgTGTAGGAGGGGCCTTCCCCAAGAGTGCTTTCGCCAAA‑ACAACATTTTTTTCTTGT‑CGGT‑CCTTTCCAAGGCCAACCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGc                                                           <  2:10459808/137‑1 (MQ=255)
                                                                                   ctttGGAGGAGTGGGCTTTCCCGAGAATGCTCCCGCCCAA‑ACACCATTTTTTTTTGTT‑TGGT‑CCTCTTCGCGGAAACCCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTaa                                                       <  2:4356746/135‑1 (MQ=255)
                                                                                    tttGTGGGATTGGGTTTGCCCGAGAGTGTGTCCGATCAG‑ACAACATTTTTTTTGGTT‑TGGT‑CCTCTTCGTGGGCAATCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTaaa                                                      <  2:12445488/136‑1 (MQ=255)
                                                                                    ttgggaGGAGTGGGCTTACCTGGGAATGCGTACGATAAG‑AAAACAGTTTTTTTTGGT‑TGGG‑CTTCCTCGCGGAAAATCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTaaa                                                      <  2:12299640/132‑1 (MQ=255)
                                                                                    ctggggggAGTGGGCTGGCCAGAGAATGCGTCAGAAAAG‑AAAACAGTTTTGTTTGGT‑TGGT‑CCTCTTCGGGGGCAATCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTaaa                                                      <  2:16489868/131‑1 (MQ=255)
                                                                                    ctgggaGGAGCGGGCTTCCCAGCAAGGACGTCAAACAAG‑AAAACAATTTTTTCTGGT‑TGGG‑CCTTTTTGCGGGAAACCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTaaa                                                      <  2:18732253/132‑1 (MQ=255)
                                                                                    ctTTTAGGAGTGTGCTTGCCCGCAGGTTCTCCTGCTCAG‑ACAACAGTTTTGTTTGGT‑CGTT‑CCTCTTCGCGGGCAACCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTaaa                                                      <  2:18617813/137‑1 (MQ=255)
                                                                                    ctTGGAGGGGGGGGTTGGCCCGTGAGTGTGTACGACAAG‑AAAACACTTTTGTTTGGT‑CGGG‑CTTTTCCGCGGACAATCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTaaa                                                      <  2:15199776/137‑1 (MQ=255)
                                                                                    ctTGGAGGAGCGTTCTTGCCCTCGACTGCTCCCAGCCGG‑ACAACGTTTTTTTTTGTT‑TGTT‑CCTCTTTGCGGGCACCCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTaaa                                                      <  1:1618572/137‑1 (MQ=255)
                                                                                       gTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAG‑ACAACAATTTTGTCTGGT‑CGGG‑CCTCTTCTCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATgg                                                   <  1:2818013/137‑1 (MQ=255)
                                                                                        ggggAGCGGGTTTACCCGCGAATGTGCCCGACCAA‑ACAACATTTTTTTTTGGT‑TGGT‑CCTCTTCGCGGAAAATCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGc                                                  <  2:15770629/135‑1 (MQ=255)
                                                                                        ggggAGCGGGCTTGCCCGTGGGTACTTCAGATCAG‑AAAACGCTTTTGT‑TGGTGTGGT‑CCTTTTCGTGGGCACCCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGc                                                  <  2:3348951/135‑1 (MQ=255)
                                                                                        gaggaggGGGCTTACCCGGGAGTGCTTCCAGTAAG‑AAAAAACTTTGGTTTGGT‑TGGT‑CCTTTCCGCGGACAATCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGc                                                  <  2:8143715/136‑1 (MQ=255)
                                                                                         gggAGTGGGTTTCCCCGAGAGTGCGTCCGATCAA‑AAAACACTTTTGTCTGGT‑CGGT‑CTCTTTCGCGGACAACCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCa                                                 <  2:16994526/136‑1 (MQ=255)
                                                                                         gggAGCGGTCTTTCCGGCAAGTGTGTCAGATAAG‑ACAACAGTTTTTTTGGTT‑TGGT‑CCTTTTCGTGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCa                                                 <  2:11089015/136‑1 (MQ=255)
                                                                                         aGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAG‑ACAACAATTTTGTCTGGT‑CGGG‑CCTCTTCGCGAGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCa                                                 <  1:3531225/137‑1 (MQ=255)
                                                                                             gggggTTTGCCCGCGAGTGCGTCAGATCGA‑ACAACAGTTTTTTTTGGT‑TGGT‑CCTTTTCGGGGACAACCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCAt                                             <  2:234076/135‑1 (MQ=255)
                                                                                             gCGGGTTTTCCCGAGAATGCGTCCGA‑CAGTAAAACATTTTTTTTTGGT‑TGGG‑CCTCTTCGTGGACAATCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCAt                                             <  2:9190507/137‑1 (MQ=255)
                                                                                                tgCTTTCCCGTGACTTCTTCCGCTCAC‑ACAACATTTTTTTTTGTT‑TGTC‑CCTCTTCGATGCCACCCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGt                                          <  1:8993422/136‑1 (MQ=255)
                                                                                                 gCTTGCCCGCGAATGCGTCCGATCAG‑ACAACAACTTTTTCTGGT‑CGGG‑CCTCTTCTCGGGCAAGACCACTCCTCCAAAGGGTTTTTGCATGGTATTACTGATTTGCCCTTTTTCGTGGCGTAAACGGCAACATAGTc                                         >  2:14590159/1‑137 (MQ=255)
                                                                                                         gaaaaTGCGTCCGATCAG‑ACAACAGTTTTTTTTGTT‑TGGT‑CCTCTTCGGGGGCAACCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTg                                 <  2:13537921/134‑1 (MQ=255)
                                                                                                             aTGCGTCCGATCAG‑ACAACAATTTTGTCTGGT‑CGGG‑CCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAg                             >  1:7672357/1‑137 (MQ=255)
                                                                                                             aTGCGTCCGATCAG‑ACAACAATTTTGTCTGGT‑CGGG‑CCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAg                             >  2:729573/1‑137 (MQ=255)
                                                                                                             aTGCGTCCGATCAG‑ACAACAATTTTGTCTGGT‑CGGG‑CCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAg                             >  1:15795619/1‑137 (MQ=255)
                                                                                                              tGCGTCCGATCAG‑ACAACAATTTTGTCTGGT‑CGGG‑CCTCTTCGCGGGCCAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGAAAGa                            >  1:4439984/1‑137 (MQ=255)
                                                                                                              tGCGTCCGATCAG‑ACAACAATTTTGTCTGGT‑CGGG‑CCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGa                            >  2:15492123/1‑137 (MQ=255)
                                                                                                              tGCGTCCGATCAG‑ACAACAATTTTGTCTGGT‑CGGG‑CCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGa                            >  1:7575344/1‑137 (MQ=255)
                                                                                                              tGCGTCCGATCAG‑ACAACAATTTTGTCTGGT‑CGGG‑CCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGa                            >  1:158543/1‑137 (MQ=255)
                                                                                                              tGCGTCCGATCAG‑ACAACAATTTTGTCTGGT‑CGGG‑CCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGa                            >  2:17230200/1‑137 (MQ=255)
                                                                                                              tGCGTCCGATCAG‑ACAACAATTTTGTCTGGT‑CGGG‑CCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGa                            >  2:17858833/1‑137 (MQ=255)
                                                                                                              tGCGTCCGATCAG‑ACAACAATTTTGTCTGGT‑CGGG‑CCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGa                            >  2:18441591/1‑137 (MQ=255)
                                                                                                              ggCGTCCGACCAG‑AAAACATTTTGTTTTGTT‑CGGG‑CTTTTTCGCGGAAAACCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGa                            <  2:4006125/136‑1 (MQ=255)
                                                                                                               ttgTCAGATCAG‑ACAGCGTTTTTGTTTGTG‑TGTT‑CTTCTCCGCGGACAATCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGaa                           <  2:4615561/135‑1 (MQ=255)
                                                                                                                cGTCCGATCAG‑ACAACAATTTTGTCTGGT‑CGGG‑CCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGaaa                          >  1:17452225/1‑137 (MQ=255)
                                                                                                                   ccaataaG‑ACAACACTTTTTTCTGTT‑TGGT‑CTTCTCCGCGGCCACCCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCt                       <  2:10618110/131‑1 (MQ=255)
                                                                                                                   ccGCTCAG‑ACACCGTTTTTTTTTGGT‑TGGT‑CTTTTTCGCGGGCAGTCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCt                       <  1:1432580/137‑1 (MQ=255)
                                                                                                                   ccGATCAG‑ACAACAATTTTGTCTGGT‑CGGG‑CCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCt                       >  1:181471/1‑137 (MQ=255)
                                                                                                                    cGATCAG‑ACAACAATTTTGTCTGGT‑CGGG‑CCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTg                      >  1:10541407/1‑137 (MQ=255)
                                                                                                                    cGACCAG‑ACACCACTTTTTTTTGTT‑TGGT‑CCTCTTCGCGGGCAATCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTg                      <  2:14872021/137‑1 (MQ=255)
                                                                                                                         aG‑ACAACAATTTTGTCTGGT‑CGGG‑CCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTa                 >  1:2459970/1‑137 (MQ=255)
                                                                                                                         aG‑ACAACAATTTTGTCTGGT‑CGGG‑CCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTa                 >  2:6644323/1‑137 (MQ=255)
                                                                                                                            acaacaGTTTGTTCTGTT‑CGGT‑CCTCTTCGCGGACAATCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTaaa               <  2:18692243/137‑1 (MQ=255)
                                                                                                                            aaaacagTTTTGTTTGTT‑GGGT‑CTTCTTCGCGGGCAACCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTaaa               <  2:18797766/135‑1 (MQ=255)
                                                                                                                             caacaaTTTTGTCTGGT‑CGGG‑CCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAAc              >  1:7443857/1‑137 (MQ=255)
                                                                                                                             caacaaTTTTGTCTGGT‑CGGG‑CCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAAc              >  1:6133705/1‑137 (MQ=255)
                                                                                                                               acaaTTTTGTCTGGT‑CGGG‑CCTCTTCGCGGGCAAGCCCGCTCCTCCAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATAGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGt            >  1:4408698/1‑137 (MQ=255)
                                                                                                                               acaaTTTTGTCTGGT‑CGGG‑CCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGt            >  1:4408651/1‑137 (MQ=255)
                                                                                                                               acaaTTTTGTCTGGT‑CGGG‑CCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGt            >  1:4408652/1‑137 (MQ=255)
                                                                                                                               acaaTTTTGTCTGGT‑CGGG‑CCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGt            >  1:5615756/1‑137 (MQ=255)
                                                                                                                               acaaTTTTGTCTGGT‑CGGG‑CCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGt            >  1:5921500/1‑137 (MQ=255)
                                                                                                                               acaaTTTTGTCTGGT‑CGGG‑CCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTAGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGt            >  1:5921471/1‑137 (MQ=255)
                                                                                                                                 gtttttGTTTGGT‑TGGT‑CCCTTTCGAGGACAACCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTaa          <  1:12813914/135‑1 (MQ=255)
                                                                                                                                 agttttttCTGTT‑CGTT‑CCTCTCCTCGGGCAACCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTaa          <  1:13762740/135‑1 (MQ=255)
                                                                                                                                   ttttGTCTGGT‑CGGG‑CCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGt        >  2:3940196/1‑137 (MQ=255)
                                                                                                                                    tttGTCTGTT‑TGGT‑CCTCTTCGCGGACAACCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGtt       <  2:17680116/137‑1 (MQ=255)
                                                                                                                                    tttGTCTGGT‑CGGG‑CCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGtt       >  2:18219252/1‑137 (MQ=255)
                                                                                                                                      tGTCTGGT‑CGGG‑CCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTc     >  2:4536263/1‑137 (MQ=255)
                                                                                                                                       ttCTGTT‑CGGT‑CCTTTTCGCGGACAATCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCa    <  2:4968141/136‑1 (MQ=255)
                                                                                                                                         cTGGT‑CGGG‑CCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGc                                                           >  2:13338980/1‑80 (MQ=255)
                                                                                                                                         cTGGT‑CGGG‑CCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGc                                                           <  1:13338980/80‑1 (MQ=255)
                                                                                                                                         cTGGT‑CGGG‑CCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCaaa  >  2:5792235/1‑137 (MQ=255)
                                                                                                                                         cTGGT‑CGGG‑CCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCaaa  >  2:14469414/1‑137 (MQ=255)
                                                                                                                                         cTGGT‑CGGG‑CCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCaaa  >  2:14469458/1‑137 (MQ=255)
                                                                                                                                         cTGGT‑CGGG‑CCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCaaa  >  2:2248334/1‑137 (MQ=255)
                                                                                                                                         cTGGT‑CGGG‑CCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCaaa  >  1:18803064/1‑137 (MQ=255)
                                                                                                                                         cTGGT‑CGGG‑CCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCaaa  >  1:17016810/1‑137 (MQ=255)
                                                                                                                                         cTGGT‑CGGG‑CCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCaaa  >  1:17016769/1‑137 (MQ=255)
                                                                                                                                         |                                                                                                                                          
CGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAG‑ACAACAATTTTGTCTGGT‑CGGG‑CCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAA  >  AE015451_phaM1/444976‑445248

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