Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451_phaM15,199,6560GT28.7% 121.6 / 30.9 130intergenic (+31/+42)PP_4579 CDS (Transcriptional regulator_ LysR family)/PP_4580 CDS (conserved protein of unknown function)–/–
Reads supporting (aligned to +/- strand):  ref base G (21/61);  new base T (22/11);  total (53/77)
Fisher's exact test for biased strand distribution p-value = 8.48e-05
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 140 total aligned reads displayed.

TGGCAGACCGGGGTGGCGTTGGTCGATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTC‑TCCGCCTGTAGCGG‑CCCTTTCGCGGGCATGCCCGCGAAAGGG‑CCGGTACAGGTACAG‑CAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGCGCCGCCCGCCCAGGATCAGCCAGCCAATCCCCAGCAACAGGCTTTCGACCACAAAG  >  AE015451_phaM1/5199521‑5199782
                                                                                                                                         |                                                                                                                                
tggcAGACCGGGGTGGCGTTGGTCGATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTC‑TCCGCCTGTAGCGG‑CCCTTTcgc                                                                                                                                 <  1:13451738/137‑1 (MQ=255)
tggcAGACCGGGGTGGCGTTGGTCGATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTC‑TCCGCCTGTAGCGG‑CCCTTTcgc                                                                                                                                 <  1:4347836/137‑1 (MQ=255)
tggcAGACCGGGGTGGCGTTGGTCGATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTA‑TCCGCCTGTAGCGG‑CCCCTTcgc                                                                                                                                 >  2:16876094/1‑137 (MQ=255)
 ggcAGACCGGGGTGGCGTTGGTCGATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTC‑TCCGCCTGTAGCGG‑CCCTTTcgcg                                                                                                                                <  1:10459431/137‑1 (MQ=255)
 ggcAGACCGGGGTGGCGTTGGTCGATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTC‑TCCGCCTGTAGCGG‑CCCTTTcgcg                                                                                                                                <  2:11683818/137‑1 (MQ=255)
 ggcAGACCGGGGTGGCGTTGGTCGATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTC‑TCCGCCTGTAGCGG‑CCCTTTcgcg                                                                                                                                <  1:18959518/137‑1 (MQ=255)
 ggcAGACCGGGGTGGCGTTGGTCGATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTC‑TCCGCCTGTAGCGG‑CCCTTTcgcg                                                                                                                                <  1:542665/137‑1 (MQ=255)
   cAGACCGGGGTGGCGTTGGTCGATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTC‑TCCGCCTGTAGCGG‑CCCTTTCGCggg                                                                                                                              >  2:11832960/1‑137 (MQ=255)
    aGACGGGGGTGGCGTTGGTAGATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGTGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTC‑TCCGCCTGTAGCGG‑TCCTTTCGCgggg                                                                                                                             >  2:3694177/1‑136 (MQ=255)
    aGACCGGGGTGGCGTTGGTCGATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGTTGGTCTGACCAACAGTC‑TCCGCCTGTAGCGG‑CCCTTTCGCGGGc                                                                                                                             <  1:11522356/137‑1 (MQ=255)
    aGACCGGGGTGGCGTTGGTCGATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTC‑TCCGCCTGTAGCGG‑CCCTTTCTCGGGc                                                                                                                             >  2:7528079/1‑137 (MQ=255)
    aGACCGGGGTGGCGTTGGTCGATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTC‑TCCGCCTGTAGCGG‑CCCTTTCGCgggg                                                                                                                             >  2:18879988/1‑136 (MQ=255)
    aGACCGGGGTGGCGTTGGTCGATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTC‑TCCGCCTGTAGCGG‑CCCTTTCGCGGGc                                                                                                                             <  2:17750770/137‑1 (MQ=255)
    aGACCGGGGTGGCGTTGGTCGATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTC‑TCCGCCTGTAGCGG‑ACCTTTCTCggga                                                                                                                             >  1:16004857/1‑136 (MQ=255)
     gACCGGTGTCGCGTTGGTCGATGAGGTGGTTGGGGTGTGCCGGTAGGTGTTGGCGCGTTAGGCTTGGGATGTGTGTATTTCGCGGATCGTGCTGTCCTGACCAATAGTC‑TCCGCCTGTAGCGT‑CCCTTTCGCGGGCa                                                                                                                            <  2:4478941/137‑1 (MQ=255)
     gACCGGGGTGGCGTTGGTCGATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTT‑TCCGCCTGTAGCGG‑CCCTTTCGTgggtg                                                                                                                            >  2:16426545/1‑135 (MQ=255)
     gACCGGGGTGGCGTTGGTCGATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTC‑TCCGCCTGTAGCGG‑CCCCTTCTCgggaa                                                                                                                            >  2:663917/1‑135 (MQ=255)
      aCCGGGGTGGCGTTGGTCGATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTC‑TCCGCCTGTAGCGG‑CCCTTTCGCGGGCAt                                                                                                                           <  2:7568147/137‑1 (MQ=255)
         ggggTGGCGTTGGTCGATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTC‑TCCGCCTGTAGCGGACCCTTTCTCgggatggc                                                                                                                         >  2:10011824/1‑132 (MQ=255)
         ggggTGGCGTTGGTCGATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTC‑TCCGCCTGTAGCGG‑GCCTTTTTCGGGGAGGcc                                                                                                                        >  2:1722709/1‑137 (MQ=255)
         ggggTGGCGTTGGTCGATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTC‑TCCGCCTGTAGCGG‑ACCTTTCGCgggtacttc                                                                                                                        >  2:2997654/1‑131 (MQ=255)
           ggTGGCGTTGGTCGATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGAGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTC‑TCCGCCTGTAGCGG‑GCCCTTTTTGGGGATGcccc                                                                                                                      >  1:4008281/1‑136 (MQ=255)
            gTGGCGTTGGTCGATGAGGGGTTTGGGTTTTGCTGGCAGTTGTTGGCGCGGTATGCGTGGGGTGTGTGTGGGCAGGGGATCGTGCTGGTCTGACCAACAGCC‑TCCGCCTGTAGTTG‑CCTTTTCGCGGGCGGGCCcgc                                                                                                                     <  2:13095556/137‑1 (MQ=255)
                   tGGTCGATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTC‑TCCGCCTGTAGCGG‑CCCCTTTTTGGGTATGTCCCCGAAAagg                                                                                                              >  2:2537404/1‑137 (MQ=255)
                   tGGTCGATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTC‑TCCGCCTGTAGCGG‑CCCCTTTCCGGGGATGCCCCCGAAAagg                                                                                                              >  2:13917272/1‑137 (MQ=255)
                     gTCGATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTC‑TCCGCCTGTAGCGG‑CCCTTTCGCGGGGGCGCCCGCAAAAagggc                                                                                                            >  1:12412556/1‑135 (MQ=255)
                     gTCGATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTC‑TCCGCCTGTAGCGG‑ACCTTTTTCGGGGATGCCCCCGaaaaagaa                                                                                                            >  1:9310653/1‑132 (MQ=255)
                     gTCGATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTC‑TCCGCCTGTAGCGG‑ACCCTTCTTGGGGATGACCACaaaaaagac                                                                                                            >  1:4836203/1‑132 (MQ=255)
                       cGATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTC‑TCCGCCTGTAGCGG‑CCCTTTCTCGGGGATGCCCACGAAAGGG‑ACgg                                                                                                         >  2:19050927/1‑137 (MQ=255)
                       cGATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGCCTGACCAATAGTC‑TCCGCCTGTAGCGG‑CCCCTTCACTGGGATGGCCACGAAAAGG‑ACgg                                                                                                         >  2:6248626/1‑137 (MQ=255)
                        gATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTC‑TCCGCCTGTAGCGG‑CCCTTTCTCGGGGATGTCCGCGAAAGGG‑GCggt                                                                                                        >  2:14447246/1‑137 (MQ=255)
                        gATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTC‑TCCGCCTGTAGCGG‑CCCTTTCGCGGGGCTGCCCACGAAAGGG‑ACggt                                                                                                        >  2:805676/1‑137 (MQ=255)
                        gATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTC‑TCCGCCTGTAGCGG‑CCCTTTCGCGGGCATGCCCCCGAAAGGG‑ACggg                                                                                                        >  2:8852993/1‑136 (MQ=255)
                        gATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTC‑TCCGCCTGTAGCGG‑CCCCTTTTCGGGGATTGCCCCGAAAagggcagt                                                                                                         >  2:18167802/1‑132 (MQ=255)
                        gATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTC‑TCCGCCTGTAGCGG‑ACCTTTTTCTGGTATGCCCCCGAAAaggaagag                                                                                                         >  1:8199706/1‑132 (MQ=255)
                        gATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTC‑TCCGCCTGTAGCGG‑ACCCTTTGCTGGGAGGCCCTCGAAAagggagag                                                                                                         >  1:17069579/1‑132 (MQ=255)
                          tgggtgggtGGGGGTTGGCCGGAGGGTGTGGGCGTGGGAGAGGTGGGAGGGGGGTGGGAAGTGGATCGTTTTGTTCGGACCAAAGGTC‑TCCGCCGGTAGTGT‑CCCTTTCGTGGGCATGCCCGCGAAAGGG‑CCggtac                                                                                                      <  2:17616390/131‑1 (MQ=255)
                          tGAGGTGTTTGGGGTTTGCCGGCATGTTTTGTCGCGGTGTGCGTGGGATGTGGGTGGGAGGTGGGTTGTGTTGTTTTGACAAACGGTC‑TCTGTTTGTAGCGG‑CCCTTTTGTGGGCGTACCCGCGAAAGGG‑CCggtac                                                                                                      <  2:18378084/137‑1 (MQ=255)
                          tGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGTCGCGGTATGCGCGGGATGTGTGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTC‑TCCGCCTGTAGCGG‑CCCTTTCGCGGGCATGCCCACGAAAGGG‑CCggtac                                                                                                      <  2:13388328/137‑1 (MQ=255)
                          tGAGGCGGTTGTGTTTTTTTGGCGGGTGTGGGTGTGGTGGGCTCGGGGTTTGTGTGGGCGGCGGAGCGTGCTGGTCTGACCAA‑AGTCTTCTGCTGGGAGCGG‑CCCTTTCGCGGGCGTGCCCGCGAAAGGG‑CCggtac                                                                                                      <  2:4445291/137‑1 (MQ=12)
                            aGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTC‑TCCGCCTGTAGCGG‑TCCTTTCGCGGGCATGCCCACAAAAGGG‑CCagtacca                                                                                                    >  2:9659793/1‑135 (MQ=255)
                            aGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTC‑TCCGCCTGTAGCGG‑CCCTTTTCGGGGGACGCCCGCGAAAAGG‑ACgggcaaa                                                                                                    >  2:5927427/1‑132 (MQ=255)
                            aGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTC‑TCCGCCTGTAGCGG‑CCCTTTCGAGGGTATGCCCGCGAAAGGG‑CCggtacag                                                                                                    >  2:9340402/1‑137 (MQ=255)
                            aGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTC‑TCCGCCTGTAGCGG‑CCCCTTCACGGGGATGCCCGCAAAAAGG‑CCagtacaa                                                                                                    >  2:4007422/1‑136 (MQ=255)
                               tggttggtTTTTCCCGGCAGGTTTTGGCGCGGGCTGCGCGGGGTGTGGGTGGGCAGCGGATCGTGCTGTTCTGACAAAAATTC‑CCCACCTGGAGCGT‑CCCTTTCGCGGGAATGCCCACGAAAGGG‑CCggtacaggta                                                                                                 <  2:11664919/137‑1 (MQ=255)
                                  ttGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTC‑TCCGCCTGTAGCGG‑CCCTTTCGCGGGCACGCCCCCAGAAAGG‑ACggcacagggaaaa                                                                                              >  2:15508417/1‑132 (MQ=255)
                                  ttGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTC‑TCCGCCTGTAGCGG‑ACCTTTTTCTGGGACGCCCCCAAAAAGG‑ACGgtgcaagtgcag                                                                                              >  1:10530884/1‑137 (MQ=255)
                                       gTTTGTCGGCGGGTGTGGGCGTGGGAGGTGCGGGATGGGGGGGGGTAGCGGGTCGTGTTGGTCTGCCCAACATTC‑TCCGCCTGTATTGG‑CCCTTTCCCGGGCATCCCCGCGAAAGGG‑CCGGTACAGGTACAG‑CAGaa                                                                                        <  2:13203607/137‑1 (MQ=255)
                                       gTTTGCTGGCGGGTGTTGGCTTGGTATGCACGAGCTGTGTTTGGACAGCGGATTGTTCTGTTTTGACCAACATTC‑TCTGCCTGTGGCGC‑CCCTTTCGCGGGCATTCCCGCGAAAGGG‑CCGGTACAGGTACAG‑CAGaa                                                                                        <  2:2154982/137‑1 (MQ=255)
                                       gTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTC‑TCCGCCTGTAGCGG‑CCCCTTCTCGGGCACTCCCCCCAAAAGG‑ACAGTACAAGTACCA‑CAGaa                                                                                        >  2:9829145/1‑137 (MQ=255)
                                       gTTTGCCGGCAGGTGTTGGCATTGTATGCGCTTGGTTTGTGTGGTTATCGGATCGTGTTGGTTCTACCAACAGTC‑TCCGACGGGAGCGG‑CACTTTCTTAGGCATGCTCGCGAAAGGG‑CCGTTACAGGTACAG‑CATaa                                                                                        <  2:14212503/137‑1 (MQ=255)
                                                cAGTTGTTTGTGCGGCGTGCGTGGGATGTGGGTGGGCAGTGGATCGTGCTGTTCTGACAAATGGTT‑TCCGCCTGTAGGGG‑CCTTTTTGTGGGCATCCCTGCGAAGGGG‑CCGGTACAGGTACAG‑CAGAACTTAGCGcc                                                                               <  2:12067795/137‑1 (MQ=255)
                                                aggtttttGGCGGGGGGTGCGGGGGGTGTGGGGGGAAAGAGGGTTGTGCTGGTCGAACAAAAATTC‑TCCGCTTGGAGCGT‑CCTTTTCGAGGGCATACCAAAGAAAGGG‑CCGGTACAGGTACAG‑CAGAACTTAGCGcc                                                                               <  2:5709127/131‑1 (MQ=255)
                                                  ggTGTTGGCGCGGGATGCGCGGGATGTGGGTGGGCGGTGGATCGTGCTGGTCTGACCAACAGTC‑TCCGCCTGTAGCGG‑CCTTTTTGCGGGCATACCCGCGAAAGGG‑CCGGTACAGGTACAG‑CAGAACTTAGCGCCAg                                                                             <  2:4646375/137‑1 (MQ=255)
                                                  ggTGTGGGCGTGGGATGTGTGGGGGGTGGGTGGGCAGGGGGTTGTGCTGGTTGGAAAACCATTC‑TCCTCTTGGGGCGG‑CCCTTTTGCGGGCATGCCCGGGAAGGGG‑CCGGTACAGGTACAG‑CAGAACTTAGCGCCAg                                                                             <  2:3927546/137‑1 (MQ=255)
                                                    ttttGTCGCGGTGTGCTTGGGATGTGGTTGGGCAGCGGATCGTGCTGGTTGGACCAACATTC‑TCCGCTTGTAGTGT‑CCTTTTCGTGGGCATCCCCGCGAAGGGG‑CCGGTACAGGTACAG‑CAGAACTTAGCGCCAGct                                                                           <  2:19475700/135‑1 (MQ=255)
                                                         cctcgGTCTGCTTGGGCTTTGTGTGGGCAGCGAATCGTGCTGTTCGGACCAACATTC‑CCCGCTTGTATCGT‑CCCTTTCTCGGACATCCCCACAAAAGGG‑CCGGTACAGGTACAG‑CAGAACTTAGCGCCAGCTCTCtt                                                                      <  2:10539562/134‑1 (MQ=255)
                                                          ggcgGGATGCGTGGGGGGTGTGGGGGCAGCGACGCGGGTTGGCCGGACCAACAGTC‑TCCGCCTGTAGCTG‑CCCGTTCGTGGGCATCCCCGAGAAAGGG‑CCGGTACAGGTACAG‑CAGAACTTAGCGCCAGCTCTCttt                                                                     <  2:12177791/136‑1 (MQ=255)
                                                             ggTGTGCGTGGGGTGGGGTTGGGAAGTGGATTGTGTTGTTCGGACCAACAGTC‑TCCGCCTGTAGTGG‑CCCTTTTGAGGGCATACCCGCGAAAGGG‑CCGGTACAGGTACAG‑CAGAACTTAGCGCCAGCTCTCTTTcac                                                                  <  2:14935077/137‑1 (MQ=255)
                                                                       ggATGTGGTGGGGCATGGGATTGTGTTGTTCGAACCAACGTTC‑CCCTCTTGGAGTGC‑CCTTTTCGTGGGCGTTCCCGAAAAAGGG‑CCGGTACAGGTACAG‑CAGAACTTAGCGCCAGCTCTCTTTCACCACCCTgcgc                                                        <  2:8377595/137‑1 (MQ=255)
                                                                       gaaTGGGGGTGGGCAGCGGATCGGGTTGGTTTGACCAACGTTT‑TCCGCTTGGAGCGG‑CCCTTTCGTGGGCAGGCCCGCGAAAGGG‑CCGGTACAGGTACAG‑CAGAACTTAGCGCCAGCTCTCTTTCACCACCCTgcgc                                                        <  2:18588517/135‑1 (MQ=255)
                                                                        gATGTGGGTGGGCGGTGGATTGTGTTGGTTTGACCAACATTC‑TCCGCTTGGAGCGT‑CCCTTTCGCGGGCATACCCGCGAAAGGG‑CCGGTACAGGTACAG‑CAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGcgcc                                                       <  2:16800190/137‑1 (MQ=255)
                                                                        gATGTGGGTGGACAGCGGATTGTGCTGGTCTGACAAACAGTC‑TCCGCCTGTAGCGG‑CCCTTTCGCGGGCATGCCCGCGAAAGGG‑CCGGTACAGGTACAG‑CAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGcgcc                                                       <  2:3798914/137‑1 (MQ=255)
                                                                            tgtgtgggCGGTGGGTCGTGCTGTTCTGACCAACATTC‑TCCGCTTGTAGCGC‑CCTTTTCGCGGGCATGCCCGCGAAGGGG‑CCGGTACAGGTACAG‑CAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGCGccgccc                                                   <  2:8885488/134‑1 (MQ=255)
                                                                            tgggtgggCAGTGGCTTGTGTTGGTCTGACCAACGGTC‑TCTGCCTGTAGTGT‑CCTTTTTGTGGGCATTCCCGCGAAAGGG‑CCGGTACAGGTACAG‑CAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGCGccgccc                                                   <  1:8413297/137‑1 (MQ=255)
                                                                            tgggtgggCAGCGGATCGTGCTGGTCTGACCAACAGTC‑TCCGCCTGTAGCGG‑CCCCTTCCTGGGTATGCCCCCGAAAAGG‑ACCGTAAAGGTACCGCCAAAACTT‑GCGCCAGACCTCTTTCACCAACCCTCGcccccc                                                   >  2:3940055/1‑137 (MQ=255)
                                                                            tgggtggaCAGCGGATCGTGCTGTTCTGACAAACATTC‑TCCTCTTGGAGCGT‑CCTTTTCGCGGGCATGCCCGCGAAAGGG‑CCGGTACAGGTACAG‑CAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGCGccgccc                                                   <  2:17992639/137‑1 (MQ=255)
                                                                            gggttgggAAGGGGGTTGTGTTGTTTTGACCAACGGTC‑TCCTCTTGTAGCGT‑CCCTTTCTCGGGCATACCTGCGAAAGGG‑CCGGTACAGGTACAG‑CAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGCGccgccc                                                   <  2:10794436/133‑1 (MQ=255)
                                                                                     aGCGGATCGTGCTGGTCTGACCAACAGTC‑TCCGCCTGTAGCGG‑CCCTTTCTCGGGGACGCCCCCAAAAGGG‑GCGGTACAAGTACAA‑CAAAACCTCGCGCCCACGCTCCTTCCCCCCCCCGcccccccccccccagac                                          >  2:14155242/1‑132 (MQ=255)
                                                                                     aGCGGATCGGGTGGGTCCGACCAACAGTC‑TCCGCCTGGAGTGG‑CCCTTTCGCGGGCATGCCAGCGAAAGGG‑CCGGTACAGGTACAG‑CAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGCGCCGCCCGCCCAGGAt                                          <  1:9241078/137‑1 (MQ=38)
                                                                                      gCGGATCGTGCTGGTCTGACCAACAGTC‑TCCGCCTGTAGCGG‑ACCCTTCGCGGGAGCGTACGAGAAAGGG‑ACGGGACAAGGACAG‑CAGAAACTCGCGCCCAATCTCTTTCCCCACACCGCGCACcccacccaaaagc                                         >  2:9257503/1‑132 (MQ=255)
                                                                                          gttgtgtTGGTCTGACCAACGGCC‑TCCTCTTGGAGTGT‑CCTTTTCGCGGGCATTCCCGCGAAGGGG‑CCGGTACAGGTACAG‑CAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGCGCCGCCCGCCCAGGATcagcc                                     <  2:2794036/134‑1 (MQ=255)
                                                                                          gcggtgtTGGTTTGCCCAACGTTC‑CCCGCTTTTGGCGG‑CCTTTTCTCGGGCATTCCTGTGAAGGGG‑CCGGTACAGGTACAG‑CAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGCGCCGCCCGCCCAGGATcagcc                                     <  2:2837464/134‑1 (MQ=255)
                                                                                          aTCGTGCTGGTCTGACCAACAGTC‑TCCGCCTGTAGCGG‑CCCTTTCGCGGGGATGCCCGCGAAAGGG‑AAGGGGCAGGTACAG‑CAGAAATTAGTGACAGCTCTCTTTTACCACCCTGCGCCGCACTCCCAGGATcagcc                                     >  2:14517222/1‑137 (MQ=255)
                                                                                            tgtgtTGGTTTGACAAACGTTC‑TCTGCTGGGAGGGT‑CCTTTTCGCGGGAATGCCTACGAAGGGG‑CCGGTACAGGTACAG‑CAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGCGCCGCCCGCCCAGGATcagccag                                   <  2:1934402/136‑1 (MQ=255)
                                                                                            tgTGCTGGTCGGACAAACAGTC‑TCCGCCTGGAAAGT‑CCTTTTCTCGGACATGCCCGCAAAAGGG‑CCGGTACAGGTACAG‑CAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGCGCCACCCGCCCAGGATcagccag                                   <  2:16947675/136‑1 (MQ=255)
                                                                                             gggCTGGCCTGACAAACAGCC‑TCCGCTTTTATCGG‑ACCTGTCATGGGCGGGCCCGAGAAAGGG‑CCGGTACAGGTACAG‑CAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGCGCCGCCCGCTCAGGATcagccagc                                  <  2:2652977/135‑1 (MQ=255)
                                                                                                 tgtccTGACCAACATCC‑TCCTCCTGGAGCGC‑CCCTTCCATGGACATCCCCGGGAAAGGG‑CCGGTACAGGTACAG‑CAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGCGCCGCCCGCCCAGGATCAGCCAGCCAAt                              <  2:18175782/133‑1 (MQ=255)
                                                                                                 tgtccTGACCAACAGTC‑TCCGCTGGTAGTGG‑CCCTTTCGCGGGCATTCCCGCGAAAGGG‑CCGGTACAGGTACAG‑CAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGCGCCGCCCGCCCAGGATCAGCCAGCCAAt                              <  2:18077620/133‑1 (MQ=255)
                                                                                                   ttttGACCAACAGTC‑TCCCCTTGTGATTT‑CCCTTTTTTGGGCATTCCCCCGAAGGGGCCCGGTACAGGTACAG‑CAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGCGCCGCCCGCCCAGGATCAGCCAGCCAATc                             <  1:19097558/134‑1 (MQ=255)
                                                                                                      tGACCAACAGTC‑TCCGCCTGTAGCGG‑CCCCTTTCCAGGTATGCCCCCGAAAAGG‑CCAAGCCAAGTACAG‑GAAAACTTCGCACCCAAGCTCTTTCCCCACCACGCGCCACCCACCCAAAAGGAGACAGACAAACCCCa                         >  2:1903682/1‑137 (MQ=255)
                                                                                                        acaaaCAGTC‑TCCGCTTGTGGCGG‑CCCTTTCGCGGGCATGCCCGCGAAAGGG‑CCGGTACAGGTACAG‑CAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGCGCCGCCCGCCCAGGATCAGCCAGCCAATCCCCAGc                       <  1:3957815/134‑1 (MQ=255)
                                                                                                        aCCAACATTC‑TCCGCCTGTAGCGG‑CCCTTGCGCGGGCATGCCTGCGAAAGGG‑CCGGTACAGGTACAG‑CAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGCGCCGCCCGCCCAGGATCAGCCAGCCAATCCCCAGc                       <  2:3654137/137‑1 (MQ=255)
                                                                                                        aCCAACAGTC‑TCCGCTTGTAGCGT‑CCCTTTAGCGGGTATACCCGCGAAAGGG‑CCGGTACAGGTACAG‑CAGAACTTAGCGCCAGCTTTCTTTCACCACCCTGCGCCGCCCGCCCAGGATCAGCCAGCCAATCCCCAGc                       <  2:18849662/137‑1 (MQ=255)
                                                                                                        aCCAACAGTC‑TCCGCCTGGAGTGG‑CCCTTCCGCGGGCATGCCCGCGAAAGGG‑CCGGTACAGGTACAG‑CAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGCGCCGCCCGCCCAGGATCAGCCAGCCAATCCCCAGc                       <  2:10012752/137‑1 (MQ=255)
                                                                                                                tc‑tcCGTTTGTAGCGG‑CCTTTTCGCGGGCATGCCGGCGAAAGGG‑CCGGTACAGGTACAG‑CAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGCGCCGCCCGCCCAGGATCAGCCAGCCAATCCCCAGCAACAGGCt               <  2:13530780/137‑1 (MQ=255)
                                                                                                                tc‑tcCGCCTGTAGCGG‑CCCTTTCGTGGGCATGTCCGTTAAAGGG‑GCGGTGCAGGTGCCG‑CAGAACCTCGCGCAAGGTGTCTTTCCCCACCCCTGGCCGCCCGCACAGGATTAGACGGGCAATCTACAACAAAAgggt               >  2:3689494/1‑135 (MQ=255)
                                                                                                                   tcCGCTTGTAGCGT‑CCCTTTCGTGGGCATGCCCGCGAAAGGG‑CCGGTACAGGTACAG‑CAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGCGCCGCCCGCCCAGGATCAGCCAGCCAATCCCCAGCAACAGGCttt             <  2:15605667/137‑1 (MQ=255)
                                                                                                                   tcCGCCTGTAGCGT‑CCTTTTCTTGGCCATCCCCGCAAAAGGG‑CCGGTACAGGTACAG‑CAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGCGCCGCCCGCCCAGGATCAGCCAGCCAATCCCCAGCAACAGGCttt             <  1:9007145/137‑1 (MQ=255)
                                                                                                                   tcCGCCTGTAGCGC‑CCCTTTCGCGGGCATACCCACGAAAGGG‑CCGGTACAGGTACAG‑CAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGCGCCGCCCGCCCAGGATCAGCCAGCCAATCCCCAGCAACAGGCttt             <  2:10784858/137‑1 (MQ=255)
                                                                                                                   cccGCCTTGATCGG‑CCCTTTCACGGGCACTCCCGCGAAAGGT‑CCGGTACAGGTACAG‑CAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGCGCCGCCCGCCCAGGATCAGCCAGCCAATCCCCAGCAACAGGCttt             <  1:17191424/136‑1 (MQ=255)
                                                                                                                      gcttGTAGTGT‑CCCTTTCGTGGGCGTGCCCGCGAAAGGG‑CCGGTACAGGTACAG‑CAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGCGCCGCCCGCCCAGGATCAGCCAGCCAATCCCCAGCAACAGGCTTTCGa          <  2:3301786/134‑1 (MQ=255)
                                                                                                                        ttgcggCGT‑CCCTTTCGCGGGCATGCCCGCGAAAGGG‑CCGGTACAGGTACAG‑CAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGCGCCGCCCGCCCAGGATCAGCCAGCCAATCCCCAGCAACAGGCTTTCGAcc        <  2:4251737/132‑1 (MQ=255)
                                                                                                                        cTGTAGTGT‑CCTTTTCGCGGGCATGCCCGCGAAAGGG‑CCGGTACAGGTACAG‑CAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGCGCCGCCCGCCCAGGATCAGCCAGCCAATCCCCAGCAACAGGCTTTCGAcc        <  2:16385845/137‑1 (MQ=255)
                                                                                                                        cTGGAGCGG‑CCCTTTCGCGGGCATACCCTCGAAAGGG‑CCGGTACAGGTACAG‑CAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGCGCCGCCCGCCCAGGATCAGCCAGCCAATCCCCAGCAACAGGCTTTCGAcc        <  2:19141237/137‑1 (MQ=255)
                                                                                                                             gCGT‑CCCTTTCGCGGTCATGCCCTCGAAAGGG‑CCGGTACAGGTACAG‑CAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGCGCCGCCCGCCCAGGATCAGCCAGCCAATCCCCAGCAACAGGCTTTCGACCACaaa   <  2:9462757/137‑1 (MQ=255)
                                                                                                                             gCGG‑CCTTTTCTCGGGCATGCCCGTGAAGGGG‑CCGGTACAGGTACAG‑CAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGCGCCGCCCGCCCAGGATCAGCCAGCCAATCCCCAGCAACAGGCTTTCGACCACaaa   <  2:18660928/137‑1 (MQ=255)
                                                                                                                              cGG‑CCTTTTTGGGGACATTCCCGCGAAAGGG‑CCGGTACAGGTACAG‑CAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGCGCCGCCCGCCCAGGATCAGCCAGCCAATCCCCAGCAACAGGCTTTCGACCACAAAg  <  1:18915182/137‑1 (MQ=255)
                                                                                                                               tgtccttttCGCGGGCAGGCCCGCGAAGGGG‑CCGGTACAGGTACAG‑CAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGCGCCGCCCGCCCAGGATCAGCCAGCCAATCCCCAGCAACAGGCTTTCGACCACAAAg  <  1:18088170/131‑1 (MQ=255)
                                                                                                                               tgccccTTTCGGGGGCGTGCCCGCGAAAGGG‑CCGGTACAGGTACAG‑CAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGCGCCGCCCGCCCAGGATCAGCCAGCCAATCCCCAGCAACAGGCTTTCGACCACAAAg  <  1:18088211/134‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                
TGGCAGACCGGGGTGGCGTTGGTCGATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTC‑TCCGCCTGTAGCGG‑CCCTTTCGCGGGCATGCCCGCGAAAGGG‑CCGGTACAGGTACAG‑CAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGCGCCGCCCGCCCAGGATCAGCCAGCCAATCCCCAGCAACAGGCTTTCGACCACAAAG  >  AE015451_phaM1/5199521‑5199782

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