Predicted mutation
evidence seq id position mutation annotation gene description
RA AE015451_phaM1 4,586,057 +C intergenic (+140/+75) PP_4061 CDS (conserved protein of unknown function) → / ← PP_4063 CDS (putative Long‑chain‑fatty‑acid‑CoA ligase) –/–

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451_phaM14,586,0561.C98.2% 418.1 / 0.6 110intergenic (+139/+76)PP_4061 CDS (conserved protein of unknown function)/PP_4063 CDS (putative Long‑chain‑fatty‑acid‑CoA ligase)–/–
Reads supporting (aligned to +/- strand):  ref base . (2/0);  new base C (23/85);  total (25/85)
Fisher's exact test for biased strand distribution p-value = 5.00e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 126 total aligned reads displayed.

TAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTCTGC‑GG‑TAAACCCGCTCCTACCAGGG‑A‑CGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACT  >  AE015451_phaM1/4585921‑4586187
                                                                                                                                           |                                                                                                                                   
  gACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGG‑A‑cc                                                                                                                                   <  1:10443037/137‑2 (MQ=39)
  gACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGG‑A‑cc                                                                                                                                   <  2:2444805/137‑2 (MQ=39)
  gACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGG‑A‑cc                                                                                                                                   <  2:1222766/137‑2 (MQ=39)
  gACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGG‑A‑cc                                                                                                                                   <  1:6993812/137‑2 (MQ=39)
  gACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGG‑A‑cc                                                                                                                                   <  1:5408701/137‑2 (MQ=39)
  gACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGG‑A‑cc                                                                                                                                   <  1:5408656/137‑2 (MQ=39)
  gACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATGGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGG‑A‑cc                                                                                                                                   <  1:11622341/137‑2 (MQ=38)
   aCCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGG‑A‑ccg                                                                                                                                  <  1:8940968/137‑3 (MQ=38)
   aCCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGG‑A‑ccg                                                                                                                                  <  1:4298101/137‑3 (MQ=38)
   aCCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGG‑A‑ccg                                                                                                                                  <  1:17612920/137‑3 (MQ=38)
   aCCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGG‑A‑ccg                                                                                                                                  <  1:13697312/137‑3 (MQ=38)
     cACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGTCTCTTCGCGGG‑TAAAACCGCTCCTCCCAGGG‑A‑ccgtg                                                                                                                                >  2:9057176/1‑133 (MQ=255)
     cACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGG‑A‑ccgtg                                                                                                                                <  1:7993382/137‑5 (MQ=38)
         gCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGG‑ACCGTGCCgg                                                                                                                             <  2:16691234/137‑1 (MQ=38)
               gCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCATGG‑ACCGTGCCGGCCTTGa                                                                                                                       <  2:4159186/137‑1 (MQ=38)
               gCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCATGG‑ACCGTGCCGGCCTTGa                                                                                                                       <  1:9129932/137‑1 (MQ=38)
               gCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGG‑ACCGTGCCGGCCTTGa                                                                                                                       <  1:1183988/137‑1 (MQ=38)
               gCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGG‑ACCGTGCCGGCCTTGa                                                                                                                       <  1:104403/137‑1 (MQ=38)
               gCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGG‑ACCGTGCCGGCCTTGa                                                                                                                       <  2:1269767/137‑1 (MQ=38)
                cGTGAACAGCTGATTTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGAA‑CATTCCGGCCTTTaa                                                                                                                      >  2:13271552/1‑137 (MQ=37)
                    aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCATGG‑ACCGTGCCGGCCTTGAATgcg                                                                                                                  <  2:3804049/137‑1 (MQ=38)
                    aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCATGG‑ACCGTGCCGGCCTTGAATgcg                                                                                                                  <  1:10646698/137‑1 (MQ=38)
                    aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCATGG‑ACCGTGCCGGCCTTGAATgcg                                                                                                                  <  2:35579/137‑1 (MQ=38)
                    aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCATGG‑ACCGTGCCGGCCTTGAATgcg                                                                                                                  <  2:2657347/137‑1 (MQ=38)
                    aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCATGG‑ACCGTGCCGGCCTTGAATgcg                                                                                                                  <  1:16896849/137‑1 (MQ=38)
                    aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGG‑ACCGTGCCGGCCTTGAATgcg                                                                                                                  <  2:10830778/137‑1 (MQ=38)
                    aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGG‑ACCGTGCCGGCCTTGAATgcg                                                                                                                  <  2:169648/137‑1 (MQ=38)
                    aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGG‑ACCGTGCCGGCCTTGAATgcg                                                                                                                  <  1:4578707/137‑1 (MQ=38)
                    aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGG‑ACCGTGCCGGCCTTGAATgcg                                                                                                                  <  2:2365724/137‑1 (MQ=38)
                    aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGG‑ACCGTGCCGGCCTTGAATgcg                                                                                                                  <  1:6396178/137‑1 (MQ=38)
                    aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGG‑ACCGTGCCGGCCTTGAATgcg                                                                                                                  <  1:4330067/137‑1 (MQ=38)
                    aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGG‑ACCGTGCCGGCCTTGAATgcg                                                                                                                  <  1:14356160/137‑1 (MQ=38)
                    aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGG‑ACCGTGCCGGCCTTGAATgcg                                                                                                                  <  2:11352480/137‑1 (MQ=38)
                    aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGG‑ACCGTGCCGGCCTTGAATgcg                                                                                                                  <  1:11820122/137‑1 (MQ=38)
                    aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGG‑ACCGTGCCGGCCTTGAATgcg                                                                                                                  <  1:17360352/137‑1 (MQ=38)
                    aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGG‑ACCGTGCCGGCCTTGAATgcg                                                                                                                  <  1:17346293/137‑1 (MQ=38)
                    aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGG‑ACCGTGCCGGCCTTGAATgcg                                                                                                                  <  1:17346263/137‑1 (MQ=38)
                    aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGG‑ACCGTGCCGGCCTTGAATgcg                                                                                                                  <  1:16896892/137‑1 (MQ=38)
                    aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGG‑ACCGTGCCGGCCTTGAATgcg                                                                                                                  <  2:14398359/137‑1 (MQ=38)
                    aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGG‑ACCGTGCCGGCCTTGAATgcg                                                                                                                  <  2:5383967/137‑1 (MQ=38)
                    aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGG‑ACCGTGCCGGCCTTGAATgcg                                                                                                                  <  1:1569783/137‑1 (MQ=38)
                    aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGG‑ACCGTGCCGGCCTTGAATgcg                                                                                                                  <  2:5761172/137‑1 (MQ=38)
                    aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGG‑ACCGTGCCGGCCTTGAATgcg                                                                                                                  <  2:6163328/137‑1 (MQ=38)
                    aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGG‑ACCGTGCCGACCTTGAATgcg                                                                                                                  <  1:9311555/137‑1 (MQ=38)
                     aCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGG‑ACCGTGCCGGCCTTGAATgcgg                                                                                                                 <  1:11280243/137‑1 (MQ=38)
                      cAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGG‑ACCGTGCCGGCCTTGAATgcggc                                                                                                                <  1:5686893/137‑1 (MQ=38)
                      cAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGG‑ACCGTGCCGGCCTTGAATgcggc                                                                                                                <  2:4127095/137‑1 (MQ=38)
                      cAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGG‑ACCGTGCCGGCCTTGAATgcggc                                                                                                                <  2:4127138/137‑1 (MQ=38)
                       ggCTGATCTGCGCTGCTCCATTGAAAGGTTGGGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGG‑ACCGTGCCGGCCTTGAATgcggcg                                                                                                               <  2:14021107/136‑1 (MQ=38)
                       aGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCATGG‑ACCGTGCCGGCCTTGAATgcggcg                                                                                                               <  1:4801406/137‑1 (MQ=38)
                       aGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGG‑ACCGTGCCGGCCTTGAATgcggcg                                                                                                               <  1:7108199/137‑1 (MQ=38)
                       aGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGG‑ACCGTGCCGGCCTTGAATgcggcg                                                                                                               <  1:7108153/137‑1 (MQ=38)
                       aGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGG‑ACCGTGCCGGCCTTGAATgcggcg                                                                                                               <  1:5887575/137‑1 (MQ=38)
                       aGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGG‑ACCGTGCCGGCCTTGAATgcggcg                                                                                                               <  1:5629539/137‑1 (MQ=38)
                       aGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGG‑ACCGTGCCGGCCTTGAATgcggcg                                                                                                               <  1:5629538/137‑1 (MQ=38)
                       aGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGG‑ACCGTGCCGGCCTTGAATgcggcg                                                                                                               <  2:8619136/137‑1 (MQ=38)
                       aGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGG‑ACCGTGCCGGCCTTGAATgcggcg                                                                                                               <  2:4408124/137‑1 (MQ=38)
                       aGCTGATCTGCGCTGCTCCATTGAAAGGGTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGG‑ACCGTGCCGGCCTTGAATgcggcg                                                                                                               <  2:13164239/137‑1 (MQ=38)
                         cTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGG‑ACCGTGCCGGCCTTGAATGCGGCGCt                                                                                                             <  1:7604799/137‑1 (MQ=38)
                          tGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGG‑ACCGTGCCGGCCTTGAATGCGGCGCTg                                                                                                            <  2:13401140/137‑1 (MQ=38)
                          tGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGG‑ACCGTGCCGGCCTTGAATGCGGCGCTg                                                                                                            <  2:2534968/137‑1 (MQ=38)
                              cTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGcccc                                                                                                        <  1:1633147/137‑1 (MQ=38)
                              cTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGcccc                                                                                                        <  1:10447191/137‑1 (MQ=38)
                               tGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCCt                                                                                                       <  2:11032899/137‑1 (MQ=38)
                                   ctgctCCATTGAAAGTTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTa                                                                                                   <  2:13090952/137‑1 (MQ=38)
                                   ctgctCCATTGAAAGGTTGAGGGTGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTa                                                                                                   <  1:1717479/137‑1 (MQ=38)
                                   ctgctCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTa                                                                                                   <  2:13092572/137‑1 (MQ=38)
                                   ctgctCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTa                                                                                                   <  1:14450539/137‑1 (MQ=38)
                                   ctgctCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTa                                                                                                   <  2:10266975/137‑1 (MQ=38)
                                   ctgctCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTa                                                                                                   <  1:14647153/137‑1 (MQ=38)
                                   ctgctCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTa                                                                                                   <  1:8079579/137‑1 (MQ=38)
                                   ctgctCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTa                                                                                                   <  2:6103187/137‑1 (MQ=38)
                                   ctgctCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTa                                                                                                   <  1:11740570/137‑1 (MQ=38)
                                   ctgctCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGC‑GGATAAACCCGCTCCTACCAGGG‑ACCGTGCCTGCCTTTAATGCGGCGCTGCCCCTGTTa                                                                                                   <  2:4525905/137‑1 (MQ=37)
                                   ctgctCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCGTCGC‑GGATAAATCCGCTCCTACTAGGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTa                                                                                                   <  2:4525906/137‑1 (MQ=255)
                                    tgctCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAg                                                                                                  <  1:5547242/137‑1 (MQ=38)
                                    tgctCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAg                                                                                                  <  1:15597573/137‑1 (MQ=38)
                                              aaaGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGG‑ACCGTGCCGGCCTTGAATGCGGCGCt                                                                                                             >  2:5330625/1‑116 (MQ=37)
                                              aaaGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGG‑ACCGTGCCGGCCTTGAATGCGGCGCt                                                                                                             <  1:5330625/116‑1 (MQ=37)
                                                        gATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TGAAACCACTCCCCCCAGGGAA‑CGTTCCGGCCCTTAATGCGGAGCTTCCCCTGGTAGGAGCGGGTTTccccccaaaa                                                                              >  2:14078240/1‑131 (MQ=255)
                                                                                                                                 tACCAGGG‑ACCGTGCCGGCCTTGCATGCGGCGCTGCCACTAGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCACGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGGTCTCACGCATGCgg      >  2:4823434/1‑137 (MQ=38)
                                                                                                                                 tACCAGGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCgg      >  1:15845893/1‑137 (MQ=255)
                                                                                                                                 tACCAGGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCgg      >  1:13536238/1‑137 (MQ=255)
                                                                                                                                 tACCAGGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCgg      >  2:4869363/1‑137 (MQ=255)
                                                                                                                                 tACCAGGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCgg      >  2:7241560/1‑137 (MQ=255)
                                                                                                                                 tACCAGGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCgg      >  2:8990408/1‑137 (MQ=255)
                                                                                                                                 tACCAGGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCgg      >  1:1776514/1‑137 (MQ=255)
                                                                                                                                 tACCAGGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCgg      >  1:223173/1‑137 (MQ=255)
                                                                                                                                 tACCAGGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCgg      >  1:223204/1‑137 (MQ=255)
                                                                                                                                 tACCAGGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCgg      >  1:11573817/1‑137 (MQ=255)
                                                                                                                                 tACCAGGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCgg      >  1:12552963/1‑137 (MQ=255)
                                                                                                                                 tACCAGGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCgg      >  2:15734396/1‑137 (MQ=255)
                                                                                                                                 tACCAGGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCgg      >  2:14420655/1‑137 (MQ=255)
                                                                                                                                 tACCAGGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAAACCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCgg      >  2:1761542/1‑137 (MQ=255)
                                                                                                                                     aGGG‑ACCGTGCCGGCCTTGAATGCGGCGTTGCCCATGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACt  >  2:3269413/1‑137 (MQ=255)
                                                                                                                                     aGGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACt  >  2:2907564/1‑137 (MQ=255)
                                                                                                                                     aGGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACt  >  1:10944903/1‑137 (MQ=255)
                                                                                                                                     aGGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACt  >  2:820258/1‑137 (MQ=255)
                                                                                                                                     aGGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACt  >  2:8947896/1‑137 (MQ=255)
                                                                                                                                     aGGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACt  >  2:15211575/1‑137 (MQ=255)
                                                                                                                                           |                                                                                                                                   
TAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTCTGC‑GG‑TAAACCCGCTCCTACCAGGG‑A‑CGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACT  >  AE015451_phaM1/4585921‑4586187

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