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.

CGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAA  >  AE015451_phaM1/444976‑445248
                                                                                                                                         |                                                                                                                                          
cgtcAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTc                                                                                                                                            <  1:8177158/137‑1 (MQ=255)
cgtcAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTc                                                                                                                                            <  1:8177205/137‑1 (MQ=255)
 gtcAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCt                                                                                                                                           <  1:3970228/137‑1 (MQ=255)
 gtcAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCt                                                                                                                                           <  1:3970267/137‑1 (MQ=255)
 gtcAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCt                                                                                                                                           <  2:2195801/137‑1 (MQ=255)
  tcAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTg                                                                                                                                          <  2:6067758/137‑1 (MQ=255)
  tcAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTg                                                                                                                                          <  2:6067804/137‑1 (MQ=255)
    aaGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGGTTGGCCGCGAAGGCGTACGATCCGAAAAAAACTTTGTTTGGt                                                                                                                                        >  2:13189842/1‑137 (MQ=255)
      ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCg                                                                                                                                     <  2:12288370/137‑1 (MQ=255)
       gCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCgg                                                                                                                                    <  2:16126604/137‑1 (MQ=255)
        cATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCggg                                                                                                                                   <  2:17618769/137‑1 (MQ=255)
           ccGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGGTTGTCCCTGAACGCGTCCAATCCAACAACAAATTTTTTTGGTCtggacc                                                                                                                                >  2:3954878/1‑134 (MQ=255)
            ctacgaAATGCTTTGCCGGATCACCGACACCTCTTCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCctc                                                                                                                              <  2:12452641/135‑1 (MQ=255)
            cgacgaAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCctc                                                                                                                              <  1:5154403/137‑1 (MQ=255)
            cgacgaAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCctc                                                                                                                              <  2:10385170/137‑1 (MQ=255)
             gacgaAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCctct                                                                                                                             <  2:129612/137‑1 (MQ=255)
             gacgaAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCctct                                                                                                                             <  1:17668359/137‑1 (MQ=255)
             gacgaAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCctct                                                                                                                             <  1:2063131/137‑1 (MQ=255)
             gacgaAATGCTTGGCCGGATCAGCGACACCTGCGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTATGAGTGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCctct                                                                                                                             <  2:16326420/137‑1 (MQ=255)
              acgaAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGATTTTCCTCTAACGCGTCCGCCCAGAAAAAAAATTTGTTCGGTCtggtacctt                                                                                                                             >  2:4816065/1‑131 (MQ=255)
                  agggtttGGCTGGATCATCGACACCTACGCCTGAATATTTGCGTCATCCCCTGTTTCAGCGCCTGCTTTGGAGGAGCGGGCTTGCCCGTGAGTGCTCCCGGTAAGACAACAGTTTTGTTTGGTGGGGCCTTTTcgcg                                                                                                                        <  2:18837935/132‑1 (MQ=255)
                     gCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCt                                                                                                                                           >  1:8182349/1‑117 (MQ=255)
                     gCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCt                                                                                                                                           <  2:8182349/117‑1 (MQ=255)
                     gCTGGGCCGGACCACAAACACCTACCCCTGAATACTTGCTCCCCCCCCTATTTCAACTCCTGCTTTGTAGGAATGGTCTTGCCTGCGAATTCGTCCGATCAAACAACACTTTTGTCTGGTCGGGCCTCTTCGCggga                                                                                                                     <  2:14093294/137‑2 (MQ=255)
                        tGGCCGGACCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGAGTGCCCACGAACGTGTCTAACCCGAAAAAAACTTTGTTTGTTTGGGCCTCCTCTCGGACAAg                                                                                                                  >  2:9790552/1‑137 (MQ=255)
                          gCCGGATCATCGACACCTACGCCTGAATACTTTCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCATACAACAATTTTGCCTGGTCGGGCCTCTTCGCGGGCAAGcc                                                                                                                <  2:8504750/137‑1 (MQ=255)
                          gCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTTCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGcc                                                                                                                <  2:7076484/137‑1 (MQ=255)
                            cGGATCATCGACCCCGACGCCGAAATACTTGCTTCCTCCCCTCTTCCATCTCTTGCTTTGTAGGAGTGGGCTTGCCCGCGAATGCTCCCGCCCAGACAACACTTTTTTTTGGTCGGGCCTCTTCGCGGACAACCCCg                                                                                                              <  2:12932661/137‑1 (MQ=255)
                                                       cTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCACGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTggcatggcat                                                                                   >  1:15003563/1‑137 (MQ=255)
                                                             tCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAAACAACAATTTTGTCTGGTCGGGCCTCCTCGCGGGCAAAACCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGa                                                                             >  1:10512709/1‑137 (MQ=255)
                                                                  cccTATTTCAGCGCCTGCCTTGTAGGAGCGGGGTTTTCCCCGAATGCGTCCTACCCAACAACAAATCTTTCTTGTTGGGTCTCTTCTCTAGCAAACCCCCCCCCCCCCAAGGTTGTGGCATGGAATGACTGAtcttc                                                                        >  2:8424201/1‑133 (MQ=255)
                                                                        ttCATCACCTGCTTTGGGGGGGTGGGCTTCCCCGCGAGTGCTTCTGCTCAGACAACGTTTTTTTTTGTTCTTTCTTCTTCGCGGGCAACCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCAtttt                                                                  <  1:17633657/137‑1 (MQ=255)
                                                                               ctttctttGGAGGGGCGGTCTTACCCAAGAGTGTTTCCAACCAAAAAACACTTTTTTTGGTTGGGTCCTTTTCGCGAACAACCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGc                                                           <  2:6380874/131‑1 (MQ=255)
                                                                               cctgctgTGTAGGAGGGGCCTTCCCCAAGAGTGCTTTCGCCAAAACAACATTTTTTTCTTGTCGGTCCTTTCCAAGGCCAACCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGc                                                           <  2:10459808/137‑1 (MQ=255)
                                                                                   ctttGGAGGAGTGGGCTTTCCCGAGAATGCTCCCGCCCAAACACCATTTTTTTTTGTTTGGTCCTCTTCGCGGAAACCCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTaa                                                       <  2:4356746/135‑1 (MQ=255)
                                                                                    tttGTGGGATTGGGTTTGCCCGAGAGTGTGTCCGATCAGACAACATTTTTTTTGGTTTGGTCCTCTTCGTGGGCAATCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTaaa                                                      <  2:12445488/136‑1 (MQ=255)
                                                                                    ttgggaGGAGTGGGCTTACCTGGGAATGCGTACGATAAGAAAACAGTTTTTTTTGGTTGGGCTTCCTCGCGGAAAATCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTaaa                                                      <  2:12299640/132‑1 (MQ=255)
                                                                                    ctggggggAGTGGGCTGGCCAGAGAATGCGTCAGAAAAGAAAACAGTTTTGTTTGGTTGGTCCTCTTCGGGGGCAATCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTaaa                                                      <  2:16489868/131‑1 (MQ=255)
                                                                                    ctgggaGGAGCGGGCTTCCCAGCAAGGACGTCAAACAAGAAAACAATTTTTTCTGGTTGGGCCTTTTTGCGGGAAACCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTaaa                                                      <  2:18732253/132‑1 (MQ=255)
                                                                                    ctTTTAGGAGTGTGCTTGCCCGCAGGTTCTCCTGCTCAGACAACAGTTTTGTTTGGTCGTTCCTCTTCGCGGGCAACCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTaaa                                                      <  2:18617813/137‑1 (MQ=255)
                                                                                    ctTGGAGGGGGGGGTTGGCCCGTGAGTGTGTACGACAAGAAAACACTTTTGTTTGGTCGGGCTTTTCCGCGGACAATCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTaaa                                                      <  2:15199776/137‑1 (MQ=255)
                                                                                    ctTGGAGGAGCGTTCTTGCCCTCGACTGCTCCCAGCCGGACAACGTTTTTTTTTGTTTGTTCCTCTTTGCGGGCACCCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTaaa                                                      <  1:1618572/137‑1 (MQ=255)
                                                                                       gTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCTCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATgg                                                   <  1:2818013/137‑1 (MQ=255)
                                                                                        ggggAGCGGGTTTACCCGCGAATGTGCCCGACCAAACAACATTTTTTTTTGGTTGGTCCTCTTCGCGGAAAATCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGc                                                  <  2:15770629/135‑1 (MQ=255)
                                                                                        ggggAGCGGGCTTGCCCGTGGGTACTTCAGATCAGAAAACGCTTTTGTTGGTGTGGTCCTTTTCGTGGGCACCCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGc                                                  <  2:3348951/135‑1 (MQ=255)
                                                                                        gaggaggGGGCTTACCCGGGAGTGCTTCCAGTAAGAAAAAACTTTGGTTTGGTTGGTCCTTTCCGCGGACAATCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGc                                                  <  2:8143715/136‑1 (MQ=255)
                                                                                         gggAGTGGGTTTCCCCGAGAGTGCGTCCGATCAAAAAACACTTTTGTCTGGTCGGTCTCTTTCGCGGACAACCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCa                                                 <  2:16994526/136‑1 (MQ=255)
                                                                                         gggAGCGGTCTTTCCGGCAAGTGTGTCAGATAAGACAACAGTTTTTTTGGTTTGGTCCTTTTCGTGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCa                                                 <  2:11089015/136‑1 (MQ=255)
                                                                                         aGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGAGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCa                                                 <  1:3531225/137‑1 (MQ=255)
                                                                                             gggggTTTGCCCGCGAGTGCGTCAGATCGAACAACAGTTTTTTTTGGTTGGTCCTTTTCGGGGACAACCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCAt                                             <  2:234076/135‑1 (MQ=255)
                                                                                             gCGGGTTTTCCCGAGAATGCGTCCGACAGTAAAACATTTTTTTTTGGTTGGGCCTCTTCGTGGACAATCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCAt                                             <  2:9190507/137‑1 (MQ=255)
                                                                                                tgCTTTCCCGTGACTTCTTCCGCTCACACAACATTTTTTTTTGTTTGTCCCTCTTCGATGCCACCCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGt                                          <  1:8993422/136‑1 (MQ=255)
                                                                                                 gCTTGCCCGCGAATGCGTCCGATCAGACAACAACTTTTTCTGGTCGGGCCTCTTCTCGGGCAAGACCACTCCTCCAAAGGGTTTTTGCATGGTATTACTGATTTGCCCTTTTTCGTGGCGTAAACGGCAACATAGTc                                         >  2:14590159/1‑137 (MQ=255)
                                                                                                         gaaaaTGCGTCCGATCAGACAACAGTTTTTTTTGTTTGGTCCTCTTCGGGGGCAACCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTg                                 <  2:13537921/134‑1 (MQ=255)
                                                                                                             aTGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAg                             >  1:7672357/1‑137 (MQ=255)
                                                                                                             aTGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAg                             >  2:729573/1‑137 (MQ=255)
                                                                                                             aTGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAg                             >  1:15795619/1‑137 (MQ=255)
                                                                                                              tGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCCAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGAAAGa                            >  1:4439984/1‑137 (MQ=255)
                                                                                                              tGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGa                            >  2:15492123/1‑137 (MQ=255)
                                                                                                              tGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGa                            >  1:7575344/1‑137 (MQ=255)
                                                                                                              tGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGa                            >  1:158543/1‑137 (MQ=255)
                                                                                                              tGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGa                            >  2:17230200/1‑137 (MQ=255)
                                                                                                              tGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGa                            >  2:17858833/1‑137 (MQ=255)
                                                                                                              tGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGa                            >  2:18441591/1‑137 (MQ=255)
                                                                                                              ggCGTCCGACCAGAAAACATTTTGTTTTGTTCGGGCTTTTTCGCGGAAAACCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGa                            <  2:4006125/136‑1 (MQ=255)
                                                                                                               ttgTCAGATCAGACAGCGTTTTTGTTTGTGTGTTCTTCTCCGCGGACAATCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGaa                           <  2:4615561/135‑1 (MQ=255)
                                                                                                                cGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGaaa                          >  1:17452225/1‑137 (MQ=255)
                                                                                                                   ccaataaGACAACACTTTTTTCTGTTTGGTCTTCTCCGCGGCCACCCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCt                       <  2:10618110/131‑1 (MQ=255)
                                                                                                                   ccGCTCAGACACCGTTTTTTTTTGGTTGGTCTTTTTCGCGGGCAGTCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCt                       <  1:1432580/137‑1 (MQ=255)
                                                                                                                   ccGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCt                       >  1:181471/1‑137 (MQ=255)
                                                                                                                    cGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTg                      >  1:10541407/1‑137 (MQ=255)
                                                                                                                    cGACCAGACACCACTTTTTTTTGTTTGGTCCTCTTCGCGGGCAATCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTg                      <  2:14872021/137‑1 (MQ=255)
                                                                                                                         aGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTa                 >  1:2459970/1‑137 (MQ=255)
                                                                                                                         aGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTa                 >  2:6644323/1‑137 (MQ=255)
                                                                                                                            acaacaGTTTGTTCTGTTCGGTCCTCTTCGCGGACAATCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTaaa               <  2:18692243/137‑1 (MQ=255)
                                                                                                                            aaaacagTTTTGTTTGTTGGGTCTTCTTCGCGGGCAACCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTaaa               <  2:18797766/135‑1 (MQ=255)
                                                                                                                             caacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAAc              >  1:7443857/1‑137 (MQ=255)
                                                                                                                             caacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAAc              >  1:6133705/1‑137 (MQ=255)
                                                                                                                               acaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTCCAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATAGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGt            >  1:4408698/1‑137 (MQ=255)
                                                                                                                               acaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGt            >  1:4408651/1‑137 (MQ=255)
                                                                                                                               acaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGt            >  1:4408652/1‑137 (MQ=255)
                                                                                                                               acaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGt            >  1:5615756/1‑137 (MQ=255)
                                                                                                                               acaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGt            >  1:5921500/1‑137 (MQ=255)
                                                                                                                               acaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTAGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGt            >  1:5921471/1‑137 (MQ=255)
                                                                                                                                 gtttttGTTTGGTTGGTCCCTTTCGAGGACAACCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTaa          <  1:12813914/135‑1 (MQ=255)
                                                                                                                                 agttttttCTGTTCGTTCCTCTCCTCGGGCAACCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTaa          <  1:13762740/135‑1 (MQ=255)
                                                                                                                                   ttttGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGt        >  2:3940196/1‑137 (MQ=255)
                                                                                                                                    tttGTCTGTTTGGTCCTCTTCGCGGACAACCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGtt       <  2:17680116/137‑1 (MQ=255)
                                                                                                                                    tttGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGtt       >  2:18219252/1‑137 (MQ=255)
                                                                                                                                      tGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTc     >  2:4536263/1‑137 (MQ=255)
                                                                                                                                       ttCTGTTCGGTCCTTTTCGCGGACAATCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCa    <  2:4968141/136‑1 (MQ=255)
                                                                                                                                         cTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGc                                                           >  2:13338980/1‑80 (MQ=255)
                                                                                                                                         cTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGc                                                           <  1:13338980/80‑1 (MQ=255)
                                                                                                                                         cTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCaaa  >  2:5792235/1‑137 (MQ=255)
                                                                                                                                         cTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCaaa  >  2:14469414/1‑137 (MQ=255)
                                                                                                                                         cTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCaaa  >  2:14469458/1‑137 (MQ=255)
                                                                                                                                         cTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCaaa  >  2:2248334/1‑137 (MQ=255)
                                                                                                                                         cTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCaaa  >  1:18803064/1‑137 (MQ=255)
                                                                                                                                         cTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCaaa  >  1:17016810/1‑137 (MQ=255)
                                                                                                                                         cTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCaaa  >  1:17016769/1‑137 (MQ=255)
                                                                                                                                         |                                                                                                                                          
CGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAA  >  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: