Predicted mutation
evidence seq id position mutation annotation gene description
RA AE015451_phaM1 4,740,816 T→G intergenic (+53/+700) 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,8160TG94.2% 1188.5 / 3.1 343intergenic (+53/+700)gltA CDS (citrate synthase)/yeiW CDS (putative oxidoreductase)–/–
Reads supporting (aligned to +/- strand):  ref base T (1/19);  new base G (194/129);  total (195/148)
Fisher's exact test for biased strand distribution p-value = 7.17e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 343 total aligned reads displayed.

AGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTC‑GCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTG  >  AE015451_phaM1/4740690‑4740944
                                                                                                                              |                                                                                                                                 
aGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGCCCGCAAGTAAGCTGCAGTAGCCGAACGCAAGCAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCgg                                                                                                                        >  1:7348495/1‑137 (MQ=255)
aGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAATGCTGCCCTCTGGCAGCC‑TTTTTTATTTGTCCTGTGCCgg                                                                                                                        <  1:11616935/137‑1 (MQ=255)
    cGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCctc                                                                                                                    <  1:6539566/137‑1 (MQ=255)
    cGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCctc                                                                                                                    <  1:4056620/137‑1 (MQ=255)
    cGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCctc                                                                                                                    <  2:6149811/137‑1 (MQ=255)
       ccGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTc                                                                                                                 >  2:11018479/1‑137 (MQ=255)
         gTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGGAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTc‑gc                                                                                                              >  2:3258091/1‑137 (MQ=255)
         gTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTctt                                                                                                               >  1:17357742/1‑135 (MQ=255)
         gTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTctc                                                                                                               >  1:6079626/1‑135 (MQ=255)
         gTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTc‑gc                                                                                                              >  2:9401773/1‑137 (MQ=255)
         gTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTc‑gc                                                                                                              >  1:3708854/1‑137 (MQ=255)
         gTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTc‑gc                                                                                                              >  2:439883/1‑137 (MQ=255)
         gTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCCCTTc‑gc                                                                                                              >  1:1035162/1‑137 (MQ=255)
           ccGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGNCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTC‑GCgg                                                                                                            <  1:16269606/137‑1 (MQ=255)
           ccGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTC‑GCgg                                                                                                            <  1:4323886/137‑1 (MQ=255)
           ccGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTC‑GCgg                                                                                                            <  2:17515121/137‑1 (MQ=255)
              cgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCAGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTC‑GCGGCCt                                                                                                         >  1:11057275/1‑137 (MQ=255)
                     tGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGTCCTCTTC‑GCGGCCTTGCCcgc                                                                                                  <  2:2689572/137‑1 (MQ=255)
                         cacCGGCGAGCAGAAGCTCGACGTCGTTTCCCTGAAAGACCTCAATTAAGCGCCAGTAGCCAAACGCAAGAAGTCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCTTGCCCCTTC‑GCGGCCTTACCCGCGAag                                                                                              <  2:10516049/137‑1 (MQ=255)
                         cacCGGCGAGCAGAAGCGTGACATCGTTACCTTGAAAGACCTCAAGTAAGCTGCAGTAGCCAAATGCAAGAAGGCGGCCCTCGGGAATCC‑TTTTTTATTTGTCTTGTGCCGGCCTCTTT‑GCGGCCTTGCCCGCGAag                                                                                              <  2:2298609/137‑1 (MQ=255)
                         cacCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTC‑GCGGCCTTGCCCGCGAag                                                                                              <  2:2394144/137‑1 (MQ=255)
                                gAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTC‑GCGGCCTTGCCCGCGAAGAGGCCgg                                                                                       >  2:6316884/1‑137 (MQ=39)
                                 aGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTC‑GCGGCCTTGCCCTCGAAGAGGACGGc                                                                                      >  1:14635130/1‑137 (MQ=255)
                                 aGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTC‑GCGGCCTTGCCCGCGGAGAGGCCGGc                                                                                      >  1:11771127/1‑137 (MQ=255)
                                 aGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTC‑GCGGCCTTGCCCGCGAAGAGGCCGGc                                                                                      >  2:11945438/1‑137 (MQ=255)
                                 aGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTC‑GCGGCCTTGCCCGCGAAGAGGCCGGc                                                                                      >  1:4022216/1‑137 (MQ=39)
                                 aGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTC‑GCGGCATTGCCCGCGAAGAGGCCGGc                                                                                      >  2:2061497/1‑137 (MQ=255)
                                  gCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTC‑GCGGCCTTGCCCGCGAAGAGGCCGGca                                                                                     >  1:15486865/1‑137 (MQ=39)
                                  gCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTC‑GCGGCCTTGCCCGCGAAGAGGCCGGca                                                                                     >  2:7204507/1‑137 (MQ=255)
                                   cAGAAGTGCGACATCGTTGCCTTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCTTTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCTGGCCTCTTC‑GTGGCCTTGCCCGCGAAGAGGCCGGcac                                                                                    <  1:6747382/137‑1 (MQ=255)
                                    aGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGTCGGCCTCTTCATCGGCCTT‑CCCACGAAGAGGCCAGcaaa                                                                                   >  1:14419572/1‑135 (MQ=255)
                                    aGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTC‑TTGGCCTTGCCCGCAAAGAGGACAGcaca                                                                                   >  1:8827886/1‑137 (MQ=255)
                                    aGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTC‑TCGGCCCTTCCCGCGAAGAGGCCGGcaca                                                                                   >  1:1071993/1‑137 (MQ=255)
                                    aGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTC‑GCGGCCTTGCCCGCGAAGAGGCCGGcaca                                                                                   >  1:12775985/1‑137 (MQ=255)
                                    aGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTC‑GCGGCCTTGCCCGCGAAGAGGCCGGcaca                                                                                   >  1:1787057/1‑137 (MQ=39)
                                    aGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTC‑GCGGCCTTGCCCGCGAAGAGGCCGGcaca                                                                                   >  1:5562433/1‑137 (MQ=39)
                                    aGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTC‑GCGGCCTTGCCCGCGAAAAGGCCAGaaca                                                                                   >  1:15000678/1‑137 (MQ=255)
                                    aGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTC‑GCGGCCTTGCCCACGAAGAGGCCGGcaca                                                                                   >  1:598431/1‑137 (MQ=255)
                                    aGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTC‑GCGGCATTGCCCGCGAAGAGGCCGGcaca                                                                                   >  2:7537339/1‑137 (MQ=255)
                                    aGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCAGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTC‑GCGGCCTTGCCCGCGAAGAGGCCGGcaca                                                                                   >  1:8452192/1‑137 (MQ=255)
                                      aaGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTACAGTAGCCGAACGCCAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGACTCTTC‑ATGGACCTTCCCAAGAAAAGGCCGGCACAGa                                                                                 >  1:2343932/1‑137 (MQ=255)
                                              cATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTC‑GCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACaa                                                                         >  2:1623327/1‑137 (MQ=255)
                                              cATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTC‑GCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACaa                                                                         >  2:5288193/1‑137 (MQ=255)
                                                tCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTC‑GCGGCCTTTCCCCCGAAGAGGCCAGCACAGAAACCACaaaa                                                                       >  1:6969422/1‑137 (MQ=255)
                                                tCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTC‑GCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACaaaa                                                                       <  2:17039483/137‑1 (MQ=255)
                                                      cccTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTA‑GCGGCCTTGCCCGCGAAGAGGACGGCACAGACACCACAAAAACTCAg                                                                 >  1:4469927/1‑137 (MQ=255)
                                                        cTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTC‑GCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGtt                                                               <  2:9789276/137‑1 (MQ=255)
                                                               aCCGCAAGGAAGCTGCAGTAGCCGAACGAAAGAAGGCTGCCCTAGGGCAGCCTTTTTTGTTTGTCCT‑GTGCCGGCTTCTTC‑GCGGCCTTGCCCGCGAAGNGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTcc                                                        <  1:17032054/137‑1 (MQ=255)
                                                                    aaGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTC‑TCGGCCTTTCCCCCAAAAAGGCCGGCACAGACACCACAAAACCTCAATTTTTCTccgcccg                                                   >  1:6763326/1‑137 (MQ=255)
                                                                    aaGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTC‑GCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACCAAACCTCCGTTTTTCTccgcccg                                                   >  2:8299446/1‑137 (MQ=255)
                                                                    aaGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTC‑GCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTccgcccg                                                   >  1:8685044/1‑137 (MQ=39)
                                                                    aaGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTC‑GCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAAACTCAGTTTTTCTccgcccg                                                   >  2:9397761/1‑137 (MQ=255)
                                                                    aaGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTC‑GCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAAACTCAGTTTTTCTccgcccg                                                   >  2:4135809/1‑137 (MQ=255)
                                                                    aaGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTC‑GCGGCCTTGCCCGCGAAGAGGCCGGCACAGAAACCACAAAAACTCAGTTTTTCTccgcccg                                                   >  2:7811655/1‑137 (MQ=255)
                                                                    aaGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTC‑GCGGACTTGCCCGCAAAGAGGCAGGCACAGAAACCACAAAACCTCCGTTTTTCTccgcccg                                                   >  1:1381773/1‑137 (MQ=255)
                                                                     aGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTC‑GCGGCCTTGCCCGCGAAGAGGCCGGCACAGGCACCACAAAAACTCAGTTTTTCTccacccgc                                                  >  1:12258348/1‑137 (MQ=38)
                                                                        aaGTTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTC‑GCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCAAAAGACATCAGTTTTTCTCCGCCCACTCa                                               >  2:11441420/1‑137 (MQ=255)
                                                                        aaGCTGCAGTAGCCGAACGCAAGAAAGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTC‑GCGGGATTGCCCGGGAAGAGGCCGGCACAGACACCAAAAAAACTCAATTTTTCTCCGCCCACTCa                                               >  1:9068321/1‑137 (MQ=255)
                                                                          gttGCAGTAGCCAAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTGTTTGTCTTGTGCCGTCCTCTTC‑GTGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACg                                             <  1:10338278/135‑1 (MQ=255)
                                                                          gCTGCAGTAGCCGAACGCAAGAATGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTC‑GCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTTAGTTTTTCTCCGCCCGCTCACg                                             <  2:13103494/137‑1 (MQ=38)
                                                                          gCTGCAGTAGCCGAACGCAAGAAGTCTGCCCTCGGGCAGCC‑TTTTTTATTTTTCCTGTGCCGGCCTCTTC‑GCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACg                                             <  1:5105530/137‑1 (MQ=255)
                                                                          gCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTGTTTGTCTTGTGCCGGCCTTTTC‑GCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACg                                             <  1:3595960/137‑1 (MQ=255)
                                                                          gCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTC‑GCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACg                                             <  2:3708854/137‑1 (MQ=255)
                                                                          gCTGCAGTAGCCGAACGAAAGAAGGCTGCCCTCGGGCAG‑CTTTTTTTATTTGTCCTGTGCCGGCCTCTTC‑GCGGCTTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACg                                             <  1:13428938/137‑1 (MQ=255)
                                                                          gCTGCAGGAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTGTTTGTCCTGTGCCGTCCTCTTC‑GTGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACg                                             <  1:12660561/137‑1 (MQ=255)
                                                                           cTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTC‑GCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGc                                            <  2:12412041/137‑1 (MQ=255)
                                                                           cTGCAGTAGCCGAACACAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTC‑GCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGc                                            <  2:15046589/137‑1 (MQ=255)
                                                                                gTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTC‑GCGGCCTTGCCCACGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGccc                                       <  2:9136326/137‑1 (MQ=255)
                                                                                 tAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCctc                                                                                                                    <  1:11318914/60‑1 (MQ=255)
                                                                                 tAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCctc                                                                                                                    >  2:11318914/1‑60 (MQ=255)
                                                                                  ggCCGAACGCAAGAATGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCTTGTGCCGGCCTCTTC‑GCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCtt                                     <  2:13185488/136‑1 (MQ=255)
                                                                                  atccGAACGCAAGAAGGTTGCCCTCGGGAAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTC‑GAGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCtt                                     <  1:10054007/135‑1 (MQ=255)
                                                                                  aGCCGAACGCAAGAAGGCTGCCCTTGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTC‑GCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCtt                                     <  1:8133702/137‑1 (MQ=255)
                                                                                  aGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTC‑GTGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCtt                                     <  2:7058880/137‑1 (MQ=255)
                                                                                  aGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTC‑GCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCtt                                     <  1:2670577/137‑1 (MQ=255)
                                                                                  aGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTC‑GCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCtt                                     <  2:10244367/137‑1 (MQ=255)
                                                                                  aGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTC‑GCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCtt                                     <  2:12109085/137‑1 (MQ=255)
                                                                                  aGCCAAACGCAAGAAGGCTGCCCTCGGGCAACC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTC‑GTGGCCTTGCCCGTGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCtt                                     <  1:15631422/137‑1 (MQ=255)
                                                                                    ccGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTC‑GCGGCATTGCCCGCGAAGAGGCCGGCACCGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCCCGCAGCCCTTTc                                   >  2:5694970/1‑137 (MQ=255)
                                                                                    ccGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATGTGTCCTGTGCCGGCCTCTTC‑GCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTc                                   <  2:14448758/137‑1 (MQ=255)
                                                                                       aaCGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTC‑GCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAgg                                <  2:8811056/137‑1 (MQ=255)
                                                                                           cAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCTTGTGCTGGCCTCTTC‑GCGGCCTTGCCCGCGAGGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTAt                            <  1:1902540/137‑1 (MQ=255)
                                                                                           cAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTC‑TCGGCCTTTCCCCCGAAGAGGCCGGCACAAAAACCCCCAAACCTCCGTTTTTTTCCTCCCGCTCACGCAGCCCTTTCAGGGTAt                            >  1:3269969/1‑137 (MQ=255)
                                                                                             agaagGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTC‑GCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTg                          <  2:8646170/137‑1 (MQ=255)
                                                                                              gagtgCTGCCCTCGTGCAGCC‑TTTTTTATTTGCCCTTTGCCGGCCTCTTC‑GCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTCTTCTCCGCCCGCTCACGCATCCCTTTCAGGGTATTGa                         <  2:1517478/133‑1 (MQ=37)
                                                                                              gaagGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCCCCTC‑GCGGCCCTGCCCACAAAAAAGACAGCACAGAAACCACAAAAACCCAATTTTTTTCCGCCCACCCACACAACCCCTTCAGGGTAtttg                         >  1:544344/1‑135 (MQ=255)
                                                                                                agGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGACTCTTC‑GCGGCCTTGCCCGCGAAGAGGCCGGAACCGAAAACCCAAAAACCCCGTTTTTTTCCACCCCCTCACGCAGCCCTTTCAGGGTATTTaaa                       >  1:14342498/1‑136 (MQ=255)
                                                                                                 ggCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTC‑GCGGCCTTGCCCGCAAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACg                      >  2:6644599/1‑137 (MQ=38)
                                                                                                         tCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTC‑GCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCCGCCCTTTCAGGGTATTGAACGGCGCATcc              >  2:11540979/1‑137 (MQ=255)
                                                                                                         tCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTC‑GCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATcc              >  2:9520557/1‑137 (MQ=255)
                                                                                                           gggCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTC‑GCGGCCTTTCCCGCGAAGAAGCCGGCACAGAAACCACAAAAACTCCGTTTTTCTTCGCCCGCTCAAACAGCCCTTTTAGGGTATTGAACGGCGCATccac            >  1:6834453/1‑137 (MQ=255)
                                                                                                            ggCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTC‑GCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATccacc           <  2:11648688/137‑1 (MQ=255)
                                                                                                                agccttttttATTTGTCCTGTGCCGGCCTCTTC‑GCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCAAAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCCTTAAGGGTATTGAACGGCGCATCCACCACg        >  2:13596163/5‑137 (MQ=38)
                                                                                                                 gcTTTTTTTGTTTTTCCTGTGCTGGCCTCTTC‑GCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGa       <  1:10888339/136‑1 (MQ=39)
                                                                                                                   cttttttATTTGTCCTGTGCCGGCCTCTTC‑GCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAAc     <  1:11987906/136‑1 (MQ=255)
                                                                                                                   cttttttATTTGTCCTGTGCCGGCCTCTTC‑GCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAAc     <  1:8801735/136‑1 (MQ=255)
                                                                                                                    ttttttATTTGTCCTGTGCCGGCCTCTTC‑GCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAAACTCAGTTTTTCTCCGCCCGCTCAAGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACt    >  2:1274906/1‑137 (MQ=39)
                                                                                                                    ttttttATTTGTCCTGTGCCGGCCTCTTC‑GCGGCCTTGCCCGCGAAAAGGCCGGCACAAACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAACCCCTTCAGGGTATTGAACGGCGCATCCCCCACAAACt    >  2:14254857/1‑137 (MQ=37)
                                                                                                                     tttttATTTGTCCTGTGCCGGCCTCTTC‑GCGGCCTTGCCCGCGAAGAAGCCTGCCCAGACACCACAAAACCTCCGTTTTTCTCCGCCCGCCCACGCAGCACTTTCAAGGTATTGAACGGTGCACCCACCACGAACtt   >  2:16710369/1‑137 (MQ=255)
                                                                                                                     tttttATTTGTCCTGTGCCGGCCTCTTC‑GCAGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAAACTCAGTTTTTCTCCGCCCGCTCACGCAGCCCCTTCAGGGTATTGAACGGCGCATCCACCACGAACtt   >  2:4963945/1‑137 (MQ=255)
                                                                                                                              |                                                                                                                                 
AGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTC‑GCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTG  >  AE015451_phaM1/4740690‑4740944

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