Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451_phaM14,889,6230AG39.7% 130.4 / 50.3 143intergenic (+69/‑84)hyi CDS (hydroxypyruvate isomerase)/glxR CDS (tartronate semialdehyde reductase)–/–
Reads supporting (aligned to +/- strand):  ref base A (44/41);  new base G (0/56);  total (44/99)
Fisher's exact test for biased strand distribution p-value = 7.22e-13
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 143 total aligned reads displayed.

AGTACAAGCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGGAACAACTCGGGACCCTGTGGGAGCGGGCTCGCCCGCGAACACCGGCAAAGCCGGTGCCATCCACCGA‑GTTG‑TCTGCTTCGCGGGCGAGCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGCCCATGGC  >  AE015451_phaM1/4889490‑4889756
                                                                                                                                     |                                                                                                                                          
       gCCGCTGACCACCACCGAAGCGTGCCTGGGCTTGCCGAAAACCCACAACGCAATCTGAGGAACAACTCGGGACCCTGTGGGAGCGGGCTCGCCCGCGAACACCGGCAAAGCCGGTGCCATCCACCGA‑GTTG‑TCTGCt                                                                                                                                <  1:14781095/137‑1 (MQ=255)
       gCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGGAACAACTCGGGACCCTGTGGGAGCGGGCTCGCCCGCGAACACCGGCAAAGCCGGTGCCATCCACCGA‑GTTG‑TCTGCt                                                                                                                                <  1:7212302/137‑1 (MQ=255)
               ccaccaccGAGGTGGACCGGGGCTGGTGGAAAACCCACACCGCAACCTGAGAAACAACTCGGGACCCTGTGGGAGCGGGCTCTCCTGCGAACACCGGCAAAGCCGGTGCCATCCCCCGA‑GTTG‑TCTGCTTCGCGGGc                                                                                                                        <  2:15940890/137‑1 (MQ=25)
               ccaccaccGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGGAACAACTCGGGACCCTGTGGGAGCGGGGTTTCCCGCGAAAAACAGCACAGACAGGGACATCCCCCGA‑GTTT‑TGTTCTTCTCgggg                                                                                                                        >  2:16407975/1‑136 (MQ=255)
                  ccaccGAAGTGGGCCTGGGTTGGCGGAAAACCCACACCGCATCCTGGGAAACAACCTGGGACCCTTTTGGGCCGGGCCTTCCCGCGAACACCGGCAAATCCGGGGCCACCCCCCGG‑TTTG‑TTGGCTTCGCGGGCggg                                                                                                                     <  2:15983761/137‑3 (MQ=255)
                      cGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGGAACAACTCGGGACCCTGTGGGAGCGGGATCGGCCGCGAACCACAGAAAAGCCAGTGTCACCCCCCGA‑AGTG‑TCTTCCTCGCGGGGGAGgc‑‑cc                                                                                                               >  2:14959588/1‑136 (MQ=255)
                      cGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGGAACAACTCGGGACCCTGTGGGAGAGGGGTGGCTCAGGAAAAAAGGAAAAAAAGGTGCCAGCCACCAA‑AATG‑GCTGCTTCCTGGGCgggggcc                                                                                                                 >  1:15545692/1‑131 (MQ=255)
                             ggctgggggTGGCTGAAAGCCCACAACGCAGCCTGAGAAAAACTTGGGGACCTTGTGGGAGAGGTCCCTCCGTCGAACACTGGCAAATTCGGTGCCGTCCACCGA‑GTTG‑TCTGCTTCGCGGGCGACCC‑‑CGCTCCcac                                                                                                        <  2:13415246/132‑1 (MQ=255)
                                cTGGGCTGGCTGAAAACCCACAACGCAATCTGAGGAACAACTCGGGACCCTGTGGGAGCGGGATCGCCCCCCCACACCCACAAAACCAGTGTCACCCCCCCA‑GGTG‑TCTTTTTCTCCGGCGCGCC‑‑CCCTCCccc‑acg                                                                                                    >  2:4169174/1‑135 (MQ=255)
                                      gggcggAAAACCCACAATGAACCTTGGGGAACACCTTTGGACCCTGTGGGAGCGTGCTCGCCAGCGAACACCGACAAATCCGTTCCCACCCACCGG‑TTTT‑TCGGTTTCGTGGGCGACCC‑‑CGCTCCCAC‑AGGGGGCCa                                                                                              <  2:16683452/136‑1 (MQ=255)
                                       ggtgaaaacCCCACAACGCACTCGGAGAAACAACTCGGGACCTTTTGGGAGTGGCCTCCCCTGCAAACACCGAAAAATCTGTTCCCCCCCCCCGG‑TTTT‑TCTTCTTCGTGGGCAACCC‑‑CGCTCCCAC‑AGGGGGCCaa                                                                                             <  2:8622144/132‑1 (MQ=255)
                                       ggcgaaaaaCCCACAACGCACTCGGAGGAACACCTTGGGACCTTGGGGGGGTGGGCTCCCCAGAGAACACCGGAAAATCCGGTCCCATCCACTGA‑GTTG‑TTTGTTTTGGGGGCGACCC‑‑CGCTCCCAC‑AGGGGGCCaa                                                                                             <  2:12291778/132‑1 (MQ=255)
                                       ggagaaaacCCCAAAACGCAATTGGAGAAAAAATTTGGGACCTGGGGGGAGGGGGCTCACCAGAAAACACAGAAAAAGCTGGTGCCATCCACGGG‑GTTT‑TTTGTTTCGGGGACAACCC‑‑CGCTCCCAC‑AGGGGGCCaa                                                                                             <  2:9356593/132‑1 (MQ=255)
                                       ggaaaaaaaCCAACAACACAACCGGAGAAAAAACTCGGGACCTGGGGGGAGCGGGCTCCCCCGCAAACACGGGCAAATTGGGTGCCATCCACGGG‑GTTT‑TTGGTTTCGGGGGCGACCC‑‑CGCTCCCAC‑AGGGGGCCaa                                                                                             <  2:10819541/132‑1 (MQ=255)
                                           gAAAACCCACAACGCAATCTGAGGAACAACTCGGGACCCTGTGGGAGCGGGGTCTCCCTCGGACCCCCGCAAAAACAGTGTCCCCCCCCAA‑AGTG‑TTTTCTTCTCCGGCGGGCC‑‑CCCTCCCCCAAGGGGACCaaaac                                                                                          >  2:11888839/1‑135 (MQ=255)
                                           aaaaaCCAAAAACGCAACCGAAGAAACACCTAGGGACCTTGGGGGAGCGGCCTCACCAAAAAAAACAGAAAAAGCTGGTACCACCCACCGG‑GTTG‑TCGGCTTCGGGGACAACCC‑‑CGCTCCCAC‑AGGGGGCCAAAgct                                                                                         <  2:2045633/136‑1 (MQ=255)
                                            aaaccccaaaccccACTCCGAGAAACACCTTGGGACCTTGTGGGAGCGGTCTCTCCCGCAAACACCGGAAAAGCCGGTGCCCCCCACCGG‑TTTT‑TTTGTTTCGTGGGCGACCC‑‑CGCTCCCAC‑AGGGGGCCAAAgctt                                                                                        <  2:15821525/137‑1 (MQ=255)
                                            aaacCCCACAACGAAGGTGGAGAAAAAACTTGGGACCTTGGGGGGGCGGGCCCGCCAGAAAACACCGGAAAAGCTGGTACCACCCA‑CGG‑GTTGTTCGGTTTCGTGGGCGACCC‑‑CGCTCCCAC‑AGGGGGCCAAAgctt                                                                                        <  2:9252755/137‑1 (MQ=255)
                                                  caaaaCACACTCTAAGGAACAACTTGGAACCTGGGGGGAAGGGGCTCGCCCGAGAACACGGACAAAGCGGGTGCACTCCACCGA‑GTGG‑TCTGCTTCGCGGGCGATCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTga                                                                                  <  2:4557323/134‑1 (MQ=255)
                                                    ccacccAATCTGAGAAAAGGCTCGGGACCTGGGGGGAGCGGTCTCCCCAGAGAACATCGGAAACGCGGGTCCCATCCACCGG‑GTGT‑TCTGTTTCGGGGTAGACCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTgaga                                                                                <  2:1858479/132‑1 (MQ=255)
                                                    caACGCCCTCTGAGAAACACCTCGGAACCCTGTGGGAGTGGTTTCTCCCGAAAACACCGGAAAACCTGGTACCACCCACCGA‑TTGG‑TTTGCTTCGCGGGCAACCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTgaga                                                                                <  2:2984556/137‑1 (MQ=255)
                                                    caACGCACTCGGAGGAACAATTTGGGACCTTGGGGGGGTGGGTTCTCCCGAGAACACCGAAAAATCTGGTGCCACCCCCCGG‑GTTG‑TTTTCTTCGCGGGCAACCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTgaga                                                                                <  2:12376114/137‑1 (MQ=255)
                                                    aactgaaaTCTGAGAAACAACTCGGAACCTTGGGGGATCGGACTCGCCAGGAAACACC‑GAAAAACCGGGGCCCGCCAACTAGGTGG‑TTGGCTTTGTGGGGGATCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTgaga                                                                                <  2:3709651/131‑1 (MQ=255)
                                                    aaaCGCAATCTGAGAAAAAACTCGGAACCTTTTGGGAGCGGGCTTGCCCGCAAACACCGGAAAACCTGGTGCCACCCACCGG‑TTTG‑TTTTCTTCGCGGGCAACCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTgaga                                                                                <  2:13120628/136‑1 (MQ=18)
                                                    aaaCGCAATCGGAGGAACAACTTGGGACTTTGGGGGAGTGGGCTCGCCCGAGAACCCAAACAAAGTTGTTCCCCTCCCCTGA‑TTTT‑TTGGTTTCGCGGGCGACCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTgaga                                                                                <  2:13470055/136‑1 (MQ=255)
                                                      ccGCACTCTGAGAAACACTTCGGGACCTTGTGGGAGTGGGCTCGCCCGCGAACACCGGCAAATCTGTTGCCATCCACTGG‑GTTT‑TCTGTTTCGCGGGCGACCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGata                                                                              <  2:11974373/136‑1 (MQ=25)
                                                          aaTCTGAAGAACAACTTGGGACCTGGTGGGAGGGGTCTTACCTACGAACACCGAAAAATCCGGTGCCACCCACCGG‑TTTG‑TCTGCTTAGTGGGAAAACC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATcac                                                                          <  2:14202402/137‑1 (MQ=255)
                                                              ggAGGAACACTTCGGGACCTTTTGGGAGTGGGCTGGCCAGCGAAAACTGGAAAAGCTGGGGCCATCCACCGG‑GTTG‑TCTGTTTCGTGGGTGAGCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAAtt                                                                      <  2:15405234/136‑1 (MQ=255)
                                                                ggaaacaccTCGGGACCTTGTGGGGGTGGGCTCTCCCGAGAACACGGGCAAAGCTGGTTCCCTCCCCTGG‑TTTT‑TTTTTTTTGTGGGCGATCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTAc                                                                    <  2:15111742/134‑1 (MQ=255)
                                                                aGGAACACCTGGGGACCTTGGGGGAGCGTGTTCTCCCGCAAACACTGAAAAAGCGGGTGCCACCCACCGG‑TTTG‑TTTGTTTTGCGGGAGATCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTAc                                                                    <  2:11213010/137‑1 (MQ=255)
                                                                            ggACCTTTTGGGAGCGTTCTCTCCCGAGAACATCGGCAAGTCTGGTTCCGTCCCCCGA‑TTTG‑TCTGCTCCGTGGGCGATCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAAtt                                                        <  2:5237750/137‑1 (MQ=255)
                                                                             gACCTTGTGGGGGGGGTTCCGCCCGAGAACACTGGAAAGGCCGGTGCCCCCCCCGGG‑TTTG‑TTTTTTTCGCGGGCGACCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAAttt                                                       <  2:15429122/137‑1 (MQ=255)
                                                                             aaCCTTGGGGGAGGGGCCTCACCAAAGAAAACCAGAAAAGCTGGGGCCATCCACTGA‑GTTG‑TTTGTTTCGAGGTCAACCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAAttt                                                       <  2:3889396/136‑1 (MQ=255)
                                                                             aaCCCTGTGGGAGCGGTCCCGCCTGCGAAAACCGGCAAAGCCGGTGCCATCCACCGA‑GTTT‑TCTGCTTCGTGGGCGGTCC‑‑CGCTCCCAC‑AGTGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAAttt                                                       <  2:12130579/136‑1 (MQ=37)
                                                                                        gggggtcTCGCCCAAAAAAACAGGAAAATCGGGTTCCATCCCCCGG‑TTTG‑TTTGCTTCGTGGGCGACCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTaa                                            <  2:5656304/131‑1 (MQ=255)
                                                                                        agtggccTCCCCCGAAAACACCGGAAAATCCGTTTCCCCCCCCCGA‑GTTT‑TCTCCTCTACGTCCGAACC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTaa                                            <  1:11305931/131‑1 (MQ=255)
                                                                                            ggCTCGCCTGAAAACACAGGCAAAACCGGGCCCATCCACTGA‑GGTG‑TCTGCTTCGCGGGCGGACC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTGGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATc                                        <  2:9694325/137‑1 (MQ=255)
                                                                                                 gCCCGGAAACACCGGAAAAGCCGGTGCCATCCACCGA‑GTTG‑TCTGCTTGGTGGGTGAGCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGttt                                   <  2:16042328/137‑1 (MQ=255)
                                                                                                  cccgaaaaCACCGACAAAACCGGTGCCATCCCCGGG‑GTTG‑TTTGCTTTGCGGGCGAGCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTc                                  <  2:869272/131‑1 (MQ=255)
                                                                                                  cccGAGAACACCGAAAAATCTGGTGCCACCCACAGA‑GTTG‑TCTGTTTCGCGGGCGACCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTc                                  <  2:2145220/137‑1 (MQ=255)
                                                                                                   ccgaaaaCACCGGAAAAGCCGGTACCACCCACTGG‑GTTT‑TCTTCTTCGGGGGTGACCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCa                                 <  2:3682033/132‑1 (MQ=255)
                                                                                                    tgcgAACACTGAAAAACCCGGTGCCGTCCACTGA‑GTTG‑TCTGTTTTCCGGGCGAGCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGCTTCAt                                <  2:8425042/136‑1 (MQ=255)
                                                                                                         acacTGGAAAGGCCGGTGCCCCCCACCGA‑GTTG‑TCGGTTTCGTGGACGAGCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCa                           <  2:14192994/137‑1 (MQ=255)
                                                                                                         acacTGACAAAGCGGGTGCCACCCACTGG‑TTTT‑TCTTCTTCTCGGGCGCCCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCa                           <  2:4271353/137‑1 (MQ=255)
                                                                                                         acacCGGAAAATCTGGTGCCATCCACTGG‑GTTG‑TTTGCTTCGTGGGCGACCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCa                           <  2:13891690/137‑1 (MQ=255)
                                                                                                          cacCGGCAAAGCCGGTGCCATCCACCGA‑GTTG‑TCTGCTTCGCGGGCGAGCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCAc                          >  1:10022970/1‑137 (MQ=255)
                                                                                                           acCAGCAAAGACGGTGCCATCCACCAA‑GGTG‑TCTGCTTCGCGGGCGGACC‑‑CGCTCCCAC‑AAGGGGCCAAAACTTGCTTGAGATATCACAATTGCAAAGGGGTAACTTCTCATGGCCAAAACCGGTTTCATTGGCAcc                         >  2:4507868/1‑137 (MQ=255)
                                                                                                              ggCAAATCTGGTGCCATCCCCTGG‑TTTG‑TCTGCTTCGTGGGCGACCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGc                      <  2:2421094/137‑1 (MQ=255)
                                                                                                              ggCAAAGCGGGTGCCATCCACCGA‑GTTG‑TCTGCTTCGCGGGCGGTCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGc                      <  2:15939624/137‑1 (MQ=255)
                                                                                                              ggCAAAGCCGGTGCCATCCACCGA‑GTTG‑TCTGCTTCGCGGGCGAGCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTACAATCGGTTTCATCGGCACCGGc                      >  1:3293683/1‑137 (MQ=255)
                                                                                                              gacAAATCCGTTCCCGTCCACCGG‑GTTT‑TTTGCTTCGCGGGCGATCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGc                      <  2:11389189/135‑1 (MQ=255)
                                                                                                              gacAAAGCTGTTGCCATCCACCGA‑GTTG‑TTTGCTTCGCGGGCGACCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGc                      <  2:8903263/135‑1 (MQ=255)
                                                                                                              gacAAACCTGGTCCCCTCCCCCGT‑GTTT‑TTTTTTTCGCGGGCGACCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGc                      <  1:4560915/135‑1 (MQ=255)
                                                                                                              gaaaaaGCCGTTACCCTCCACAGG‑TTTG‑TCGGTTTAGCGGGCAATCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGc                      <  2:7344855/134‑1 (MQ=255)
                                                                                                               gaaaaGCGGGTGCCATCCCCTGG‑GTTG‑TTTGTTTGGGGGGCGACCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGca                     <  2:6306735/135‑1 (MQ=255)
                                                                                                                cAAAGCCGGTGCCATCCACCGA‑GTTG‑TCTGCTTCGCGGGCGAGCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGcat                    >  1:3817755/1‑137 (MQ=255)
                                                                                                                cAAAGCCGGTGCCATCCACCGA‑GTTG‑TCTGCTTCGCGGGCGAGCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGcat                    >  2:10582112/1‑137 (MQ=255)
                                                                                                                cAAAGCCGGTGCCATCCACCGA‑GTTG‑TCTGCTTCGCGGGCGAGCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGcat                    >  2:2753671/1‑137 (MQ=255)
                                                                                                                cAAAGCCGGTGCCATCCACCGA‑GTTG‑TCTGCTTCGCGGGCGAGCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGcat                    >  2:2753630/1‑137 (MQ=255)
                                                                                                                cAAAGCCGGTGCCATCCACCGA‑GTTG‑TCTGCTTCGCGGGCGAGCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGcat                    >  1:13151232/1‑137 (MQ=255)
                                                                                                                cAAAGCCGGTGCCATCCACCGA‑GTTG‑TCTGCTTCGCGGGCGAGCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGcat                    >  2:16982297/1‑137 (MQ=255)
                                                                                                                cAAAGCCGGTGCCATCCACCGA‑GTTG‑TCTGCTTCGCGGGCGAGCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGcat                    >  1:1367315/1‑137 (MQ=255)
                                                                                                                cAAAGCCGGTGCCATCCACCGA‑GTTG‑TCTGCTTCGCGGGCGAGCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGcat                    >  1:16144719/1‑137 (MQ=255)
                                                                                                                cAAAGCCGGTGCCATCCACCGA‑GTTG‑TCTGCTTCGCGGGCGAGCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGcat                    >  1:17132215/1‑137 (MQ=255)
                                                                                                                cAAAGCCGGTGCCATCCACCGA‑GTTG‑TCTGCTTCGCGGGCGAGCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGcat                    >  1:5157756/1‑137 (MQ=255)
                                                                                                                cAAAGCCGGTGCCATCCACCGA‑GTTG‑TCTGCTTCGCGGGCGAGCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGcat                    >  1:3817727/1‑137 (MQ=255)
                                                                                                                cAAAGCCGGTGCCATCCACCGA‑GTTG‑TATGCTTCGCGGGCGAGCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGcat                    >  1:14085839/1‑137 (MQ=255)
                                                                                                                cAAAGCCGGTGCCACCCACCGA‑GTTG‑TCTGCTTCGCGGGCGAGCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAACTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGcat                    >  1:11325186/1‑137 (MQ=255)
                                                                                                                  acaccGGTCCCATCCCCCGG‑TTTT‑CCTGCTCATCGGATGAGCC‑‑CACTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGcatca                  <  1:2737675/134‑1 (MQ=255)
                                                                                                                   gtctggTGCCATCCCCTGG‑TTTT‑TCTGTTTCGTGGGCGCCCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGcatcat                 <  1:606636/133‑1 (MQ=255)
                                                                                                                   atctgttGCCCTCCCCGGA‑GTTG‑TTTGCTTCAGGGGCAACCC‑‑CTCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAACGACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGcatcat                 <  2:11356578/131‑1 (MQ=255)
                                                                                                                    tctgttGCCATCCCCCGG‑TTTT‑TTTGCTTAGCGTGCGGACC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTACTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATg                <  2:8653852/132‑1 (MQ=255)
                                                                                                                    tctggTGCAATCCACTGA‑GTTT‑TCTGCTCCGTGGGCAAGCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATg                <  2:8653814/134‑1 (MQ=255)
                                                                                                                    gttggTGCCATCAATTGG‑GTTG‑TTTGTTTCGTGGGAGATCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATg                <  2:14441501/134‑1 (MQ=255)
                                                                                                                    ccggttGCCATCCCCCGG‑TTTT‑TTTTTTTCGCGGACACCCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATg                <  1:10268236/132‑1 (MQ=255)
                                                                                                                    actggTTCCATCTCCCAA‑TTTG‑TTTTCTTAGCGGGCGAGCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATg                <  2:12693508/134‑1 (MQ=255)
                                                                                                                      cGGTGCCATCCACCGG‑TTTT‑CCTGCTTCGAGTGTGAGCCCGCGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATg                <  2:9872558/137‑1 (MQ=255)
                                                                                                                      cGGTGCCATCCACCGA‑GTTG‑TCTGCTTCGCGGGCGAGCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATggg              >  1:6425490/1‑137 (MQ=255)
                                                                                                                      cGGTGCCATCCACCGA‑GTTG‑TCTGCTTCGCGGGCGAGCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATggg              >  1:9374620/1‑137 (MQ=255)
                                                                                                                       ggTGCCATCCACCGA‑GTTG‑TCTGCTTCGCGGGCGAGCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCg                              <  2:5507994/120‑1 (MQ=255)
                                                                                                                        ttGCCATCCACCGA‑GTTG‑TCTGTTTCGTGGGCAACCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCa            <  1:13900318/136‑1 (MQ=255)
                                                                                                                        gTGCCATCCACCGA‑GTTG‑TCTGCTTCGCGGGCGAGCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCa            >  2:10147335/1‑137 (MQ=255)
                                                                                                                          gCCATCCACCGG‑TTTT‑TTTGTTTCGTGGGCGATCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGGGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAg          <  2:16828139/137‑1 (MQ=255)
                                                                                                                           ccgtCCACCGG‑TTTT‑TTTGTTTCGCGGGCGACCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGc         <  2:7471916/134‑1 (MQ=255)
                                                                                                                           ccacccacCGG‑GTTT‑TCTGTTTCGGGGGAAACCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGc         <  1:2806675/137‑1 (MQ=255)
                                                                                                                           ccacccacCGA‑TTTT‑TCTGTTTCGTGGGCGACCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGc         <  2:11755947/137‑1 (MQ=255)
                                                                                                                           ccATCCCCGGG‑GTTG‑TCTGTTTCGTGGGCGACCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGc         <  1:10085962/137‑1 (MQ=255)
                                                                                                                           ccATCCCCCGA‑TTTT‑TTTGCTCCGCGGGAGACCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACACTTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATAGGCACCGGCATCATGGGCAAGc         <  2:1167024/137‑1 (MQ=255)
                                                                                                                           ccATCCCCCGA‑TTTG‑TTTGCTTCGCGGGCGATCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGc         <  2:15565399/137‑1 (MQ=255)
                                                                                                                           ccATCCCCCGA‑GTTT‑TTTTCTTCGCGGGCGAGCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGc         <  2:8199727/137‑1 (MQ=255)
                                                                                                                           ccATCCACTGG‑GTTT‑TTTGTTTCGTGGGCGAGCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGc         <  2:10598874/137‑1 (MQ=255)
                                                                                                                           ccATCCACGGG‑GTTT‑TCTTTTTTGTGGTCGACCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGc         <  2:6509812/137‑1 (MQ=255)
                                                                                                                           ccATCCACCGG‑TTTT‑TTTTTTTCGGGGGCGACCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGATGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAgt         <  1:16645905/137‑2 (MQ=255)
                                                                                                                           ccATCCACCGG‑GTTG‑TCTGCTTCGCGGGCGACCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGc         <  2:8199728/137‑1 (MQ=255)
                                                                                                                           ccATCCACCGA‑GTTT‑TCGGTTTCGGGGGAAAGCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGc         <  2:16398829/137‑1 (MQ=255)
                                                                                                                            cATCCACCGA‑GTTG‑TCTGCTTCGCGGGCGAGCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGcc        >  1:7453873/1‑137 (MQ=255)
                                                                                                                            aaTCCAACGA‑GTTG‑TATGCTTCGCGGGCGAGCC‑‑CGCTCCCAA‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTGAAGTCGGTTTCATCGGCACCGGCATCATGGGCAAGcc        >  2:1100256/2‑137 (MQ=255)
                                                                                                                                 aCCGA‑GTTG‑TCTGCTTCGCGGGCGAGCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGCCCATgg   >  1:5743194/1‑137 (MQ=255)
                                                                                                                                 aCCGA‑GTTG‑TCTGCTTCGCGGGCGAGCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGCCCATgg   >  2:9097302/1‑137 (MQ=255)
                                                                                                                                 aCCGA‑GTTG‑TCTGCTTCGCGGGCGAGCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGCCCATgg   >  2:15010016/1‑137 (MQ=255)
                                                                                                                                     |                                                                                                                                          
AGTACAAGCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCTGAGGAACAACTCGGGACCCTGTGGGAGCGGGCTCGCCCGCGAACACCGGCAAAGCCGGTGCCATCCACCGA‑GTTG‑TCTGCTTCGCGGGCGAGCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGCCCATGGC  >  AE015451_phaM1/4889490‑4889756

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