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)
                 aGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCgg                                                                                                                >  2:8531471/1‑137 (MQ=255)
                     cGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCctc                                                                                                            <  2:12979809/137‑1 (MQ=255)
                       gCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCtt                                                                                                          <  2:11945550/137‑1 (MQ=255)
                        ccGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTc                                                                                                         >  2:2879385/1‑137 (MQ=255)
                          gTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTcgc                                                                                                       >  2:1827414/1‑137 (MQ=255)
                          gTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTcgc                                                                                                       >  1:15605673/1‑137 (MQ=255)
                          gTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTcgc                                                                                                       >  2:3988338/1‑137 (MQ=255)
                          gTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTcgc                                                                                                       >  2:8216793/1‑137 (MQ=255)
                          gTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTcgc                                                                                                       >  2:14874999/1‑137 (MQ=255)
                          gTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTcgc                                                                                                       >  2:12889342/1‑137 (MQ=255)
                          gTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCAGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCtttgg                                                                                                       >  1:10304577/1‑134 (MQ=255)
                                      tGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTTGGGCAGCCTTTTTTATTTGTCCTGT‑‑‑‑GCCTCTTCGCGGCCTTGCCCGCGAag                                                                                       <  1:3521340/137‑1 (MQ=255)
                                          cacCGGCGAGCAGAATCGCGACATCGTTACCCTGAAAGACCGCAAGGAAGCTGCAGTAGCAGAACGAAAGAAGGCTGCCCTCGGGAAGCCTTTTTTATTTGTCCTGTGCCGTCCTCTTCGCGGCCTTGCCCGCGAag                                                                                       <  2:8074125/137‑1 (MQ=255)
                                          cacCGGCGAGCAGAAGCTCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAGGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCCCTTCGCGGCCTTGCCCGCGAag                                                                                       <  2:9740599/137‑1 (MQ=255)
                                           acCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAaga                                                                                      <  2:903848/137‑1 (MQ=255)
                                           acCGGCGAGCAGAAACGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTACCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAaga                                                                                      <  2:13769786/137‑1 (MQ=255)
                                             cGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACAGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAgg                                                                                    >  1:19344322/1‑137 (MQ=255)
                                                  aGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGANGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCAGc                                                                               >  1:17960773/1‑137 (MQ=38)
                                                  aGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGc                                                                               >  2:9210522/1‑137 (MQ=255)
                                                  aGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGc                                                                               >  2:5611254/1‑137 (MQ=39)
                                                  aGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCACGGCCCTGCCCACGAAGAGGACgga                                                                               >  1:12307412/1‑136 (MQ=255)
                                                   gCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGca                                                                              >  1:4667325/1‑137 (MQ=255)
                                                   gCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCAGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCAGCCTTTCCCGCTCAGAAGCCAGca                                                                              >  1:4244508/1‑137 (MQ=255)
                                                   gCAGAAGCGCGACATCGCTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCATCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGca                                                                              >  2:4514620/1‑137 (MQ=38)
                                                    cAGAAGCTCAACATCTTTGCCCTGAAAGACCGCAGGGAAGCTGAAGTAGCAAAACACAAGAAGGCTACCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGcac                                                                             <  1:16505578/137‑1 (MQ=255)
                                                    cAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAGGGCTGCCCTTGGGCAGCTTTTTTTATTTGTCCTGTTCCGTCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGcac                                                                             <  1:1580569/137‑1 (MQ=255)
                                                    cAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAGGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTTTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGcac                                                                             <  2:6797966/137‑1 (MQ=38)
                                                    aagagGCGCGACGTCGTTGCCTTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAGGAGGGTTGCCCTCGGGCAGCCTTTTTTGTTTTTCCTGTGCCGTCCTTTTCGTGGCCTTGCCCGCGAAGAGGCCGGcac                                                                             <  1:3765178/136‑1 (MQ=255)
                                                     aGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGcaca                                                                            >  1:7314779/1‑137 (MQ=39)
                                                     aGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGcaca                                                                            >  2:3200191/1‑137 (MQ=255)
                                                          cgcgACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAAAGGCCGGCACAAacac                                                                       >  1:18746379/1‑137 (MQ=38)
                                                            cgACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTCTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGAcacca                                                                     <  2:15141292/137‑1 (MQ=255)
                                                               cATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGCAGAGGCCGGCACAGACACCACaa                                                                  >  2:4863262/1‑137 (MQ=255)
                                                                 ttgTTGCCCTGAAAGACCGCAATTAAGCGGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACaaaa                                                                <  1:13270626/135‑1 (MQ=38)
                                                                 tctttGCCCGGAAAGACCGCAGGTAAGCTGCAGTAGCCGAACACAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTTTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACaaaa                                                                <  1:8550469/134‑1 (MQ=255)
                                                                 tCGTTGCCTTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCTTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACaaaa                                                                <  2:15653912/137‑1 (MQ=255)
                                                                 tCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACaaaa                                                                <  2:7934280/137‑1 (MQ=255)
                                                                     tGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTc                                                            <  2:7255601/137‑1 (MQ=255)
                                                                     tGCCCGGAAAGACCGCAAGGAATCTGCAGGAGCCGAACGAAAGAAGGCTGCCCTTGGGAAGCTTTTTTTGTTTGTCTTGTACCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTc                                                            <  1:200192/137‑1 (MQ=37)
                                                                     tGCCCGGAAAGACCAAAGGGAAGTTGCAGTAGCAGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTc                                                            <  1:7144448/137‑1 (MQ=255)
                                                                       cccTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTTCCCGCGAAAAGGCCGGCACAGACACCACAAAACCTCaa                                                          >  1:14316404/1‑136 (MQ=255)
                                                                       cccTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGACCTGCCCGCGAAAAGGCCGGCAAAGAAAACACAAAACCCCaa                                                          >  1:17404399/1‑136 (MQ=255)
                                                                              aGACCGCAAGTAGGCGGCAGTAGCCGAACGAAAGAAGGCTGCCTTCGGGCAGCTTTTTTTGTTTTTCTTGTGCCGGCCTTTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTtct                                                   <  1:3590430/137‑1 (MQ=255)
                                                                              aGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTtct                                                   <  2:1255998/137‑1 (MQ=255)
                                                                               gACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAAACTCAGTTTTtctc                                                  >  2:10335965/1‑137 (MQ=255)
                                                                               gACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCCTCCCGGCCCTGCCCACCAAAAAGCCAGCACAAAAACCCCACAAACCCAATTTTtctc                                                  >  1:2437215/1‑137 (MQ=255)
                                                                               gACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCATTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGGCGGCACAGACACCACAAAACCTCCGTTTTtctc                                                  >  1:9995844/1‑137 (MQ=255)
                                                                                aCCGCAAGTAAGCTGCGGTAGCAGAAAGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTGGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTcc                                                 <  2:11043551/137‑1 (MQ=255)
                                                                                  cGCAAGTAAGCGGCAGTAGCCAAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTccgc                                               <  1:2330057/137‑1 (MQ=255)
                                                                                    cAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGAAACCACCAAAACCCAGTTTTTCTccaccc                                             >  1:19151973/1‑137 (MQ=37)
                                                                                     aaGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCTCGGCCTTGCCCTCGAAAAGGACGGCACAGACACCACAAAACCCCAATTTTTCTcctccct                                            >  1:5057515/1‑136 (MQ=255)
                                                                                     aaGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCCCAAACACCACAAAACCTCAGTTTTTCTccgcccg                                            >  2:13366322/1‑137 (MQ=255)
                                                                                     aaGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTccgcccg                                            >  2:6398751/1‑137 (MQ=255)
                                                                                     aaGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCAGGCACAGACAACCCAAAACCTCAGTTTTTCTccgcccg                                            >  2:1898918/1‑137 (MQ=255)
                                                                                     aaGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCCCTTCGCGGCCTTGCCCGCTAAGAGGACGGCACAGACACCACAAAAACTCAGTTTTTCTccgcccg                                            >  1:2761908/1‑137 (MQ=255)
                                                                                       gTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCt                                          <  2:8827825/137‑1 (MQ=255)
                                                                                        tAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCTGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTc                                         <  1:12951447/137‑1 (MQ=38)
                                                                                           gCTGCGGTAGCAGAACGAAAGAAGGTTGCCCTTGGGCAGCCTTTTTTGTTTGTCTTGTGCCGGCCTTTTTGTGGCTTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACg                                      <  1:1059292/137‑1 (MQ=255)
                                                                                           gCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCTGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACg                                      <  1:13527048/137‑1 (MQ=255)
                                                                                           gCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACg                                      <  2:9462691/137‑1 (MQ=255)
                                                                                           gCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACg                                      <  1:12028882/137‑1 (MQ=255)
                                                                                           gCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACg                                      <  1:4051892/137‑1 (MQ=255)
                                                                                            cTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGc                                     <  2:2303000/137‑1 (MQ=255)
                                                                                            cTGCAGTAGCCGAAAACAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTTTCCTGTGCCTGCCTCTTCGCGGCTTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGc                                     <  1:17003419/137‑1 (MQ=255)
                                                                                              gCAGTAGCAGAACGCAAGAAGGCTGCCCTCGGGCAGCCTGGTGTATTTGTCCTGTGCCGGCCTCTTCAGTGTCTTGCCCGTGAAGAGGACGGCACAGCCAACACCAAACGTCAGTTTTTCTCCGCCCGCTCcagcag                                   >  2:14609246/1‑137 (MQ=255)
                                                                                                aGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGcc                                 <  1:13393343/137‑1 (MQ=255)
                                                                                                  tAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCCCTTCGCGGCCTTTTCCGCGAAGAGACCAGAACAAACACCCCAAAAACCCAGTTTTTCTTCTCCCACTCACACAACCCt                               >  1:5825709/1‑137 (MQ=255)
                                                                                                   aGCCGAACGCAAGAAGGTTGCCCTCGGGTAGCCTTTTTGATTTGTCTTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCTGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCtt                              <  2:6132808/137‑1 (MQ=255)
                                                                                                   aGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCtt                              <  2:10648025/137‑1 (MQ=39)
                                                                                                   aGCCGAACACAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCTCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCtt                              <  1:19455625/137‑1 (MQ=255)
                                                                                                       gAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTTCCCACGAAGAGGCCAGCACAGACACCCCAAAACCTCAGTTTTTCTCCACCCGCTCACACAGCCCTTTcca                          >  1:15215663/1‑135 (MQ=255)
                                                                                                        aaCGCAAGAATGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAgg                         <  1:8859914/137‑1 (MQ=255)
                                                                                                        aaCGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAACCCTTTCAgg                         >  2:7163770/1‑137 (MQ=255)
                                                                                                            aaagaagGCTTCCCTCGGGCATCTTTTTTTGTTTTTCCTGTGCTGTCCTTTTTGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTAt                     <  1:8247797/136‑1 (MQ=255)
                                                                                                             aaaaagGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGTCCGCGAAGAGGCCGGAACAGACACCACAAAAACTCCATTTTTCTCCGCCCGCTCCCGCCGCCCTTTCAGGGTAtt                    >  1:112029/4‑137 (MQ=255)
                                                                                                              agaagGCTGCCCTCGTGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTg                   <  2:17796160/137‑1 (MQ=255)
                                                                                                                  ggCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACg               >  2:16415500/1‑137 (MQ=39)
                                                                                                                  ggCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCCTGCCCACGAAGAGGCCGGCACAAACACCACAAAAACCCAGTTTTTCTCCGCCCGCTCACGCAACCCCTTCAGGGTATTGAACg               >  1:14567275/1‑137 (MQ=255)
                                                                                                                         cTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGAAACCACAAAAACCCAGTTTTTCTCCGCCCGCCCACGCAGCCCTTTCAAGGTATTTAACGGCGCATc        >  1:7621969/1‑137 (MQ=37)
                                                                                                                          tCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCTCGGCCTTGCCCGCGAAGAGGCCGGAACAAACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATcc       >  1:3905822/1‑137 (MQ=255)
                                                                                                                             ggCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATccacc    <  2:14773490/137‑1 (MQ=255)
                                                                                                                              gCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATccacca   <  2:12250558/137‑1 (MQ=255)
                                                                                                                               cAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAAACACCACAAAACCCCCGTTTTTCTCCTCCCGCTCCCGCAACCCCTTCAGGGTATTGAACGGTGTATCcaccac  >  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: