Predicted mutation
evidence seq id position mutation annotation gene description
RA AE015451_phaM1 4,741,232 C→T intergenic (+469/+284) 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,741,2320CT85.0% 272.3 / 29.9 100intergenic (+469/+284)gltA CDS (citrate synthase)/yeiW CDS (putative oxidoreductase)–/–
Reads supporting (aligned to +/- strand):  ref base C (9/6);  new base T (46/39);  total (55/45)
Fisher's exact test for biased strand distribution p-value = 7.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑AACC‑TCC‑GCCAAC‑TGCAATTCGTGAA‑CCAAGCAAAAGGCCCCGAAAAGGCCCTCC‑GCAGCCCCCCCAACCTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGT  >  AE015451_phaM1/4741187‑4741343
                                              |                                                                                                                    
gAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACgg                       >  2:9168592/1‑137 (MQ=255)
gAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACgg                       <  1:1495419/137‑1 (MQ=255)
 aaCTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGt                      >  1:5805458/1‑137 (MQ=255)
 aaCTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGt                      >  2:6605224/1‑137 (MQ=255)
 aaCTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGt                      >  2:10478250/1‑137 (MQ=255)
 aaCTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGt                      >  2:10478291/1‑137 (MQ=255)
 aaCTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGt                      >  2:14358006/1‑137 (MQ=255)
 aaCTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGt                      >  1:13363157/1‑137 (MQ=255)
 aaCTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGt                      >  1:18602202/1‑137 (MQ=255)
 aaCTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGt                      >  2:7960132/1‑137 (MQ=255)
  aCTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTc                     <  2:5517157/137‑1 (MQ=255)
  aCTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTc                     <  1:7137337/137‑1 (MQ=255)
  aCTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTc                     <  1:18942677/137‑1 (MQ=255)
  aCTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTc                     <  1:16988682/137‑1 (MQ=255)
  aCTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGC‑TTCTGCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTc                     <  1:3306998/137‑1 (MQ=255)
  aCTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGC‑TTC‑TGCGACTTGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACATGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTc                     <  2:12568127/137‑1 (MQ=255)
   cTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTca                    <  2:18878761/137‑1 (MQ=255)
   cTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTca                    <  1:1145707/137‑1 (MQ=255)
   cTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTca                    <  1:17722279/137‑1 (MQ=255)
   cTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAA‑TCCAGGTACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTca                    <  1:18961090/137‑1 (MQ=255)
      ggTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCacac                 >  1:6338841/1‑137 (MQ=255)
      ggTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCacac                 >  2:5195707/1‑137 (MQ=255)
      ggTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCacac                 >  2:10675537/1‑137 (MQ=255)
      ggTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCacac                 >  2:5870818/1‑137 (MQ=255)
      ggTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCacac                 >  2:14428294/1‑137 (MQ=255)
      ggTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCacac                 >  2:6745464/1‑137 (MQ=255)
      ggTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCacac                 >  1:15859094/1‑137 (MQ=255)
      ggTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGC‑TTC‑TGCGAC‑TTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCacac                 >  2:5870778/1‑137 (MQ=255)
       gTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACAcc                >  1:10263982/1‑137 (MQ=255)
         atatttCCAGCTCTGGTCGCCTTCGTGCTTGAG‑CAGC‑TTC‑TGCGAC‑TTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGATACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgc              <  2:10356196/132‑1 (MQ=255)
         ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTTGCCAGCGGCGCCTTGATCACGGTCACACcgc              <  1:2315996/137‑1 (MQ=255)
         ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgc              >  1:7528140/1‑137 (MQ=255)
         ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgc              <  1:8414392/137‑1 (MQ=255)
         ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgc              >  2:7294476/1‑137 (MQ=255)
         ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgc              <  1:7987193/137‑1 (MQ=255)
         ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgc              <  1:7987160/137‑1 (MQ=255)
         ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgc              <  1:9131656/137‑1 (MQ=255)
         ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgc              <  2:4978183/137‑1 (MQ=255)
         ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgc              >  2:10741043/1‑137 (MQ=255)
         ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgc              <  2:14052537/137‑1 (MQ=255)
         ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgc              <  2:4306392/137‑1 (MQ=255)
         ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgc              <  2:2642543/137‑1 (MQ=255)
         ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgc              >  2:17438991/1‑137 (MQ=255)
         ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgc              >  2:17438992/1‑137 (MQ=255)
         ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgc              >  2:18526243/1‑137 (MQ=255)
         ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgc              >  2:18526287/1‑137 (MQ=255)
         ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgc              >  2:7913346/1‑137 (MQ=255)
         ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgc              >  1:17341335/1‑137 (MQ=255)
         ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgc              >  1:17452356/1‑137 (MQ=255)
         ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgc              <  1:17050563/137‑1 (MQ=255)
         ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgc              >  2:7993254/1‑137 (MQ=255)
         ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgc              <  1:14678053/137‑1 (MQ=255)
         ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgc              >  1:13970889/1‑137 (MQ=255)
         ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgc              <  1:13271468/137‑1 (MQ=255)
         ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgc              <  1:13191722/137‑1 (MQ=255)
         ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgc              >  1:2225172/1‑137 (MQ=255)
         ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgc              >  1:3497055/1‑137 (MQ=255)
         ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgc              >  1:4658548/1‑137 (MQ=255)
         ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgc              >  1:1196630/1‑137 (MQ=255)
         ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgc              >  1:7340470/1‑137 (MQ=255)
         ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgc              >  1:6580696/1‑137 (MQ=255)
         ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgc              <  1:4187319/137‑1 (MQ=255)
         ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgc              <  1:1865541/137‑1 (MQ=255)
         ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CCGC‑TTC‑TGCGAC‑TTGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCAAGCGGCGCCTTGATCACGGCCACACcgc              >  1:7528179/1‑137 (MQ=255)
         ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGC‑TTC‑TGCGACTTGC‑ATTCGTGAATCCATGC‑ACATGCCCCGGACACGTCCTCCTACAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgc              <  2:7730138/137‑1 (MQ=255)
         ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGC‑TTC‑TGCGAC‑TTGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgc              >  1:9511285/1‑137 (MQ=255)
         ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGC‑TTC‑TGCGAC‑TTGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgc              >  2:17438954/1‑137 (MQ=255)
          taGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgcg             >  2:3172020/1‑137 (MQ=255)
          taGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgcg             <  2:1354999/137‑1 (MQ=255)
            ggTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGAt           >  2:4827788/1‑137 (MQ=255)
              tCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCTTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAg         <  2:4377187/137‑1 (MQ=255)
              tCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAg         >  1:13177999/1‑137 (MQ=255)
              tCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAg         >  1:11269786/1‑137 (MQ=255)
              tCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAg         >  1:7646517/1‑137 (MQ=255)
              tCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAg         <  2:6321729/137‑1 (MQ=255)
              tCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAg         <  1:17715009/137‑1 (MQ=255)
              tCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAg         >  1:18666551/1‑137 (MQ=255)
              tCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAg         <  2:1266666/137‑1 (MQ=255)
              tCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAg         <  2:1091589/137‑1 (MQ=255)
              tCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAg         >  1:5516723/1‑137 (MQ=255)
              tCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAg         >  1:7646467/1‑137 (MQ=255)
              tCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAA‑TCCAGGAACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAg         <  1:14713601/137‑1 (MQ=255)
              tCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAA‑TCCAGGAACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAg         <  1:9266605/137‑1 (MQ=255)
              tCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGC‑TTC‑TGCGAC‑TTGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAg         >  2:1551267/1‑137 (MQ=255)
              tCCAGCTCTGGTCACCTTCGTACTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGGCCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAg         <  2:15013714/137‑1 (MQ=255)
               ccAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAgg        <  1:200713/137‑1 (MQ=255)
                cAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGc       <  2:9759664/137‑1 (MQ=255)
                cAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGc       <  1:7931931/137‑1 (MQ=255)
                cAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGc       <  1:2078294/137‑1 (MQ=255)
                cAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGc       <  1:13872157/137‑1 (MQ=255)
                cAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGTGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGc       >  2:8016632/1‑137 (MQ=255)
                cAGCTCTGGTCACCTTCGTGCTTGAG‑CAGC‑TTC‑TGCGAC‑TTGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGc       <  2:4693737/137‑1 (MQ=255)
                cAGCTCTGGTCACCTTCGTGCTTGAG‑CAGC‑TTC‑TGCGAC‑TTGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGc       <  1:8363108/137‑1 (MQ=255)
                cAGCTCTGGTCACCTTCGTGCTTGAACAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGc       <  2:10015644/137‑1 (MQ=255)
                   ctctGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCtt    >  1:4833220/1‑137 (MQ=255)
                   ctctGGTCACCTTCGTGCTTGAG‑CAGC‑TTC‑TGCGAC‑TTGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTAACACCGCGATAGGCCtt    >  2:18892969/1‑137 (MQ=255)
                   ctctGGTCACCTTCGTGCTTGAG‑CAGC‑TTC‑TGCGAC‑TNGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCtt    >  1:18302307/1‑137 (MQ=255)
                     ctGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGt  >  1:12702846/1‑137 (MQ=255)
                     ctGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGTGCCTTGATCACGGTCACACCGCGATAGGCCTTGt  >  2:16923628/1‑137 (MQ=255)
                     ctGGTCACCTTCGTGCTTGAG‑CAGC‑TTC‑TGCGAC‑TTGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGt  >  1:15946509/1‑137 (MQ=255)
                                              |                                                                                                                    
GAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑AACC‑TCC‑GCCAAC‑TGCAATTCGTGAA‑CCAAGCAAAAGGCCCCGAAAAGGCCCTCC‑GCAGCCCCCCCAACCTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGT  >  AE015451_phaM1/4741187‑4741343

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑