Predicted mutation
evidence seq id position mutation annotation gene description
RA AE015451_phaM1 4,740,819 2 bp→CT intergenic (+56/+696) 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,8190TC96.0% 1317.1 / ‑5.1 380intergenic (+56/+697)gltA CDS (citrate synthase)/yeiW CDS (putative oxidoreductase)–/–
Reads supporting (aligned to +/- strand):  ref base T (0/15);  new base C (177/188);  total (177/203)
Fisher's exact test for biased strand distribution p-value = 7.17e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*AE015451_phaM14,740,8200GT99.2% 1280.7 / NA 380intergenic (+57/+696)gltA CDS (citrate synthase)/yeiW CDS (putative oxidoreductase)–/–
Reads supporting (aligned to +/- strand):  ref base G (0/3);  new base T (177/197);  total (177/203)
Fisher's exact test for biased strand distribution p-value = 2.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 380 total aligned reads displayed.

AAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCAC  >  AE015451_phaM1/4740689‑4740952
                                                                                                                                  ||                                                                                                                                    
 aGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCgg                                                                                                                               >  2:8200209/1‑137 (MQ=255)
     cGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAACTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCctc                                                                                                                           <  2:5303617/137‑1 (MQ=255)
      ggCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCctct                                                                                                                          <  1:8894677/137‑1 (MQ=255)
      ggCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCctct                                                                                                                          <  1:11485559/137‑1 (MQ=255)
          gTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTctc                                                                                                                      >  1:16377830/1‑135 (MQ=255)
          gTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTcgc                                                                                                                      >  2:6733616/1‑137 (MQ=255)
          gTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTcgc                                                                                                                      >  2:1739270/1‑137 (MQ=255)
          gTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCCTcgc                                                                                                                      >  1:9734782/1‑137 (MQ=255)
            ccGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCgg                                                                                                                    <  1:6610789/137‑1 (MQ=255)
            ccGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCgg                                                                                                                    <  2:2941961/137‑1 (MQ=255)
              gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGTCGGACTCTTTGCGGcc                                                                                                                  >  1:15413048/1‑137 (MQ=255)
               cgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGGCCCTTCGCGGACt                                                                                                                 >  1:1170653/1‑137 (MQ=255)
               cgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCt                                                                                                                 >  1:15297448/1‑137 (MQ=255)
                 ccAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTg                                                                                                               <  2:4110902/137‑1 (MQ=255)
                 ccAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTg                                                                                                               <  2:1415575/137‑1 (MQ=255)
                          cacCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTGTTTGTCCT‑GTGCCGGCCTCTTCGCGGCCTTGCCCGCGAag                                                                                                      <  2:6150983/137‑1 (MQ=255)
                          cacCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAag                                                                                                      <  2:13492827/137‑1 (MQ=255)
                          cacCGGCAAGCAGAAGCGCGACGTCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTTGGGCAGCC‑TTTTTGATTTGTCCTGTGCCGGCCTCTTTGTGGCCTTGCCCGCGAag                                                                                                      <  2:11599590/137‑1 (MQ=255)
                            ccGGCGAGCAGAAGCGCGACATCGTGGCCCTGAAAGACGGCAAGCAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTTTGCCGGCCTCTTCGCGGCCTTTCCCGCGAagag                                                                                                    <  2:6865963/137‑1 (MQ=255)
                             cGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCCCTTCGCGGCCTTGCCCACAAAAAgg                                                                                                   >  1:4143082/1‑137 (MQ=255)
                                  aGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCgga                                                                                              >  2:2609464/1‑136 (MQ=255)
                                  aGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGc                                                                                              >  2:15889376/1‑137 (MQ=39)
                                  aGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGc                                                                                              >  2:9266423/1‑137 (MQ=255)
                                  aGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCCCAGAGAGGCCGGc                                                                                              >  1:10394349/1‑137 (MQ=255)
                                  aGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACACAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCATCTTCGCGGCCTTGCCCGCGAAGAGGCCGGc                                                                                              >  2:11982861/1‑137 (MQ=255)
                                  aGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCAGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGc                                                                                              >  2:9604022/1‑137 (MQ=255)
                                    cAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGTCCTTGCCCGCGAAGAGGCCGGcac                                                                                            <  1:13144201/137‑1 (MQ=255)
                                     aGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTTCCCGCGAAGAGGCCGGcaca                                                                                           >  1:13283323/1‑137 (MQ=255)
                                     aGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGcaca                                                                                           >  2:1967558/1‑137 (MQ=255)
                                     aGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGcaca                                                                                           >  1:9013958/1‑137 (MQ=255)
                                     aGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCAGCCTTGCCCGCGACGAGGCCGGcaca                                                                                           >  2:10281756/1‑137 (MQ=255)
                                     aGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCCTCGCGGCCTTGCCCGCGAAGAGGCCAGcaca                                                                                           >  1:7513357/1‑137 (MQ=255)
                                     aGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGCAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTTTCCTGTGCCGGCCTCTTCGTGGCCTTGCCCGCGAAGAGGCCGGcaca                                                                                           >  2:18779796/1‑137 (MQ=255)
                                     aGAAACGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGGAGCCGAACGCAAGAAGGCTGCCTTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGTGAAGAGGCCGGcaca                                                                                           <  1:12122333/137‑1 (MQ=255)
                                       aaGCGCGACATTGTTGCCCTGAAAGACCGCAAGGAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGa                                                                                         <  1:18767402/137‑1 (MQ=255)
                                            cgACATTTTTGCCCTGAAAGACCGCAAGGAATCTGCGGTAGCCGAACGCAAGAAGGCTCCCCTTGGGCAGCC‑TTTTTTGTTTGTCCTGTGCTGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGAcacca                                                                                    <  1:10953088/137‑1 (MQ=25)
                                               cATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACaa                                                                                 >  2:16470485/1‑137 (MQ=255)
                                               cATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGATGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACaa                                                                                 >  1:3292054/1‑137 (MQ=38)
                                                 tCGTTGCCCTGAAAGACCGCAAGTAAGCTGCNGTAGCCGAACGCAATAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTTGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACaaaa                                                                               <  2:3772391/137‑1 (MQ=255)
                                                 tCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACaaaa                                                                               <  2:13283323/137‑1 (MQ=255)
                                                 tCGTTGCCATGAAAGACCACAATTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCTTTGCCCGCGAAGAGGCCGGCACAGACACCACaaaa                                                                               <  2:9191884/137‑1 (MQ=255)
                                                     ttcccGGAAAGACCGAAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTGTTTGTCCTGTGCCGGCCTCTTCGTGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTc                                                                           <  1:14244468/135‑1 (MQ=255)
                                                               gACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCCCAGACACCACAAAACCTCAGTTTTtctc                                                                 >  1:12731015/1‑137 (MQ=39)
                                                               gACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTtctc                                                                 >  2:10703451/1‑137 (MQ=255)
                                                               gACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGACGGCACAAACACCACAAAAACCCAGTTTTtctc                                                                 >  1:17548420/1‑137 (MQ=255)
                                                               gACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGACTTGCCCGCGAAGAGGCCGGCACAGACACCACCAAACCTCAGTTTTtctc                                                                 >  1:3117481/1‑137 (MQ=38)
                                                                  cGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTccgc                                                              <  2:15051738/137‑1 (MQ=255)
                                                                    cAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGACTTGCCCGCGAAGAGGCCGGCACAGACACCACAACAACCCAGTTTTTCTccgccc                                                            >  1:8183369/1‑137 (MQ=38)
                                                                     aaGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCTCGGTCTTTCCCCCTAAGAAGCAGGCACAGACACCACAAAACCCCCATTTTTTTccccccc                                                           >  1:5127741/1‑136 (MQ=255)
                                                                     aaGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGACTTGCCCGCGAAGAGGCCAGCACAGACACCACAAAAAACCAATTTTTCTccgcccg                                                           >  1:5327201/1‑137 (MQ=37)
                                                                       gTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTGTTTGTCCTGTTCCGGCTTCTTCGCGGCCTTTCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCt                                                         <  1:7039070/137‑1 (MQ=255)
                                                                         aaGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGTGGACTTGCCCGCGAAGAAGACGGCACAGACACCCCAAAACCTCAATTTTTCTCCGCCCGCTcc                                                       >  1:4925581/1‑136 (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)
                                                                           gCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCTGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTCTTCTCCGCCCGCCCAcc                                                     <  2:12682200/137‑2 (MQ=255)
                                                                           gCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑CTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACg                                                     <  2:1040675/137‑1 (MQ=255)
                                                                           gCTGCAGTAGCCGAACACAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTTCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACg                                                     <  1:12888518/137‑1 (MQ=255)
                                                                            cTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGc                                                    <  2:5517925/137‑1 (MQ=255)
                                                                            cTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGc                                                    <  2:13062221/137‑1 (MQ=255)
                                                                            cTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGGGCCGGCCTTTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGc                                                    <  2:11412455/137‑1 (MQ=255)
                                                                               cAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGACTTGCCCGCGAAGAAGCCGGTACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAga                                                 >  1:1590726/1‑136 (MQ=255)
                                                                                aGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCCCTTCGCGGCCTTGCCCGCGAAAAGGACGGCACAGAAACCACAAAACCTCAATTTTTCTCCGCCCGCTCACACAAcc                                                >  1:6082338/1‑137 (MQ=255)
                                                                                 gTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCCCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGccc                                               <  1:18534759/137‑1 (MQ=255)
                                                                                  tAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAAACTCAGTTTTTCTCCGCCCGCTTACGCAGCCCt                                              >  2:1282865/1‑137 (MQ=38)
                                                                                  tAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCCTGCCCACGAAGAAGCCGGCACAGACACCACAAAACATCAGTTTTTCTTCGCCCGCTCACGCAGCCCt                                              >  2:5055988/1‑137 (MQ=255)
                                                                                   ggacGAACGCAAGAGGGCTGCCTTCGGGCAGCC‑TTTTTTGTTTGTCCTGTTCTTGCCTCTTCGCGGCCTTTCCCGCGAGGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCtt                                             <  1:1000671/134‑1 (MQ=255)
                                                                                   aGCCGAACGCAAGAAGTCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCtt                                             <  1:16712556/137‑1 (MQ=255)
                                                                                   aGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTGTTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCtt                                             <  2:146969/137‑1 (MQ=255)
                                                                                   aGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCtt                                             <  1:2590141/137‑1 (MQ=255)
                                                                                   aGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCtt                                             <  2:17846437/137‑1 (MQ=255)
                                                                                   aGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCtt                                             <  2:8819195/137‑1 (MQ=255)
                                                                                   aGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCtt                                             <  2:12144982/137‑1 (MQ=255)
                                                                                   aGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCtt                                             <  1:14548394/137‑1 (MQ=255)
                                                                                   aGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCtt                                             <  2:6094451/137‑1 (MQ=255)
                                                                                   aGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCtt                                             <  2:11848915/137‑1 (MQ=255)
                                                                                   aGCCGAACGCAAGAAGGCTGCCCTCGGGCAG‑CTTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCtt                                             <  2:11205431/137‑1 (MQ=255)
                                                                                   aGCCGAACGAAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTGTTTGTCCTGTGCCGTCCTTTTCGCGGCCTTTCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCtt                                             <  1:1931448/137‑1 (MQ=255)
                                                                                       gAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCATTGCCCGCGAAGAGGCCGGCACAAGCACCACCAAACCCCAGTTTTTCTCAGACCGCTCAAGCCACCCTTTCAg                                         >  2:12297560/1‑137 (MQ=255)
                                                                                        aaCGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGACGGAACAGACACCACAAAAACTCAGTTTTTTTCCGCCCGCTCACACAACCCTTTCAgg                                        >  1:12528035/1‑137 (MQ=255)
                                                                                          cGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGACTCTTCGCGGCCTTGCCCACGAAGAGGCCGGCACAGACACCACCAAAAATCAATTTTTCTCCGCCCACTCAAACCGCCCCTTCAAGGt                                      >  1:8410259/1‑137 (MQ=255)
                                                                                            cAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGTCTCTTCGCGGCCTTGCCCGCGAAGAGGCCAGCACAGAAACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCCTTCAGGGTAt                                    >  1:15585102/1‑137 (MQ=255)
                                                                                            cAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCCCAAAACCTAAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTAt                                    >  2:6415075/1‑137 (MQ=255)
                                                                                              agaagGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTg                                  <  2:16241314/137‑1 (MQ=255)
                                                                                                aatgCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGCCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTCCTCCGCCCGCTCACGCAGCCCTCTCAGGGTATTGaa                                <  2:8082246/134‑1 (MQ=255)
                                                                                                  ggCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCCTGCCCGCGAAGAGGCCCGCACAGACACCCCAAAACCTCAGTTTTTCTCCGCCCGCTCACACAGCACTTTCAGGGTATTGAACg                              >  1:2338573/1‑137 (MQ=255)
                                                                                                  ggCTGCCCTCGGGCAG‑CTTTTTTTATTTGTCCTTTGCCGGCCTCTTCGCGGCCTGGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACg                              <  1:3715360/137‑1 (MQ=255)
                                                                                                     tGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGgcg                           <  1:2842765/137‑1 (MQ=255)
                                                                                                      gCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGgcgc                          <  2:7292406/137‑1 (MQ=255)
                                                                                                        ccTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGGAGAGGCCGGCACAGACACCACAAAAACTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCGAGGTATTGAACGGCGCAt                        >  2:5798913/1‑137 (MQ=255)
                                                                                                          tCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTTACGCAGCCCTTTCAGGGTATTGAACGGCGCAttc                      >  2:8515867/1‑135 (MQ=255)
                                                                                                          tCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATcc                      >  2:12409080/1‑137 (MQ=255)
                                                                                                          tCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGAAACCACAAAACCTCCATTTTTCTCCACCCGCTCCCACAGCCCTTTCAGGGTATTGGACAGCGCATcc                      >  1:18271363/1‑137 (MQ=255)
                                                                                                          ccGGGCAGCC‑TTTTTTCTTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCCGGGTATTGAACGGTGCATcc                      >  2:17068948/2‑137 (MQ=255)
                                                                                                           cGGGCAG‑CTTTTTTTATTTTTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATcca                     <  1:17062969/137‑1 (MQ=255)
                                                                                                            gggCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGTCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCGCACGCAGACCTTTCAGGGTATTGAACTGTGCATccac                    >  2:344977/1‑137 (MQ=255)
                                                                                                               cAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCAGCACAGACACCACAAAAACTCAGTTTTTCTCCGCCCACTCACGCAGCCCTTTCAGGGTATTGAACGGGGCATCcaccac                 >  2:5684793/1‑137 (MQ=255)
                                                                                                                    cttttttATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAAc             <  1:17347929/136‑1 (MQ=255)
                                                                                                                        tttATTTGTCCTTTGCTGGCCTTTTCGTGGCTTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACttgt         <  1:3590148/137‑1 (MQ=255)
                                                                                                                                  ||                                                                                                                                    
AAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCAC  >  AE015451_phaM1/4740689‑4740952

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