Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451_phaM14,675,3600TG44.2% 157.4 / 68.4 204intergenic (+59/+45)PP_4137 CDS (putative Outer membrane siderophore receptor)/chrR CDS (Chromate reductase)–/–
Reads supporting (aligned to +/- strand):  ref base T (61/45);  new base G (1/83);  total (63/141)
Fisher's exact test for biased strand distribution p-value = 2.34e-19
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 205 total aligned reads displayed.

CGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCATTGAGCTTGCCCGAA  >  AE015451_phaM1/4675226‑4675495
                                                                                                                                      |                                                                                                                                       
tgtccAGGACAGTGCCTTTTCGACCCTGGGCGAAGGATTTGGTGGGTATGCCGGGGTTTCCGACAGCTTCCGCTGACGGGGGGGCGCGGGTTTACCAGCGAGTGCGTCGTTGAATTTGCCGACGCGTTCGCGGGGaa                                                                                                                                       <  2:8447473/134‑1 (MQ=255)
tggCCAGGACAATGCCCGTTCGACACTGGGTGAAGGGTTTGGTGGGTATGCTGGGGTTTCCTACGGCTTCCGCTGACGGGGGGGCGCGGGTTTACCCGCGAATGTGTCGTGGAATTTTCCGCCGCATTCGCGGGGaa                                                                                                                                       <  2:1301941/136‑1 (MQ=255)
cgacagggACAGGGCCAGTTCGACACTGGGCGAAGAGTTGGGCGTGGGGGCTGGGTTTTCCTACACCTTCCGCGGATGGGTGGGCGAGGGTTTACCCGCGAATGCGTCATGGACTTTCCCGACCCGTCCGCGGGTaa                                                                                                                                       <  2:7259586/131‑1 (MQ=255)
cGGCCGGGACAATTCCTGTTCGACCCGGGGCGAAGGGTTGGGCGTGTATGCTGGGGTTTCCTACAGCTTCCGTTGATGGGTGGGTGGGGGTTTACCCACGAATTTTTCATTGAATTTGCCGACGCATTCGTGGGTaa                                                                                                                                       <  2:10446729/137‑1 (MQ=255)
cGGCCGGGACAATGCCTGTTCGACACCGGGCGAAGGGTTTGGCGGGTGTGCTGGGGTTTCCTACAGCTTCCGCTGATGGGTGGGCGTGGTTTTACCCGCGAATGCGTCGTTGAGTTTGCCGACGCGTTCGTGGGGaa                                                                                                                                       <  2:4883902/137‑1 (MQ=255)
cGGCCAGGACAGTTCCTGTTCGACACTGGGCGAAGGGTTTGGCGTGTATGCTGGGGTTCCCTACAGCTTCCGCTGACGGGTGGGTGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCTGACGCGTTCGGGGGGaa                                                                                                                                       <  2:5255185/137‑1 (MQ=255)
cGGCCAGGACAGTACCTGTTCGACACCGGGCGAAGGATTGGGCGGGGACGCTGGGGTTCCTTACAACTCCCGCGGATGGGTGGGCGTGGGTTTACCCGAGAATGCTTCGTGGAGTTTTCCGACTCATCCGCGGGGaa                                                                                                                                       <  2:3926938/137‑1 (MQ=255)
cGGCCAGGACAATGCCTGTTCGACACCGGGTGAAGGTTTTGGCGTGGGTGCTGGGGTTTCTTACGTTTTCCTCTGATGGGTGGGCGCGGGTTTACCCGCGAGTGCTTCATTAAATTTTCCAACGCGTTCGCGGGGaa                                                                                                                                       <  2:11797535/137‑1 (MQ=255)
cGGCCAGGACAATCCCCGTTCGACACCGGGAGAAGGGTTTGGTGGGGGGGCTGGGGTTTCCCACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATTCGTCGTGGAATTTGCCGACGCATTCGTGGGGaa                                                                                                                                       <  2:1249724/137‑1 (MQ=255)
   ccAGAACAATGCCTGTTCGACACCGGGCGAAGGGTTGGGCGTGGATGCTGGGGTTTCCTACAGCTTCCGCGGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGGAAAcc                                                                                                                                    <  2:1860558/137‑1 (MQ=37)
           aTGCCCGTTCGCCACCTGCCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACATCTTCCGCTGCCGTGTGGGCGCGGGTTTACCCGCGAATACGTCATGAAATTTACCAACGCATTCGCGGGTAAACCCGCTCCTa                                                                                                                            <  2:9711336/137‑1 (MQ=37)
               ccGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGGTTGCCCGCTCCCGCGTTCTTTGATTTTCCTACGCATTTTCGGGTAAACCCCCTccccccca                                                                                                                        >  1:16982898/1‑131 (MQ=255)
                     gACCCCTGGCGAGGGATTTGGCTTTTATTCCGGGGTTTCCTACAGCTTCCTCTGATGGGTGGGAACGGGTTTACCAGCGAATGCGTCATTGAATTTGCCGACCCATTCGCGGGGAAACCCGCTCCTACACAGGGGCg                                                                                                                  <  2:2314463/137‑1 (MQ=255)
                     gACACCGGGAGAAGGATTTGGTGTGTGTGCTGGGTTTTCCTACAGCTTCCGCTGATGGGTGGGCGCGGGTTTACCCGCGAGTGCGTCGTTGAGTTGGCCGACGCGTTCGCGGGGAAACCCGCTCCTACACAGGGGCg                                                                                                                  <  2:256882/137‑1 (MQ=255)
                          cGGGCGAAGGATTTGTCGTGGAGGCCGGGGTTTCCTACAGCTTCTGCGGATGGGTGGGCGTGGGTTTACCCGCGAATGTTTCATGGAGTTTGCCGACGCATTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGca                                                                                                             <  1:11774165/137‑1 (MQ=37)
                           ggtcGAAGGATTTGGTGTGGATTCCGGGGTTCCCTACACCTTCCGCGAACGGGGGTGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCGTTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGcac                                                                                                            <  2:8801252/134‑1 (MQ=25)
                           gggCGAAGGATTTGGCGTGTGTGCTGGGTTTTCCTACAGCTTCCGCTGACGGGTGGACGCGGGTTTACCCACGAATACGTCATCGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGcac                                                                                                            <  2:13876786/137‑1 (MQ=38)
                           gggCGAAGGATTTGGCGTGTGTGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATGGAATTTGCCGACGCATTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGcac                                                                                                            <  1:7575283/137‑1 (MQ=255)
                           gggCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTCCGGCTTCCTTGGGCGGGTGGGCGCGGGTTTACCCGCGAGGGCGTCGGGGAGTTTGCCGACGCGTTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGcac                                                                                                            <  1:10602176/137‑1 (MQ=255)
                           gggAGAAGGGTTTGGTGTGGGTGCAGGGGTTTCCTACATCTTCTGTTGATGGGTGGGCGTGGGTTTACCTGCGAGTGTGTCGTGGAGTTTTCCGACGCATTCGCGGGGAACCCCGCTCCTACACAGGGGCGTTGcac                                                                                                            <  1:5176336/137‑1 (MQ=255)
                                  gggtttGGTGTGTATGCTGGGGTTCCCTACAGCTTCCGCGGATGGGGGGGCGGGGGTTGACCCGAGAAGGCGTAGGTGAATTTGCCGACGCGTTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGa                                                                                                     <  2:1752194/134‑1 (MQ=255)
                                  ggATTGGGCGTTTATGCCGGGTTCTCCTACATCTTCCGCCAATGGTTGGGCGCGTTTTTACCCGCGGATGCTTCATTGAGTTTACCGCCGCGTTCGCGGGCAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGa                                                                                                     <  1:13288795/137‑1 (MQ=255)
                                   ggtttGGTGTGTCTGCGGGGTTTTCTCACCTCTTCCTCTGATGGGGGGGCGTGGGTTTCCCCGAGACTTCGCAATTAACTTTGCAGACGCATTCCCGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGAt                                                                                                    <  2:2434937/135‑1 (MQ=255)
                                   gATTTGGCGTGTGTTCCGGGTTTTCCTCCATCTCCCCCTGACGGGTGGGCGCGGGTTTACCCGCGAATACTTCGTTGAATTTGCCGACGCGTTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGAt                                                                                                    <  2:11285039/137‑1 (MQ=255)
                                       tGGTGTGTATACCGGGGTTTCCTACACCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGCCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAgg                                                                                                <  2:6941908/137‑1 (MQ=255)
                                            tgtATGCCGGGGTTTCTTACATCTTCTGCTGACGGGTGGGCGCGGTTTTACCCGCGAATGCTTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTc                                                                                           <  2:15524692/137‑1 (MQ=38)
                                            tgtAGGCTGGGGTTTCCTACATCTTCCGTGAAAGGGTGGGCGCGGGTTCACCCGCGAGTGTTTCGTTGAATTTTCCGACGCGTTCGCGGGAAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTc                                                                                           <  1:16649595/137‑1 (MQ=255)
                                            tggaggCCGGGGTTTCCTACATCTTCCGCGGACGGGTGGGCGCGGGTTTACCCGCAAGTGTGTCATTGAATTTCCCGCTGCATTCGCGGGCAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTc                                                                                           <  1:6485597/132‑1 (MQ=255)
                                               aTGCTGGGTTTTCCTACAGCTCCCTTTGATGGGTGGGCGCGGGTTTACCCGCAAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGa                                                                                        <  2:1734045/137‑1 (MQ=37)
                                                 tccGGGGTTTCCTCAAGCTTCCACTGACGGGTGGGCGCGGTTTTACCCGCAAATGCTTCATTGAATTTGCCGACGCATTCGCGGGGAAACCCGCTCCTACACAGGGGCGGTGCACAGCAGGATTAGGGGTTCAGAcc                                                                                      <  2:6445726/136‑1 (MQ=38)
                                                 gctggggTTTCCTACATCTTCCGCGGATGGGTGGGCGTGGGTTTACCCGCGAATGCTTCGTTGAGTTTGCCGACGCGTTCGTGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGAcc                                                                                      <  1:12479317/134‑1 (MQ=37)
                                                 gctggggTTTCCTACATCTTCCGCGGACGGGTGGGCGTGGTTTTACCCGCGAATGCTTCATTGAGTTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGAcc                                                                                      <  1:12908514/134‑1 (MQ=38)
                                                 gCCGGGGTTTCCTAAAGCTCCCTCTGATGGGGGGGCGCGGGTTGACCAGAGAATGTGTCATGGAATTTTCCGACGCGTTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGAcc                                                                                      <  1:13094559/137‑1 (MQ=255)
                                                 gCCGGGGTTCCCTACAGCTCCCGCCGACGGGGGGGCGTGGTTTTACCCACGAATTTGTCATTGAGTTTGCCGACGCATTTGCGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGGTTAGGGGTTCAGAcc                                                                                      <  1:12031982/137‑1 (MQ=255)
                                                      gttttCCAACAGCTTCCACTGACGGGTGGCTGCGGGTTGACCCAGGAGGCCGTCATGGAGTTGGCCGATGCGTTAGAGGGCAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCt                                                                                 <  1:14930534/135‑1 (MQ=255)
                                                      gtttcccTACATCCTCCTCCGACGGTTGGTCACGCTCTCACCCCCGAATCCTTCATTAACTTTTCTGCCGCGTTCGCGCGCAAACCCACCCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCt                                                                                 <  2:1076269/132‑1 (MQ=255)
                                                      ggTTTCCTACAGCTTCCGCCGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCTGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCt                                                                                 <  2:16228729/137‑1 (MQ=255)
                                                                  tttCCGTTAACGGGTGGGCGCGGGTTGACCCGCGAGTGTGTAATTGAATTTTCCGACTAATTCGAGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCAc                                                                     <  1:1828670/136‑1 (MQ=255)
                                                                  cTTCTGCTGATGGTTGGGCGTGGTTTTACCCGCGAATGCGTCGTGGAGTTTGCCGACACGTTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCAc                                                                     <  1:9297729/137‑1 (MQ=255)
                                                                  cTTCCGCTGATGGGGGGTCGTGGGTTTACTCGAGAGTGCGGCGGTGAGTTTGCCGACGCATTTGAGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCAc                                                                     <  2:16915745/137‑1 (MQ=255)
                                                                  cTTCCGCTGACGGGTGGGCGTGGGTTTACCCGCGAATGTGTCATTGAGTTTGCCGACACATTCGTGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCAc                                                                     <  1:15383925/137‑1 (MQ=255)
                                                                  cTTCCGCTGACGGGTGGGCGAGGGTTTACCCGCGAATGCTTCATGGAATTTGCCGACGCGTTCGTGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCAc                                                                     <  2:17085022/137‑1 (MQ=255)
                                                                  cTTCCCCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCAc                                                                     <  2:8066420/137‑1 (MQ=255)
                                                                     ccGCTGACGGGTGGGCGCGGGGTTTCCCCCGAACGCATCCTTGGATTTTCCGACGCACTCTCTGGTAAAACCACCCCCCCACAAAGGAGTTGCCCAACAAGAGTCAGGGTTCCAACCGCCCCCTTCCAATTCCCCCa                                                                  >  2:4147875/1‑137 (MQ=255)
                                                                         ggagggggggGTGAGGGTTTACCCGAGAAGGCGTCATGGATTTTACCGACACATTCGCGGGGAACCCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGAc                                                              <  1:6156924/133‑1 (MQ=255)
                                                                         ggaagtgTGGGCGTGGGTTCACCTGCAACTCCGTCATTGAGTTTTCCGACTCATTCGTGGTGAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGAc                                                              <  2:1766238/131‑1 (MQ=255)
                                                                            cGGGTGGGCGCGGGTTTACCCGCGAGTAAGTCATTGAATTTGCCAACACGTTCGGGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCa                                                           <  2:9021366/137‑1 (MQ=255)
                                                                            cGGGTGGGCGCGGGTTTACCCAAGAATCTGTCATTGAATTTGCCAACGCGTTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCa                                                           <  2:13489835/137‑1 (MQ=255)
                                                                              cgcattcgcgGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCaaa                                                         >  2:10783648/7‑137 (MQ=255)
                                                                                tGGGCGCGGGTTTACCCGCGAATGCGTCGTTGAGTTTTCCGACGCATTTGTGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAAcc                                                       <  2:3386346/137‑1 (MQ=38)
                                                                                     gtgggTTTACCCGCGAGTGTGTCGTTGGGTTTGCTGACGCGTTCGGGGGGAACCCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCg                                                  <  1:13700927/135‑1 (MQ=255)
                                                                                      cgGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGa                                                 >  1:8718514/1‑137 (MQ=255)
                                                                                       gggTTTACCCGCGAGTGCGTCATGGAATTTGCCGACGCATTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGAc                                                <  1:12744161/137‑1 (MQ=255)
                                                                                       gggTTTACCCGCGAATGCTTCGTTGAATTTGCCGACGCATTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGAc                                                <  1:10282566/137‑1 (MQ=255)
                                                                                       gggTTTACCCGCGAATGCGTCATTGAGTTTGCAGACGCATGCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGAc                                                <  2:7521389/137‑1 (MQ=255)
                                                                                       gggTTGACCGGAGAAGTCGTCATTGAGTTTGCCGCCGCATTCCCGGGGAACCCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAGACCTTTCGAc                                                <  2:6592139/137‑1 (MQ=255)
                                                                                        ggTTTACCCGCGAAGGCGTCATGGAATTGGCTGACGCATTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACa                                               <  2:11055600/137‑1 (MQ=255)
                                                                                                cgcgAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTg                                       >  2:14789971/1‑137 (MQ=255)
                                                                                                cgcgAATGCGTCACTGAATTTTCCGAAGCAGTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGTACAGGAGGATTCGGGGTTCCGACCGACCTGTTTAACTTCACCCATGACGAAACCCTGTCGACAAACGCCTg                                       >  2:2958415/1‑137 (MQ=255)
                                                                                                   gAATGCGTCGTTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCaa                                    <  2:153264/137‑1 (MQ=255)
                                                                                                      tGCGTTATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCCGGATTAGGGGTTCAGACCGCACTGTTCCACTTCACCCATGACGCAAACCTGTCGACAAGCGCACGCAAGaa                                 >  2:1216080/1‑137 (MQ=38)
                                                                                                      tGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACCCAGGGGCGTTGCACCGCAGGATTCGGGGTTCCGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCACGaa                                 >  1:4819397/1‑137 (MQ=39)
                                                                                                      tGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGaa                                 >  1:8411111/1‑137 (MQ=255)
                                                                                                      tGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGaa                                 >  1:2781739/1‑137 (MQ=255)
                                                                                                      tGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCATGCAAGaa                                 >  2:9563350/1‑137 (MQ=255)
                                                                                                       gCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAAc                                >  2:5802334/1‑137 (MQ=255)
                                                                                                       gCGTCATTGAATTTGCCGACACATTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAAc                                <  2:3336121/137‑1 (MQ=255)
                                                                                                         gtcgtggAATTTGCCGACGCATTCGTGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACgg                              <  1:11286231/131‑1 (MQ=255)
                                                                                                           cATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACATGTCGACAAACGCCTGCAAGAACGGTc                            >  2:9512185/1‑137 (MQ=255)
                                                                                                           cATTAAATTTGCCTACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTc                            <  2:3709873/137‑1 (MQ=255)
                                                                                                             ttGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTcgc                          <  1:12142754/137‑1 (MQ=255)
                                                                                                             ttGAATTTGCCGACGCATTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTcgc                          <  2:8469543/137‑1 (MQ=255)
                                                                                                             tggagtttGCCGACGCGTTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTcgc                          <  1:16473777/132‑1 (MQ=255)
                                                                                                             tggagtttGCCGACGCGTTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTcgc                          <  1:570866/132‑1 (MQ=255)
                                                                                                             tggAATTTGCCGACGCGTTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTcgc                          <  1:8061165/135‑1 (MQ=255)
                                                                                                              tGAATTTGCCGACGCATTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTcgcg                         <  2:16982898/137‑1 (MQ=255)
                                                                                                               gagtttCCCGACGCGTTCGTGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGt                        <  1:3409741/134‑1 (MQ=255)
                                                                                                               gAATTTTCCAACGCATTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTcgagt                        <  1:9981238/137‑4 (MQ=255)
                                                                                                               gAATTTGCCGACGCATTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGt                        <  2:14296875/137‑1 (MQ=255)
                                                                                                                aaTTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTc                       >  2:12690152/1‑137 (MQ=255)
                                                                                                                aaTTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTc                       >  2:8541314/1‑137 (MQ=255)
                                                                                                                   ttGCCGACGCGTTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTg                    <  1:11894721/137‑1 (MQ=255)
                                                                                                                    tGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGt                   >  2:6240224/1‑137 (MQ=255)
                                                                                                                     gCCGACGAATTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTc                  <  2:4191736/137‑1 (MQ=255)
                                                                                                                        gACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCAtt               >  2:13858125/1‑137 (MQ=255)
                                                                                                                        gACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCAtt               >  2:15057776/1‑137 (MQ=255)
                                                                                                                        gACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCAtt               >  1:14443339/1‑137 (MQ=255)
                                                                                                                        gACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCAtt               >  1:3085605/1‑137 (MQ=255)
                                                                                                                        gACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCAtt               >  2:1674644/1‑137 (MQ=255)
                                                                                                                        gACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCAACCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCAtt               >  2:7061996/1‑137 (MQ=255)
                                                                                                                          cGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCACTGTTCACCTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCATTGa             >  1:4060276/1‑137 (MQ=255)
                                                                                                                          cGCATTCGCGGGGAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGNATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCATTGa             <  1:14254200/137‑1 (MQ=255)
                                                                                                                             aTTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCATTGAGCt          >  1:11175612/1‑137 (MQ=255)
                                                                                                                             aTTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCATTGAGCt          >  1:11038431/1‑137 (MQ=255)
                                                                                                                             aTTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCATTGAGCt          >  1:9663194/1‑137 (MQ=255)
                                                                                                                             aTTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCATTGAGCt          >  2:9190068/1‑137 (MQ=255)
                                                                                                                             aTTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCATTGAGCt          >  1:12171606/1‑137 (MQ=255)
                                                                                                                             aTTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCATTGAGCt          >  1:544628/1‑137 (MQ=255)
                                                                                                                             aTTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCATTGAGCt          >  1:15684271/1‑137 (MQ=255)
                                                                                                                                cgcgGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTTTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCATTGAGCTTGc       <  2:15940920/137‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                       
CGGCCAGGACAATGCCCGTTCGACACCGGGCGAAGGATTTGGCGTGTATGCCGGGGTTTCCTACAGCTTCCGCTGACGGGTGGGCGCGGGTTTACCCGCGAATGCGTCATTGAATTTGCCGACGCATTCGCGGGTAAACCCGCTCCTACACAGGGGCGTTGCACAGCAGGATTAGGGGTTCAGACCGCCCTGTTCAACTTCACCCATGACGCAAACCTGTCGACAAACGCCTGCAAGAACGGTCGCGTCTTGTCATTGAGCTTGCCCGAA  >  AE015451_phaM1/4675226‑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: ‑