Predicted mutation
evidence seq id position mutation annotation gene description
RA AE015451_phaM1 4,740,811 (T)7→6 intergenic (+48/+705) gltA CDS (citrate synthase) → / ← yeiW CDS (putative oxidoreductase) –/–

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451_phaM14,740,8050T.95.8% 1660.9 / 46.4 361intergenic (+42/+711)gltA CDS (citrate synthase)/yeiW CDS (putative oxidoreductase)–/–
Reads supporting (aligned to +/- strand):  ref base T (0/15);  new base . (163/177);  total (169/192)
Fisher's exact test for biased strand distribution p-value = 8.22e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 439 total aligned reads displayed.

AAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCAC  >  AE015451_phaM1/4740669‑4740937
                                                                                                                                        |                                                                                                                                    
    gCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttt                                                                                                                                  <  2:9581396/137‑1 (MQ=255)
    gCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttt                                                                                                                                  <  1:363268/137‑1 (MQ=255)
    gCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttt                                                                                                                                  <  1:11037004/137‑1 (MQ=255)
    gCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttt                                                                                                                                  <  1:4926400/137‑1 (MQ=255)
    gCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttt                                                                                                                                  <  2:17530958/137‑1 (MQ=255)
    gCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttt                                                                                                                                  <  2:6095606/137‑1 (MQ=255)
    gCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttt                                                                                                                                  <  1:14719546/137‑1 (MQ=255)
    gCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCGGTACACCGGCGAGCAGAAGCGCGGCATCGTTGCCCTGAAAGATCGCAAGTAAGCTGTAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttt                                                                                                                                  <  2:4986958/137‑1 (MQ=255)
     ctctCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCtttttt                                                                                                                                 <  1:18772005/137‑1 (MQ=255)
     ctctCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCtttttt                                                                                                                                 <  1:5482679/137‑1 (MQ=255)
        tcCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCNGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCANTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttatt                                                                                                                              >  1:15608146/1‑134 (MQ=255)
        tcCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCANAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttatt                                                                                                                              >  1:15891386/1‑134 (MQ=255)
        tcCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttatt                                                                                                                              >  1:18674/1‑134 (MQ=255)
        tcCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttatt                                                                                                                              >  1:18098848/1‑134 (MQ=255)
        tcCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttatt                                                                                                                              >  2:12736662/1‑134 (MQ=255)
         ccAGCCCGTACAAGATCTGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGGAGGCTGCCCTCTGGCAGCCttttttcttt                                                                                                                             <  2:14149872/137‑5 (MQ=255)
         ccAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAGAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttattt                                                                                                                             <  1:16450895/137‑5 (MQ=255)
         ccAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttattt                                                                                                                             <  2:6953756/137‑5 (MQ=255)
         ccAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttattt                                                                                                                             <  1:17874886/137‑5 (MQ=255)
         ccAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttattt                                                                                                                             <  1:3019431/137‑5 (MQ=255)
         ccAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttattt                                                                                                                             <  1:2900787/137‑5 (MQ=255)
                          ggCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCctct                                                                                                           <  1:8894677/137‑1 (MQ=255)
                          ggCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCctct                                                                                                           <  1:11485559/137‑1 (MQ=255)
                              gTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTTAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGACAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTcgc                                                                                                       >  2:11484399/1‑137 (MQ=255)
                              gTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCCCTTcgc                                                                                                       >  1:10634849/1‑137 (MQ=255)
                              gTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAG‑CATTTTTTATTTGTCCTGTGCCGGCCTCTTcgc                                                                                                       >  2:8140762/1‑137 (MQ=255)
                              gTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCAGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTcgc                                                                                                       >  2:4153899/1‑137 (MQ=255)
                                ccGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCgg                                                                                                     <  2:11376322/137‑1 (MQ=255)
                                ccGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCgg                                                                                                     <  1:446731/137‑1 (MQ=255)
                                  gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGTCGGACTCTTTGCGGcc                                                                                                   >  1:15413048/1‑137 (MQ=255)
                                   cgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGGCCCTTCGCGGACt                                                                                                  >  1:1170653/1‑137 (MQ=255)
                                   cgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCt                                                                                                  <  2:5174447/137‑1 (MQ=255)
                                   cgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCt                                                                                                  >  1:15297448/1‑137 (MQ=255)
                                        gCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGccc                                                                                             >  2:1935553/1‑137 (MQ=255)
                                          tGTACACCGGCAATCAGACGCTCGACATCGTTCCCCCGAACGCCCGCAATTAATCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTTTCCTGTGCCGGCCTCTCCGCGACCTTCCCcgc                                                                                           <  2:7766399/137‑1 (MQ=255)
                                              cacCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTGTTTGTCCT‑GTGCCGGCCTCTTCGCGGCCTTGCCCGCGAag                                                                                       <  2:6150983/137‑1 (MQ=255)
                                              cacCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAag                                                                                       <  2:13492827/137‑1 (MQ=255)
                                                      aGCAGACGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGACCGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGGCCCTTCTCTGCCCTGGCCTCGAAAAGACCAGc                                                                               >  1:7046953/1‑137 (MQ=255)
                                                      aGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCgga                                                                               >  2:2609464/1‑136 (MQ=255)
                                                      aGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCgga                                                                               >  2:10996132/1‑136 (MQ=255)
                                                      aGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGc                                                                               >  2:5466882/1‑137 (MQ=255)
                                                      aGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCAGCATTGCCCACGAAGAGGACggg                                                                               >  1:239338/1‑136 (MQ=255)
                                                        cAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGTCCTTGCCCGCGAAGAGGCCGGcac                                                                             <  1:13144201/137‑1 (MQ=255)
                                                        aaGAAGCGCGACATCGTTCCCCTGAAAGACCGCAAGTAAGCGGCGGTAGCCGAACGCAAGAAGGCTGCCCTTGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGcac                                                                             <  1:14305567/136‑1 (MQ=255)
                                                         aGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGGACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGcaca                                                                            >  1:2071524/1‑137 (MQ=38)
                                                         aGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTTCCCGCGAAGAGGCCGGcaca                                                                            >  1:13283323/1‑137 (MQ=255)
                                                         aGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGcaca                                                                            >  2:15871525/1‑137 (MQ=255)
                                                         aGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGcaca                                                                            >  2:3514009/1‑137 (MQ=255)
                                                         aGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGcaca                                                                            >  1:9013958/1‑137 (MQ=255)
                                                         aGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAAAGGCCGGcaca                                                                            >  2:12988893/1‑137 (MQ=38)
                                                         aGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGCAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTTTCCTGTGCCGGCCTCTTCGTGGCCTTGCCCGCGAAGAGGCCGGcaca                                                                            >  2:18779796/1‑137 (MQ=255)
                                                         aGAAACGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGGAGCCGAACGCAAGAAGGCTGCCTTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGTGAAGAGGCCGGcaca                                                                            <  1:12122333/137‑1 (MQ=255)
                                                           aaGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGa                                                                          >  2:16116651/1‑137 (MQ=39)
                                                                cgACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCCCTTCTCTGCCCTTCCCCCTAAGAGGACGGCCCAAAaaccc                                                                     >  1:6679274/1‑136 (MQ=255)
                                                                   cATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGATGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACaa                                                                  >  1:3292054/1‑137 (MQ=38)
                                                                     tCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACaaaa                                                                <  2:11737075/137‑1 (MQ=255)
                                                                     tCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACaaaa                                                                <  2:13283323/137‑1 (MQ=255)
                                                                     tCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACACAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACaaaa                                                                <  2:12569344/137‑1 (MQ=39)
                                                                             cTGAAAGACCGCAAGTAAGTTGCAGTAGCCTAACGCAAGAAGGTTGCCCGCGGGCAGCC‑TTTTTTATTTGTCCGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGtt                                                        <  2:12151177/137‑1 (MQ=255)
                                                                                   gACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCCCAGACACCACAAAACCTCAGTTTTtctc                                                  >  1:12731015/1‑137 (MQ=39)
                                                                                   gACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTtctc                                                  >  2:6864439/1‑137 (MQ=255)
                                                                                      cGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCTGCCTTTCCCGCGAAGAGGCCGGCACAGACACCACAAAACCCCCGTTTTTCTccgc                                               >  1:17100196/1‑137 (MQ=255)
                                                                                        cAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCTGCCTTTCCCGCGAAGAGGACGGCACAGACACCACCAAACCTCAGTTTTTCTccgccc                                             >  1:16113607/1‑137 (MQ=255)
                                                                                         aaGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCTCGGCCTTTTCCCCTAAGAGGCCGGTACCAAAACCCCAAAACCTCCGTTTTTCTcccccca                                            >  1:5327152/1‑136 (MQ=255)
                                                                                         aaGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGACAGCC‑TTCTTTATTTGTCCTGTGCCGGCCCCTTCTCGGCCTTGCCCGCGCAGAGGCCGGCACAGCCACCACAAAACCTCAGGTTTTCTccgcccg                                            >  2:562118/1‑137 (MQ=255)
                                                                                            tAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAG‑CTTTTTTTGTTTGTCTTGTGCCGGCCTCTTCGTGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTc                                         <  1:1842512/137‑1 (MQ=255)
                                                                                            gaaGCTGCAGGAGCAGAACGCGAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTc                                         <  1:17493138/136‑1 (MQ=255)
                                                                                             aaGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCCCAGACACCACAAAACCTCAGTTTTTCTCCGCACGCTCa                                        >  2:18172993/1‑137 (MQ=255)
                                                                                              aGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCATTTTTTCTCCGCCCGCTCAc                                       <  2:2359871/137‑1 (MQ=255)
                                                                                              aGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCAc                                       <  2:15605306/137‑1 (MQ=255)
                                                                                               gcggCAGTAGCCGAACACAAGAAGGCGGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTTGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACg                                      <  1:15013526/134‑1 (MQ=255)
                                                                                               gCTGCATGAGCCGAACGCAAGAATGCTGCCCTCGTGTAG‑CTTTTTTTATGTTTCTTGTTCCGGCCCCTTCGCGGCTTTGCCCGTGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACg                                      <  1:9796784/137‑1 (MQ=255)
                                                                                               gCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACg                                      <  2:8667780/137‑1 (MQ=255)
                                                                                               gCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACg                                      <  2:10377002/137‑1 (MQ=255)
                                                                                               gCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACg                                      <  2:647587/137‑1 (MQ=255)
                                                                                               gCTGCAGTAGCCGAACACAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTTCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACg                                      <  1:12888518/137‑1 (MQ=255)
                                                                                                cTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCActc                                     <  2:8426099/137‑3 (MQ=255)
                                                                                                cTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGc                                     <  2:1231447/137‑1 (MQ=255)
                                                                                                cTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGc                                     <  2:3870366/137‑1 (MQ=255)
                                                                                                    aGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCCCTTCGCGGCCTTGCCCGCGAAAAGGACGGCACAGAAACCACAAAACCTCAATTTTTCTCCGCCCGCTCACACAAcc                                 >  1:6082338/1‑137 (MQ=255)
                                                                                                       ggacGAACGCAAGAGGGCTGCCTTCGGGCAGCC‑TTTTTTGTTTGTCCTGTTCTTGCCTCTTCGCGGCCTTTCCCGCGAGGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCtt                              <  1:1000671/134‑1 (MQ=255)
                                                                                                       aGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTGTTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCtt                              <  2:146969/137‑1 (MQ=255)
                                                                                                       aGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCtt                              <  2:17846437/137‑1 (MQ=255)
                                                                                                       aGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTGGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCtt                              <  2:9237248/137‑1 (MQ=255)
                                                                                                       aGCCAAACACAAGAAGGCTGCCCTCGGGAAGCC‑TTTTTTGTTTGTCCTGTGCCGGCCTTTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCtt                              <  1:10266813/137‑1 (MQ=255)
                                                                                                       aGCAGAACGAAAGAGGGCTGCCCTTGGGCAGCC‑TTTTTTGTTTGTCTTGTGCCGTCCTCTTTGCGGCCTTTCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCtt                              <  1:2776944/137‑1 (MQ=255)
                                                                                                         ccGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTc                            <  2:1498764/137‑1 (MQ=255)
                                                                                                           gAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAAAGGCCGGCACAGACACCACAAAACCCCAGTTTTTCTCCGCCCACTCACGCAGCCCTTTCAg                          >  1:8024540/1‑137 (MQ=38)
                                                                                                            aaCGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAgg                         <  1:4054368/137‑1 (MQ=255)
                                                                                                            aaCGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGACGGAACAGACACCACAAAAACTCAGTTTTTTTCCGCCCGCTCACACAACCCTTTCAgg                         >  1:12528035/1‑137 (MQ=255)
                                                                                                              cGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGACTCTTCGCGGCCTTGCCCACGAAGAGGCCGGCACAGACACCACCAAAAATCAATTTTTCTCCGCCCACTCAAACCGCCCCTTCAAGGt                       >  1:8410259/1‑137 (MQ=255)
                                                                                                                cAAGAAGGGTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCATTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTAt                     <  2:8982435/137‑1 (MQ=255)
                                                                                                                cAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTAt                     <  1:8379598/137‑1 (MQ=255)
                                                                                                                            ccTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCTCGGCCTTTCCCTCGAAGAGGACGGAACAGACACCACAAAACCCCAATTTTTTTCCGCCCCCTCCCGCCACCCCTTCCGGGGATTTAACGGCGCAt         >  1:16760839/1‑137 (MQ=255)
                                                                                                                            ccTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGGAGAGGCCGGCACAGACACCACAAAAACTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCGAGGTATTGAACGGCGCAt         >  2:5798913/1‑137 (MQ=255)
                                                                                                                            ccTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAATTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCAt         >  2:11924473/1‑137 (MQ=39)
                                                                                                                             cTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCAGCACAGACACCACAAAAACTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATc        >  2:1640065/1‑137 (MQ=255)
                                                                                                                              tCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATcc       >  2:10624663/1‑137 (MQ=255)
                                                                                                                              tCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATcc       >  2:3142907/1‑137 (MQ=255)
                                                                                                                              ccGGGCAGCC‑TTTTTTCTTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCCGGGTATTGAACGGTGCATcc       >  2:17068948/2‑137 (MQ=255)
                                                                                                                                        |                                                                                                                                    
AAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCAC  >  AE015451_phaM1/4740669‑4740937

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