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.94.3% 1764.1 / 74.5 393intergenic (+42/+711)gltA CDS (citrate synthase)/yeiW CDS (putative oxidoreductase)–/–
Reads supporting (aligned to +/- strand):  ref base T (0/22);  new base . (202/161);  total (210/183)
Fisher's exact test for biased strand distribution p-value = 3.88e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 496 total aligned reads displayed.

GCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCAC  >  AE015451_phaM1/4740673‑4740937
                                                                                                                                    |                                                                                                                                    
gCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAATGCTGCCCTCGTGCAGCCttttt                                                                                                                                  <  2:2761074/137‑1 (MQ=255)
gCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttt                                                                                                                                  <  1:6797636/137‑1 (MQ=255)
gCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttt                                                                                                                                  <  2:11527329/137‑1 (MQ=255)
gCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttt                                                                                                                                  <  2:2177433/137‑1 (MQ=255)
    tcCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttatt                                                                                                                              >  1:650544/1‑134 (MQ=255)
    tcCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCAGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttatt                                                                                                                              >  2:17254087/1‑134 (MQ=255)
     ccAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTTCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttattt                                                                                                                             <  2:16581922/137‑5 (MQ=255)
     ccAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttattt                                                                                                                             <  1:5163500/137‑5 (MQ=255)
     ccAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttattt                                                                                                                             <  1:18463103/137‑5 (MQ=255)
     ccAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttattt                                                                                                                             <  2:14274819/137‑5 (MQ=255)
     ccAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttattt                                                                                                                             <  2:16680281/137‑5 (MQ=255)
     ccAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttattt                                                                                                                             <  2:18640767/137‑5 (MQ=255)
     ccAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttattt                                                                                                                             <  1:13944191/137‑5 (MQ=255)
     ccAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttattt                                                                                                                             <  2:3566749/137‑5 (MQ=255)
     ccAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttattt                                                                                                                             <  1:10017416/137‑5 (MQ=255)
      cAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttatttg                                                                                                                            <  1:17753443/137‑6 (MQ=255)
       aGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttatttgt                                                                                                                           <  1:10784953/137‑7 (MQ=255)
                 aGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCgg                                                                                                                >  2:8531471/1‑137 (MQ=255)
                     cGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCctc                                                                                                            <  2:12979809/137‑1 (MQ=255)
                       gCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCtt                                                                                                          <  2:11945550/137‑1 (MQ=255)
                        ccGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTc                                                                                                         >  2:2879385/1‑137 (MQ=255)
                          gTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTcgc                                                                                                       >  2:1827414/1‑137 (MQ=255)
                          gTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTcgc                                                                                                       >  1:15605673/1‑137 (MQ=255)
                          gTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTcgc                                                                                                       >  2:3988338/1‑137 (MQ=255)
                          gTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTcgc                                                                                                       >  2:8216793/1‑137 (MQ=255)
                          gTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTcgc                                                                                                       >  2:14874999/1‑137 (MQ=255)
                          gTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTcgc                                                                                                       >  2:12889342/1‑137 (MQ=255)
                          gTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCAGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCtttgg                                                                                                       >  1:10304577/1‑134 (MQ=255)
                                      tGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTTGGGCAGCC‑TTTTTTATTTGTCCTGT‑‑‑‑GCCTCTTCGCGGCCTTGCCCGCGAag                                                                                       <  1:3521340/137‑1 (MQ=255)
                                          cacCGGCGAGCAGAATCGCGACATCGTTACCCTGAAAGACCGCAAGGAAGCTGCAGTAGCAGAACGAAAGAAGGCTGCCCTCGGGAAGCC‑TTTTTTATTTGTCCTGTGCCGTCCTCTTCGCGGCCTTGCCCGCGAag                                                                                       <  2:8074125/137‑1 (MQ=255)
                                          cacCGGCGAGCAGAAGCTCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAGGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCCCTTCGCGGCCTTGCCCGCGAag                                                                                       <  2:9740599/137‑1 (MQ=255)
                                           acCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAaga                                                                                      <  2:903848/137‑1 (MQ=255)
                                           acCGGCGAGCAGAAACGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTACCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAaga                                                                                      <  2:13769786/137‑1 (MQ=255)
                                             cGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACAGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAgg                                                                                    >  1:19344322/1‑137 (MQ=255)
                                                  aGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGANGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCAGc                                                                               >  1:17960773/1‑137 (MQ=38)
                                                  aGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGc                                                                               >  2:9210522/1‑137 (MQ=255)
                                                  aGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGc                                                                               >  2:5611254/1‑137 (MQ=39)
                                                  aGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCACGGCCCTGCCCACGAAGAGGACgga                                                                               >  1:12307412/1‑136 (MQ=255)
                                                   gCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGca                                                                              >  1:4667325/1‑137 (MQ=255)
                                                   gCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCAGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCAGCCTTTCCCGCTCAGAAGCCAGca                                                                              >  1:4244508/1‑137 (MQ=255)
                                                   gCAGAAGCGCGACATCGCTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCATCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGca                                                                              >  2:4514620/1‑137 (MQ=38)
                                                    cAGAAGCTCAACATCTTTGCCCTGAAAGACCGCAGGGAAGCTGAAGTAGCAAAACACAAGAAGGCTACCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGcac                                                                             <  1:16505578/137‑1 (MQ=255)
                                                    cAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAGGGCTGCCCTTGGGCAG‑CTTTTTTTATTTGTCCTGTTCCGTCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGcac                                                                             <  1:1580569/137‑1 (MQ=255)
                                                    cAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAGGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTTTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGcac                                                                             <  2:6797966/137‑1 (MQ=38)
                                                    aagagGCGCGACGTCGTTGCCTTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAGGAGGGTTGCCCTCGGGCAGCC‑TTTTTTGTTTTTCCTGTGCCGTCCTTTTCGTGGCCTTGCCCGCGAAGAGGCCGGcac                                                                             <  1:3765178/136‑1 (MQ=255)
                                                     aGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGcaca                                                                            >  1:7314779/1‑137 (MQ=39)
                                                     aGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGcaca                                                                            >  2:3200191/1‑137 (MQ=255)
                                                          cgcgACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAAAGGCCGGCACAAacac                                                                       >  1:18746379/1‑137 (MQ=38)
                                                            cgACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TCTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGAcacca                                                                     <  2:15141292/137‑1 (MQ=255)
                                                               cATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGCAGAGGCCGGCACAGACACCACaa                                                                  >  2:4863262/1‑137 (MQ=255)
                                                                 ttgTTGCCCTGAAAGACCGCAATTAAGCGGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACaaaa                                                                <  1:13270626/135‑1 (MQ=38)
                                                                 tctttGCCCGGAAAGACCGCAGGTAAGCTGCAGTAGCCGAACACAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTTTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACaaaa                                                                <  1:8550469/134‑1 (MQ=255)
                                                                 tCGTTGCCTTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCTTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACaaaa                                                                <  2:15653912/137‑1 (MQ=255)
                                                                 tCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACaaaa                                                                <  2:7934280/137‑1 (MQ=255)
                                                                     tGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTc                                                            <  2:7255601/137‑1 (MQ=255)
                                                                     tGCCCGGAAAGACCGCAAGGAATCTGCAGGAGCCGAACGAAAGAAGGCTGCCCTTGGGAAG‑CTTTTTTTGTTTGTCTTGTACCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTc                                                            <  1:200192/137‑1 (MQ=37)
                                                                     tGCCCGGAAAGACCAAAGGGAAGTTGCAGTAGCAGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTc                                                            <  1:7144448/137‑1 (MQ=255)
                                                                       cccTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTTCCCGCGAAAAGGCCGGCACAGACACCACAAAACCTCaa                                                          >  1:14316404/1‑136 (MQ=255)
                                                                       cccTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGACCTGCCCGCGAAAAGGCCGGCAAAGAAAACACAAAACCCCaa                                                          >  1:17404399/1‑136 (MQ=255)
                                                                              aGACCGCAAGTAGGCGGCAGTAGCCGAACGAAAGAAGGCTGCCTTCGGGCAG‑CTTTTTTTGTTTTTCTTGTGCCGGCCTTTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTtct                                                   <  1:3590430/137‑1 (MQ=255)
                                                                              aGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTtct                                                   <  2:1255998/137‑1 (MQ=255)
                                                                               gACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAAACTCAGTTTTtctc                                                  >  2:10335965/1‑137 (MQ=255)
                                                                               gACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCCTCCCGGCCCTGCCCACCAAAAAGCCAGCACAAAAACCCCACAAACCCAATTTTtctc                                                  >  1:2437215/1‑137 (MQ=255)
                                                                               gACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAG‑CATTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGGCGGCACAGACACCACAAAACCTCCGTTTTtctc                                                  >  1:9995844/1‑137 (MQ=255)
                                                                                aCCGCAAGTAAGCTGCGGTAGCAGAAAGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTGGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTcc                                                 <  2:11043551/137‑1 (MQ=255)
                                                                                  cGCAAGTAAGCGGCAGTAGCCAAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTccgc                                               <  1:2330057/137‑1 (MQ=255)
                                                                                    cAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGAAACCACCAAAACCCAGTTTTTCTccaccc                                             >  1:19151973/1‑137 (MQ=37)
                                                                                     aaGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCTCGGCCTTGCCCTCGAAAAGGACGGCACAGACACCACAAAACCCCAATTTTTCTcctccct                                            >  1:5057515/1‑136 (MQ=255)
                                                                                     aaGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCCCAAACACCACAAAACCTCAGTTTTTCTccgcccg                                            >  2:13366322/1‑137 (MQ=255)
                                                                                     aaGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTccgcccg                                            >  2:6398751/1‑137 (MQ=255)
                                                                                     aaGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCAGGCACAGACAACCCAAAACCTCAGTTTTTCTccgcccg                                            >  2:1898918/1‑137 (MQ=255)
                                                                                     aaGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCCCTTCGCGGCCTTGCCCGCTAAGAGGACGGCACAGACACCACAAAAACTCAGTTTTTCTccgcccg                                            >  1:2761908/1‑137 (MQ=255)
                                                                                       gTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCt                                          <  2:8827825/137‑1 (MQ=255)
                                                                                        tAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCTGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTc                                         <  1:12951447/137‑1 (MQ=38)
                                                                                           gCTGCGGTAGCAGAACGAAAGAAGGTTGCCCTTGGGCAGCC‑TTTTTTGTTTGTCTTGTGCCGGCCTTTTTGTGGCTTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACg                                      <  1:1059292/137‑1 (MQ=255)
                                                                                           gCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCTGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACg                                      <  1:13527048/137‑1 (MQ=255)
                                                                                           gCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACg                                      <  2:9462691/137‑1 (MQ=255)
                                                                                           gCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACg                                      <  1:12028882/137‑1 (MQ=255)
                                                                                           gCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACg                                      <  1:4051892/137‑1 (MQ=255)
                                                                                            cTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGc                                     <  2:2303000/137‑1 (MQ=255)
                                                                                            cTGCAGTAGCCGAAAACAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTTTCCTGTGCCTGCCTCTTCGCGGCTTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGc                                     <  1:17003419/137‑1 (MQ=255)
                                                                                              gCAGTAGCAGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TGGTGTATTTGTCCTGTGCCGGCCTCTTCAGTGTCTTGCCCGTGAAGAGGACGGCACAGCCAACACCAAACGTCAGTTTTTCTCCGCCCGCTCcagcag                                   >  2:14609246/1‑137 (MQ=255)
                                                                                                aGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGcc                                 <  1:13393343/137‑1 (MQ=255)
                                                                                                  tAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCCCTTCGCGGCCTTTTCCGCGAAGAGACCAGAACAAACACCCCAAAAACCCAGTTTTTCTTCTCCCACTCACACAACCCt                               >  1:5825709/1‑137 (MQ=255)
                                                                                                   aGCCGAACGCAAGAAGGTTGCCCTCGGGTAGCC‑TTTTTGATTTGTCTTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCTGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCtt                              <  2:6132808/137‑1 (MQ=255)
                                                                                                   aGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCtt                              <  2:10648025/137‑1 (MQ=39)
                                                                                                   aGCCGAACACAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCTCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCtt                              <  1:19455625/137‑1 (MQ=255)
                                                                                                       gAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTTCCCACGAAGAGGCCAGCACAGACACCCCAAAACCTCAGTTTTTCTCCACCCGCTCACACAGCCCTTTcca                          >  1:15215663/1‑135 (MQ=255)
                                                                                                        aaCGCAAGAATGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAgg                         <  1:8859914/137‑1 (MQ=255)
                                                                                                        aaCGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAACCCTTTCAgg                         >  2:7163770/1‑137 (MQ=255)
                                                                                                            aaagaagGCTTCCCTCGGGCA‑TCTTTTTTTGTTTTTCCTGTGCTGTCCTTTTTGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTAt                     <  1:8247797/136‑1 (MQ=255)
                                                                                                             aaaaagGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGTCCGCGAAGAGGCCGGAACAGACACCACAAAAACTCCATTTTTCTCCGCCCGCTCCCGCCGCCCTTTCAGGGTAtt                    >  1:112029/4‑137 (MQ=255)
                                                                                                              agaagGCTGCCCTCGTGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTg                   <  2:17796160/137‑1 (MQ=255)
                                                                                                                  ggCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACg               >  2:16415500/1‑137 (MQ=39)
                                                                                                                  ggCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCCTGCCCACGAAGAGGCCGGCACAAACACCACAAAAACCCAGTTTTTCTCCGCCCGCTCACGCAACCCCTTCAGGGTATTGAACg               >  1:14567275/1‑137 (MQ=255)
                                                                                                                         cTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGAAACCACAAAAACCCAGTTTTTCTCCGCCCGCCCACGCAGCCCTTTCAAGGTATTTAACGGCGCATc        >  1:7621969/1‑137 (MQ=37)
                                                                                                                          tCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCTCGGCCTTGCCCGCGAAGAGGCCGGAACAAACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATcc       >  1:3905822/1‑137 (MQ=255)
                                                                                                                             ggCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATccacc    <  2:14773490/137‑1 (MQ=255)
                                                                                                                              gCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATccacca   <  2:12250558/137‑1 (MQ=255)
                                                                                                                               cAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAAACACCACAAAACCCCCGTTTTTCTCCTCCCGCTCCCGCAACCCCTTCAGGGTATTGAACGGTGTATCcaccac  >  1:5393906/1‑137 (MQ=25)
                                                                                                                                    |                                                                                                                                    
GCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCAC  >  AE015451_phaM1/4740673‑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: ‑