Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451_phaM14,675,3600TG47.1% 210.7 / 76.5 246intergenic (+59/+45)PP_4137 CDS (putative Outer membrane siderophore receptor)/chrR CDS (Chromate reductase)–/–
Reads supporting (aligned to +/- strand):  ref base T (84/36);  new base G (0/107);  total (84/162)
Fisher's exact test for biased strand distribution p-value = 1.37e-33
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 247 total aligned reads displayed.

CCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCT‑‑‑‑GCAAGAACGGTCGCGTCTTGTCATTGAGCTTGCCCGAA  >  AE015451_phaM1/4675225‑4675495
                                                                                                                                       |                                                                                                                                                
 tgtccAGGACAGTGCCCGTTCGACACTGGGCGAAGGGTTTGGTGGGTGTGCTGGGGTTTCCGAAAGCTCCCTCTAATGGGTGGGAGTGGGTTTACCCGCGAATGCGTCATGGAATTTGCCAACGCATTCGCGGGGaa                                                                                                                                                <  2:9485346/134‑1 (MQ=255)
 cGGCCAGGACAGTGCCTGTTCGACACTGGGCGAAGGGTTGGGCGTGTATGCCGGGGTTTCCTACAGTTTCCGCTGACGGGTGGGCGCGGGTTTACCCGAGAATGCGTCGTGGAATTTGCCGACGCATTCGTGGGGaa                                                                                                                                                <  2:7417269/137‑1 (MQ=255)
 cGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCTGGGTTTCCTACAGCCCCCGCGGGTGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTaa                                                                                                                                                <  2:1380556/137‑1 (MQ=38)
 cGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGGTTGCCCCCAAACACGTCCCTGAACTTTCCGAAACACTTTCGGGTaa                                                                                                                                                >  2:15166329/1‑137 (MQ=255)
 cGGCCAGAACAATGCCCGTTCGACACTGGGCGAAGGGTTTTGCGTGTATGCCGGGGTTTCCTACAGCTTCATCTGATGGGGGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTaa                                                                                                                                                <  2:13106311/137‑1 (MQ=37)
 cGGCCAGAACAATGCCCGTTCGACACCGTACAAAGGATTTGGCGTGCATACCGGGGTTTCTTACAGCTCCCGTTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTaa                                                                                                                                                <  2:10754717/137‑1 (MQ=37)
 cGGCAAGGACAATGCCTGTTCGACACGGGGCGAAGGATTTGTTGTGTGTGCCGGGGTTTCCTACAGCTTCCTCTGACGGGTGGGCGCGGTTTTACCCGCGAAGGCGTCGTTGAGTTTGCCGACGCATTCGCGGGGaa                                                                                                                                                <  2:12380303/137‑1 (MQ=255)
      gggAAAATGCCCTTTCGACACCGGGCGAAGGATTTGGTGTGTATGCCGGGGTTTCCTACATCTTCCGCGGATGGGTGGGCGCGGTTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGTGGGGAAACCCg                                                                                                                                           <  2:8910732/136‑1 (MQ=37)
            attcccGTTCGACACCGTTCAAAGGATTGGGCGTTTATGCTGGGTTCTCCCACATCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCAAATTCTTCATTGACTTTTCTGACGCGTTAGCGGGGAAACCCGCTCCTa                                                                                                                                     <  2:11311994/134‑1 (MQ=255)
                ccGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTCCCCCCGAACGCGTCATTTCATTTGCCGACACACTCGCAGGTAAAACCACCCCcacaca                                                                                                                                 >  2:16004361/1‑137 (MQ=255)
                ccGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGGTTTCCCTCCAACTCCTCCCTTAACTTTCCCACGCACTCGTCAGTAAACCCCCCCCTCcaca                                                                                                                                 >  2:8709506/1‑137 (MQ=255)
                ccGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGGTTCCCCGCGAAGGTGTCAGTTAATTCGTCGACGCAATCTCAGGTAAACCCACCCCtcccca                                                                                                                                 >  2:7709355/1‑132 (MQ=255)
                 tgTTCAACACTGGGAGAAGGGTTTGGTGTGTATGCTGGGGTTTCCTACATCTTCCGCGGACGGGGGGGCGTGGGTTTACCCGCGAATGCGTAATGGAGTTTGCCGACACATTCGCGGGGAAACCCGCTCCTACACAg                                                                                                                                <  1:17091231/136‑1 (MQ=25)
                 cGTTCGACACCGGGCGAAGGATTTGGTGTGTATGCTGGGGTTTCCTACATCTCCCGCTGATGGGTGGGCGTGGGTTTACCCGCGAATGCTTCATTGAGTTTGCCGACGCATTCGCGGGGAAACCCGCTCCTACACAg                                                                                                                                <  2:2936442/137‑1 (MQ=37)
                        cacCGGGCGAAGGATTTGGCGTTTATGCCGGGTTTTCCTAAATCTTCCGCTGACGGGTGGGCACGGGTTTACCCGCGAGTACGTCGTGGAATTTGCCAACACATTCGCGGGGAAACCCGCTCCTACACAGGGGCGtt                                                                                                                         <  2:7100514/137‑1 (MQ=255)
                        cacCGGGAAAAGGATTTGACGTGTATGTCGGGGTTCCCTCCAACCTCCGCGGATGGGTGGGCGAGGGTTTACCTGCGAGTGCGTCATGGAGTTTGCCGACACGTTCGCGGGCAAACCCGCTCCTACACAGGGGCGtt                                                                                                                         <  2:5144104/137‑1 (MQ=255)
                          ccGGGCGAAGGGTTTGGCGTGTGTGCCGGGGTTTCCGACATCTTCCGCCGACGGGTGGGCGCGGGTTTACCCGAAAATGCGTCATTGAATTTGCCGACGCATTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGc                                                                                                                       <  1:11367073/137‑1 (MQ=37)
                           cGGGCGAATGGTTTGGCGTGTATGCCGGGGTTTCCTACATCTTCCTCTGACGGTTGGGCGCGGGTTTACCCGCGAATTCGTCATTGAATTTGCCCACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGca                                                                                                                      <  1:18067060/137‑1 (MQ=37)
                            tggCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACACCTTCCGCTGACGGTTGGGCGCGGGTTTACCCGCGAATGCTTCATGGAATTTTCCGACGCATTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGcac                                                                                                                     <  2:3024142/136‑1 (MQ=38)
                            gggCGAAGGGTTTGGCGTGTATGCTGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGTGGGTTTACCCACGAGTGCGTCATTGAATTTGCCGACGCATTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGcac                                                                                                                     <  2:6555356/137‑1 (MQ=38)
                            ggacGAAGGATTTGGCGTGTGTGCAGGGGTTTCCTACATCTTCCTTTGCCGGGTGGGTGCGGGTTTACCCGCGAATGCGTCATTGAATTTTCCAACGCATTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGcac                                                                                                                     <  2:2331624/134‑1 (MQ=25)
                                   ggATTTGGGGTGTATGCTGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑Ga                                                                                                         <  2:18406191/137‑1 (MQ=255)
                                   ggATTTGGCGTGGATGCCGGGGTTTCCTACAGCTTCCGCTGACGGTTGGGCGCGGGTTTACCCGCGAATGCGTCATGGAATTTGCCGACGCATTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑Ga                                                                                                         <  2:4044113/137‑1 (MQ=255)
                                   ggATTTGGCGGGTATGCTGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCGTTGAATTGGCCGACGCATTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑Ga                                                                                                         <  2:15384635/137‑1 (MQ=38)
                                    gATTTGTCGTGGATGCGGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGGGGGTTCACCCGCGAGTTCGTCGTGGATTTTTCCGACGCATTCGGGGGCAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GAt                                                                                                        <  1:1044621/137‑1 (MQ=255)
                                    gATTTGGCGTGTATGCCGGGTTTTCCTACAGCTTCCGCTGACGGGGGGTCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GAt                                                                                                        <  2:8228663/137‑1 (MQ=255)
                                           cGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGGTTCCCCCCAAACCCACCCTCGAATTTTCCAACACCCTCTCGGGTAAACCCACTCCCCCACAAAGGCGTTTTACCA‑CAA‑‑‑‑GATTCggggg                                                                                                 >  2:18712150/1‑136 (MQ=255)
                                             tgtgtgCCGGGGTTTCCTACATCTTCCGCGGATGGGGGGGCGTGGGTTTACCCGTGAGTGTGTCGTGGAATTTGCCGACGCATTCGCGGGCAACCCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTc                                                                                               <  1:16650052/137‑1 (MQ=255)
                                             tgtAGGCAGGGGTTCCCTACACTTTCCGCTGATGGGTGGGCGCGGGTTTACCCGCGAGTGCGTCGTGGAGTTTGCCGATGCATTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTc                                                                                               <  1:5070809/137‑1 (MQ=255)
                                                gtGCTGGGGTTTCCTACATCTTCCGCGGATGGGGGGGCGGGGGTTTACCCGTGAGGGCGTCATGGAGTTTACCGACGCATTCGAGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGa                                                                                            <  1:5391034/136‑1 (MQ=255)
                                                  tccGGGGTTCCTTACATCTTCCGTGGACGGGTGGGCGCGGGTTTACCCGTGAATGCGTCATTGAGTTTGCCGACGTATTTGTGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGAcc                                                                                          <  1:15605592/136‑1 (MQ=255)
                                                  gCCGGTGTTTCCGACAGCTTCCGCTGACGGTGGGGCGCGGGTTCCCCCACAAATGCGCCATTGACTTTCCCCACGCCTTCGCGGGAAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGAcc                                                                                          <  2:6404880/137‑1 (MQ=255)
                                                  gCCGGGTTTTCCTACAGCTTCCGCTGATGGGGGGGCGTGGGTTTACCCGCGAATGCGTCGTGGAGTTTGCCGACGCATTCGTGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGAcc                                                                                          <  1:3922992/137‑1 (MQ=37)
                                                  gCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGGTTCCCCTCGAAGGCATCCTTTAACTTTCCCAAACATTCTCCGGTCACCCCCCCTCCCCCCCGTGGCGGTTCACAGACAG‑‑‑‑GATCAGGGTTTCaacc                                                                                           >  2:9342467/1‑134 (MQ=255)
                                                   ctggggTTTCCGACAGCTTCTGTTAATGGGTGGGCGCGGGTTTACCCGCGAATGCTTCGTTGAGTTTTCCGACGCATTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCg                                                                                         <  1:8456763/135‑1 (MQ=37)
                                                      gggTTCCTTACATCTCCCGCGAATGGGGGGTCTTGGGTTTACCTACGAATGTTTCGTGGAGTTTGCCAACGCGTTTGTGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGccc                                                                                      <  1:13080663/137‑1 (MQ=255)
                                                       gttttCCTACACTTTCCTCTGATGGGTGGGCGCGGTTTTACCCGCGAATGTGTCGTGGAATTTTCCGACGCGTTCGTGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCt                                                                                     <  1:6367498/135‑1 (MQ=255)
                                                       ggTTTCCTACGGCCTCCGCTGACGGGTTGTCGTGGTTTTACCCGCGAATGCGTCCTTGACTTTACCGCCGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCt                                                                                     <  2:17368107/137‑1 (MQ=37)
                                                       ggTTTCCTACAGTTTCCGCTGACGGGTGGGCGCGGGTTTACCCGAGAATGCGTCATTGAGTTTTCCGACGCATCCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCt                                                                                     <  1:12295457/137‑1 (MQ=255)
                                                       ggTTTCCTACAGCTCCCTCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCTTCATTGAATTTGCCGACGCGTTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCt                                                                                     <  2:17954084/137‑1 (MQ=38)
                                                                  ccTTCCGCTAACGGGTGGGCGCGGTTTTACCCGCGAATGCGCCATGGAATTTGCAGACGCATTCGTGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCa                                                                          <  1:8485332/136‑1 (MQ=255)
                                                                   tttCCGCGGATGGGGGGGCGGGGGTTTACCCGAGAATGCGTCGTGGAGTTTGCCGACGCGTTCGTGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCAc                                                                         <  1:2200567/136‑1 (MQ=255)
                                                                   cTTCCGCTGATGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTTCCGATGCATTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCAc                                                                         <  2:1729988/137‑1 (MQ=255)
                                                                          tGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATGGAGTTTGCCGACGCGTTCGCGGGGAAACCCGCTCCTACACAGGGGCGTCGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGAc                                                                  <  1:15851409/137‑1 (MQ=38)
                                                                          tGACGGGGGGGCGTGGGTTTACCCGCGAATATGTCATTGAATTTTCCGACGCATTCGTGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGAc                                                                  <  1:14784103/137‑1 (MQ=255)
                                                                             tggggggGCGTGGGTTTACCCGCGAATGTGTCGTGGAGTTTGCCGACGCATTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCa                                                               <  1:11964961/136‑1 (MQ=38)
                                                                             tgggTGGGCGCGGGTTTACCCGCGAGTGCGTCGGTGAATTTGCCGACGCGTTCGTGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCa                                                               <  2:646824/136‑1 (MQ=255)
                                                                               ggtgggtgGGGGTTTACCCGCGAGTGTGTCGTGGAGTTTGCCGACGCGTTTGCGGGGAACCCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCaaa                                                             <  1:3296426/137‑1 (MQ=255)
                                                                                    tcgcgGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGt                                                        >  1:10561996/2‑137 (MQ=255)
                                                                                    tcgcgGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGt                                                        >  1:7538725/2‑137 (MQ=255)
                                                                                       cgGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGa                                                     >  1:4461587/1‑137 (MQ=255)
                                                                                        gggTTTACCCGCGAGTGCGTCGTGGAGTTTGCAGACGCATTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGAc                                                    <  1:4194827/137‑1 (MQ=255)
                                                                                         ggTTTACCCGCGAATGCTCCATTGAATTTGCCGACACATTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACa                                                   <  2:6154905/137‑1 (MQ=255)
                                                                                             tACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACg                                               >  2:9188612/1‑137 (MQ=255)
                                                                                                ccGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCCACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCt                                            >  1:6525215/1‑137 (MQ=255)
                                                                                                 cgcgAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCT‑‑‑‑g                                       >  2:3601573/1‑137 (MQ=255)
                                                                                                  gagAATGCGTCGTTGAGTTTGCCGACGCATTCGCGGGGAACCCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCT‑‑‑‑Gc                                      <  1:10698588/135‑1 (MQ=255)
                                                                                                   cgAATGCGTCATGGAGTTTGCCGACACATTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCT‑‑‑‑GCa                                     <  1:11164423/137‑1 (MQ=255)
                                                                                                      aTGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCT‑‑‑‑GCAAGa                                  >  1:5934262/1‑137 (MQ=255)
                                                                                                      aTGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCCGTCGACAAACGCCT‑‑‑‑GCAAGa                                  >  2:6651636/1‑137 (MQ=255)
                                                                                                       tGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCT‑‑‑‑GCAAgca                                 >  1:4811959/1‑135 (MQ=255)
                                                                                                       tGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCT‑‑‑‑GCAAGaa                                 >  1:16897297/1‑137 (MQ=255)
                                                                                                       tGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCT‑‑‑‑GCAAGaa                                 >  2:16600660/1‑137 (MQ=255)
                                                                                                       tGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCT‑‑‑‑GCAAGaa                                 >  1:18528420/1‑137 (MQ=255)
                                                                                                       tGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCT‑‑‑‑GCAAGaa                                 >  1:5670412/1‑137 (MQ=255)
                                                                                                       tGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCT‑‑‑‑GCAAGaa                                 >  1:668098/1‑137 (MQ=255)
                                                                                                       tGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCT‑‑‑‑GCAAGaa                                 >  1:9351370/1‑137 (MQ=255)
                                                                                                       tGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCACTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCT‑‑‑‑GCAAGaa                                 >  2:3812360/1‑137 (MQ=255)
                                                                                                        gtgTCGTTGACTTTGCCGACGCATTCGAGGGCAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCT‑‑‑‑GCAAGAAc                                <  1:7680368/135‑1 (MQ=255)
                                                                                                        gCGTCGTTGAGTTTGCCGACGCATTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCT‑‑‑‑GCAAGAAc                                <  1:15283406/137‑1 (MQ=255)
                                                                                                        gCGTCGTTGAGTTTGCAGAAGCGTTAGCGGGAAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCT‑‑‑‑GCAAGAAc                                <  1:4942020/137‑1 (MQ=255)
                                                                                                        gCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCT‑‑‑‑GCAAGAAc                                >  1:1238650/1‑137 (MQ=255)
                                                                                                            cATTGAATTTGCCGACGCATTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCT‑‑‑‑GCAAGAACGGTc                            <  2:13195995/137‑1 (MQ=255)
                                                                                                              ttGAATTTGCCGACGAATTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCT‑‑‑‑GCAAGAACGGTcgc                          <  1:1966350/137‑1 (MQ=255)
                                                                                                               tGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAGGATTGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCCTGACGCAAACCTGTCGACAAACGCCT‑‑‑‑GCAAGAACGGt                             >  2:13289900/1‑137 (MQ=38)
                                                                                                               tGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAGGATTGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAACCGCCT‑‑‑‑GCAAGCACGGt                             >  2:13289869/1‑137 (MQ=38)
                                                                                                               tGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCT‑‑‑‑GCAAGAACGGTcgcg                         <  2:13957107/137‑1 (MQ=255)
                                                                                                               ggAATTTGCCGACGCGTTCGTGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCT‑‑‑‑GCAAGAACGGTcgcg                         <  1:17455894/136‑1 (MQ=255)
                                                                                                                 aaTTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCT‑‑‑‑GCAAGAACGGTCGCGTc                       >  1:1243229/1‑137 (MQ=255)
                                                                                                                    ttGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCT‑‑‑‑GCAAGAACGGTCGCGTCTTg                    <  1:8920989/137‑1 (MQ=255)
                                                                                                                    ttGCCGACGCATTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCT‑‑‑‑GCAAGAACGGTCGCGTCTTg                    <  1:12607877/137‑1 (MQ=255)
                                                                                                                    ttGCCGACACGTTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCT‑‑‑‑GCAAGAACGGTCGCGTCTTg                    <  2:10790171/137‑1 (MQ=255)
                                                                                                                         gacacgttCGCGGGCAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCT‑‑‑‑GCAAGAACGGTCGCGTCTTGTCAtt               <  1:17222285/131‑1 (MQ=255)
                                                                                                                         gACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGCAAGAACGGTCGCGTCTTGt                   >  1:17741696/1‑137 (MQ=255)
                                                                                                                         gACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCT‑‑‑‑GCAAGAACGGTCGCGTCTTGTCAtt               >  2:11908658/1‑137 (MQ=255)
                                                                                                                         gACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCT‑‑‑‑GCAAGAACGGTCGCGTCTTGTCAtt               >  1:15413/1‑137 (MQ=255)
                                                                                                                         gACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCT‑‑‑‑GCAAGAACGGTCGCGTCTTGTCAtt               >  2:10129380/1‑137 (MQ=255)
                                                                                                                         gACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCACTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCT‑‑‑‑GCAAGAACGGTCGCGTCTTGTCAtt               >  2:6832257/1‑137 (MQ=255)
                                                                                                                          aCGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCT‑‑‑‑GCAAGAACGGTCGCGTCTTGTCATTg              <  2:16448543/137‑1 (MQ=255)
                                                                                                                              aTTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCACCTTCACCCATGACGCAAACCTGTCGACAAACGCCT‑‑‑‑GCAAGAACGGTCGCGTCTTGTCATTGAGCt          >  1:3616908/1‑137 (MQ=255)
                                                                                                                              aTTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCT‑‑‑‑GCAAGAACGGTCGCGTCTTGTCATTGAGCt          >  2:9596416/1‑137 (MQ=255)
                                                                                                                              aTTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCT‑‑‑‑GCAAGAACGGTCGCGTCTTGTCATTGAGCt          >  2:12582202/1‑137 (MQ=255)
                                                                                                                              aTTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCT‑‑‑‑GCAAGAACGGTCGCGTCTTGTCATTGAGCt          >  2:9637681/1‑137 (MQ=255)
                                                                                                                              aTTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCT‑‑‑‑GCAAGAACGGTCGCGTCTTGTCATTGAGCt          >  1:13303223/1‑137 (MQ=255)
                                                                                                                              aTTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCT‑‑‑‑GCAAGAACGGTCGCGTCTTGTCATTGAGCt          >  2:16403310/1‑137 (MQ=255)
                                                                                                                              aTTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCT‑‑‑‑GCAAGAACGGTCGCGTCTTGTCATTGAGCt          >  1:3873752/1‑137 (MQ=255)
                                                                                                                                 cgcgGTTAAACCCACTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCT‑‑‑‑GCAAGAACGGTCGCGTCTTGTCATTGAGCTTGc       <  2:7994564/137‑1 (MQ=255)
                                                                                                                                 cgcgGGTAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCT‑‑‑‑GCAAGAACGGTCGCGTCTTGTCATTGAGCTTGc       <  2:5797057/137‑1 (MQ=255)
                                                                                                                                 agcgGGCAAACCCGCTCCTACACAGGGGCGGTGCACAG‑CAG‑‑‑‑GATTGGGGGTTCAGACCGCCCTGTTCAACTTCAAACAGGAAGCAGGCCTGTCGACAGACGCCG‑‑‑‑GCAGGAACGGTCGTGTCTGGTCATTGAGCTTGc       <  2:5568763/136‑1 (MQ=255)
                                                                                                                                      gTAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCT‑‑‑‑GCAAGAACGGTCGCGTCTTGTCATTGAGCTTGCCCGaa  <  1:18696857/137‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                                
CCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAG‑CAG‑‑‑‑GATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCT‑‑‑‑GCAAGAACGGTCGCGTCTTGTCATTGAGCTTGCCCGAA  >  AE015451_phaM1/4675225‑4675495

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