Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451_phaM14,344,3230GA29.1% 242.3 / 17.5 147intergenic (+54/‑102)PP_3817 CDS (Polyamine ABC transporter_ ATP‑binding protein)/PP_3818 CDS (Periplamic phosphate‑binding protein)–/–
Reads supporting (aligned to +/- strand):  ref base G (92/8);  new base A (40/1);  total (136/11)
Fisher's exact test for biased strand distribution p-value = 4.48e-01
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 147 total aligned reads displayed.

CGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGC‑CGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTAT‑CATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTGCTGCTCACCCT  >  AE015451_phaM1/4344192‑4344456
                                                                                                                                    |                                                                                                                                      
cGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGC‑CGGGGCGCATTAACGCTCTTGTACGTTGCAGCGGCTGTCAGATCGAGCACCGCCCCCGCGGCACACGAGCTg                                                                                                                                   >  2:2140775/1‑137 (MQ=255)
   tGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGC‑CGGGGCGCATTAACGCTCTTGTACGTTGCAGCGTCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCa                                                                                                                                >  2:8358276/1‑137 (MQ=255)
   tGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGC‑CGGGGCGCATTAACGCTCTTGTACGTTGCAGCGTCTGCCAGATCGAATACCGCCCCCGTGGCACACGACCTGTCa                                                                                                                                >  2:16480769/1‑137 (MQ=255)
   tGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGC‑CGGGGCGCATTAACGCTCTTGTACGTTGCAGCGGCTTCCCGATCGAACACCACCCACACGGCACACGAACTGGCa                                                                                                                                >  1:7563703/1‑137 (MQ=255)
   tGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGC‑CGGGGCGCATTAACGCTCTTGTACGTTGCAGCGGCTGTCCAACCGAGAACCACCCACACGGCACACAATccagca                                                                                                                                >  1:7826109/1‑132 (MQ=255)
   tGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGC‑CGGGGCGCATTAACGCTCTTGTACGTTGCAGCGGCTGTCAGACCGAGCACCCCCCACACAGCACAAAAACTGGCa                                                                                                                                >  1:6484021/1‑137 (MQ=255)
   tGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGC‑CGGGGCGCATTAACGCTCTTGTACGTTGCAGCGGCTGCCCGACCGAGCACCGCCCCCACGGAACACAAACTTGCa                                                                                                                                >  1:7837227/1‑137 (MQ=255)
   tGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGC‑CGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCCGATCGAGCACCACCCACGCGGCACACAAGCCGGCa                                                                                                                                >  1:9517996/1‑137 (MQ=255)
   tGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGC‑CGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCCGCCAGAGCGAACACCGCCCCCACGGCAAAAGACCTGGCa                                                                                                                                >  2:16318347/1‑137 (MQ=255)
     ggCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGC‑CGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTTCCCGATCGAGCACCACCCACGCGGCACACGAGCTGGTAgg                                                                                                                              >  2:303724/1‑137 (MQ=255)
        aGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGC‑CGGGGCGCATTAACGCTCTTGTACGTTGCAGCGACCGCCCGATCTAAAACCCACCACGCGGCACACAAGCCGGCaaggga                                                                                                                           >  1:9021120/1‑132 (MQ=255)
          tGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGC‑CGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGTCAAACCGAGCACCGCCCGCGCGGCACACGATCTGGCAggggatg                                                                                                                         >  2:4112747/1‑132 (MQ=255)
          tGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGC‑CGGGGCGCATTAACGCTCTTGTACGTTGCAGCGACCGCCCAATTTAGCACCACCCACGCAAAACACAATCTTTCCGGCGCtt                                                                                                                         >  2:104442/1‑136 (MQ=255)
           ggCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGC‑CGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCACCCACACGGCACACGATCTGGCAggagatgt                                                                                                                        >  1:13766072/1‑131 (MQ=255)
             cAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGC‑CGGGGCGCATTAACGCTCTTGTACGTTGCAGCGGATGCCAGATCAAAAACCACCCCCACAGCACACAATACTGCAGGGGATGCaa                                                                                                                      >  1:16038858/1‑136 (MQ=255)
             cAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGC‑CGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCCAACCGAACACCCCCCACGCGGCACACAACCTGGCCGGCGATGCAc                                                                                                                      >  1:17790624/1‑137 (MQ=255)
             cAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGC‑CGGGGCGCATTAACGCTCTTGTACGTTGCAGCGACTGCCCAACCGAACACCCCCCCCACGGCACACAACCCTGCAGGAGCttccc                                                                                                                      >  1:3826910/1‑133 (MQ=255)
             cAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGAC‑CGGGGCGCATTAACGCTCTTGTACGTTGCAGCGGCCGCCAGATCGAGCACCACCCACGCAGCACACGATCCGGCAAGGGATGCaa                                                                                                                      >  1:11125988/1‑136 (MQ=255)
              aGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGC‑CGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCCGCCCAGCCGAGAACCGCCCGCGCGGCACAAGCTCTGGCAGGAGCTGCACt                                                                                                                     >  2:13185545/1‑137 (MQ=255)
              aGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAATCGGC‑CGGGGCGCATTAACGCTCCTGTACGTTGCAGCGCCTTCCAGGTCTAACACCGCCCGCGCGGCACACGGCCTGGCAGGCACTTCACt                                                                                                                     >  2:13185573/1‑137 (MQ=255)
               gACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGC‑CGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGAGCGAACACCGACCACGCAGCACACAACCTCGCCGGCGATgtaatc                                                                                                                    >  2:565840/1‑132 (MQ=255)
                  gCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGC‑CGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTTCCAAATCGAGAACCACCCACACGGCACACTATTTGTCAGgcactgcactcctg                                                                                                                 >  2:15526776/1‑133 (MQ=255)
                  gCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGC‑CGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCCAATCGAGCACCCCCCACACGGCCCACAACCTTGCAGGCGATGCCCCATTg                                                                                                                 >  2:14956201/1‑137 (MQ=255)
                  gCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGC‑CGGGGCGCATTAACGCTCTTGTACGTTGCAGCGACTGCCCGACCGAACACCGCCCACACGGAACACAACCTGGGAGGGGATGCACTAGTg                                                                                                                 >  1:14741521/1‑137 (MQ=255)
                       cATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGC‑CGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCCTCCCAATCGAGCACCACCCCCCCGGCACACAACCTGGCAGGAGCCGCACCACTTGAgga                                                                                                            >  2:11571356/1‑136 (MQ=255)
                         tttATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGC‑CGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCCGATCGAGGACCACCCGCGCAGCACAAGAGCTGGCAAGCAATGCACTATTGGAGGCGt                                                                                                          >  2:7064461/1‑137 (MQ=255)
                         tttATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGC‑CGGGGCGCATTAACGCTCTTGTACGTTGCAGCGACCTCCCGATCGAGCACCGCCCGCGCGGCACACAATCTGGCAGGAGCTGCACTACTGGaaggat                                                                                                          >  1:11568267/1‑132 (MQ=255)
                               ccGAATGGGCCGAGAGCGACCTGAGCAAGTCGGC‑CGGGGCGCATTAACGCTCTTGTACGTTGCAGCGGATGCCCGAGCGAGAAACACCCCCACGGCCCACAACCTTGCAAACGCTGGACTCTTTGGGGAGTACGccc                                                                                                    >  1:17561223/1‑137 (MQ=255)
                               ccGAATGGGCCGAGAGCGACCTGAGCAAGTCGGC‑CGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACCCGACCTGGCAGGCGCTGACCCCCTGGAGGCGGCCGccc                                                                                                    >  2:11996971/1‑137 (MQ=255)
                                 gAATGGGCCGAGAGCGACCTGAGCAAGTCGGC‑CGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCCGGCAAGCCTAGCCCCGCCCACGCTGCACACGACCTGGCAAGATCCGCACTCTTGGAGGACTACGcccca                                                                                                  >  2:16109668/1‑136 (MQ=255)
                                    tggttagagaGCTACCTGAG‑AAGTCTGCTCGGGGAGCATGAACTCTCTTGTGTGTTAAAGCGCCTGCCAGACAAAACACCCCCCGCGCGACACCCCCTCTGACAGCCGCTGCACTATTGGAGGCGTACGCCCCGTca                                                                                               <  2:14183485/131‑1 (MQ=255)
                                    tGGGCCGAGAGCGACCTGAGCAAGTCGGC‑CGGGGCGCATTAACGCTCTTGTACGTTGCAGCGGCTTCCAAATCGAGAACCACCCACCCGGCCCCCAACCTTGCCAGCGCCGCACCACTTGAGGCGTCCAccccaccc                                                                                               >  2:15017641/1‑133 (MQ=255)
                                     gggCCGAGAGCGACCTGAGCAAGTCGGC‑CGGGGCGCATTAACGCTCTTGTACGTTGCAGCGTCTGCCCGATCGAGCACCGCCCGTGCGGCACACGAGCTGGCAAGGGATGCACCATTTGAGGAGTCCGCCCCGTcac                                                                                              >  2:16333994/1‑137 (MQ=255)
                                     gggCCGAGAGCGACCTGAGCAAGTCGGC‑CGGGGCGCATTAACGCTCTTGTACGTTGCAGCGGCCTCCAGACCGAGCACCACCCACGCGGCACACGATCTGGCAGGAGCTGCAATCTTTGAGGCGTACACCCCCTcac                                                                                              >  2:12198916/1‑137 (MQ=255)
                                     gggCCGAGAGCGACCTGAGCAAGTCGGC‑CGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTTCCCAATCGAACACCACCCCCGCGGCACACAACCCGGCCGGGGCTGCACCCCTTGAGGAGTACTCCCCCTcac                                                                                              >  2:174183/1‑137 (MQ=255)
                                     gggCCGAGAGCGACCTGAGCAAGTCGGC‑CGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAAATCGAGCACCGCCCGCGCGGCACACAATCTGGCAGGAGCTGCACTATTTGAGGAGTACACCCCATcac                                                                                              >  2:13741019/1‑137 (MQ=255)
                                     gggCCGAGAGCGACCTGAGCAAGTCGGC‑CGGGGCGCATTAACGCTCTTGTACGTTGCAGCGACCGCCCAACCTAACACCCCCCCCACAGAACACAACCTGGCAAGCGATGAACCACTGGAGGCGTACGCCCCATcac                                                                                              >  2:5775843/1‑137 (MQ=255)
                                      ggCCGAGAGCGACCTGAGCAAGTCGGC‑CGGGGCGCATTAACGCTCTTGTACGTTGCAGCGTCTGCCCAACCGAACACCACCCCCACGGCACACCACCCTGCAGGCGCTTCCCTCTTTTAGGCGTCCAccccccccca                                                                                             >  1:10777379/1‑131 (MQ=255)
                                          gagaGCGACCTGAGCAAGTCGGC‑CGGGGCGCATTAACGCTCTTGTACGTTGCAGCGGCTGCCCAATCGAGCACCGCCCCCGCTGAACACAAGCTGGCAGGTGGTGTACCCCTTGAGGAGGACGCCCCAAcacaccat                                                                                         >  1:6175387/1‑133 (MQ=255)
                                          gagaGCGACCTGAGCAAGTCGGC‑CGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCCGATCTAGCACCACCCACGCGGCACACGATCTTGCAGGGGCTGCACTATTTGAGGCGTACACCCCATCACACTat                                                                                         >  2:15675777/1‑137 (MQ=255)
                                               cGACCTGAGCAAGTCGGC‑CGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCCGCCAGATCGAGCACCGCCCACGCGGCACACGATCTGGCAGGTGCTGGACCATTGGAGGAGTACACCCCATCAAACTat‑caccc                                                                                   >  2:4591465/1‑134 (MQ=255)
                                                    tGAGCAAGTCGGC‑CGGGGCGCATTAACGCTCTTGTACGTTGCAGCGGCTGCCCGATCGAGCACCGCCCACGCGGAACACAATCTGGCAGGCGCCGCACCCTTTGAGGCGTACGCCCCATCACATTAT‑CCCCACGCCa                                                                              >  2:12543831/1‑137 (MQ=255)
                                                    tGAGCAAGTCGGC‑CGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAAATCGAGCACCCCCCACGCGGCACACAATCTGGTAGGAGCTGTAATATTTGGGGCGTAAGCCCCATCACATTCT‑CCCCACGCCa                                                                              >  2:3293549/1‑137 (MQ=255)
                                                    tGAGCAAGTCGGC‑CGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCCGCCAGAGTGGGGACCACCCGCAAGGAACAAAAGCTTGGAGGCGATGAACTATTTGGGGCGTACGCCCCGTTACATTAT‑CCCCACGccc                                                                              >  2:6895076/1‑136 (MQ=255)
                                                    tGAGCAAGTCGGC‑CGGGGCGCATTAACGCTATTGTACGTTGCAGCGTCCTCCAGATCGAAGACCGCCCGCGCGGTACACGATCTGGTAGGCGTTGCACTACTGGAAGCGTGCGCCCCATCAAACTCt‑catcaaacac                                                                              >  2:3293522/1‑132 (MQ=255)
                                                              ggC‑CGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCCGATCGAGCCCCGCCCGCACAGCACACGATCTTGCAGGCGCTGCACCATTGGCGGCGTACGCCCCATCCCACTAT‑CCTCACACCCCCGTCGTGAt                                                                    >  2:1263641/1‑137 (MQ=255)
                                                                    gggCGCATTAACGCTCTTGTACGTTGCAGCGCCCGCCAGATCGAACACCACCCACGCAGCACACAACCTGGCAAGAGATGCACTACTGGGGGCGTACGCCCCGCCACACTAT‑CATCACACCACTGTTGTAATCTTcc                                                               >  2:2835572/1‑137 (MQ=255)
                                                                                       tACGTTGCAGCGCCTGCCAAACCGAGCACCACCCGCGCCGGACACGATCTGGCAGGGAATGTACCGTTGGAAGAGTACACCCCGTCACATTCTCCAACA‑ACAAATGTCGTAATCTGTCAGCCAAATCCCTGCcacaa                                            >  2:6645219/1‑136 (MQ=255)
                                                                                               aGCGCCTGCCAGATCGAGCACCGCCCGCGTGGCACCCGATCTGGCAGGCGCTGGACTCTTTGAGGAGTACGCCCCCTTACATTAT‑CATCACGCTACTTTCGTCACCTGTCGGCCGACCCTCTGCCACAGGTCGcccc                                    >  2:13730412/1‑136 (MQ=255)
                                                                                                 cggcTGCCAGATCGAGCACCGCCCACGCGGCACACGATCTGGCAGGCGCTGCACCACTGGAGGCGTACGCCCCATCACACTAC‑CATCCCGACACTGTCGTAATCTGCCGGCCCAACCTATGCCAAAGGTCACCCGGc                                  >  2:4465762/4‑137 (MQ=255)
                                                                                                 cgcCTTCCAGACCTAGCACCACCCCCGCGGCCCACGATCTGTTAAGCGCTGCACCCTTGGAGGCGTACGCCCCGTCCCATTTT‑CATCCCACCCCCGTCGTACTCTTCCTGCCCACCCTCTGCCACCAGTCCCCCGGc                                  >  2:14707355/1‑137 (MQ=255)
                                                                                                 cgcCTGCCAGATCGAGCACCACCCACGCGGCCCACAACCTGGCAGGCGCTGCACCACTGGAGGCGTACGCCCCGTCACACTAT‑CCCCACGCCACCTTCGTCATCTTCCCGCCAACCCCCTGTCACAGGTTACCCGGc                                  >  2:9938330/1‑137 (MQ=255)
                                                                                                 cgcCTGCCAAATCGAGCACCGCCCGCGCGGCACACGACCTTGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCCCATTGT‑CATCACGCCAATGTCTTAAGCTGTCAGCCAACTCTGTGCCACAGGTCAACCAGc                                  >  2:12786390/1‑137 (MQ=255)
                                                                                                         aGATCAAGCACCGCCCCCACGGCACACGATCTGGCAGGCGCTGCACCACTGGAGGCGTACACCCCATCACACTAC‑CCTCACGCCACTGTCGTAATCTTTCGGCCAAACCCATGCCACAAGTCAACCAGCATGCTCCg                          >  2:3094462/1‑137 (MQ=255)
                                                                                                         aGACCGAACACCGCCCACACGGCACACAACCTGGCAAGCGATACAATACTGGAGGCGTACGCCCCACCCAACTAT‑CCTCAAGCCACCTTCATCATCTTCCGGCCAACTCTATGCCACAAGGCCCCCAGCACGCTCCg                          >  2:12595388/1‑137 (MQ=255)
                                                                                                              gAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTAT‑CATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGc                     >  2:90253/1‑137 (MQ=255)
                                                                                                               aGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTAT‑CATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCt                    >  2:11258978/1‑137 (MQ=255)
                                                                                                                ccACCGCCCGCGAGGACCACCCTTTAGCATCCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTAT‑CATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCtt                   <  1:7088636/136‑1 (MQ=255)
                                                                                                                 cACCCCCCGCGCGAAAAACACTCGGGAAGCCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTAT‑CCTCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTg                  <  2:11274925/137‑1 (MQ=255)
                                                                                                                    cgcccgcGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTAT‑CATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCt               >  2:18516793/1‑137 (MQ=255)
                                                                                                                    cgcccgcGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTAT‑CATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCt               >  1:12101990/1‑137 (MQ=255)
                                                                                                                    cgcccgcGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTAT‑CATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCt               >  2:16635740/1‑137 (MQ=255)
                                                                                                                    cgcccgcGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTAACATTAT‑CATCACGCCACTGTCGTACTCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCt               >  1:3245466/1‑137 (MQ=255)
                                                                                                                        tggggggCACACGGTTTGGCATCCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTAT‑CATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtgct           <  2:5106380/132‑1 (MQ=255)
                                                                                                                        cgcgcgGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTAT‑CATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtgct           >  1:14933452/1‑137 (MQ=255)
                                                                                                                        cgcgcgGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTAT‑CATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtgct           >  1:6000464/1‑137 (MQ=255)
                                                                                                                        cgcgcgGCACACGATCTGACAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTAT‑CATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCNACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtgct           <  1:18293853/137‑1 (MQ=255)
                                                                                                                        cgcgcgGCACACAATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTAT‑CATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtgct           >  2:7452271/1‑137 (MQ=255)
                                                                                                                         gcgcgGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTAT‑CATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtgctg          >  2:7157114/1‑137 (MQ=255)
                                                                                                                         gcgcgGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTAT‑CATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtgctg          >  1:8131258/1‑137 (MQ=255)
                                                                                                                         gcgcgGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTAT‑CATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtgctg          >  1:14361698/1‑137 (MQ=255)
                                                                                                                         gcgcgGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTAT‑CATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtgctg          >  2:7640468/1‑137 (MQ=255)
                                                                                                                         gcgcgGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTAT‑CATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtgctg          >  2:8760890/1‑137 (MQ=255)
                                                                                                                         gcgcgGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTAT‑CATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtgctg          >  1:6361329/1‑137 (MQ=255)
                                                                                                                         gcgcgGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTAT‑CATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtgctg          >  1:14457114/1‑137 (MQ=255)
                                                                                                                         gcgcgGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTAT‑CATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtgctg          >  2:9467787/1‑137 (MQ=255)
                                                                                                                         gcgcgGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTAT‑CATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtgctg          >  2:13484523/1‑137 (MQ=255)
                                                                                                                         gcgcgGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTAT‑CATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtgctg          >  1:3811091/1‑137 (MQ=255)
                                                                                                                         gcgcgGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTAT‑CATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtgctg          >  1:5633851/1‑137 (MQ=255)
                                                                                                                         gcgcgGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTAT‑CATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtgctg          >  1:15904518/1‑137 (MQ=255)
                                                                                                                         gcgcgGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTAT‑CATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtgctg          >  2:18960055/1‑137 (MQ=255)
                                                                                                                         gcgcgGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTAT‑CATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtgctg          >  1:5633852/1‑137 (MQ=255)
                                                                                                                         gcgcgGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTAT‑CATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtgctg          >  2:17978780/1‑137 (MQ=255)
                                                                                                                         gcgcgGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTAT‑CATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtgctg          >  1:2605301/1‑137 (MQ=255)
                                                                                                                         gcgcgGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTAT‑CATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtgctg          >  2:16884518/1‑137 (MQ=255)
                                                                                                                         gcgcgGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTAT‑CATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtgctg          >  2:1668601/1‑137 (MQ=255)
                                                                                                                         gcgcgGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTAT‑CATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtgctg          >  1:18167092/1‑137 (MQ=255)
                                                                                                                         gcgcgGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTAT‑CATCACGCCACTGTCGTAATCTGCCAGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtgctg          >  2:1768373/1‑137 (MQ=255)
                                                                                                                         gcgcgGCACACGATCTGGCAGGCACTGCACTATTGGAGGCGTACGCCCCGTCACATTAT‑CATCACGCCACCGTCGTAATCTGCCGGCCAACTCTATGCAACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtgctg          >  2:8171862/1‑137 (MQ=255)
                                                                                                                            cgaaacacGGTTTGGCAGACGCTGCACTATTGGAGGCGTACGCCCCGTCACATTAT‑CATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTGCTGCTc       <  2:2985554/133‑1 (MQ=255)
                                                                                                                            cgGCACACGCCTGGGCAGTCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTAT‑CATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTGCTGCTc       <  2:12340136/137‑1 (MQ=255)
                                                                                                                            cgGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTcac                                                                                              >  2:10303472/1‑51 (MQ=255)
                                                                                                                            cgGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTcac                                                                                              <  1:10303472/51‑1 (MQ=255)
                                                                                                                             ggCACACGATCTGGCAGGCGCTGCACTATTGGCGGCGTACGCCCCGTCACATTAT‑CATCACGCCACTGTCGTAATCTGCCGGCCAAATCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTGCTGCTCa      >  1:13212419/1‑137 (MQ=255)
                                                                                                                             ggCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTAT‑CATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCAcc                                      <  1:4623627/105‑1 (MQ=255)
                                                                                                                             ggCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTAT‑CATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCAcc                                      >  2:4623627/1‑105 (MQ=255)
                                                                                                                             ggCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTAT‑CATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGTTTGTCTTGCTGTTCa      >  2:16299640/1‑137 (MQ=255)
                                                                                                                             ggCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTAT‑CATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTGCTGCTCa      >  2:16979329/1‑137 (MQ=255)
                                                                                                                             ggCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTAT‑CATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTGCTGCTCa      >  1:17376820/1‑137 (MQ=255)
                                                                                                                             ggCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTAT‑CATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTGCTGCTCa      >  1:5582755/1‑137 (MQ=255)
                                                                                                                                    |                                                                                                                                      
CGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGC‑CGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTAT‑CATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTGCTGCTCACCCT  >  AE015451_phaM1/4344192‑4344456

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