Predicted mutation
evidence seq id position mutation annotation gene description
RA AE015451_phaM1 1,419,549 +C intergenic (‑38/+98) PP_1242 CDS (conserved protein of unknown function) ← / ← PP_1244 CDS (conserved protein of unknown function) –/–

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451_phaM11,419,5481.C98.1% 1617.5 / 19.9 420intergenic (‑37/+99)PP_1242 CDS (conserved protein of unknown function)/PP_1244 CDS (conserved protein of unknown function)–/–
Reads supporting (aligned to +/- strand):  ref base . (0/8);  new base C (46/365);  total (46/374)
Fisher's exact test for biased strand distribution p-value = 6.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 425 total aligned reads displayed.

GGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCT‑GTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑G‑CTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCTGACCTGCTAATGCAAGGTGTGAGGCTCTGTGGGTAAA  >  AE015451_phaM1/1419416‑1419677
                                                                                                                                       |                                                                                                                                 
 ggTTCGATGCGGAACAGTGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCT‑GTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGggctg‑g‑cctct                                                                                                                              >  2:17266125/1‑133 (MQ=255)
 ggTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCT‑GTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGggctg‑g‑cctct                                                                                                                              <  1:7521659/137‑5 (MQ=255)
     cGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCT‑GTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTcgc                                                                                                                           <  1:3522658/137‑1 (MQ=255)
         gCGGAACAGTGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGTTGGGTTCT‑GTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGc                                                                                                                       <  2:18014184/137‑1 (MQ=255)
         gCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCT‑GTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGc                                                                                                                       <  2:18270912/137‑1 (MQ=255)
           ggAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCT‑GTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCAc                                                                                                                     <  1:12068675/137‑1 (MQ=255)
           ggAACAGCGGGTGAGTGTTGTTTTCGCCGTGGTAGAAGCAGGTCTTGCCGACTGGGGTGGGTGCT‑GTTGTGTGGTAGGGGGGCAGTGTTTTGTTTTCCTTGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCAc                                                                                                                     <  2:2445239/137‑1 (MQ=255)
            gAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCT‑GTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCACg                                                                                                                    <  1:8685350/137‑1 (MQ=255)
             aaCAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCT‑GTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCACGc                                                                                                                   <  2:9197837/137‑1 (MQ=255)
             aaCAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCT‑GTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCACGc                                                                                                                   <  2:11581084/137‑1 (MQ=255)
             aaCAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCT‑GTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCACGc                                                                                                                   <  1:13953223/137‑1 (MQ=255)
             aaCAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCT‑GTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCACGc                                                                                                                   <  1:11653294/137‑1 (MQ=255)
             aaCAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCT‑GTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCACGc                                                                                                                   <  1:3204389/137‑1 (MQ=255)
             aaCAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCT‑GTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCACGc                                                                                                                   <  2:1590924/137‑1 (MQ=255)
               cAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCT‑GTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCACGccc                                                                                                                 <  1:14755580/137‑1 (MQ=255)
               cAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCT‑GTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCACGccc                                                                                                                 <  2:13622571/137‑1 (MQ=255)
               cAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCT‑GTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCACGccc                                                                                                                 <  2:17653945/137‑1 (MQ=255)
               cAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCT‑GTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCACGccc                                                                                                                 <  2:3927360/137‑1 (MQ=255)
               cAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCT‑GTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCACGccc                                                                                                                 <  1:11123738/137‑1 (MQ=255)
               cAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCT‑GTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCACGccc                                                                                                                 <  1:8414123/137‑1 (MQ=255)
               cAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCT‑GTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCACGccc                                                                                                                 <  1:9665148/137‑1 (MQ=255)
               cAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCT‑GTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCACGccc                                                                                                                 <  1:9471784/137‑1 (MQ=255)
                aGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCACGccc                                                                                                                 <  2:6028870/137‑1 (MQ=255)
                   gggTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCT‑GTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCACGCCCGCTc                                                                                                             <  1:4617908/137‑1 (MQ=255)
                   gggTGAGGGTTGTTTTCTCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCT‑GTTGTGTTGTCGGTGGGCGTTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCACGCCCGCTc                                                                                                             <  2:6532125/137‑1 (MQ=255)
                    ggTGAGTGTTTTTTTCGCCTTGTTAGAAGCAGGTCTTGCTGACTGTGGGGGGTTCT‑GTTGTGTTGTCGGTGGGCATTGGTTTGTTTCCCTGGATTCTGGTTTGACTGGGCGG‑GCCTCTTCGCGGGCACGCCCGCTcc                                                                                                            <  2:8581170/137‑1 (MQ=255)
                    ggTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCT‑GTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCACGCCCGCTcc                                                                                                            <  1:9941040/137‑1 (MQ=255)
                    ggTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCT‑GTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCACGCCCGCTcc                                                                                                            <  1:6815713/137‑1 (MQ=255)
                    ggTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCT‑GTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCACGCCCGCTcc                                                                                                            <  2:14114151/137‑1 (MQ=255)
                    ggTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCT‑GTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCACGCCCGCTcc                                                                                                            <  1:2837144/137‑1 (MQ=255)
                    ggTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCT‑GTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCACGCCCGCTcc                                                                                                            <  2:10163673/137‑1 (MQ=255)
                    ggTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCT‑GTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCACGCCCGCTcc                                                                                                            <  2:10458464/137‑1 (MQ=255)
                    ggTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCT‑GTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCACGCCCGCTcc                                                                                                            <  2:10931922/137‑1 (MQ=255)
                    ggTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCT‑GTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCACGCCCGCTcc                                                                                                            <  2:627578/137‑1 (MQ=255)
                    ggTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCT‑GTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCACGCCCGCTcc                                                                                                            <  2:4416562/137‑1 (MQ=255)
                    ggTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCT‑GTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCACGCCCGCTcc                                                                                                            <  1:18338728/137‑1 (MQ=255)
                    ggTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCT‑GTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCACGCCCGCTcc                                                                                                            <  1:266507/137‑1 (MQ=255)
                    ggTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCT‑GTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCACGCCCGCTcc                                                                                                            <  1:2591156/137‑1 (MQ=255)
                    ggTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCT‑GTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCACGCCCGCTcc                                                                                                            <  1:1883528/137‑1 (MQ=255)
                    ggTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCT‑GTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCACGCCCGCTcc                                                                                                            <  2:14522398/137‑1 (MQ=255)
                    ggTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGGGGTGGGTTCT‑GTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCACGCCCGCTcc                                                                                                            <  2:14406988/137‑1 (MQ=255)
                      tGAGTGTTGTTTTCGCCTTGGTGGAAGCAGGTCTTGCCGACTGTGGTGGGTTCT‑GTTGTGTTGTCGGTGGGCATGGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCACGCCCGCTCCca                                                                                                          <  1:11413183/137‑1 (MQ=255)
                      tGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCT‑GTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCACGCCCGCTCCca                                                                                                          <  1:8930955/137‑1 (MQ=255)
                      tGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCT‑GTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCACGCCCGCTCCca                                                                                                          <  1:16535705/137‑1 (MQ=255)
                      tGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCT‑GTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCACGCCCGCTCCca                                                                                                          <  2:841500/137‑1 (MQ=255)
                      tGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCT‑GTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCACGCCCGCTCCca                                                                                                          <  2:16647770/137‑1 (MQ=255)
                      tGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCT‑GTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCACGCCCGCTCCca                                                                                                          <  1:4674066/137‑1 (MQ=255)
                      tGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCT‑GTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCACGCCCGCTCCca                                                                                                          <  2:15238062/137‑1 (MQ=255)
                         gtgtTGTTTTCGCCTTGGTAGAAGCAGTTCTTGCCGACTGTGGTGGGTTCT‑GTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCACGCCCGCTCCCACAg                                                                                                       <  1:13399387/137‑1 (MQ=255)
                         gtgtTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCT‑GTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCACGCCCGCTCCCACAg                                                                                                       <  1:15818305/137‑1 (MQ=255)
                          tgttgtTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCT‑GTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCACGCCCGCTCCCACAgg                                                                                                      <  1:15322737/137‑1 (MQ=255)
                          tgttgtTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCT‑GTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCACGCCCGCTCCCACAgg                                                                                                      <  2:9330311/137‑1 (MQ=255)
                          tgttgtCTTCGCCTTGGTAGAAGCAGGTCTTTCCGACTGTGGTGGGTTCT‑GTTGTGTTGTCGGTGGGCAGTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCACGCCCGCTCCCACAgg                                                                                                      <  2:3201363/137‑1 (MQ=39)
                              gttTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCT‑GTTGTGTTGTCGGTGTGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACt                                                                                                  <  2:2877383/137‑1 (MQ=255)
                              gttTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCT‑GTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACt                                                                                                  <  1:15144251/137‑1 (MQ=255)
                               ttttCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCT‑GTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACtg                                                                                                 <  1:16206903/137‑1 (MQ=255)
                                 ttCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCT‑GTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACtgtg                                                                                               <  1:14238890/137‑1 (MQ=255)
                                  tCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCT‑GTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGc                                                                                              >  2:5605231/1‑137 (MQ=255)
                                     ccTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCT‑GTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCACGCCCGCTcc                                                                                                            <  2:8101865/120‑1 (MQ=255)
                                      cTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCT‑GTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTg                                                                                          <  1:10628970/137‑1 (MQ=255)
                                           tAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCT‑GTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGa                                                                                     <  2:12412444/137‑1 (MQ=255)
                                           tAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCT‑GTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGa                                                                                     <  2:12631586/137‑1 (MQ=255)
                                            aGAAGCAGGTCTTGCCGACTGTGGTGGGTTCT‑GTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGAc                                                                                    <  1:9271777/137‑1 (MQ=255)
                                              aaGCAGGTCTTGCCGACTGTGGTGGGTTCT‑GTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACgg                                                                                  <  1:4514381/137‑1 (MQ=255)
                                                  aGGTCTTGCCGACTGTGGTGGGTTCT‑GTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCtgt                                                                              <  2:15587423/137‑1 (MQ=255)
                                                  aGGTCTTGCCGACTGTGGTGGGTTCT‑GTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCtgt                                                                              <  1:11284628/137‑1 (MQ=255)
                                                  aGGTCTTGCCGACTGTGGTGGGTTCT‑GTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCtgt                                                                              <  1:4012521/137‑1 (MQ=255)
                                                  aGGTCTTGCCGACTGTGGTGGGTTCT‑GTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCtgt                                                                              <  2:4901955/137‑1 (MQ=255)
                                                  aGGTCTTGCCGACTGTGGTGGGTTCT‑GTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCtgt                                                                              <  2:7351522/137‑1 (MQ=255)
                                                  aGGTCTTGCCGACTGTGGTGGGTTCT‑GTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCtgt                                                                              <  2:1334676/137‑1 (MQ=255)
                                                    gTCTTTCCGACTGTGTTGGGTTCT‑GTGGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGGGACGGCTgtgg                                                                            <  2:10790231/137‑1 (MQ=38)
                                                    gTCTTGCCGATTGTGGTGGGTTCT‑GTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTgtgg                                                                            <  2:1184320/137‑1 (MQ=255)
                                                    gTCTTGCCGACTGTGGTGGGTTCT‑GTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGTCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTgtgg                                                                            <  1:16824692/137‑1 (MQ=255)
                                                    gTCTTGCCGACTGTGGTGGGTTCT‑GTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTgtgg                                                                            <  1:5730841/137‑1 (MQ=255)
                                                    gTCTTGCCGACTGTGGTGGGTTCT‑GTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTgtgg                                                                            <  2:1930548/137‑1 (MQ=255)
                                                    gTCTTGCCGACTGTGGTGGGTTCT‑GTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTgtgg                                                                            <  1:7942308/137‑1 (MQ=255)
                                                      cTTGCCGACTGTGGTGGGTTCT‑GTTGTGTTGTCGGTGGGCATTGTTTGGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGGGCGTGACGGCTgtggtg                                                                          <  1:17549221/137‑1 (MQ=255)
                                                       ttGCCGACTGTGGTGGGTTCT‑GTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCACGCCCGCTCCCACAAGTACTGCGCAGTGCGTGACGGCTGTGGTGa                                                                         >  2:7671924/1‑137 (MQ=255)
                                                            gACTGTGGTGGGTTCT‑GTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGAAAGCCCGCTCCCACAAATGACGTGTAATGGGTGACGGAGgtggtggaaatg                                                                    >  1:10266405/1‑131 (MQ=255)
                                                                 tggtggGTTCT‑GTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCACGCCCGCTcc                                                                                                            >  1:12960439/1‑92 (MQ=255)
                                                                 tggtggGTTCT‑GTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCACGCCCGCTcc                                                                                                            <  2:12960439/92‑1 (MQ=255)
                                                                      ggTTCT‑GTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGGAGCggg                                                          >  2:4599097/1‑137 (MQ=255)
                                                                                 tgttgtCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGAccc                                                                      <  2:3351834/115‑1 (MQ=255)
                                                                                     gtCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑TCCTTTTCGTGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGATGGTTGTGGTGACCCTGTGGGAGTGGGCGTGCCCGCGAaga                                            <  2:8879847/137‑1 (MQ=37)
                                                                                      tCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACt                                                                                                  <  2:7186444/82‑1 (MQ=255)
                                                                                          ggggCATTGTTTTGTTTCCCTGGATTCTGGGTTGCTTGGGCTG‑GCCTCTTCGTGGGCACGCCCGCTCCCACAGGTATGGTGCGGTGCGTGATGGTTGTGGTGACCTTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCt                                       <  2:2140807/136‑1 (MQ=25)
                                                                                           gggCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCACGCCCACTCCCCCAAGGACTGTGTAGTGCGTGAAGGATGTGGTGGACCTGTGGGGGCGGGCGTGTCCACGAAGAGGACTg                                      >  1:8766136/1‑137 (MQ=21)
                                                                                           gggCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCACGCCCACTCCCACAGGTACTGTGCAGTGCGTTACGGATGTGGTGACCCTGTGGGAGTGGGCGTTTCCACGAAGAGGCCTg                                      >  1:2155735/1‑137 (MQ=37)
                                                                                                     ttttttCCTTGGGTTTTTGTTTTCTTGGGCTGCC‑CTCTTCGCGACCACGCCCGCCCCCACAGGGACTGGGTAGTGGGTGAGGGTTGTGGGGACTTTGTGGGGGCGGGTGTGCCCGCGAAGAGGCCTGACCTGCTAAt                            <  2:16928278/134‑1 (MQ=255)
                                                                                                     ttttttCCTTGGATTTTGGGTGGCCTGGGCGGTC‑CTCTTCGTGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGTTTTGGTGACCCTGTGGGAGCGGGCGTTCCCGCGAAGAGGCCTGACCTGCTAAt                            <  1:17261206/134‑1 (MQ=21)
                                                                                                     ttttttCCCTGGATTTTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCACGCCCGCTCCAACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCTGACCTGCTAAt                            <  2:9610244/134‑1 (MQ=37)
                                                                                                     ttgtttCCCTGGGTTCTGGTTTGCCTGGGCTG‑GCCCCTTCGTGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGGGACGGCGGTGGTGACCCTGTGGGAGTGGGCGTGCCCGCGAAGAGGCCTGACCTGCTAAt                            <  2:11124247/137‑1 (MQ=25)
                                                                                                     ttgtttCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGGAGCGGGCGTGCCTGCGAAGAGGCCTGACCTGCTAAt                            <  2:577653/137‑1 (MQ=255)
                                                                                                         ttCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGa                                                                                     >  2:5447862/1‑76 (MQ=255)
                                                                                                         ttCCCTGGATTCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGa                                                                                     <  1:5447862/76‑1 (MQ=255)
                                                                                                           tccTGGATTCTGGGTTGCCTGGGTTG‑GCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGGCGGGTGTGGTGACCCTGTGGGAGTGGGTGTGCCCTCGAAGAGGCCTGACCTGCTAATGCAAgg                      >  2:11689206/2‑137 (MQ=37)
                                                                                                                   tttGGGTTGCCTGGGTTGAC‑CTCTTAGCGGGCACGCCCGCTCCCACAGGTATTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCTGACCTGCTAATGCAAGGTGTGAGGc              <  1:5093165/135‑1 (MQ=37)
                                                                                                                   tCTGGGTTGCCTGGGCTG‑GCCTCTTCGCGGACACGCCCGCTCCCAAAGGTGCTGTGCAGTGCGTGACGGCTGTGGTGAACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCGACCTGCTCATGCAAGGTTTggggt              >  2:14773134/1‑136 (MQ=25)
                                                                                                                   tCTGGGTTGCCTGGGCTG‑GCCTCTTCGAGGGCACACCCGCTCCCACAGGTACTGTGCAGTGTGTGACGGCTGTGGCGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCTGACCTGCTAATGCAAGGTGTGAGGc              <  1:10019997/137‑1 (MQ=37)
                                                                                                                        gTTGCCTGGGCTG‑GCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGAccc                                                                      <  1:6396316/76‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                 
GGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCT‑GTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTG‑G‑CTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCTGACCTGCTAATGCAAGGTGTGAGGCTCTGTGGGTAAA  >  AE015451_phaM1/1419416‑1419677

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑