Predicted mutation
evidence seq id position mutation annotation gene description
RA AE015451_phaM1 1,499,498 1 bp→CG intergenic (+75/‑89) trpS CDS (Tryptophan‑‑tRNA ligase) → / → zapE CDS (putative enzyme with nucleoside triphosphate hydrolase domain) –/–

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451_phaM11,499,4971.C79.0% 620.1 / 264.9 308intergenic (+74/‑90)trpS CDS (Tryptophan‑‑tRNA ligase)/zapE CDS (putative enzyme with nucleoside triphosphate hydrolase domain)–/–
Reads supporting (aligned to +/- strand):  ref base . (22/42);  new base C (188/53);  total (211/97)
Fisher's exact test for biased strand distribution p-value = 1.59e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01
*AE015451_phaM11,499,4980TG84.4% 573.3 / 36.4 308intergenic (+75/‑89)trpS CDS (Tryptophan‑‑tRNA ligase)/zapE CDS (putative enzyme with nucleoside triphosphate hydrolase domain)–/–
Reads supporting (aligned to +/- strand):  ref base T (29/12);  new base G (169/53);  total (211/97)
Fisher's exact test for biased strand distribution p-value = 4.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 329 total aligned reads displayed.

TGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTG‑CCCGCGAAACAGGCGCCGCGGTGCA‑TGGCCACGGCTTTGCCGGT‑G‑TT‑T‑CGGG‑CA‑C‑G‑CC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTT  >  AE015451_phaM1/1499363‑1499633
                                                                                                                                            |||                                                                                                                                           
tGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGGGTG‑ACCACGAAAAAGGCGCCACAGTGTA‑TTGCCCCGGCTTTTTCGGT‑t‑tt‑tt                                                                                                                                              >  1:8514989/1‑136 (MQ=255)
       cGTGATGCGGCCATTCAGAGCCTGAACGAAGCATTGGCACCACAGCAGGATTGAGGCTGCCAAACGGTGGGAGAGGGAGTG‑CCCGCGAAACAGGCGCCGCGGTGCA‑TGGCACCGGCTTTGCCGGT‑G‑TTCG‑CGGG‑CA‑c                                                                                                                                     >  2:13325532/1‑137 (MQ=25)
          gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGCGGCTGCCAAACGGTGGGAGCGGGCGTG‑TCCTCGAAACAGGCGCCGCGGTGCA‑CGGCACAGGCTTTGCCAGT‑G‑TTCG‑CGGG‑CA‑a‑g‑cc                                                                                                                                >  2:12695810/1‑137 (MQ=255)
          gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTG‑CCCGCAAAACAAGCGCCACGGTGCA‑TGGCACCGGTTTTTCCGGT‑T‑TTCC‑CGGG‑CA‑c‑a‑cc                                                                                                                                >  1:8776118/1‑137 (MQ=25)
            tGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGGGTG‑CCCGCGAAACAAGAGCCGCAGTTCC‑TGGCCCCAGCTTTTTCGGT‑G‑TTCG‑CGGG‑CA‑c‑g‑ccca                                                                                                                              >  1:14031892/1‑135 (MQ=255)
            tGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTG‑CCCGCGAAACAGGCGCCGCGGTGCA‑TGGCACCGGCTTTGCCGGT‑G‑TTCG‑CGGG‑CA‑c‑g‑cccg                                                                                                                              >  2:2179692/1‑135 (MQ=25)
                 ctaTTCAGATCCTGAACGAAGCCTGGGCACCACAGCAGGATTGAGGCTGCCAAACGGTGGGAGAGGGTGGG‑CCAGCAAAACGGACGCCGCGGTGCA‑GGACACCGGCTTTGCTGGT‑G‑TTCG‑CGGA‑CA‑C‑GCCC‑GCTccc                                                                                                                         <  2:17762395/135‑1 (MQ=255)
                  cATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGGGTG‑TCCGCGAAACAAGCGCCACGGTGCA‑CGGCACCAAAGTTGCCTGT‑G‑TTCT‑CGGG‑TA‑C‑GCCC‑GCTcccc                                                                                                                        >  1:3614345/1‑136 (MQ=255)
                   aTTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTG‑CCCGCGAAAAAAGCGCCACGGTGTC‑TGGCACCGGCTTTGCCGGT‑G‑TT‑TGCGGG‑CA‑C‑GCCC‑GCTCCcac                                                                                                                       >  2:11918265/1‑137 (MQ=21)
                      cAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTG‑CCCGCGAAACAGGCGCCGCGGTGCC‑TGGCACCGGATTTGCCGGT‑G‑TTCG‑AGGG‑CA‑C‑GCCC‑GCTCCCACAGa                                                                                                                    >  1:13817184/1‑137 (MQ=21)
                       agagCCGGAACGAAGCTTTGGCCCCACAGCAGGACTGAGGTTGCCAAACGGTGGGAGCGGTCGTG‑CCCGCGAAACAGGCGCCGCGGGGCG‑TGAAACCGGTTTTGCTGGT‑G‑TTCG‑CGGG‑CA‑C‑GCCC‑GCTCCCACAGAc                                                                                                                   <  2:14932355/137‑1 (MQ=19)
                        gagCCTGAACGAATCTTTGGCACCACAGCAGGACTGGGGCTGCAAAATGTTGGGGGAGGGTGTG‑CCCGCGAAACGGGCGCCGCGGTGCG‑TGGAACTGGTTTTTCCGGT‑G‑TT‑CGTGGA‑CA‑C‑GCCC‑GCTCCCACAGACt                                                                                                                  <  2:4011239/137‑1 (MQ=255)
                            cTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTG‑CCCACGAAACAGGCGCCGCGGTGCC‑TGGCCCCAGATTTGACAGT‑G‑TTCG‑CGGG‑CA‑C‑GCCC‑ACCCCCAAAGACTGAtt                                                                                                              >  2:331963/1‑137 (MQ=21)
                             tGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTG‑ACCGCGAAACAGGCGCCACGGTGCA‑CGGCACCGGATTTGGCGGT‑G‑TTCG‑AGGG‑CA‑T‑GCCC‑GCTCCCACAGAATGATTg                                                                                                             >  2:7349217/1‑137 (MQ=21)
                              gAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTG‑CCCGCGAAACAGGCGCCGCGGTGCA‑TGGCAACGGATTTGCCGGT‑G‑TT‑CGCTGG‑CA‑C‑G‑CCAGCTCCCAAAGAATGATTGc                                                                                                            >  2:236899/1‑137 (MQ=21)
                              gAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTG‑CCCAGGAAACAGGCGCCGCGGTGCA‑TGGCAACGGTCTTGCCGGT‑G‑TTCG‑CGGG‑CA‑C‑GCCC‑GATCCCAAAGACTGATTga                                                                                                            >  2:16402286/1‑136 (MQ=25)
                                 cGAAGCTTTGGCACCACAGCAGGACTGAGGCTACCAAATGGTGGGGGCGGGCGTG‑CCCGCGAAACAGGCGCTGCGGTGCA‑TGGCACCGGTTTTGCCGGT‑G‑TT‑CGCGGG‑CA‑C‑GCCC‑GCTCCCACAGACTGATTGCCAt                                                                                                         <  2:16539234/137‑1 (MQ=18)
                                  gAAGCTTTGGCACAACGGCATGACGGAGGCTGCCAAACAGGGGAAGAGGGCGTT‑CCCGCGAACCAGGCGCCCCGGTGCC‑TGTCACCGGCTTTGCCGGT‑G‑TT‑TGTGGG‑CA‑T‑GCCC‑GCTCCCACAGACTGATTGCCAtt                                                                                                        <  2:18256861/137‑1 (MQ=12)
                                  gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAATGGTGGGAGCGGGCGTG‑TCCGCGAAAAAGGCGCCAAAGTGCA‑TGGCCACGGATTTGACAGT‑G‑TTCA‑CGGG‑GA‑C‑GCCC‑ACTTCCACAGGATGATTGCCAtt                                                                                                        >  1:18713819/1‑137 (MQ=255)
                                  gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGGGTG‑CCCACGAAACAAGCGTCGCGGTGCC‑TTGCACCGGCCTTTTCGGT‑G‑TTCG‑CGGG‑CA‑C‑GACC‑GCTCCCAAAAACAGATTTCCAtt                                                                                                        >  1:4041925/1‑137 (MQ=255)
                                  gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGGGTG‑ACCACGAAAAAAGAGCCACGGTGCC‑TGGCACCAGCTTTTCCGGT‑GTTT‑T‑CTGG‑CACC‑C‑CC‑CCCCCCCCAAACCGATTGCCacc                                                                                                        >  1:14604031/1‑135 (MQ=255)
                                  gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTG‑TCCACAAAACAGGCGCCACAGTGCC‑CTGCAACAGCCTTTCCGGT‑G‑TTCT‑CGGGACC‑C‑C‑CC‑GCCCCCCCAGAACGATTGcccct                                                                                                        >  1:3128091/1‑134 (MQ=255)
                                  gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTG‑GCCGCGAAACAAGCACCGCGGTGCA‑TGGCAACAGATTTGCCGGT‑G‑TTCT‑CGGG‑CA‑C‑GACC‑ACTCCCACAGAATGACTGCCagt                                                                                                        >  1:17589470/1‑135 (MQ=255)
                                  gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTG‑CCCGCGAAACAGGCGCCGCGGTGTA‑TGGCACCGGCTTTGCCGGT‑G‑TT‑TGCGGG‑CA‑C‑GCCC‑GCTCCCACAGACTGATTGCCAtt                                                                                                        >  2:15096648/1‑137 (MQ=21)
                                  gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTG‑CCCGCGAAAACAGCGCAGAGGTGTA‑TGGCCCCGGCCTTGCCAGTGG‑TC‑A‑CGGG‑CA‑CAC‑CC‑ACTCCCCAAGACTGATTggcctt                                                                                                        >  2:12292264/1‑132 (MQ=16)
                                      ccTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTG‑CCCGCGAAACAGGCGCCGCGGTGCA‑TGGCACCGGCTTTGCCGGT‑G‑TT‑CGCGGG‑CA‑C‑GCCC‑GCTCCCACAGAATGATTGCCATTAAGc                                                                                                    >  2:3384883/1‑137 (MQ=18)
                                        ttgtccccacATCAGGATTGAGTTTGCCGAATGGGGGGATTGGGTGTC‑CCTGCGAAACAGGCGCTGCGGTACA‑TGGCACCGGCTTTTCTGGT‑T‑TTCG‑CGGG‑CA‑C‑TCCC‑GCTCCCACAGACTGATTGCCATTAAGCgg                                                                                                  <  1:17111333/131‑1 (MQ=255)
                                        ttGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTG‑CCCGCGAAACAAGCGCCGCGGTGCA‑TGGCAACGGCTTTGCCGGT‑G‑TT‑CGAGGG‑CA‑C‑GCCC‑GATTCCAAAGAATGATTGCAATTAAGCgg                                                                                                  >  2:17827242/1‑137 (MQ=18)
                                           gCACCACAGCAGGACTGAGGCTGCCAAATGGTGGGAGTGGGCGTG‑CCTGCGAAACGGGCGCCGTGGTGCA‑TGGCACCGGCTTTGCTGGT‑G‑TTCG‑TGGG‑CA‑C‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                                               <  1:16925188/137‑1 (MQ=21)
                                           gCACCAAAGCAGGATGGAGGTTGCAAAATGGTGGGAGCGGGTGTG‑CCCGAGAAACAGGCGCCGCGGTGCA‑TGGCACCGTTTTTGCTGGT‑G‑TTCG‑CGGG‑CA‑C‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                                               <  2:16235435/137‑1 (MQ=12)
                                                acaGCAGGACTGAGGTTGCAAAATGGTGGGAGTGGTCGTG‑CCTGCGAAAAAGGCGCCGGGGGGCG‑TGGCACCGGCTTTGCTGGT‑G‑TT‑CGCGGG‑CA‑C‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCg                                                                                          <  2:2949226/137‑1 (MQ=16)
                                                   gcagGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTG‑CCCGCGAAACAGGCGCCGCGGTGCA‑TGGCACCGGCTTTGCCGGT‑G‑TTCG‑CGGG‑CA‑C‑GCCC‑TCTCCCACAGAATGATTGCCATTCAGCGGGCCTATCGCtc                                                                                       >  2:8073704/1‑136 (MQ=18)
                                                               cTGCCAAACGGTGGGAGCGGGCGTT‑CCCGCGAAACAAGCGCCGCGGTGCA‑CGGCACCGGCTTTGCCGGT‑T‑TTCG‑CGGG‑CC‑C‑GCCC‑GCTCCCCCAAAATGATTGCCATTCAGCAGGCCTATCGCTACAATGACGGccc                                                                           >  1:9345517/1‑137 (MQ=12)
                                                               cTGCCAAACGGTGGGAGCGGGCGTGCCCCGCGAAACAGGCGCCGCGGTGCA‑TGGCACCGGCTTTGCCGGT‑G‑TTCG‑CGGG‑CA‑C‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGcc                                                                            >  2:11570827/1‑137 (MQ=21)
                                                               cTGCCAAACGGTGGGAGCGGGCGTG‑CCCGCGAAACAAGCGCCGCAGTGCA‑TGGCACCAGCTTTGCCGGT‑G‑TTCG‑CGGG‑CA‑C‑GCCC‑GCTCCCACAGACTGATTGCACTTCCGCGGGCCTACCGCTACAGTGACGGcac                                                                           >  2:18149894/1‑135 (MQ=14)
                                                               cTGCCAAACGGTGGGAGCGGGCGTG‑CCCACGAAACAGGCGCCACGGTGCA‑TGGCACCGGCTTTTCCGGT‑G‑TTCG‑CGGG‑AA‑CAC‑CC‑GCTCCCACAGAATGGTTGCCATTAAACGGGAATACCGCTACAGTGACGGccc                                                                           >  2:14443830/1‑137 (MQ=16)
                                                               cTGCCAAACGGTGGGAGCGGGAGTG‑CCCGCGAAACAGGAGCCACGGTGCA‑TGGCCACAGATTTGCCGGT‑G‑TTCG‑CGGG‑CA‑CAC‑CC‑GCTCCCAAAGAATGATTGACACTAAACGGGCCCATCGCTACAGGGACGGccc                                                                           >  2:14343223/1‑137 (MQ=16)
                                                               cTGCCAAACGGTGGGAGCGGGAGTG‑CCCACAAAACAAGCGCCACGGTGCA‑CGGCCCCAGCTTTGCCGGT‑G‑TTCG‑CGGG‑CA‑C‑GCCC‑ACTCCCACAGACTGATTGCCCCTCAGCGGGCCTCCCACTCCCGTGACGGccc                                                                           >  2:731398/1‑137 (MQ=12)
                                                                    aaaCGGTGGGAGCGGGCGTG‑CCCGCGAAACAGGCGCCGCGGTGCA‑TGGCACCGGCTTTGCCGGT‑G‑TTCG‑CGGG‑CA‑C‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGTGGGCACATCGCTACAGTGACGGCCCGtttt                                                                      >  2:17239745/1‑137 (MQ=18)
                                                                    aaaCGGTGGGAGCGGGCGTG‑CCCGCGAAACAGGCGCCACGGTGCA‑TGGCACCGGCTTTGCCGGT‑G‑TTCG‑CGGG‑CA‑C‑GCCC‑GCTCCCAAAGAATGATTGCCATTCAGCGGGCCTATCGCTAAAGTGACGGCCCGtttt                                                                      >  2:15882352/1‑137 (MQ=18)
                                                                    aaaCGGTGGGAGCGGGCGTG‑CCCGCGAAACAAGTGTCGCGGTGCC‑TGGCACCGGCTTTTCCTGT‑G‑TTCG‑CGGG‑CC‑C‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGACTCTCGCTACCGTGACGGCCCGtttt                                                                      >  2:7967727/1‑137 (MQ=17)
                                                                     aaCGGTGGGAGCGGGCGTG‑CCCGCGAAACAAGCGCCACGGTGTC‑TGGCACCGGCTTTTCCGGT‑G‑TTCT‑CGGG‑CA‑CAC‑CC‑GCTCCCCCCAACCGATTGCCCTTCAGCGGGTCTCTCGCTCCAATGACGGCCCCttttt                                                                     >  2:16812437/1‑137 (MQ=9)
                                                                        ggTGGGAACGGGCGTT‑CCTGCGAAACGGGCGCCGTGTTGCG‑TGGCACCGGCTTTGCCGGT‑G‑TT‑CGCGGG‑CA‑C‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTttgt                                                                  <  2:16962199/137‑1 (MQ=21)
                                                                               gCGGGCGTG‑CCCGCGAAACAAGCACCGCGGTGCACTGCCCCCGCCTTTCCCGTT‑T‑TC‑T‑CGGG‑CA‑CAC‑CC‑CCTCCCACAGAACGACTGTCCCTCAGCGGGCCTATCGTTCCAGTGACGGCCCTTTTTTTTTTTCCTcc                                                           >  2:533098/1‑136 (MQ=12)
                                                                                       tccccGCGAAAAAGGCGCCGCGGTGCA‑TGGCACCTGCTTTGCCGGT‑G‑TTCG‑TGGG‑CA‑C‑CCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGa                                                    <  1:10522407/135‑1 (MQ=255)
                                                                                           cgtgAAACAGGCGCCGCGGGGCA‑GGGCACCTGCTTTTCTGGT‑T‑TTCG‑CGGG‑CA‑C‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                <  2:18802809/134‑1 (MQ=255)
                                                                                           cgcgAAACAGGCGCCGTGGTGCG‑TGGCACTGGCTTTGCCGGT‑G‑TTCG‑GGGG‑CA‑C‑ACCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                <  1:13295036/137‑1 (MQ=255)
                                                                                           cgcgAAACAGGCGCCGCGGTGCG‑TGGCACCGGCTTTGCTGGT‑G‑TTCG‑CGGG‑CA‑C‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                <  2:13687200/137‑1 (MQ=37)
                                                                                           cgcgAAACAGGCGCCGCGGTGCA‑TGGCACCGGTTTTTCCGGT‑G‑TT‑CGCGGG‑CA‑C‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                <  2:13788660/137‑1 (MQ=25)
                                                                                           cgcgAAACAGGCGCCGCGGTACA‑TGACACTGGCTTTACTGGT‑G‑TTCG‑CGGG‑CA‑C‑CCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                <  1:18942267/137‑1 (MQ=255)
                                                                                           cgcaAAACAGGCGCCGCGGTGCA‑GGACACCGGTTTTGCGGGT‑G‑TT‑TGCGGG‑CA‑T‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                <  1:12463548/137‑1 (MQ=37)
                                                                                           agaaaaaCGGGCGCCGCGGGAAA‑TGGCACTGGTTTTGCTGGT‑G‑TTCG‑TGGG‑CA‑C‑ACCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                <  1:6544651/133‑1 (MQ=255)
                                                                                            gcgAAACATGCGCTACGGTGCA‑TGTCACAGGCTTTGCCGGT‑G‑TTCG‑CGGG‑CA‑C‑ACCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCAt                                               <  2:3911523/137‑1 (MQ=25)
                                                                                            gcgAAACAGGCTCCGCGGTGCA‑TGGCACCTGCTTTGCCGGT‑G‑TT‑CGCGGG‑CA‑C‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCAt                                               <  2:7213849/137‑1 (MQ=25)
                                                                                            acgAAACAGGCGCCGCGGTGCA‑TGGCACCGGCTTTGCCGGT‑G‑TT‑CGCGGG‑CA‑C‑GCCC‑GCTCCCACTGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCAt                                               <  2:973168/136‑1 (MQ=25)
                                                                                               aaaCAGGCGCCGCGGTGCA‑TGGCACCGGCTTTGCCGGT‑G‑TTCG‑CGGG‑CA‑C‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGAc                                            >  1:16188829/1‑137 (MQ=37)
                                                                                                 aCAGGCGCCGCGGTGCA‑TGGCACCGGCTTTGCCGGT‑G‑TTCG‑CGGG‑CA‑C‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTc                                          >  2:12520566/1‑137 (MQ=37)
                                                                                                 aCAGGCGCCGCGGTGCA‑TGGCACCGGCTTTGCCGGT‑G‑TTCG‑CGGG‑CA‑C‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTc                                          >  2:8713361/1‑137 (MQ=37)
                                                                                                 aCAGGCGCCGCGGTGCA‑TGGCACCGGCTTTGCCGGT‑G‑TTCG‑CGGG‑CA‑C‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTc                                          >  2:14101033/1‑137 (MQ=37)
                                                                                                   gggCGCCACGGTGCA‑GGGAACTGGTTTTTCCGGT‑G‑TT‑CGCGGG‑CA‑C‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTccc                                        <  2:643311/136‑1 (MQ=25)
                                                                                                      cgctgtggTGCG‑TGGCACCGGTTTTGCTGGT‑G‑TTCT‑CGGA‑CA‑C‑CCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa                                     <  1:13180200/131‑1 (MQ=255)
                                                                                                      cgctgtggTGCA‑TGGCACTGTTTTTGCTGGT‑G‑TTCG‑TGGG‑CA‑C‑TCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa                                     <  1:712035/131‑1 (MQ=255)
                                                                                                      cgccgcGGTGCG‑GGGCACCGGCTTTTCTGGT‑G‑TT‑CGTGGA‑CA‑C‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa                                     <  2:10943708/137‑1 (MQ=255)
                                                                                                      cgccgcGGTGCA‑TGGCACCGGCTTTGCCGGT‑G‑TT‑CGTGGG‑CA‑C‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa                                     <  2:2655807/137‑1 (MQ=37)
                                                                                                      cgccgcGGGGCA‑GGGCACCGGCTTTGCCGGT‑G‑TT‑CGCGGG‑CA‑C‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa                                     <  2:10290123/137‑1 (MQ=37)
                                                                                                      caccgcGGGGCA‑TGGCACCGGCTTTGCTGGT‑G‑TT‑CGCGGG‑CA‑CTC‑CC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa                                     <  2:5912540/135‑1 (MQ=255)
                                                                                                           cgGTGCA‑TGGCACCGGCTTTGCCGGT‑G‑TTCG‑AGGG‑CA‑C‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACg                                <  2:16595578/137‑1 (MQ=37)
                                                                                                                cA‑TGGCACCGGTTTTGTCGGT‑G‑TTCG‑CGGG‑CA‑C‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATc                           <  2:4853143/137‑1 (MQ=37)
                                                                                                                   tggcaccGGCTTTGCCGGT‑G‑TTCG‑CGGG‑CA‑C‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCaa                         >  1:10094774/7‑137 (MQ=37)
                                                                                                                       accGGCTTTGCCGGT‑G‑TTCG‑CGGG‑CA‑C‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCgg                                                                                                  <  1:18054088/58‑1 (MQ=11)
                                                                                                                       accGGCTTTGCCGGT‑G‑TTCG‑CGGG‑CA‑C‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAg                     >  1:5428500/3‑137 (MQ=38)
                                                                                                                       accGGCTTTGCCGGT‑G‑TTCG‑CGGG‑CA‑C‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAg                     >  1:8411261/3‑137 (MQ=38)
                                                                                                                            cTTTGCCGGT‑G‑TTCG‑CGGG‑CA‑C‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCgg                                                                                                  <  1:11756771/55‑1 (MQ=11)
                                                                                                                            cTTTGCCGGT‑G‑TTCG‑CGGG‑CA‑C‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAGTCATGACTCCCTTAGACCGCTATCAAGCAGATCTg                >  1:10173760/1‑137 (MQ=38)
                                                                                                                            cTTTGCCGGT‑G‑TTCG‑CGGG‑CA‑C‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAGTCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                >  1:15361591/1‑137 (MQ=38)
                                                                                                                            cTTTGCCGGT‑G‑TTCG‑CGGG‑CA‑C‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                >  1:2188346/1‑137 (MQ=39)
                                                                                                                            cTTTGCCGGT‑G‑TTCG‑CGGG‑CA‑C‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                >  1:10986664/1‑137 (MQ=39)
                                                                                                                            cTTTGCCGGT‑G‑TTCG‑CGGG‑CA‑C‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                >  2:892859/1‑137 (MQ=39)
                                                                                                                            cTTTGCCGGT‑G‑TTCG‑CGGG‑CA‑C‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                >  1:12975891/1‑137 (MQ=39)
                                                                                                                            cTTTGCCGGT‑G‑TTCG‑CGGG‑CA‑C‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                >  1:14954298/1‑137 (MQ=39)
                                                                                                                            cTTTGCCGGT‑G‑TTCG‑CGGG‑CA‑C‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                >  1:15687564/1‑137 (MQ=39)
                                                                                                                            cTTTGCCGGT‑G‑TTCG‑CGGG‑CA‑C‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                >  2:5771074/1‑137 (MQ=39)
                                                                                                                            cTTTGCCGGT‑G‑TTCG‑CGGG‑CA‑C‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                >  2:4311125/1‑137 (MQ=39)
                                                                                                                            cTTTGCCGGT‑G‑TTCG‑CGGG‑CA‑C‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                >  1:7423523/1‑137 (MQ=39)
                                                                                                                            cTTTGCCGGT‑G‑TTCG‑CGGG‑CA‑C‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                >  2:3524887/1‑137 (MQ=39)
                                                                                                                            cTTTGCCGGT‑G‑TTCG‑CGGG‑CA‑C‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                >  2:12946056/1‑137 (MQ=39)
                                                                                                                            cTTTGCCGGT‑G‑TTCG‑CGGG‑CA‑C‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                >  2:2549993/1‑137 (MQ=39)
                                                                                                                              ttGCCGGT‑G‑TTCG‑CGGG‑CA‑C‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGaa              >  2:203104/1‑137 (MQ=39)
                                                                                                                              ttGCCGGT‑G‑TTCG‑CGGG‑CA‑C‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGaa              >  1:9697136/1‑137 (MQ=39)
                                                                                                                              ttGCCGGT‑G‑TTCG‑CGGG‑CA‑C‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGaa              >  2:11065676/1‑137 (MQ=39)
                                                                                                                              ttGCCGGT‑G‑TTCG‑CGGG‑CA‑C‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGaa              >  2:77210/1‑137 (MQ=39)
                                                                                                                              ttGCCGGT‑G‑TTCG‑CGGG‑CA‑C‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGaa              >  2:13896738/1‑137 (MQ=39)
                                                                                                                              ttGCCGGT‑G‑TTCG‑CGGG‑CA‑C‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGaa              >  2:15418953/1‑137 (MQ=39)
                                                                                                                                        tgttcgcGGG‑CA‑C‑TCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTccc       <  1:17374385/131‑1 (MQ=255)
                                                                                                                                        tgttcgcGGG‑CA‑C‑TCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTccc       <  1:16812437/131‑1 (MQ=255)
                                                                                                                                         gttcgcGGG‑CA‑C‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCg      >  1:8097101/6‑137 (MQ=255)
                                                                                                                                         gttcgcGGG‑CA‑C‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCg      >  1:11073181/6‑137 (MQ=255)
                                                                                                                                          ttcgcGGG‑CA‑C‑CCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGa     <  1:11134095/133‑1 (MQ=255)
                                                                                                                                            cgcGGG‑CA‑C‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGAct   >  1:3747760/3‑137 (MQ=255)
                                                                                                                                             gcGGG‑CA‑C‑GCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGActt  >  2:13133562/2‑137 (MQ=255)
                                                                                                                                            |||                                                                                                                                           
TGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTG‑CCCGCGAAACAGGCGCCGCGGTGCA‑TGGCCACGGCTTTGCCGGT‑G‑TT‑T‑CGGG‑CA‑C‑G‑CC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTT  >  AE015451_phaM1/1499363‑1499633

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