Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451_phaM15,485,1530AG31.7% 136.5 / 14.7 88intergenic (+59/‑88)purH CDS (phosphoribosylaminoimidazolecarboxamide formyltransferase_IMP cyclohydrolase)/purD CDS (Phosphoribosylamine‑‑glycine ligase)–/–
Reads supporting (aligned to +/- strand):  ref base A (20/36);  new base G (1/25);  total (21/67)
Fisher's exact test for biased strand distribution p-value = 2.13e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTG  >  AE015451_phaM1/5485018‑5485285
                                                                                                                                       |                                                                                                                                     
aTGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGAt                                                                                                                                      <  2:17482131/137‑1 (MQ=255)
aTGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGAt                                                                                                                                      <  2:14758676/137‑1 (MQ=255)
aTGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGAt                                                                                                                                      <  2:9126474/137‑1 (MQ=255)
aTGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGAt                                                                                                                                      <  2:18895951/137‑1 (MQ=255)
aTGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGAt                                                                                                                                      <  2:9126439/137‑1 (MQ=255)
aTGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGAt                                                                                                                                      <  2:17352466/137‑1 (MQ=255)
        gTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGATGCAAGTCTAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGAtt                                                                                                                              <  1:13777929/137‑1 (MQ=255)
          cATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCa                                                                                                                            <  2:3128974/137‑1 (MQ=255)
             cGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCAccg                                                                                                                         <  2:17697209/137‑1 (MQ=255)
               ctgctgCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCAccgcc                                                                                                                       <  1:11512147/137‑1 (MQ=255)
               ctgctgCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCAccgcc                                                                                                                       <  1:1596009/137‑1 (MQ=255)
               ctgctgCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCAccgcc                                                                                                                       <  1:14507988/137‑1 (MQ=255)
                  ctgcCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCa                                                                                                                    <  1:3109143/137‑1 (MQ=255)
                  ctgcCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCa                                                                                                                    <  1:15996299/137‑1 (MQ=255)
                  ctgcCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCa                                                                                                                    <  1:15996985/137‑1 (MQ=255)
                  ctgcCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCa                                                                                                                    <  1:3153030/137‑1 (MQ=255)
                  ctgcCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCa                                                                                                                    <  2:2195297/137‑1 (MQ=255)
                  ctgcCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCa                                                                                                                    <  1:14373235/137‑1 (MQ=255)
                  ctgcCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCa                                                                                                                    <  1:13898669/137‑1 (MQ=255)
                   tgcCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCAt                                                                                                                   <  2:10291264/137‑1 (MQ=255)
                   tgcCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCAt                                                                                                                   <  1:16426481/137‑1 (MQ=255)
                   tgcCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCAt                                                                                                                   <  1:12903948/137‑1 (MQ=255)
                    gcCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCAtt                                                                                                                  <  1:11594261/137‑1 (MQ=255)
                             gCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCaa                                                                                                        <  2:6157119/137‑1 (MQ=255)
                                     cgcgATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCTCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTc                                                                                                <  1:3151076/137‑1 (MQ=255)
                                               ttCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGGTTGCCCTCGAGTACAGCGGTTGGTTCCCCTTCTCATTCTCTGGCAGGCCCTCTCCCACCATGtt                                                                                      >  2:8953684/1‑137 (MQ=255)
                                                 cgcTGGCGTGCTCTACTTCCGCCACGAGTTCTTCGGATCCAGAGAAAGTTAGAGCCTTTGTGGGAGCGGGTTTACCCGAGAAGACGAGGTGGGGTTCACCACTGCATTCGTGGGCAATCCCGCTCCCACAGTGTTCg                                                                                    <  2:13734276/135‑1 (MQ=255)
                                                       cgtGCTCCCCTCCCCCCACTACTTACGCGAATCCTGAGTCAGGCCGAGCCCTTGGGGGAGCGGTCTTGCCCGAGAATCCGTTGGTGGGTTCCCCTCGGCCTTCGTGGGAAACCCCACTCCCACAGTGTTCGGCGCaa                                                                              <  2:10742311/135‑1 (MQ=255)
                                                       cgtGCACCACTCCCCCCCCTAATCACGCGTATCCCGAAACAATTCGAGATTTTGTGGGGGCGGCCTTCCCCCCGAACACGCTGGTAGCTCCACCTCCACGTCCGTGGCCACTCCCGCTCCCACAGTGTTCGGCGCaa                                                                              <  2:47619/135‑1 (MQ=255)
                                                         tGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGGTTGTCCACTAAGACAACGGGGGGGTCCCCCCCCCCCCCCCTGGCAAGCCCCCCCCCCCAACGTTCTGCGCaaac                                                                            >  2:13287623/1‑135 (MQ=255)
                                                         tGCGCAACTCCCGCCACTAATCCCGCGAATCACGAGACAAGTCGAGATCCTGTGGGAGCGGGCTTGCACGCGACGACGATGGTGGATTCACCGCCGCATTCGCAGGCAAGAACGCTCCCACAGTGTTCGGCGCAAGc                                                                            >  2:8991168/1‑137 (MQ=255)
                                                          gcgcCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGGTTGTCCCCGAAAACGACGGTGGAGTCCCCCCCCAACTCGCGGGGAAGCACCCCCCCACCATGGTTGGCGCaaaac                                                                           >  2:4632702/1‑134 (MQ=255)
                                                          gcgcCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGATTGCCCGCAAACAAGAGAGTGGATTCACTACCACACCCGCGAGCAAACCCGCCCCCCCAACGTTCTGCGCAAAcc                                                                           >  1:15774062/1‑137 (MQ=255)
                                                            tccATTTCCACCATAAGTTACTCGGGTCCCGGGGCAGGTTGCGATCCTGTGGAATTGGCCTTGCCTGCAACTACGCTTGTGGCTCCACCACCCCCTTTGCGGCAAACCCCGCTCCCACAGTGTTCGGCGCAAGCCAc                                                                         <  2:4708776/136‑1 (MQ=255)
                                                            tccACTTCCGCCACTAGTTACGCGGATCCCGAGGAAGGTCGAGATCCTGTGGGAGCGGGCTTGCCTGTGAGGACGGTGGGGGATTCACCTCCGCATTCGCGGACAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCAc                                                                         <  2:14103082/136‑1 (MQ=37)
                                                            gcCACTTCCTCCACTACTTACTCGGATCCCGGGGCAAGTCGAGGTCTTGTGGGAGCGGGTTTGCCCGAGAAGACGGTGGTGGATTCACCGCCGCATTCGCGGGAAATCCCGCTCCCACAGTGTTCGGCGCAAGCCAc                                                                         <  2:3186506/137‑1 (MQ=255)
                                                            gcCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGGTTTTCCACGAACACAACGGTGGAGTCACCACCACCCTCGCGGGACAAACCCCCCCCCCAACGTTCGGCGCAAACCAc                                                                         >  2:13446282/1‑137 (MQ=255)
                                                            gcCACTTCCACCACTAATTACGCAGACCCTGGGACAAGTCGAGCTCTGGGGGGGGCGGGCTTCCCCGCGAATATGGTGGTGGCTCCACTACCCCGTTTGGGGGAAATCCCGCTCCCACAGTGTTCGGCGCAAGCCAc                                                                         <  2:5168060/137‑1 (MQ=255)
                                                            actaCTTCCGCAACTGATTACGCGGCCCTCGGGGAAGGTAAGGCCCCGGGGGGAGCGGGCTGGCCCGAAAAAACGGGGGGGGAGCCACCGCAGCATTGGGGGAAAATCCCGCTCCCACAGTGTTCGGCGCAAGCCAc                                                                         <  2:3186538/134‑1 (MQ=255)
                                                            accACTTCCCCCGTGAATTACGCGGATCCTGAGGAGATTCGAGGTCCGGTGGGGGCGGTTTTATCCGCAAATACGCTGGTGGATTCACCCCCGCGTTCGTGGGAAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCAc                                                                         <  2:12326120/136‑1 (MQ=255)
                                                              cACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGGTTGCCCGCGAAGACAACGGTGGATTCACCACCACATTTTCGGGCAAGCCCGCTCCCACAACGTTCGGCGCAAACCAcag                                                                       >  1:774036/1‑135 (MQ=37)
                                                              cACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTTTCGGGCAAACCCACCCCCACAATGTTCGGCGCAAGCCACgg                                                                       >  1:774085/1‑137 (MQ=39)
                                                                 ttCCGCCACTAAGTACGTGGATCCAGAGACAGGTAGAGATCCTGGGGGAGCGGGTTTGCCCGAGAATACGAGGGTGGATTCACCGCCGCATTCGCGGACAACCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGAt                                                                    <  2:9173216/137‑1 (MQ=255)
                                                                         cTAATTACGTGCATCCTGAGAAGAGTTGAGGTCTTGGGGGGGCGTGCCTGCCAAGGAATCCGGGGGTGGGTTCACCTCTGCGTTCGTGGGAAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGaa                                                            <  2:13272808/137‑1 (MQ=255)
                                                                                cgcgGCTCCTGAGACAGGTAGAGCTCTTGTGGGGGCGGTCTTCCCCGAAAAGACGGTGGTGGCTCCACCGCCCCGTTCGCGGACAACCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGt                                                     <  2:12307303/137‑1 (MQ=255)
                                                                                       ccagAGGAAATTCGAGGCCCTGGGGAAGCGGTCTTCCCTACGAGGAAGGGGGGGGGTTCACGCCCCCATCCATGGGAAATCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACa                                              <  2:17004524/134‑1 (MQ=255)
                                                                                             gCAAGTGGAGCTTTTGTGGGAGTGGGCTTGCCCGAGAAGATGGTGGGGGTTCCCCCGCCCCGTTCGCGGACAATCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAg                                        <  2:10070027/137‑1 (MQ=255)
                                                                                                  tCGAGCTCCTGTGGAGGCGGTCTTACCCGAAAGTACGATGGTGGATTCACCTCCACATTCGTGGGCAATCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTg                                   <  2:15652910/137‑1 (MQ=255)
                                                                                                      gaTCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCACGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGatca                               <  1:2023408/137‑1 (MQ=255)
                                                                                                       aTCTTGGGGGGGCTGTCTTCCCCGCGAATATGGTGGGGGATTCACCGCCGCGTTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGatcat                              <  2:13634811/137‑1 (MQ=255)
                                                                                                          ttttgGGGGTGGGTTTGCCCACGAATACGGTGTTGGCTTCACCCCCGCGTTCGCGGCAAACCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCgg                           <  2:16601889/134‑1 (MQ=255)
                                                                                                            gtgGGAGCGTGCGTCCCCTAAAGTCCGGTGGTGGGTTCCCCCCCACATTCGTAGACACCCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCa                         <  1:15759250/137‑1 (MQ=255)
                                                                                                            gggggAGCGGTCTTGCCTGCGAATACGGTGGTGGATCCACCTCCTCATCTGCGGGCAACCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCa                         <  1:15759292/135‑1 (MQ=255)
                                                                                                              gggAGCGGTCTTCCCCGCGAGTCCGCTGGTGGCTCCACCCCCGCCTTCGTGGACAAACCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGc                       <  2:17274113/137‑1 (MQ=255)
                                                                                                                  gCGGTCTTGCCCGGGAGTATGGGGGTGGCTCCACCTCCGCGTTCCCGGGAACCCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTg                   <  2:18536800/137‑1 (MQ=255)
                                                                                                                  gCGGGTTTGCCTGAGAGTACGATGGTGGGTTCGCCTCCTCGTTTGTGGGCAACCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTg                   <  2:8954928/137‑1 (MQ=255)
                                                                                                                  gCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTg                   >  1:17060273/1‑137 (MQ=255)
                                                                                                                   cGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGACAAGCCCGCTCCCACAGTGGTCGGCGCAAGTCATGGGATTACGTGAATTGAGGTGTTGACATGGAAGTTTTGATTATCGGCAGCGGTgg                  <  2:9128287/137‑1 (MQ=255)
                                                                                                                       tttGCCCGCGAAGACGGTGTTGACTTCACCGCCCCCTCCGCGGCAAACCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGt              <  2:2148408/136‑1 (MQ=255)
                                                                                                                       tttGCCCGAGAATACGGGGGGGGCTCCACCACCGCATTCGTGGACACCCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGt              <  2:11446587/136‑1 (MQ=255)
                                                                                                                       cTTGCCCGCGAGTATGGTGGTGGATTAACCGCCCCATTTGCGGGCAACCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGt              <  2:18544166/137‑1 (MQ=255)
                                                                                                                       cTTCCCCGCGAAGACGATGTGGGCCTCACCGCCGCCTTCGCGGACAATCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGt              <  2:1744680/137‑1 (MQ=255)
                                                                                                                        ttGCCCGAGAAGATGATGTTGGTTTCCCCCCCACCTTCGTGGACAACCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTg             <  2:9094859/137‑1 (MQ=255)
                                                                                                                          ggccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAATTTTGATCATCGGCAGCGGTGGCCGTGAg           >  1:17211076/3‑137 (MQ=255)
                                                                                                                             tgcgAGTGCGGTGTTGGTTTCACTGCCTTGTTTGCGGGAAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCAc        <  1:2239511/136‑1 (MQ=255)
                                                                                                                             cgcgGATACGGTGGTGGATCCCCCCCCGCCTTCGCGGGCAACCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCAc        <  2:8188185/137‑1 (MQ=255)
                                                                                                                             cgcgAGTATGGTGGTGGATTCCCCTCCTCTTTCTGGGGCAATCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCAc        <  2:824880/137‑1 (MQ=255)
                                                                                                                             cgcgAGTAGGATGGGGGATCCACCTCCACACCCGCGGACACCCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCAc        <  2:9378501/137‑1 (MQ=255)
                                                                                                                             cgcgAGTACGGTGGTTGCTTCACCGCAGCCTTCGCGGAAAACCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCAc        <  2:16920706/137‑1 (MQ=255)
                                                                                                                             cgcgAGTACGGTGGTGGCTTCCCCTCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCAc        <  2:14048827/137‑1 (MQ=255)
                                                                                                                             cgcgACGACGGTGGGGGCTTCACCGCCAAATCACCGGCCAACCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTGGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCAc        <  2:6480931/137‑1 (MQ=255)
                                                                                                                             cgcgAAGATGATGGTGGCTTCACCACCACGTTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCAc        <  2:9489961/137‑1 (MQ=255)
                                                                                                                             cgagAATATGGTGGTGGATCCCCCTCCGCCTTCGCGGACAACCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCAc        <  2:7232065/134‑1 (MQ=255)
                                                                                                                             cgaaaaTAAGATGGTGGATCCCCCTCCGCCTTCGCGGCCAACCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCAc        <  2:11493625/133‑1 (MQ=255)
                                                                                                                             cccaaaTATGTTGGTGGATCCGCCTCCGCGTTCGCGGAAAACCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCAc        <  2:10219450/133‑1 (MQ=255)
                                                                                                                             cacgAATACGGTGGTGGCTCCACCTCCGCCTTCGCGGAAAACCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCAc        <  2:2119989/135‑1 (MQ=255)
                                                                                                                              gcgAAGACGGTGGGGGATTCACCTCCGCATTTACGGACAACCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACg       <  1:13960368/137‑1 (MQ=255)
                                                                                                                               cgAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGc      >  1:6790095/1‑137 (MQ=255)
                                                                                                                                  aTACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcct   >  1:17390141/1‑137 (MQ=255)
                                                                                                                                  aTACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcct   >  2:1813823/1‑137 (MQ=255)
                                                                                                                                  aTACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcct   >  1:17390142/1‑137 (MQ=255)
                                                                                                                                  aTACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcct   >  2:14472991/1‑137 (MQ=255)
                                                                                                                                   tACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCACGCCACGGGATCACTTGCATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCCGCGGTGGCCGTGAGCACGCcctg  >  2:9792736/1‑137 (MQ=255)
                                                                                                                                   tACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctg  >  2:17903209/1‑137 (MQ=255)
                                                                                                                                   tACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctg  >  2:14428259/1‑137 (MQ=255)
                                                                                                                                   tACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctg  >  2:12365840/1‑137 (MQ=255)
                                                                                                                                   tACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctg  >  2:9792697/1‑137 (MQ=255)
                                                                                                                                   tACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctg  >  2:11729042/1‑137 (MQ=255)
                                                                                                                                       |                                                                                                                                     
ATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTG  >  AE015451_phaM1/5485018‑5485285

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: