Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451_phaM14,675,3600TG46.5% 150.8 / 114.2 283intergenic (+59/+45)PP_4137 CDS (putative Outer membrane siderophore receptor)/chrR CDS (Chromate reductase)–/–
Reads supporting (aligned to +/- strand):  ref base T (74/71);  new base G (3/123);  total (77/206)
Fisher's exact test for biased strand distribution p-value = 1.73e-21
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 287 total aligned reads displayed.

GCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCATTGAGCTTGC  >  AE015451_phaM1/4675224‑4675490
                                                                                                                                        |                                                                                                                                  
  tggCCAGGACAATGCCCGTTCGACACCGGGAGAAGGGTTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGTGGGTTTACCCGCGAATGCTTACTGGAATTTTCCGACGCATTCGCGGGTaa                                                                                                                                  <  2:8548347/136‑1 (MQ=38)
  cGGCCGGAACAATCCCTGTTCGACACCGGGCAAAGGATTGGGCGGGCATTCTGGGGTTTCCTACATCTTCCGCTGACGGGTGGTCGCGGGTTTACCCGCGAATCCTTCGTTGAATTTGCCGACGCATTCGCGGGCaa                                                                                                                                  <  2:3151216/137‑1 (MQ=25)
  cGGCCAGTCCAGTTCCCGTTTGCCACGTGGCGAGGGATTTGGCGTGGATTCTGGGTTTTCCTCCATTTTCCGCTGATGGGTGGGCGCGGGTTTACCCGCGAGTGCTCCATTGAGTTTCCCAACACATTCACGGTTaa                                                                                                                                  <  2:16785427/137‑1 (MQ=255)
  cGGCCAGGACAATTCCCTTTCGCCACCGTACGAAGGATTTGGCGTGTATGCAGGGGTTTCCTACACCTTCCGTTGACGGGTGGGCGCGGGTTTCCCCGCGAATGTGTCATTGAATTTGCCAACGCATTCGCGGGTaa                                                                                                                                  <  2:19282936/137‑1 (MQ=255)
  cGGCCAGGACAATGCCTGTTCGACACTGGGCGAAGGATTTGTCTTGTATGTCGGGTTTTCCTACGCCTTCCGTTGACGGGTGGGCGCGGGTTTCCCCGCGAGTGCGTCGTTGAGTTTGCCGACGCATTCGCGGGGaa                                                                                                                                  <  2:18103446/137‑1 (MQ=255)
    gCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGTGTGCATGCCGGGGTTTCCTACAGCTTCCTCTGATGGGTGGGCGCGGGTTTACCTGCGAATGCTTCATTGAATTTGCCAACGCATTCGCGGGTAAAc                                                                                                                                <  2:16349823/137‑1 (MQ=38)
    accAGGCCAGTGCTCTTTCAACACCGGGCGAAGGATTTGGCGTGTATCCTGGGGTTTCCCACAGCTCCCACTGACGGTTGGGCGCGGTTTTACCCGCGAATGCGTCGTTGAGTTTCCCGAAACGTTCGCGGGGAAAc                                                                                                                                <  2:2416601/136‑1 (MQ=255)
       agaaCAATGCCCGTTCGACCCCGTGATAAGGATTTGGCGTGGAGTCTGGGGTTTCCTACGGCTTCCGTTGACGGTTGGGCGTGGGTTTACCCACGAAGGCGTCATGGAATTTGCCGACGCATTCGCGGGTAAACCCg                                                                                                                             <  2:12990453/134‑1 (MQ=25)
       agaaCAATCCCCGTTCGACACGGGGCGAGGGAGTGGGGGTGTGTGCTGGGGTTTCCCACATCTTCCGCGGATGGGGGGGCGCGGGTTTACCCGCGAGTGTTTCATGGAGTTTGCCGACGCCTTTGAGGGTAAACCCg                                                                                                                             <  2:5949352/134‑1 (MQ=255)
       aGGACAATTCCTGTTCGACACCGGGAGAAGGATTTGGCGTGGAGGCGGGGGTTTCCTACATCTTCCGTGGATGGGGGGGCGGGGGTTTACCCGCGAATGCGTCGTGGAGTTTGCCGACACGTTCGCGGGGAAACCCg                                                                                                                             <  2:8261362/137‑1 (MQ=255)
            aaTACCCGCTCAACACCGGGCAAAGGATTTGGCGTGTATGCCGGGTTTTCCTACATCTTCCGCTGACGGGTGGGCGGGGGTTTACCCGCGAATGCGTCATTGAATTTTCCGACGCATTCGCGGGGAAACCCGCTCCt                                                                                                                        <  2:228114/137‑1 (MQ=37)
               tcccGTTCGACACTGGGCGAAGGATTTGGCGTTTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACACATTCGCGGGTAAACCCGCTCCTaca                                                                                                                     <  2:13745544/136‑1 (MQ=255)
               ccccGTCCGACACCGGTCGAAGGGTTTGTCTTGTATGCCGGGGTTTCTTACATCTTCCGCGGATGTGTGTGCGCGGGTTTACCCGCAAATGCGTCATGGAATTTGCCGACGCATTCGTGGGTAAACCCGCTCCTaca                                                                                                                     <  2:9462169/136‑1 (MQ=25)
                 ccGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTTCCCCCGAACCCGTCACTGAAGTTTCCAAAGCACCCGCGAGGAAAACCACCCCTaccaa                                                                                                                   >  2:14476950/1‑134 (MQ=255)
                 ccGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGGTTCCCCGCGAACGCGTCATTGAAGTTGCCGACGAATTCGCGGGTAAACCCACTCCTCcaca                                                                                                                   >  2:2929003/1‑137 (MQ=255)
                  cGTTCGACACTGGGCGAGGGGTTGGGTGTGGGTGCTGGGTTTTCCTACAGCTCCCGTTGATGGGTGGGTGTGGGTTTACCCGCGAGTGTTTCATTGAATTTGCCGACGCATTCGCGGGGAAACCCGCTCCTACACAg                                                                                                                  <  2:19398455/137‑1 (MQ=19)
                  cGTTCGACACTGGGCAAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAg                                                                                                                  <  2:11511323/137‑1 (MQ=255)
                    ttCGACACCGGGCGAAGGATTTGGTGTGTATGCCGGGGTTTCCTAAAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGTGTCGTTGAGTTTGCCGACGCGTTTGCGGGGAAACCCGCTCCTACACAggg                                                                                                                <  2:15092077/137‑1 (MQ=38)
                      cGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGGTTGCCCCCGAACGCGGCAGTGGATTTTCCAACGCACTCTCGGGGGAAACCACTCCCCCACAagggg                                                                                                              >  2:6869868/1‑136 (MQ=255)
                            tgggCGAAGGGTTTGGCGTGTATGCCGGGGTTTCCGACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATGGAATTGGCCGACGCATTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGca                                                                                                        <  1:19052368/136‑1 (MQ=38)
                             gggCGAAGGATTTGGCGTGTATGCTGGGGTTTCCTACAGCTTCCGTTGATGGGTGGGCGTGGGTTTACCCGCGAATGCGTCGTTGAATTTGCCGACGCGTTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGcac                                                                                                       <  1:19378836/137‑1 (MQ=255)
                             gggAGAAGGGTTTGGCGGGTATGCCGGGGTTCCCTACAGCTTCCGCGGACGGGTGGGAGGGGGTTGACCCGAGAAGGCGTAAGGGAATTTGCCGACGCGTTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGcac                                                                                                       <  1:12386936/137‑1 (MQ=255)
                                    gggtttGGCGTGTATGCCGGGGTTTCCTACATTTTCCGTGGAGGGGGGGGCGCGGGTTGACCAGAGAGGGCGTAGTGGAGTTTGCCGAAGCATTAGTGGGAAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGa                                                                                                <  1:16098545/134‑1 (MQ=255)
                                     ggttgggCGTGGGTGCCGGGGTTTCCTCCAGCTTCCGCTGATGGGGGGGCGTGGTTTGACCTGCGAAGGTTTAATGGAATTGGCCAACGCATTCGAGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGAt                                                                                               <  2:4072216/132‑1 (MQ=255)
                                            cGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGGTTCCCCCCAAACGCTTCACCAAACTCTCCCACACACTCCCAAGTAAACCCACCCCCCCCCAAGGGAGTTGTACAACAGGAGTAggggg                                                                                        >  2:1558764/1‑136 (MQ=255)
                                              tgtATGCTGGGTTTTCCTCCATTTTCCGCGGATGGGTGGGCGTGGTTTTACCCGCGAGTGCGTCATGGAGTTTGCCGACGCCTTCGTGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTc                                                                                      <  1:7647860/137‑1 (MQ=255)
                                              tggaTGCTGGGTTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCAACGAATGCTTCGTGGAATTTCCCGACGCGTTTGTGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTc                                                                                      <  1:18131147/134‑1 (MQ=255)
                                                 gtGCTGGGTTTTCCTACATTTTCTGTGGAGGGGGGGGGGGGGGTTTACCCGGGAGTGTTTCAGGGAGTTTTCCGACGCGTTCGCGGGAAACCCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGa                                                                                   <  1:6757770/136‑1 (MQ=255)
                                                 gggcaggggTTTCCCACAGCTTCCGTTGACGGGGGGGCGCGGGTTTACCCGCGAATACGTCATTGAATTGGCCGACGCGTTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTATGGGTTCAGa                                                                                   <  2:16193761/132‑1 (MQ=38)
                                                 aTGCCGGGTTTTCCTACAGCTTCCGCTGATGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAGTTTGCCGACGCATTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGa                                                                                   <  2:11985019/137‑1 (MQ=255)
                                                 aTGCCGGGGTTCCCTACAGCTTCCTCTGAGGGGTGGGCGTGGGTTTACCCGCGAATGCGTCGTGGAGTTTGCCGACTCATTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGa                                                                                   <  1:6901389/137‑1 (MQ=37)
                                                   gctggggTTTCCTCCGTCTCCCGTTGGTGGGGGGGCGTGGGTTTACCAGCGGGTGTGTGGTGGAGTGTGCTGACGCGTTTGGGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGAcc                                                                                 <  1:9518241/134‑1 (MQ=255)
                                                   gCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAGTTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGAcc                                                                                 <  2:18252133/137‑1 (MQ=255)
                                                   gCCGGGGTTCCCTACATCTCCCGCTGACGGGTGGTCGCGGGTTTACCCGCGAATGTGTCGTGGAATTTTCCAACGCATTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGAcc                                                                                 <  2:750032/137‑1 (MQ=255)
                                                   gCCGGGGTTCCCTACACCTTCCATTGACGGGTGGGCGCGGGTTTACCCGCGAGTGCGTCATGGAATTTGCCGACGCATTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGAcc                                                                                 <  2:7279949/137‑1 (MQ=38)
                                                   cccGGGTTTTCCGACAGCTTCCGCTGATGGGTGGGCGTGGGTTTACCCGTGAATGCGTTGTTGAGTTTGTCGACGCGTTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGAcc                                                                                 <  1:14436412/136‑1 (MQ=255)
                                                              cTACATCTTCCGTGGACGGGTGGGTGTGGTTTTCCCCGCGAATGTGTCGTTGAGTTTGCGGACGCGTTTGTGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCaa                                                                      <  2:1881562/137‑1 (MQ=255)
                                                              cTACAGCTTCCGCTGATGGGGGGGCGTGGGTTTACCCGCGAATGCGTCATGGAATTTTCCGACGCATTCGTGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCaa                                                                      <  1:13741126/137‑1 (MQ=255)
                                                                    cTTCCGCTGATGGGTGGGCGCGGGTTTACCCGAGAGTGTGTCGGGGAGTTTGCCGACGAATTCGGGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCAc                                                                <  1:8103683/137‑1 (MQ=255)
                                                                    cTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCGTTGAGTTTGCCAACGCATTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCAc                                                                <  2:11054139/137‑1 (MQ=255)
                                                                           tggcGGGGGGGAGTGGGTTGACCTGAGAATGCGTCGTGGAATTTGCCGACGCATTCGTGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGAc                                                         <  1:1777530/134‑1 (MQ=255)
                                                                           tGACGGGTGGGCGGGGGTTGACCGGTGAGTGCGTCGTTGAATTTGCTGACGCGTTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGAc                                                         <  2:4922602/137‑1 (MQ=255)
                                                                           tGACGGGTGGGCGCGGGTTTACCCACGAATGCGCCATGGAGTTTTCCACCGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGACCTGTTCAACTTCACCCATGAc                                                         <  1:3733561/137‑1 (MQ=38)
                                                                              gggggggGCGGGGGTTTACCCGCGAGTGTGTCGTTGAGTTTGCCGACGCGTTCGTGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCa                                                      <  1:15983058/136‑1 (MQ=255)
                                                                              cGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCGTTCGCGGGTAAACCCGCTCCTACATAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCa                                                      <  2:4534100/137‑1 (MQ=255)
                                                                              cGGGGGGGCGAGGGTTTACCCGCGAATGCGTCATTGACTTTTCCGACGCGTTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCa                                                      <  1:18840119/137‑1 (MQ=38)
                                                                                 gccttcgcgGGTTCACCCACGAATGCGTCATCGAACCCGCCGACGCATCCACGAGTACACCCGCTCCCACACCGGGGCGTTGCACAGCAAGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCAGGACGCaaaa                                                   >  2:11673716/6‑136 (MQ=255)
                                                                                    ggcacgGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCCCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCCGACCGCCCTGTTCACCTTCCCCCATGACGCAAACCTg                                                >  1:13929709/5‑137 (MQ=38)
                                                                                      cgcgGGTTTACCCGCGAATGCGTAATGGAGTTGGCCGACGCATTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTc                                              <  2:17374717/137‑1 (MQ=39)
                                                                                           ttttACCAGCGAGTGCGTCATGGAATTTGCCGACGCGTTCGTGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACaa                                         <  1:12847815/136‑1 (MQ=255)
                                                                                             ttACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACCTCACCCATGACGCAAACCTGTCGAAAAAc                                       >  2:572428/1‑137 (MQ=255)
                                                                                                cccGCGAATGCGTCATTGAATTTGCCGACGCATTCTCGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGcc                                    <  2:5878734/137‑1 (MQ=255)
                                                                                                 ccGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCt                                   >  2:3401805/1‑137 (MQ=255)
                                                                                                     gagttcgtcgtTGAGTTTGCCGACGCATTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCaa                               <  1:11099776/132‑1 (MQ=255)
                                                                                                     gAATGCGTCATTGAATTTGCCGACGCGTTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAACAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCaa                               <  2:15392737/137‑1 (MQ=255)
                                                                                                        tGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGaa                            >  1:10256439/1‑137 (MQ=255)
                                                                                                        tGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGaa                            >  1:16670669/1‑137 (MQ=255)
                                                                                                        tGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGaa                            >  1:13571418/1‑137 (MQ=255)
                                                                                                        tGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGaa                            >  2:17762197/1‑137 (MQ=255)
                                                                                                        tGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGaa                            >  1:3534147/1‑137 (MQ=255)
                                                                                                        tGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCATACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGaa                            >  1:4477176/1‑137 (MQ=255)
                                                                                                         gCGTCATTGAGTTTGCCAACGCATTAGGGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAAc                           <  1:8777832/137‑1 (MQ=255)
                                                                                                         gCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAAc                           <  1:13778568/137‑1 (MQ=255)
                                                                                                         gCGTCATGGAATTTGCCGACGCATTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAAc                           <  2:4401918/137‑1 (MQ=255)
                                                                                                         gCGTAATGGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAAc                           <  2:17046507/137‑1 (MQ=255)
                                                                                                            tCATTGAATTTGCCGACGCGTTCGTGGGTAAGCCCGCTCCTACACAGGGGCGGTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAACCCTGTCGACAAACGCCTGCAAGAACGGt                        <  2:15685089/137‑1 (MQ=255)
                                                                                                             cATTGGATTTGCCGACGCATTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTc                       <  2:14310719/137‑1 (MQ=255)
                                                                                                             cATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTc                       >  1:2989991/1‑137 (MQ=255)
                                                                                                             cATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTc                       >  1:17477066/1‑137 (MQ=255)
                                                                                                             cATTGAATTTCCAGACGAATTCGCGGGCAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTc                       <  1:16870349/137‑1 (MQ=255)
                                                                                                              aTTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTcg                      >  1:3984098/1‑137 (MQ=255)
                                                                                                               ttGAGTTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTcgc                     <  2:17633155/137‑1 (MQ=255)
                                                                                                               ttGACTTTGCCGGCTCATTTGCGGGTAAACCCGCCCCTACACAGGGGCGTGGCACAGCAGGATTAGGGGTTCAGACCGCCCTATTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTcgc                     <  2:762762/137‑1 (MQ=255)
                                                                                                               ttGAATTTGCCGACGCGTTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTcgc                     <  1:13116380/137‑1 (MQ=255)
                                                                                                               ttGAATTGGCCGACGCGTTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTcgc                     <  2:10245977/137‑1 (MQ=255)
                                                                                                               tggagtttGCCGACGCGTTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTcgc                     <  1:11659746/132‑1 (MQ=255)
                                                                                                               tggAATTTGCCGACGCGTTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTcgc                     <  2:62304/135‑1 (MQ=255)
                                                                                                                  aaTTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTTGACAAACGCCTGCAAGAACGGTCGCGTc                  >  1:9719505/1‑137 (MQ=255)
                                                                                                                  aaTTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTc                  >  2:13511047/1‑137 (MQ=255)
                                                                                                                  aaTTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACATGTCGACAAACGCATGCAAGAACGGTCGCGTc                  >  2:12408108/1‑137 (MQ=255)
                                                                                                                     tttccGACGCGTTCGGGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTg               <  1:5099468/134‑1 (MQ=255)
                                                                                                                     ttGCCGACGCGTTTGTGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTg               <  1:2289979/137‑1 (MQ=255)
                                                                                                                     ttGCCGACGCGTTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTg               <  1:465650/137‑1 (MQ=255)
                                                                                                                     ttGCCGACGCATTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTg               <  2:9240728/137‑1 (MQ=255)
                                                                                                                      ttccGACGCATTTGCGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGt              <  1:19491431/135‑1 (MQ=255)
                                                                                                                         cGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCAt           >  1:14922283/1‑137 (MQ=255)
                                                                                                                          gACGCGTTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCAtt          <  2:16582973/137‑1 (MQ=255)
                                                                                                                          gACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCAtt          >  2:11344501/1‑137 (MQ=255)
                                                                                                                          gACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCAtt          >  1:17113936/1‑137 (MQ=255)
                                                                                                                          gACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCAtt          >  1:13021621/1‑137 (MQ=255)
                                                                                                                          gACGCATTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCAtt          <  2:3630700/137‑1 (MQ=255)
                                                                                                                          gACGCATTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTCCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCAtt          <  2:18743143/137‑1 (MQ=255)
                                                                                                                           aCGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCATTg         <  2:7080315/137‑1 (MQ=255)
                                                                                                                               aTTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCATTGAGCt     >  1:2674089/1‑137 (MQ=255)
                                                                                                                               aTTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCATTGAGCt     >  1:6543319/1‑137 (MQ=255)
                                                                                                                               aTTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCATTGAGCt     >  1:15701217/1‑137 (MQ=255)
                                                                                                                               aTTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCATTGAGCt     >  1:15387798/1‑137 (MQ=255)
                                                                                                                                  tgcgGGGAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCGGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCATTGAGCTTGc  <  2:2727450/136‑1 (MQ=255)
                                                                                                                                  cggggggAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCATTGAGCTTGc  <  1:5624452/134‑1 (MQ=255)
                                                                                                                                     gggTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCATTGa        <  2:11200312/128‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                  
GCCGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCATTGAGCTTGC  >  AE015451_phaM1/4675224‑4675490

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