Predicted mutation
evidence seq id position mutation annotation gene description
RA AE015451_phaM1 4,741,228 +C intergenic (+465/+288) 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,2281.C85.0% 340.9 / 59.5 100intergenic (+465/+288)gltA CDS (citrate synthase)/yeiW CDS (putative oxidoreductase)–/–
Reads supporting (aligned to +/- strand):  ref base . (9/6);  new base C (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

Only 100 of 115 total aligned reads displayed.

TCTTCTTTCAGGTTACGGGTCAGGCTGCCATCCGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑AACC‑TCC‑GCCAAC‑TGCAATTCGTGAA‑CCAAGCAAAAGGCCCCGAAAAGGCCCTCC‑GCAGCCCCCCCAACCTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGT  >  AE015451_phaM1/4741092‑4741343
                                                                                                                                         |                                                                                                                        
tcttctTTCAGGTTACGGGTCAGGCTGCCATCCGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTgag                                                                                                                           <  1:14350071/137‑1 (MQ=255)
  ttcttTCAGGTTACGGGTCAGGCTGCCATCCGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTgagca                                                                                                                         >  2:10405460/1‑135 (MQ=255)
  ttcttTCAGGTTACGGGTCAGGCTGCCATCCGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTgagca                                                                                                                         >  1:12275390/1‑135 (MQ=255)
   tcttTCAGGTTACGGGTCAGGCTGCCATCCGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTgagcag                                                                                                                        >  2:10090243/1‑134 (MQ=255)
   tcttTCAGGTTACGGGTCAGGCTGCCATCCGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTgagcag                                                                                                                        >  2:11934495/1‑134 (MQ=255)
   tcttTCAGGTTACGGGTCAGGCTGCCATCCGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTgagcag                                                                                                                        >  1:7373266/1‑134 (MQ=255)
   tcttTCAGGTTACGGGTCAGGCTGCCATCCGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTgagcag                                                                                                                        >  2:16739376/1‑134 (MQ=255)
   tcttTCAGGTTACGGGTCAGGCTGCCATCCGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTgagcag                                                                                                                        >  1:3073530/1‑134 (MQ=255)
   tcttTCAGGTTACGGGTCAGGCTGCCATCCGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTgagcag                                                                                                                        >  2:6534555/1‑134 (MQ=255)
   tcttTCAGGTTACGGGTCAGGCTGCCATCCGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTgagcag                                                                                                                        >  2:6346679/1‑134 (MQ=255)
   tcttTCAGGTTACGGGTCAGGCTGCCATCCGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTgagcag                                                                                                                        >  2:6346678/1‑134 (MQ=255)
   tcttTCAGGTTACGGGTCAGGCTGCCATCCGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTgagcag                                                                                                                        >  1:2416578/1‑134 (MQ=255)
   tcttTCAGGTTAAGGGTCCGGCTGCCATCCGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCCGCTCTGGTCACATTCGTGCTTgagcag                                                                                                                        >  2:11934450/1‑134 (MQ=255)
    cttTCAGGTTACGGGTCAGGCTGCCATCCGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTgagcagc                                                                                                                       >  1:6276294/1‑133 (MQ=255)
                                                                                               gAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACgg                       >  2:9168592/1‑137 (MQ=255)
                                                                                                aaCTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGt                      >  2:7960132/1‑137 (MQ=255)
                                                                                                aaCTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGt                      >  1:5805458/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                      >  2:10478291/1‑137 (MQ=255)
                                                                                                aaCTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGt                      >  1:18602202/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)
                                                                                                 aCTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTc                     <  1:16988682/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                     <  2:5517157/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                    <  1:1145707/137‑1 (MQ=255)
                                                                                                  cTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTca                    <  2:18878761/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                 >  2:14428294/1‑137 (MQ=255)
                                                                                                     ggTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCacac                 >  2:5870818/1‑137 (MQ=255)
                                                                                                     ggTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCacac                 >  2:5195707/1‑137 (MQ=255)
                                                                                                     ggTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCacac                 >  2:6745464/1‑137 (MQ=255)
                                                                                                     ggTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCacac                 >  1:6338841/1‑137 (MQ=255)
                                                                                                     ggTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCacac                 >  1:15859094/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              >  2:10741043/1‑137 (MQ=255)
                                                                                                        ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgc              <  1:8414392/137‑1 (MQ=255)
                                                                                                        ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgc              >  2:7993254/1‑137 (MQ=255)
                                                                                                        ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgc              >  2:7913346/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:2642543/137‑1 (MQ=255)
                                                                                                        ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgc              <  2:4306392/137‑1 (MQ=255)
                                                                                                        ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgc              <  2:4978183/137‑1 (MQ=255)
                                                                                                        ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgc              >  1:2225172/1‑137 (MQ=255)
                                                                                                        ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgc              >  1:1196630/1‑137 (MQ=255)
                                                                                                        ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgc              <  1:9131656/137‑1 (MQ=255)
                                                                                                        ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgc              >  1:3497055/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:4658548/1‑137 (MQ=255)
                                                                                                        ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgc              <  1:13271468/137‑1 (MQ=255)
                                                                                                        ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgc              >  1:13970889/1‑137 (MQ=255)
                                                                                                        ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgc              >  1:17452356/1‑137 (MQ=255)
                                                                                                        ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgc              >  1:7528140/1‑137 (MQ=255)
                                                                                                        ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgc              >  1:7340470/1‑137 (MQ=255)
                                                                                                        ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgc              <  1:14678053/137‑1 (MQ=255)
                                                                                                        ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgc              >  1:17341335/1‑137 (MQ=255)
                                                                                                        ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgc              <  1:17050563/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              >  2:17438954/1‑137 (MQ=255)
                                                                                                        ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGC‑TTC‑TGCGAC‑TTGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgc              >  1:9511285/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         <  2:6321729/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:18666551/1‑137 (MQ=255)
                                                                                                             tCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAg         >  1:7646517/1‑137 (MQ=255)
                                                                                                             tCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAg         >  1:7646467/1‑137 (MQ=255)
                                                                                                             tCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAg         <  2:1091589/137‑1 (MQ=255)
                                                                                                             tCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAg         <  2:1266666/137‑1 (MQ=255)
                                                                                                             tCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAg         >  1:5516723/1‑137 (MQ=255)
                                                                                                             tCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAA‑TCCAGGAACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAg         <  1:9266605/137‑1 (MQ=255)
                                                                                                             tCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAA‑TCCAGGAACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAg         <  1:14713601/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)
                                                                                                               cAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGc       <  1:13872157/137‑1 (MQ=255)
                                                                                                               cAGCTCTGGTCACCTTCGTGCTTGAGCAGCT‑TCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGc       <  1:2078294/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)
                                                                                                                                         |                                                                                                                        
TCTTCTTTCAGGTTACGGGTCAGGCTGCCATCCGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑AACC‑TCC‑GCCAAC‑TGCAATTCGTGAA‑CCAAGCAAAAGGCCCCGAAAAGGCCCTCC‑GCAGCCCCCCCAACCTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGT  >  AE015451_phaM1/4741092‑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: ‑