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.C99.0% 396.8 / ‑3.1 103intergenic (+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 . (0/1);  new base C (15/87);  total (15/88)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 122 total aligned reads displayed.

GACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TCTGC‑GGTAAACCCGCTCCTACCAGGGA‑C‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACT  >  AE015451_phaM1/4585923‑4586187
                                                                                                                                        |                                                                                                                                    
gACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCATGGA‑cc                                                                                                                                    <  2:3049154/137‑2 (MQ=39)
gACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGA‑cc                                                                                                                                    <  1:16052116/137‑2 (MQ=39)
gACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGA‑cc                                                                                                                                    <  2:14328740/137‑2 (MQ=39)
gACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGA‑cc                                                                                                                                    <  1:5436404/137‑2 (MQ=39)
gACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGA‑cc                                                                                                                                    <  2:6426239/137‑2 (MQ=39)
gACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGA‑cc                                                                                                                                    <  2:3468347/137‑2 (MQ=39)
 aCCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGTAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCATGGA‑ccg                                                                                                                                   <  2:12258082/137‑3 (MQ=38)
 aCCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGA‑ccg                                                                                                                                   <  1:2118092/137‑3 (MQ=38)
 aCCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGA‑ccg                                                                                                                                   <  1:378889/137‑3 (MQ=38)
 aCCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGA‑ccg                                                                                                                                   <  2:2153982/137‑3 (MQ=38)
 aCCACGGCACCGGCGTGAACAGATGATCTGCACTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGA‑ccg                                                                                                                                   <  1:2366770/137‑3 (MQ=38)
  ccACGGCACCGGCGTGAACAGCGGATTTGCACGGCTCCATTGAAGGGTGGAGGATGA‑CTGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGA‑ccgtg                                                                                                                                 <  1:14997911/137‑5 (MQ=25)
   cACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGA‑ccgtg                                                                                                                                 <  1:10175231/137‑5 (MQ=38)
      ggCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCCTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCg                                                                                                                              <  2:11453755/137‑1 (MQ=38)
      ggCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGACGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCATGGACC‑GTGCCg                                                                                                                              <  2:11453721/137‑1 (MQ=38)
                  aaCAGTTGATCTGCGCGGCTCCAGTGAAAGTTTTAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATgcg                                                                                                                  <  1:5353033/137‑1 (MQ=37)
                  aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGTATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATgcg                                                                                                                  <  1:3479831/137‑1 (MQ=38)
                  aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATgcg                                                                                                                  <  1:8792608/137‑1 (MQ=38)
                  aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATgcg                                                                                                                  <  2:5468486/137‑1 (MQ=38)
                  aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATgcg                                                                                                                  <  1:2901678/137‑1 (MQ=38)
                  aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATgcg                                                                                                                  <  2:5726968/137‑1 (MQ=38)
                  aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATgcg                                                                                                                  <  1:3479867/137‑1 (MQ=38)
                  aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATgcg                                                                                                                  <  2:12384353/137‑1 (MQ=38)
                  aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATgcg                                                                                                                  <  2:14720375/137‑1 (MQ=38)
                  aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATgcg                                                                                                                  <  1:6897914/137‑1 (MQ=38)
                  aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATgcg                                                                                                                  <  1:9148184/137‑1 (MQ=38)
                  aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATgcg                                                                                                                  <  2:12030057/137‑1 (MQ=38)
                  aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATgcg                                                                                                                  <  1:10816443/137‑1 (MQ=38)
                  aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATgcg                                                                                                                  <  2:866389/137‑1 (MQ=38)
                  aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATgcg                                                                                                                  <  1:11201592/137‑1 (MQ=38)
                  aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATgcg                                                                                                                  <  2:8490875/137‑1 (MQ=38)
                  aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATgcg                                                                                                                  <  2:7484436/137‑1 (MQ=38)
                  aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATgcg                                                                                                                  <  1:9938772/137‑1 (MQ=38)
                  aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATgcg                                                                                                                  <  2:1415154/137‑1 (MQ=38)
                  aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATgcg                                                                                                                  <  1:9104124/137‑1 (MQ=38)
                  aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATgcg                                                                                                                  <  1:16581300/137‑1 (MQ=38)
                  aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTCCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATgcg                                                                                                                  <  1:13667774/137‑1 (MQ=38)
                  aaCAGCTGAGCTGCGCTGATCCATTGAAAGGTCGAGGATAACCGGGTGCTTGGCAAGGTGAAAACCTTGCATGGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATgcg                                                                                                                  <  2:5206304/137‑1 (MQ=25)
                     aGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTTGCAAGGTGAAGATCTTGCGTTGATCGGGCCGGTCGCGTCGCG‑GGTAAACCCGCTTCTACCAGGGA‑CTGTGCCGGCCTTGAATgcggcg                                                                                                               <  2:2149135/137‑1 (MQ=255)
                     aGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCATGGACC‑GTGCCGGCCTTGAATgcggcg                                                                                                               <  1:16547880/137‑1 (MQ=38)
                     aGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATgcggcg                                                                                                               <  2:13016728/137‑1 (MQ=38)
                     aGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATgcggcg                                                                                                               <  1:8257398/137‑1 (MQ=38)
                     aGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATgcggcg                                                                                                               <  1:8030727/137‑1 (MQ=38)
                     aGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATgcggcg                                                                                                               <  2:11275872/137‑1 (MQ=38)
                     aGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATgcggcg                                                                                                               <  2:16057466/137‑1 (MQ=38)
                     aGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATgcggcg                                                                                                               <  2:10951385/137‑1 (MQ=38)
                     aGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATgcggcg                                                                                                               <  1:948703/137‑1 (MQ=38)
                     aGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATgcggcg                                                                                                               <  1:16547920/137‑1 (MQ=38)
                     aGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATgcggcg                                                                                                               <  1:10387600/137‑1 (MQ=38)
                     aGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATgcggcg                                                                                                               <  2:8516368/137‑1 (MQ=38)
                     aGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATgcggcg                                                                                                               <  1:11173003/137‑1 (MQ=38)
                     aGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATgcggcg                                                                                                               <  1:13350844/137‑1 (MQ=38)
                     aGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATgcggcg                                                                                                               <  1:14277771/137‑1 (MQ=38)
                     aGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATgcggcg                                                                                                               <  1:10271670/137‑1 (MQ=38)
                     aGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATgcggcg                                                                                                               <  2:7319464/137‑1 (MQ=38)
                     aGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATgcggcg                                                                                                               <  1:1786245/137‑1 (MQ=38)
                     aGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATgcggcg                                                                                                               <  1:1786277/137‑1 (MQ=38)
                     aGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATgcggcg                                                                                                               <  1:2140111/137‑1 (MQ=38)
                     aGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATgcggcg                                                                                                               <  2:5129658/137‑1 (MQ=38)
                     aGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATgcggcg                                                                                                               <  1:7423039/137‑1 (MQ=38)
                     aGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATgcggcg                                                                                                               <  1:5417129/137‑1 (MQ=38)
                     aGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATgcggcg                                                                                                               <  1:6024659/137‑1 (MQ=38)
                     aGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATgcggcg                                                                                                               <  1:6024697/137‑1 (MQ=38)
                     aGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGAACGACCGGGTGCTTGGCAAGTTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCTGGTAAACCCGCTCCTACCATGGACC‑GTGCCGGCCTTGAATgcggcg                                                                                                               <  2:9365169/137‑1 (MQ=255)
                       cTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCt                                                                                                             <  1:8283085/137‑1 (MQ=38)
                       cTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCt                                                                                                             <  2:7482454/137‑1 (MQ=38)
                       cTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCt                                                                                                             <  2:7080326/137‑1 (MQ=38)
                       cTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCt                                                                                                             <  2:4186018/137‑1 (MQ=38)
                       cTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCt                                                                                                             <  1:98543/137‑1 (MQ=38)
                       cTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCt                                                                                                             <  1:10246407/137‑1 (MQ=38)
                        tGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTg                                                                                                            <  1:8524261/137‑1 (MQ=38)
                        tGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTg                                                                                                            <  2:3230527/137‑1 (MQ=38)
                        tGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTg                                                                                                            <  1:5870022/137‑1 (MQ=38)
                         gATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGc                                                                                                           <  1:917887/137‑1 (MQ=38)
                         gATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGc                                                                                                           <  1:11062010/137‑1 (MQ=38)
                         gATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGc                                                                                                           <  1:16004065/137‑1 (MQ=38)
                         gATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGc                                                                                                           <  1:16004095/137‑1 (MQ=38)
                         gATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGc                                                                                                           <  1:8442098/137‑1 (MQ=38)
                             tGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCt                                                                                                       <  1:819811/137‑1 (MQ=38)
                                 ctgctCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTa                                                                                                   <  1:5165294/137‑1 (MQ=38)
                                 ctgctCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTa                                                                                                   <  1:4994798/137‑1 (MQ=38)
                                 ctgctCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTa                                                                                                   <  1:4927534/137‑1 (MQ=38)
                                 ctgctCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTa                                                                                                   <  1:1715444/137‑1 (MQ=38)
                                  tgctCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATgcg                                                                                                                  <  2:4824816/121‑1 (MQ=37)
                                  tgctCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATgcg                                                                                                                  >  1:4824816/1‑121 (MQ=37)
                                  tgctCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAg                                                                                                  <  1:11132982/137‑1 (MQ=38)
                                       cATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGA‑ccg                                                                                                                                   >  1:15962527/1‑97 (MQ=37)
                                       cATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCAGGGA‑ccg                                                                                                                                   <  2:15962527/99‑3 (MQ=37)
                                                                                                                              cTACCAGGGACC‑GTGCCGGCCTTGACTGCGGCGCTGCCCCTGGTAGGAGCGGGTTTCCCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACACGCAGATGCAGCCGCAATCCCCGCCACACTGATCTCACGCATGCg       >  2:5436462/1‑137 (MQ=255)
                                                                                                                              cTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCg       >  2:5436508/1‑137 (MQ=255)
                                                                                                                              cTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCg       >  1:5926199/1‑137 (MQ=255)
                                                                                                                               tACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAANCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCgg      >  1:13650778/1‑137 (MQ=255)
                                                                                                                               tACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCgg      >  2:5306362/1‑137 (MQ=255)
                                                                                                                               tACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCgg      >  1:11354621/1‑137 (MQ=255)
                                                                                                                               tACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCgg      >  2:4377200/1‑137 (MQ=255)
                                                                                                                               tACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCgg      >  2:14547497/1‑137 (MQ=255)
                                                                                                                               tACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACACATGCgg      >  1:3166063/1‑137 (MQ=255)
                                                                                                                               tACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCAGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCgg      >  2:14547529/1‑137 (MQ=255)
                                                                                                                               tACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCATGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTAATCTCACGCATGCgg      >  2:16011147/1‑137 (MQ=255)
                                                                                                                                   aGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACt  >  1:4590263/1‑137 (MQ=255)
                                                                                                                                        |                                                                                                                                    
GACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TCTGC‑GGTAAACCCGCTCCTACCAGGGA‑C‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACT  >  AE015451_phaM1/4585923‑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: ‑