Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451_phaM14,344,3080GA34.5% 125.3 / 23.7 123intergenic (+39/‑117)PP_3817 CDS (Polyamine ABC transporter_ ATP‑binding protein)/PP_3818 CDS (Periplamic phosphate‑binding protein)–/–
Reads supporting (aligned to +/- strand):  ref base G (43/29);  new base A (38/0);  total (92/31)
Fisher's exact test for biased strand distribution p-value = 4.84e-07
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 162 total aligned reads displayed.

CGGGCGACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCT‑‑‑‑GGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGT‑CACATTAT‑CATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCT  >  AE015451_phaM1/4344172‑4344443
                                                                                                                                        |                                                                                                                                             
  ggCGACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGCGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGACTTCCAGCTCTAACaccacc                                                                                                                                             >  1:2004469/1‑137 (MQ=255)
  ggCGACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGTCTGCCAGATCGAGCAccgcc                                                                                                                                             >  2:4762463/1‑137 (MQ=255)
  ggCGACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGGCTGCCAGACCAAGCaccacc                                                                                                                                             >  2:11779563/1‑137 (MQ=255)
  ggCGACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCaccacc                                                                                                                                             >  2:12183653/1‑137 (MQ=255)
  ggCGACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCAccgcc                                                                                                                                             >  2:18293853/1‑137 (MQ=255)
  ggCGACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGACCTCCAAAGCGAGAaccacc                                                                                                                                             >  2:16810221/1‑137 (MQ=255)
  ggCGACAGCGTGAAGGTGCGTTGGACAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGTCTGTCAGATCGAGAAccgcc                                                                                                                                             >  2:3179613/1‑137 (MQ=255)
   gCGACGGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCAccgccc                                                                                                                                            <  1:11481425/137‑1 (MQ=255)
   gCGACAGCGTGAAGGTTCGTGGGGCAGTGGAGGACGCTTGCATTTATACCGAATGGGCCGAGAGAGACCTGAGCAGGTCGGCCGGGGCGCGTTAGCGCTCTTGTACGTGGCAGTGCCTGCCAGATCGAGCAtcgccc                                                                                                                                            <  2:14912439/137‑1 (MQ=255)
   gCGACAGCGTGAAGGTGCGTTGTGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCAccgccc                                                                                                                                            <  2:9063793/137‑1 (MQ=255)
   gCGACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCAccgccc                                                                                                                                            <  2:4152846/137‑1 (MQ=255)
   gCGACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCAccgccc                                                                                                                                            <  2:11411740/137‑1 (MQ=255)
   gCGACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCAccgccc                                                                                                                                            <  1:14925167/137‑1 (MQ=255)
   gCGACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCAccgccc                                                                                                                                            <  1:12151688/137‑1 (MQ=255)
   gCGACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCAccgccc                                                                                                                                            <  2:3077600/137‑1 (MQ=255)
   gCGACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCAccgccc                                                                                                                                            <  2:3077599/137‑1 (MQ=255)
   gCGACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCAccgccc                                                                                                                                            <  2:7227543/137‑1 (MQ=255)
   gCGACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCAccgccc                                                                                                                                            <  2:8481522/137‑1 (MQ=255)
   gCGACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCAccgccc                                                                                                                                            <  2:10304949/137‑1 (MQ=255)
   gCGACAGCGTGAAGGTGCGTTGGGCAGTGACAGACGCTTACATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCTGCCGGGGCGCAGTAACGCTCTTGTACGTTGCAGCGCCTGCCGGATCGAGCAccgccc                                                                                                                                            <  2:12348573/137‑1 (MQ=255)
    cGACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCAccgcccg                                                                                                                                           <  1:15534480/137‑1 (MQ=255)
    cGACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCAccgcccg                                                                                                                                           <  1:15534507/137‑1 (MQ=255)
        aGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCcgcgcg                                                                                                                                       <  1:6192955/137‑1 (MQ=255)
        aGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCcgcgcg                                                                                                                                       <  1:3970483/137‑1 (MQ=255)
        aGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCcgcgcg                                                                                                                                       <  1:3970455/137‑1 (MQ=255)
        aGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCcgcgcg                                                                                                                                       <  2:13058253/137‑1 (MQ=255)
        aGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCcgcgcg                                                                                                                                       <  2:15498288/137‑1 (MQ=255)
        aGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCcgcgcg                                                                                                                                       <  1:13557995/137‑1 (MQ=255)
        aGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCcgcgcg                                                                                                                                       <  1:10652580/137‑1 (MQ=255)
         gCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCgg                                                                                                                                      <  2:12608821/137‑1 (MQ=255)
         gCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCgg                                                                                                                                      <  1:2130905/137‑1 (MQ=255)
          cGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGGCTTCCAAAGTGAGCACCACCCCCCCgga                                                                                                                                     >  1:3388945/1‑136 (MQ=255)
          cGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGc                                                                                                                                     <  1:2620467/137‑1 (MQ=255)
          cGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGc                                                                                                                                     <  2:61869/137‑1 (MQ=255)
          cGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGc                                                                                                                                     <  2:4863933/137‑1 (MQ=255)
          cGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGc                                                                                                                                     <  1:10193797/137‑1 (MQ=255)
          cGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGc                                                                                                                                     <  1:4614154/137‑1 (MQ=255)
           gTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTTCCAGCTCGAGCCCCACCCGTGCGGca                                                                                                                                    >  2:18254496/1‑137 (MQ=255)
              aaGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCacac                                                                                                                                 <  2:16737995/137‑1 (MQ=255)
              aaGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCacac                                                                                                                                 <  2:637200/137‑1 (MQ=255)
               aGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCCGATCGAGCACCGCCCCCTCGGTacaca                                                                                                                                >  2:11299815/1‑136 (MQ=255)
                    cGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGGCTGTCAGATCGAGCACCGCCCCCGCGGCACACGAGCT‑‑‑‑g                                                                                                                       >  2:2140775/1‑137 (MQ=255)
                       tGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGTCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCT‑‑‑‑GGCa                                                                                                                    >  2:8358276/1‑137 (MQ=255)
                       tGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGTCTGCCAGATCGAATACCGCCCCCGTGGCACACGACCT‑‑‑‑GTCa                                                                                                                    >  2:16480769/1‑137 (MQ=255)
                       tGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGGCTTCCCGACCGAGCACCGCCCACGCGGAACACAATCT‑‑‑‑GGCa                                                                                                                    >  1:16838943/1‑137 (MQ=255)
                       tGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGGCTGTCAGACCGAGCACCCCCCACACAGCACAAAAACT‑‑‑‑GGCa                                                                                                                    >  1:6484021/1‑137 (MQ=255)
                       tGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGGCTGCCCGACCGAGCACCGCCCCCACGGAACACAAACT‑‑‑‑TGCa                                                                                                                    >  1:7837227/1‑137 (MQ=255)
                       tGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCC‑‑‑‑GGCa                                                                                                                    >  2:12598303/1‑137 (MQ=255)
                       tGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAAACCTAACACCGCCCCCGCAGTACACTATCtttca                                                                                                                        >  2:17125082/1‑133 (MQ=255)
                       tGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCCTCCAGATCGAGCACCACCCACGCGGCCCACGATCT‑‑‑‑GGCa                                                                                                                    >  2:12529720/1‑137 (MQ=255)
                       tGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCCGTCCCATCAAGCACCACCCCCGCGGCACACAACCC‑‑‑‑GGCa                                                                                                                    >  1:7637301/1‑137 (MQ=255)
                       tGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCCGCCAGAGCGAACACCGCCCCCACGGCAAAAGACCT‑‑‑‑GGCa                                                                                                                    >  2:16318347/1‑137 (MQ=255)
                       tGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGACCGCCCGACCGAGCACCCCCCACGCGGCACACAAGCT‑‑‑‑GGCa                                                                                                                    >  2:16455281/1‑137 (MQ=255)
                         ggCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTTCCCAATCTAGCACCACCCCCACGGCACACAACCT‑‑‑‑GGCaag                                                                                                                  >  2:303723/1‑135 (MQ=255)
                           cAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTTAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCcgcgcg                                                                                                                                       >  1:5484288/1‑118 (MQ=255)
                           cAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTTAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCcgcgcg                                                                                                                                       <  2:5484288/118‑1 (MQ=255)
                            aGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGACCGCCCGATCTAAAACCCACCACGCGGCACACAAGCC‑‑‑‑GGCaaggga                                                                                                               >  1:9021120/1‑132 (MQ=255)
                              tGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGTCAAACCGAGCACCGCCCGCGCGGCACACGATCT‑‑‑‑GGCAggggatg                                                                                                             >  2:4112747/1‑132 (MQ=255)
                              tGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCCCCACCCACGCGGCACACGGTCT‑‑‑‑GGCCGGCGCTg                                                                                                             >  2:15710151/1‑137 (MQ=255)
                              tGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGACCGCCCAATTTAGCACCACCCACGCAAAACACAATCT‑‑‑‑TTCCGGCGCtt                                                                                                             >  2:104442/1‑136 (MQ=255)
                                 cAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGGCTGCCCAAGCTAACACCCCCCCCACTGCACACAACCT‑‑‑‑TGCCGGCGATGCAc                                                                                                          >  2:12093884/1‑137 (MQ=255)
                                 cAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGGATGCCAGATCAAAAACCACCCCCACAGCACACAATAC‑‑‑‑TGCAGGGGATGCaa                                                                                                          >  1:16038858/1‑136 (MQ=255)
                                 cAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCCAATCGAGCACCGCCCGCGCGGCACACGATCT‑‑‑‑GGCAGGCGATGCAc                                                                                                          >  2:11228095/1‑137 (MQ=255)
                                 cAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGACCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGGCCGCCAGATCGAGCACCACCCACGCAGCACACGATCC‑‑‑‑GGCAAGGGATGCaa                                                                                                          >  1:11125988/1‑136 (MQ=255)
                                  aGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAATCGGCCGGGGCGCATTAACGCTCCTGTACGTTGCAGCGCCTTCCAGGTCTAACACCGCCCGCGCGGCACACGGCCT‑‑‑‑GGCAGGCACTTCACt                                                                                                         >  2:13185573/1‑137 (MQ=255)
                                   gACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGAGCGAACACCGACCACGCAGCACACAACCT‑‑‑‑CGCCGGCGATgtaatc                                                                                                        >  2:565840/1‑132 (MQ=255)
                                      gCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTTCCAAATCGAGAACCACCCACACGGCACACTATTT‑‑‑‑GTCAGgcactgcactcctg                                                                                                     >  2:15526776/1‑133 (MQ=255)
                                      gCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGACTGCCCGACCGAACACCGCCCACACGGAACACAACCT‑‑‑‑GGGAGGGGATGCACTAGTg                                                                                                     >  1:14741521/1‑137 (MQ=255)
                                           cATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCCTCCCAATCGAGCACCACCCCCCCGGCACACAACCT‑‑‑‑GGCAGGAGCCGCACCACTTGAgga                                                                                                >  2:11571356/1‑136 (MQ=255)
                                             tttATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGGCTTCCAAAGCGAGCACCGCCCGCGCGGCACACGAGCT‑‑‑‑GGTAGGGGATGGAACATTGGAggagg                                                                                              >  1:3445159/1‑134 (MQ=255)
                                                   ccGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGGATGCCCGAGCGAGAAACACCCCCACGGCCCACAACCT‑‑‑‑TGCAAACGCTGGACTCTTTGGGGAGTACGccc                                                                                        >  1:17561223/1‑137 (MQ=255)
                                                     gAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCAGGCGCTGCACTATTGGAGGCGTACGc                                                                                          >  2:8101383/1‑137 (MQ=255)
                                                        tGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGGCTTCCAAATCGAGAACCACCCACCCGGCCCCCAACCT‑‑‑‑TGCCAGCGCCGCACCACTTGAGGCGTCCAccccaccc                                                                                   >  2:15017641/1‑133 (MQ=255)
                                                         gggCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGTCTGCCCGATCGAGCACCGCCCGTGCGGCACACGAGCT‑‑‑‑GGCAAGGGATGCACCATTTGAGGAGTCCGCCCCGT‑cac                                                                                 >  2:16333994/1‑137 (MQ=255)
                                                         gggCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTTCCCAATCGAACACCACCCCCGCGGCACACAACCC‑‑‑‑GGCCGGGGCTGCACCCCTTGAGGAGTACTCCCCCT‑cac                                                                                 >  2:174183/1‑137 (MQ=255)
                                                         gggCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAAATCGAGCACCGCCCGCGCGGCACACAATCT‑‑‑‑GGCAGGAGCTGCACTATTTGAGGAGTACACCCCAT‑cac                                                                                 >  2:13741019/1‑137 (MQ=255)
                                                         gggCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCCTCCCAATCAAACCCCACCCACACGGCCCACAACCT‑‑‑‑TGCCAGCACTGCACCACTTGAGGAGTACGCCCCCT‑cac                                                                                 >  2:7420052/1‑137 (MQ=255)
                                                         gggCCGAGAGCGACCCGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGGCTGCCATACCTAGCACCGCCCCCTTGGCACACGATCT‑‑‑‑TGCCAGCGCCTCAGCCCTTGAGGCGTACTCCCAGT‑cac                                                                                 >  1:13920042/1‑137 (MQ=255)
                                                          ggCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGTCTGCCCAACCGAACACCACCCCCACGGCACACCACCC‑‑‑‑TGCAGGCGCTTCCCTCTTTTAGGCGTCCAccccccccca                                                                                 >  1:10777379/1‑131 (MQ=255)
                                                          ggCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGGCTTCCCAATCCAGCACCACCCACACAGCCCCCAACCT‑‑‑‑GGCCGGCACTGCACTCCTTGAGGGGTACCCCCCAT‑ccca                                                                                >  1:2380834/1‑137 (MQ=255)
                                                          ggCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGGATGTCAGATCGAGAACCACCCCCACGGCACACAAGCT‑‑‑‑CGAAAGAAATGCACTATTTGAGGTGGCCGccccaccccc                                                                                 >  1:5631028/1‑131 (MQ=255)
                                                             cGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGTCTGCCCAATCTAACACCACCCACGAGGTACACGACCT‑‑‑‑GGCCGGTGGTTGACTCTTGGAAGAGTCCGCCCCCT‑cccactc                                                                             >  2:12786412/1‑134 (MQ=255)
                                                              gagaGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGTCTGCCAAATCGCGCACCACCCCCGCGGCAAACAAGCT‑‑‑‑GGAAGGAGCTGCACTAGTTGGGGCGGACGCCCCAG‑CACATTat                                                                            >  1:7062605/1‑137 (MQ=255)
                                                              gagaGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGGCTGCCCAATCGAGCACCGCCCCCGCTGAACACAAGCT‑‑‑‑GGCAGGTGGTGTACCCCTTGAGGAGGACGCCCCAA‑cacaccat                                                                            >  1:6175387/1‑133 (MQ=255)
                                                                        tGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGGCTGCCCGATCGAGCACCGCCCACGCGGAACACAATCT‑‑‑‑GGCAGGCGCCGCACCCTTTGAGGCGTACGCCCCAT‑CACATTAT‑CCCCACGCCa                                                                 >  2:12543831/1‑137 (MQ=255)
                                                                        tGAGCAAGTCGGCCGGGGCGCATTAACGCTATTGTACGTTGCAGCGTCCTCCAGATCGAAGACCGCCCGCGCGGTACACGATCT‑‑‑‑GGTAGGCGTTGCACTACTGGAAGCGTGCGCCCCAT‑CAAACTCt‑catcaaacac                                                                 >  2:3293522/1‑132 (MQ=255)
                                                                                  ggCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCCGATCGAGCCCCGCCCGCACAGCACACGATCT‑‑‑‑TGCAGGCGCTGCACCATTGGCGGCGTACGCCCCAT‑CCCACTAT‑CCTCACACCCCCGTCGTGAt                                                       >  2:1263641/1‑137 (MQ=255)
                                                                                       gggCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCTGTACAAGATCT‑‑‑‑GGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGT‑CCCATTAT‑CATCACGCCAATGTCGTAATCTGcc                                                  >  2:3818817/1‑137 (MQ=255)
                                                                                       gggCGCATTAACGCTCTTGTACGTTGCAGCGACCGCCAGAGCGAGTACCGCCCGCGTGGCACACGATCT‑‑‑‑TGCAGGCGCTGTACTACTTGAGGGGTACACCCCATCCACA‑TAT‑CATCACGCCAACGTCGGCATCTGcc                                                  >  2:2835533/1‑137 (MQ=255)
                                                                                                          tACGTTGCAGCGCCTGCCAAACCGAGCACCACCCGCGCCGGACACGATCT‑‑‑‑GGCAGGGAATGTACCGTTGGAAGAGTACACCCCGT‑CACATTCTCCAACA‑ACAAATGTCGTAATCTGTCAGCCAAATCCCTGCcacaa                               >  2:6645219/1‑136 (MQ=255)
                                                                                                                  aGCGCCTTCCCAATCGAGCACCACCCGCACGGCACACGACCT‑‑‑‑TGCAGGCGCCGCACCACTTGAAGCGTACACCCCAC‑CACACTAT‑CATCAAACCACCGTCGTCATCTGCCGGCCAACTCTATGCCACAGGTCACCCg                       >  2:4478018/1‑137 (MQ=255)
                                                                                                                  aGCGCCTGCCAGATCGAGCACCGCCCGCGTGGCACCCGATCT‑‑‑‑GGCAGGCGCTGGACTCTTTGAGGAGTACGCCCCCT‑TACATTAT‑CATCACGCTACTTTCGTCACCTGTCGGCCGACCCTCTGCCACAGGTCGcccc                       >  2:13730412/1‑136 (MQ=255)
                                                                                                                    cgcCTTCCAGACCTAGCACCACCCCCGCGGCCCACGATCT‑‑‑‑GTTAAGCGCTGCACCCTTGGAGGCGTACGCCCCGT‑CCCATTTT‑CATCCCACCCCCGTCGTACTCTTCCTGCCCACCCTCTGCCACCAGTCCCCCGGc                     >  2:14707355/1‑137 (MQ=255)
                                                                                                                            cgATCGAGCACCGCCCGCGAGGCACACGATCT‑‑‑‑GGCAGGCGCTGCAACCTTGGAGGCGTCCGCCCCGT‑CCCATTAT‑TATTCCACCAACGTCGGGATCTGCCGGCCAACTCTATGCNCCAGGTCACCCGGCATTCTCCg             >  2:9343591/2‑137 (MQ=255)
                                                                                                                                 gAGCACCGCCCGCGCGGCACACGATCT‑‑‑‑GGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGT‑CACATTAT‑CATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGc        >  2:90253/1‑137 (MQ=255)
                                                                                                                                   ccACCGCCCGCGAGGACCACCCTTT‑‑‑‑AGCATCCGCTGCACTATTGGAGGCGTACGCCCCGT‑CACATTAT‑CATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCtt      <  1:7088636/136‑1 (MQ=255)
                                                                                                                                    cACCCCCCGCGCGAAAAACACTCG‑‑‑‑GGAAGCCGCTGCACTATTGGAGGCGTACGCCCCGT‑CACATTAT‑CCTCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTg     <  2:11274925/137‑1 (MQ=255)
                                                                                                                                       cgcccgcGCGGCACACGATCT‑‑‑‑GGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGT‑CACATTAT‑CATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCt  >  2:7053741/1‑137 (MQ=255)
                                                                                                                                       cgcccgcGCGGCACACGATCT‑‑‑‑GGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGT‑CACATTAT‑CATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCt  >  1:12101990/1‑137 (MQ=255)
                                                                                                                                       cgcccgcGCGGCACACGATCT‑‑‑‑GGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGT‑CACATTAT‑CATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCt  >  1:12102015/1‑137 (MQ=255)
                                                                                                                                        |                                                                                                                                             
CGGGCGACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCT‑‑‑‑GGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGT‑CACATTAT‑CATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCT  >  AE015451_phaM1/4344172‑4344443

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