Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451_phaM14,889,6240GT32.1% 210.5 / 15.9 143intergenic (+70/‑83)hyi CDS (hydroxypyruvate isomerase)/glxR CDS (tartronate semialdehyde reductase)–/–
Reads supporting (aligned to +/- strand):  ref base G (41/54);  new base T (0/45);  total (44/99)
Fisher's exact test for biased strand distribution p-value = 3.99e-09
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 146 total aligned reads displayed.

AGTACAAGCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCT‑GAGGAACAACTCGGGACCCTGTGGGAGCGGGCTCGCCCGCGAACACCGGCAAAGCCGGTGCCATCCACCGA‑GTTG‑TCTGCTTCGCGGGCGAGCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGCCCATGGCTCAG  >  AE015451_phaM1/4889490‑4889760
                                                                                                                                        |                                                                                                                                            
aGTACAAGCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCT‑GAGGAACAACTCGGGACCCTGTGGGAGCGGGCTCGCCCGCGAACACCGGCAAAGCCGGTGCCATCCACCGA‑Gtt                                                                                                                                            <  1:3912794/137‑1 (MQ=255)
      aGCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCT‑GAGGAACAACTCGGGACCCTGTGGGAGCGGGCTCGCCCGCGAACACCGGCAAAGCCGGTGCCATCCACCGA‑GTTG‑TCTGc                                                                                                                                     <  1:15041223/137‑1 (MQ=255)
       gCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCT‑GAGGAACAACTCGGGACCCTGTGGGAGCGGGCTCGCCCGCGAACACCGGCAAAGCCGGTGCCATCCACCGA‑GTTG‑TCTGCt                                                                                                                                    <  1:7212302/137‑1 (MQ=255)
               ccaccaccGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCT‑GAGGAACAACTCGGGACCCTGTGGGAGCGGGGTTTCCCGCGAAAAACAGCACAGACAGGGACATCCCCCGA‑GTTT‑TGTTCTTCTCgggg                                                                                                                            >  2:16407975/1‑136 (MQ=255)
                  ccaccGAAGTGGGCCTGGGTTGGCGGAAAACCCACACCGCATCCT‑GGGAAACAACCTGGGACCCTTTTGGGCCGGGCCTTCCCGCGAACACCGGCAAATCCGGGGCCACCCCCCGG‑TTTG‑TTGGCTTCGCGGGCggg                                                                                                                         <  2:15983761/137‑3 (MQ=255)
                      cGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCT‑GAGGAACAACTCGGGACCCTGTGGGAGCGGGATCGGCCGCGAACCACAGAAAAGCCAGTGTCACCCCCCGA‑AGTG‑TCTTCCTCGCGGGGGAGgc‑‑cc                                                                                                                   >  2:14959588/1‑136 (MQ=255)
                      cGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCT‑GAGGAACAACTCGGGACCCTGTGGGAGAGGGGTGGCTCAGGAAAAAAGGAAAAAAAGGTGCCAGCCACCAA‑AATG‑GCTGCTTCCTGGGCgggggcc                                                                                                                     >  1:15545692/1‑131 (MQ=255)
                             ggctgggggTGGCTGAAAGCCCACAACGCAGCCT‑GAGAAAAACTTGGGGACCTTGTGGGAGAGGTCCCTCCGTCGAACACTGGCAAATTCGGTGCCGTCCACCGA‑GTTG‑TCTGCTTCGCGGGCGACCC‑‑CGCTCCcac                                                                                                            <  2:13415246/132‑1 (MQ=255)
                                cTGGGCTGGCTGAAAACCCACAACGCAATCT‑GAGGAACAACTCGGGACCCTGTGGGAGCGGGGTTGCCCTCTAACACCGGCCAAAACAGTTTCCTCCCCCAA‑GGTG‑TCTTCTTCCCTGGCGAACC‑‑CCCTCCccc‑aag                                                                                                        >  2:1595521/1‑135 (MQ=255)
                                      gggcggAAAACCCACAATGAACCTT‑GGGGAACACCTTTGGACCCTGTGGGAGCGTGCTCGCCAGCGAACACCGACAAATCCGTTCCCACCCACCGG‑TTTT‑TCGGTTTCGTGGGCGACCC‑‑CGCTCCCAC‑AGGGGGCCa                                                                                                  <  2:16683452/136‑1 (MQ=255)
                                       ggtggAAAACCCACAAC‑CAAACTGGAGAAAAACCTCGGGACCTTGGGGGGGTGGGCTTGCCCGAAAACACGGGAAAGTTTGGTGCCACCCACCGG‑TTGG‑TTTGTTTTGTGGGCGATCC‑‑CGCTCCCAC‑AGGGGGCCaa                                                                                                 <  2:10841930/133‑1 (MQ=255)
                                       ggtgaaaacCCCACAACGCACTCG‑GAGAAACAACTCGGGACCTTTTGGGAGTGGCCTCCCCTGCAAACACCGAAAAATCTGTTCCCCCCCCCCGG‑TTTT‑TCTTCTTCGTGGGCAACCC‑‑CGCTCCCAC‑AGGGGGCCaa                                                                                                 <  2:8622144/132‑1 (MQ=255)
                                       ggcgaaaaaCCCACAACGCACTCG‑GAGGAACACCTTGGGACCTTGGGGGGGTGGGCTCCCCAGAGAACACCGGAAAATCCGGTCCCATCCACTGA‑GTTG‑TTTGTTTTGGGGGCGACCC‑‑CGCTCCCAC‑AGGGGGCCaa                                                                                                 <  2:12291778/132‑1 (MQ=255)
                                       ggagaaaacCCCAAAACGCAATTG‑GAGAAAAAATTTGGGACCTGGGGGGAGGGGGCTCACCAGAAAACACAGAAAAAGCTGGTGCCATCCACGGG‑GTTT‑TTTGTTTCGGGGACAACCC‑‑CGCTCCCAC‑AGGGGGCCaa                                                                                                 <  2:9356593/132‑1 (MQ=255)
                                           gAAAACCCACAACGCAATCT‑GAGGAACAACTCGGGACCCTGTGGGAGCGGGGTCTCCCTCGGACCCCCGCAAAAACAGTGTCCCCCCCCAA‑AGTG‑TTTTCTTCTCCGGCGGGCC‑‑CCCTCCCCCAAGGGGACCaaaac                                                                                              >  2:11888839/1‑135 (MQ=255)
                                            aaaccccaaaccccACTCC‑GAGAAACACCTTGGGACCTTGTGGGAGCGGTCTCTCCCGCAAACACCGGAAAAGCCGGTGCCCCCCACCGG‑TTTT‑TTTGTTTCGTGGGCGACCC‑‑CGCTCCCAC‑AGGGGGCCAAAgctt                                                                                            <  2:15821525/137‑1 (MQ=255)
                                            aaacCCCACAACGAAGGTG‑GAGAAAAAACTTGGGACCTTGGGGGGGCGGGCCCGCCAGAAAACACCGGAAAAGCTGGTACCACCCA‑CGG‑GTTGTTCGGTTTCGTGGGCGACCC‑‑CGCTCCCAC‑AGGGGGCCAAAgctt                                                                                            <  2:9252755/137‑1 (MQ=255)
                                            aaaaCCCACAACGCAATCT‑GAGGAACAACTCGGGACCCTGTGGGAGCGGGGTTGCCCATGCACACCCGCCAAACCCGTTTCCCCCCCCCA‑GTTG‑TTTTTCTCTCGGGCGCGCC‑‑CTCCCCCAC‑CAGGGGAAAAAgact                                                                                            >  2:5121978/1‑134 (MQ=255)
                                            aaaaCCCACAACGCAATCT‑GAGGAACAACTCGGGACCCTGTGGGAGCGGGATCGCCCCCGAACCCCAGCAAAAAAAGTGTCACCCCCAAA‑GTTG‑TCTTTTTCTCTAGCGAGGC‑‑CTCTCCCACAAGGGGGACaaaact                                                                                             >  2:3780348/1‑134 (MQ=255)
                                            aaaaCACACAACGCACTCT‑GAGAAACAACTGGGGACCTTGGGGGAGCGGGCTCGCCCGAAAACACGGGAAAAGCCGGGGCAACCCACTGA‑GTTG‑TCTGCTTCGCGGGCGACCC‑‑CGCTCCCAC‑AGGGGGCCAAAgctt                                                                                            <  2:8642873/137‑1 (MQ=255)
                                                  caaaaCACACTCT‑AAGGAACAACTTGGAACCTGGGGGGAAGGGGCTCGCCCGAGAACACGGACAAAGCGGGTGCACTCCACCGA‑GTGG‑TCTGCTTCGCGGGCGATCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTga                                                                                      <  2:4557323/134‑1 (MQ=255)
                                                    ccacccAATCT‑GAGAAAAGGCTCGGGACCTGGGGGGAGCGGTCTCCCCAGAGAACATCGGAAACGCGGGTCCCATCCACCGG‑GTGT‑TCTGTTTCGGGGTAGACCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTgaga                                                                                    <  2:1858479/132‑1 (MQ=255)
                                                    caACGCCCTCT‑GAGAAACACCTCGGAACCCTGTGGGAGTGGTTTCTCCCGAAAACACCGGAAAACCTGGTACCACCCACCGA‑TTGG‑TTTGCTTCGCGGGCAACCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTgaga                                                                                    <  2:2984556/137‑1 (MQ=255)
                                                    aactgaaaTCT‑GAGAAACAACTCGGAACCTTGGGGGATCGGACTCGCCAGGAAACACC‑GAAAAACCGGGGCCCGCCAACTAGGTGG‑TTGGCTTTGTGGGGGATCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTgaga                                                                                    <  2:3709651/131‑1 (MQ=255)
                                                    aaaCGCAATCT‑GAGAAAAAACTCGGAACCTTTTGGGAGCGGGCTTGCCCGCAAACACCGGAAAACCTGGTGCCACCCACCGG‑TTTG‑TTTTCTTCGCGGGCAACCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTgaga                                                                                    <  2:13120628/136‑1 (MQ=18)
                                                    aaaCGCAATCG‑GAGGAACAACTTGGGACTTTGGGGGAGTGGGCTCGCCCGAGAACCCAAACAAAGTTGTTCCCCTCCCCTGA‑TTTT‑TTGGTTTCGCGGGCGACCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTgaga                                                                                    <  2:13470055/136‑1 (MQ=255)
                                                      ccGCACTCT‑GAGAAACACTTCGGGACCTTGTGGGAGTGGGCTCGCCCGCGAACACCGGCAAATCTGTTGCCATCCACTGG‑GTTT‑TCTGTTTCGCGGGCGACCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGata                                                                                  <  2:11974373/136‑1 (MQ=25)
                                                          aaTCT‑GAAGAACAACTTGGGACCTGGTGGGAGGGGTCTTACCTACGAACACCGAAAAATCCGGTGCCACCCACCGG‑TTTG‑TCTGCTTAGTGGGAAAACC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATcac                                                                              <  2:14202402/137‑1 (MQ=255)
                                                               ggAGGAACACTTCGGGACCTTTTGGGAGTGGGCTGGCCAGCGAAAACTGGAAAAGCTGGGGCCATCCACCGG‑GTTG‑TCTGTTTCGTGGGTGAGCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAAtt                                                                          <  2:15405234/136‑1 (MQ=255)
                                                                 gggAAAAATTCGGGACCCTGTGTGAGCGGGCTTGCCCGCAAACACTGACAAATCTGGTGCCATCCCCCGA‑TTTG‑TCTTCTTCGCGGCCGCTCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTAc                                                                        <  2:16390188/136‑1 (MQ=255)
                                                                 ggaaacaccTCGGGACCTTGTGGGGGTGGGCTCTCCCGAGAACACGGGCAAAGCTGGTTCCCTCCCCTGG‑TTTT‑TTTTTTTTGTGGGCGATCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTAc                                                                        <  2:15111742/134‑1 (MQ=255)
                                                                 aGGCACCACTCGGGACCCTTTGGGAGTGCTCTCGCCCGAGAACACCGGCAAAGCGGGTCCCCCCCACAGG‑TTTG‑TCTGTTTTGCGGGCGACCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTAc                                                                        <  2:11264442/137‑1 (MQ=255)
                                                                              gACCTTGTGGGGGGGGTTCCGCCCGAGAACACTGGAAAGGCCGGTGCCCCCCCCGGG‑TTTG‑TTTTTTTCGCGGGCGACCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAAttt                                                           <  2:15429122/137‑1 (MQ=255)
                                                                              gACCCTGTGGGAGCGGTTTCTCCCGAAAACCCCGACAAATCCGGTTCCCTCCCCCGG‑TTTT‑TCTTTTTCGGGGGCGATCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAAttt                                                           <  2:2506441/137‑1 (MQ=255)
                                                                              aaCCTTGGGGGAGGGGCCTCACCAAAGAAAACCAGAAAAGCTGGGGCCATCCACTGA‑GTTG‑TTTGTTTCGAGGTCAACCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAAttt                                                           <  2:3889396/136‑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)
                                                                                             gccTCGCCCGCGAACACCGGAAACCCTGTTGCCCCCCCCCGA‑TTTG‑TTTTCTTTTTGGGCGAGCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATc                                            <  2:12405726/135‑1 (MQ=255)
                                                                                                   cccgaaaaCACCGACAAAACCGGTGCCATCCCCGGG‑GTTG‑TTTGCTTTGCGGGCGAGCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTc                                      <  2:869272/131‑1 (MQ=255)
                                                                                                   cccGCGAACACTGACAACGTCGTTTCCGTCCACCGG‑TTTT‑TCTGTTTCGGGGGCGACCC‑‑CGCCCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTc                                      <  2:3382575/137‑1 (MQ=255)
                                                                                                   cccGAGAACACCGAAAAATCTGGTGCCACCCACAGA‑GTTG‑TCTGTTTCGCGGGCGACCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTc                                      <  2:2145220/137‑1 (MQ=255)
                                                                                                    ccGAGAACACCGACAAAGCTGGTGCCATCCCCTGG‑TTTG‑TCTGCTTCGCGGGCGACCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCa                                     <  2:12154882/137‑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)
                                                                                                               ggaaaaCCCGGTGCCCTCCCCCGA‑TTTT‑TCCTCTTTGCGGTTGACCC‑‑CGCTCCCAC‑ATGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGc                          <  2:213548/134‑1 (MQ=255)
                                                                                                               ggCAAAGCGGGTGCCATCCACCGA‑GTTG‑TCTGCTTCGCGGGCGGTCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGc                          <  2:15939624/137‑1 (MQ=255)
                                                                                                               ggCAAAGCCGGTGCCATCCACCGA‑GTTG‑TTTGCTTCGCGGGCGATCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCGTGCGTGAGATATCACAATTACAAAGAGGTAATTTCTCAGGGCTAAAATCGGTTTCATCGGCACCGGc                          <  2:11972466/137‑1 (MQ=255)
                                                                                                               gacAAAGCTGTTGCCATCCACCGA‑GTTG‑TTTGCTTCGCGGGCGACCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGc                          <  2:8903263/135‑1 (MQ=255)
                                                                                                               gaaaaaGCCGTTACCCTCCACAGG‑TTTG‑TCGGTTTAGCGGGCAATCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGc                          <  2:7344855/134‑1 (MQ=255)
                                                                                                                gaaaaGCGGGTTCCACCCACGGG‑TTTG‑TTTGTTTCGTGGGCAACCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGca                         <  2:17119296/135‑1 (MQ=255)
                                                                                                                gaaaaGCGGGTGCCATCCCCTGG‑GTTG‑TTTGTTTGGGGGGCGACCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGca                         <  2:6306735/135‑1 (MQ=255)
                                                                                                                 cAAAGCCGGTGCCATCCACCGA‑GTTG‑TCTGCTTCGCGGGCGAGCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGcat                        >  1:17132215/1‑137 (MQ=255)
                                                                                                                 cAAAGCCGGTGCCATCCACCGA‑GTTG‑TCTGCTTCGCGGGCGAGCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGcat                        >  1:1367315/1‑137 (MQ=255)
                                                                                                                 cAAAGCCGGTGCCATCCACCGA‑GTTG‑TCTGCTTCGCGGGCGAGCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGcat                        >  1:3817727/1‑137 (MQ=255)
                                                                                                                 cAAAGCCGGTGCCATCCACCGA‑GTTG‑TCTGCTTCGCGGGCGAGCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGcat                        >  1:5157756/1‑137 (MQ=255)
                                                                                                                 cAAAGCCGGTGCCATCCACCGA‑GTTG‑TCTGCTTCGCGGGCGAGCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGcat                        >  1:13151232/1‑137 (MQ=255)
                                                                                                                 cAAAGCCGGTGCCATCCACCGA‑GTTG‑TCTGCTTCGCGGGCGAGCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGcat                        >  2:2753630/1‑137 (MQ=255)
                                                                                                                 cAAAGCCGGTGCCATCCACCGA‑GTTG‑TCTGCTTCGCGGGCGAGCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGcat                        >  2:5560568/1‑137 (MQ=255)
                                                                                                                 cAAAGCCGGTGCCATCCACCGA‑GTTG‑TCTGCTTCGCGGGCGAGCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGcat                        >  2:10582112/1‑137 (MQ=255)
                                                                                                                 cAAAGCCGGTGCCATCCACCGA‑GTTG‑TCTGCTTCGCGGGCGAGCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAAGTCATTTCTCATGGCTACACTCGGTTTCATCGGCACCGGcat                        >  1:3656190/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)
                                                                                                                     gttggTGCCATCAATTGG‑GTTG‑TTTGTTTCGTGGGAGATCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATg                    <  2:14441501/134‑1 (MQ=255)
                                                                                                                     ccggttGCCATCCCCCGG‑TTTT‑TTTTTTTCGCGGACACCCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATg                    <  1:10268236/132‑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)
                                                                                                                        ggTGCCGTCCCCTGG‑TTTG‑TCGGTTTGGGGGGCGACCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGc                 <  1:1655353/137‑1 (MQ=255)
                                                                                                                        ggTGCCATCCACCGA‑GTTG‑TCTGCTTCGCGGGCGAGCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCg                                  >  1:5507994/1‑120 (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)
                                                                                                                         gTGCCATCCACCGA‑GTTG‑TCTGCTTCGCGGGCGAGCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCa                >  2:4226236/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)
                                                                                                                            ccATCCCCGGG‑GTTG‑TCTGTTTCGTGGGCGACCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGc             <  1:10085962/137‑1 (MQ=255)
                                                                                                                            ccATCCCCCGA‑TTTT‑TTTGCTCCGCGGGAGACCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACACTTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATAGGCACCGGCATCATGGGCAAGc             <  2:1167024/137‑1 (MQ=255)
                                                                                                                            ccATCCCCCGA‑TTTT‑TTTATTTCGCGGGCGACCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGc             <  2:6509761/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)
                                                                                                                            ccATCCACCGG‑TTTT‑TTTTTTTCGGGGGCGACCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGATGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAgt             <  1:16645905/137‑2 (MQ=255)
                                                                                                                            ccATCCACCGG‑GTTG‑TTGGTTTTGTGGGAGAGCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGc             <  2:8421932/137‑1 (MQ=255)
                                                                                                                            ccATCCACCGG‑GTTG‑TCTGCTTCGCGGGCGACCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGc             <  2:8199728/137‑1 (MQ=255)
                                                                                                                            ccATCCACCGG‑GTTG‑TCGGCTTCGTGGACGACCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATAGGCACCGGCATCATGGGCAAGc             <  2:1166980/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       >  2:9097302/1‑137 (MQ=255)
                                                                                                                                  aCCGA‑GTTG‑TCTGCTTCGCGGGCGAGCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGCCCATgg       >  2:15010016/1‑137 (MQ=255)
                                                                                                                                        gTTT‑TCTGCTCCGCGGACAACCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGCCCATGGCTCAg  <  1:3924673/137‑1 (MQ=255)
                                                                                                                                        gTTG‑TCTGCTTCGCGGGCGAGCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGCCCATGGCTCAg  >  1:6969533/1‑137 (MQ=255)
                                                                                                                                        gTTG‑TCTGCTTCGCGGGCGAGCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGCCCATGGCTCAg  >  2:12678571/1‑137 (MQ=255)
                                                                                                                                        |                                                                                                                                            
AGTACAAGCCGCTGACCACCACCGAAGCGGGCCTGGGCTGGCTGAAAACCCACAACGCAATCT‑GAGGAACAACTCGGGACCCTGTGGGAGCGGGCTCGCCCGCGAACACCGGCAAAGCCGGTGCCATCCACCGA‑GTTG‑TCTGCTTCGCGGGCGAGCC‑‑CGCTCCCAC‑AGGGGGCCAAAGCTTGCTTGAGATATCACAATTACAAAGAGGTAATTTCTCATGGCTAAAATCGGTTTCATCGGCACCGGCATCATGGGCAAGCCCATGGCTCAG  >  AE015451_phaM1/4889490‑4889760

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