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.7% 567.3 / 218.9 266intergenic (+74/‑90)trpS CDS (Tryptophan‑‑tRNA ligase)/zapE CDS (putative enzyme with nucleoside triphosphate hydrolase domain)–/–
Reads supporting (aligned to +/- strand):  ref base . (26/28);  new base C (155/57);  total (181/85)
Fisher's exact test for biased strand distribution p-value = 9.22e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*AE015451_phaM11,499,4980TG85.1% 532.2 / 21.8 266intergenic (+75/‑89)trpS CDS (Tryptophan‑‑tRNA ligase)/zapE CDS (putative enzyme with nucleoside triphosphate hydrolase domain)–/–
Reads supporting (aligned to +/- strand):  ref base T (25/10);  new base G (145/55);  total (181/85)
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 283 total aligned reads displayed.

TGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGC‑CACGGCTTTGCCGGTG‑T‑T‑T‑C‑GGGCAC‑G‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGAC  >  AE015451_phaM1/1499363‑1499631
                                                                                                                                           |||                                                                                                                                       
  aaGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGGAGGACTGAGGCTGCCAAACGGTGGGAGCGGGGGTGTCCGCGAAAAAGGCGCCACGGTGCATGGC‑ACCAGCCTTTTCGGTG‑t‑t‑tgtg                                                                                                                                        >  1:15865049/1‑134 (MQ=255)
     cccGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCAGCGGTGTATGGC‑ACCGGGTTTGCCGGTG‑t‑t‑tgcggga                                                                                                                                     >  2:4583512/1‑131 (MQ=37)
          gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGTGTGCCCGCGAAACAGGCGCCGCGGTGCATGGC‑ACCGGCTTTTTCGGTG‑T‑TCG‑C‑GGGCAc‑g‑cc                                                                                                                             >  2:2721418/1‑137 (MQ=25)
          gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGAGTGCCCGCGAAACAAGCGCCGCGGTGCACGGC‑ACCAGATTTTCCAGTG‑T‑TCT‑C‑GGGCAc‑g‑cc                                                                                                                             >  2:14381565/1‑137 (MQ=25)
           aTGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCAAAACAGGCGCCGAGGTGCATGGC‑ACCGGCTTTGCCGGTG‑T‑TCG‑C‑GGGCAc‑g‑ccc                                                                                                                            >  2:8668969/1‑136 (MQ=25)
            tGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGGGTGTCCGCGAAACAAACGCCACGGTGCCCGGC‑ACCCGCTTTTTCTGTGTT‑T‑T‑C‑GGGCAc‑a‑cccc                                                                                                                           >  1:7210195/1‑135 (MQ=255)
               ggCCATTAAGATCCGGAACGAAGCTTTGGCACCACAACAGGACTGAGGCTTCCAAACGGTGGGAGCGGGCGTGCCCACGAAACGGGCACTGCGGTGCATGGC‑ACTGTTTTTGCTGGTG‑T‑T‑TGG‑GGGCAc‑g‑cccgctc                                                                                                                        <  2:17977249/137‑6 (MQ=12)
                 ccATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGC‑ACCGGCTTTGCCGGTG‑T‑TCG‑C‑GGGCAC‑GCCCGCTccc                                                                                                                       >  2:2146854/1‑137 (MQ=21)
                   aTTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCCTGGC‑ACCGGATTTGCCGGTG‑T‑TCG‑C‑GGGCAC‑GCCCGCTCCcac                                                                                                                     >  2:10770619/1‑137 (MQ=21)
                   aTTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGC‑ACCGGCTTTGCCGGTG‑T‑TCG‑C‑GGGCAC‑GCCCGCTCCcac                                                                                                                     >  2:8386896/1‑137 (MQ=21)
                   aTTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGAGTGTCCGCGAAACAAGAGCCACGGGGCATGGC‑ACCGGCTTTGCCGGTG‑T‑TCG‑C‑GGGTACAC‑CCGCTCCcac                                                                                                                     >  1:17694237/1‑137 (MQ=21)
                      cAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCAAGAAACAGGAGCCGCGGGGCATGGC‑ACCGGATTTGCCAGTG‑T‑TCG‑C‑GGGCAC‑GCCCGCTCCCACAGa                                                                                                                  >  1:2601955/1‑137 (MQ=21)
                              gAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGAGAAACAGGCGCCGCGGTGCATGGC‑ACCGGCTTTGTCGGTG‑T‑TCG‑C‑GGGCACAC‑CCGCTCCCACAGACTGATTGc                                                                                                          >  2:5387511/1‑137 (MQ=21)
                              gAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCACGAAACAGGCGCCGCGGTGCATGGC‑ACCGGCTTTGTCGGTG‑T‑TCG‑C‑AGGCACAC‑CCGCTCCCACAGCCTGACTGc                                                                                                          >  2:7622155/1‑137 (MQ=18)
                              gAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGACCGCGAAAAAGGCACCGCGGTGCATGGC‑ACCAGCTTTGGCGGTT‑T‑TCA‑C‑GGGCCC‑GCCCGATCCCACAGGCTGATTGc                                                                                                          >  2:11411147/1‑137 (MQ=255)
                                 cGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGGGTGCCCGCGAAACAGGCGCCGCGGTGCCTGGCACCCGACTTT‑CCGGTG‑T‑TCG‑C‑GGGCCCAC‑CCACTCCCACAGACTGATTGCCAt                                                                                                       >  2:16131614/1‑137 (MQ=255)
                                  gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGGGTGGCCGCGAAACAAGAGCCGCGGTGTATTGC‑ACCAGCTTTGCCGGTG‑T‑T‑TGC‑GGGACC‑GCCCACTCCCACAGAATGAGTGACAtt                                                                                                      >  2:14493285/1‑137 (MQ=255)
                                  gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAAGCGCCACGGTGCACTGC‑ACCAGCTTTTCCGGTG‑T‑TCG‑C‑GGGCAC‑GCCCCCTCCCACAAAATGATTGCCAtt                                                                                                      >  2:8098455/1‑137 (MQ=18)
                                  gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCACGAAACAGGCGCCGCGGTGCATGGC‑ACCGGCTTTGCCGGTG‑T‑T‑CGC‑GGGCAC‑GCCCGCTCCCACAGACTGATTGCCact                                                                                                      >  2:15842965/1‑135 (MQ=21)
                                  gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGAGTGTCCACGAAACAGGCGCCGCAGTGTATTGT‑ACCAGCTTTTTCGGTG‑T‑TCT‑C‑GGGCACAC‑CCACTCCCACCAACCGCTTGCCAtc                                                                                                      >  1:3504729/1‑136 (MQ=255)
                                  gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGAGTGCCCGCGAAACAAGAGCCACAGTGCACGGC‑ACCAGATTTGCCAGTG‑T‑T‑TGC‑GGGCAC‑GCCCACTCCCACAAAATGATTGCCact                                                                                                      >  2:6169870/1‑135 (MQ=255)
                                  gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGAGTGCCCACGAAACAAGCACCACGGTGTCTGGC‑ACCAACTTTTCCGGTGTT‑T‑T‑T‑GGGCACAC‑CCACTCCCCCAAAATGATTGCCAtt                                                                                                      >  1:5430246/1‑137 (MQ=255)
                                  gAAGCCTTGGCACCACAGCAGGACTGAGACTGCCAAACGGTGGGAGCGGGCGTGTCCGCGAAACAGGCGCCGCGGTGCCTTGC‑ACCGGCTTTGCCGGTG‑T‑TCG‑C‑GGGCAC‑GCCCGATCCCACAGAATGATTGCCAtt                                                                                                      >  2:15030849/1‑137 (MQ=21)
                                      ttttGGCACCACGGGGGGATTGGGGTTGCAAAACGGGGGGAGTGGGTGTGCCCGAGAAACGGGCCCTGTGGTGCGGGGC‑ACTGGTTTTGCGGGTG‑T‑T‑TGT‑GGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGc                                                                                                  <  1:4187958/135‑1 (MQ=255)
                                      ccTTGACCCCACAACAGGCCTGAGGCTGCCAAACGGTGGGGAGGGGCGTACCCGCAAAAAGTACCCCGCGGTACGTAAC‑ACCGGCTTGGCTGGGG‑TCA‑G‑T‑GGGCAC‑GCCCGCTCCCACAGACTGATTACCATTAAGc                                                                                                  <  2:3107489/137‑1 (MQ=255)
                                        ttGGCACCACAGCGGGACTGAGGCTGCAAAAAGGTGGGAGCGGGTGTGCCAGAGAAACGGGCGCTGTGGTGCGGGGC‑ACCGGTTTTGCTGGTG‑T‑T‑CGC‑GGGAAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCgg                                                                                                <  2:10918848/137‑1 (MQ=255)
                                        tgggaaccacAGAAGGAGTGGGGCGACAAAAGGGGGGGAGAGGGTGTTCCCGCGAAACAGGCTCCGTGGGGCATGGC‑ACCGGTTTTGCTGGTG‑T‑TCG‑T‑GGGCAC‑CCCCGCTCCCACAGACTGATTGCCATTAAGCgg                                                                                                <  1:9022621/132‑1 (MQ=255)
                                         tGGCACCACAGCGGGATTGAGGTTGCCAAATGGTGGGAGCGGTTGTTCCCGCGAAACGGGCGCCGTGGGGCAGGGA‑ACTGGCTTTGCTGGTG‑T‑TCG‑C‑GGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCggg                                                                                               <  2:1833063/137‑1 (MQ=16)
                                         tGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGTGGGCTTGCCCGCGAAACGGGCGCCGCGGTGCATGGC‑ACTTGCTTTGCCGGTG‑T‑TCG‑C‑GGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCggg                                                                                               <  2:17134834/137‑1 (MQ=21)
                                           gaaccacAGAAGGATGGGGGCTACAAAATGGGGGGAGCGGGTGTTCCCGAGAAACAGGCGCCGTGGGGCGTGGC‑ACTGGTTTTTCCGGTT‑T‑TCG‑C‑GGGCAC‑CCCCGCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                                             <  1:9341229/135‑1 (MQ=255)
                                           gCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGC‑ACCGGCTTTGCCGGTG‑T‑T‑CGA‑GGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                                             <  2:16713815/137‑1 (MQ=21)
                                           acaccccAGCAGGACTGGGTTTGCCAAACGTGGGGATCGGTTTTCCCTGCGAAACGGGCGCCGTGGTGCGTGCC‑ACCGTTTTTACCGGTG‑T‑TCG‑C‑GGGCAC‑ACCCGCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                                             <  1:6619821/136‑1 (MQ=255)
                                              ccacAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGTCCTCGAAACAGGCGCCCCGGTGCCTGGC‑CCCGGCTTTTCCGGTG‑T‑TCG‑C‑GGGTCC‑GCCCGCTCCCACCAAATTATTGGCCTTCAACGGGCCCAt                                                                                          >  1:18222235/1‑137 (MQ=255)
                                                     agGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCATGAAACAGGCGCCGCGGTGCATGGC‑ACCAGCTTTTCCAGTG‑T‑TCG‑C‑GGGCAC‑GCCCGCTCCCACAGAATGATTGCCACTAAGGGGGACTATCGCTaaa                                                                                   >  2:2312529/1‑135 (MQ=18)
                                                           gAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCAGTTCACGGC‑ACCGGCTTTGCCGGTT‑T‑TCG‑C‑GGGCCC‑GCCCGCTCCCACCGACTGATTGCCATTAAGCGGGCCTATCGCCACAGTGAca                                                                             >  2:9880844/1‑136 (MQ=18)
                                                               cTGCCAAACGGTGGGAGCGGGCGTTCCCGCGAAACAGACGCCGCGGTGCATGGC‑ACCGGCTTTGCCGGTG‑T‑TCG‑T‑GGGCAC‑GCCCGCTCCCAAAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGcac                                                                         >  2:7252958/1‑135 (MQ=21)
                                                               cTGCCAAACGGTGGGAGCGGGCGTGTCCGCGAAACAGGCGCCGCGGTGCCTGGC‑CCCGGCTTTGCCGGTG‑T‑T‑TGC‑GGGGAC‑GCCCGCTTCCACAAACTGATTGCCATTAAGCGGGCCTATTGCTACAGCGACGGccc                                                                         >  2:6722324/1‑137 (MQ=37)
                                                               cTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGC‑ACCGGCTTTGCCGGTG‑T‑T‑TGC‑GGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTCAGCGGGCCTATCGCTACAGTGACGGccc                                                                         >  2:12508449/1‑137 (MQ=21)
                                                               cTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAAAAGGCGCCACGGTGCATGGC‑ACCGGATTTGTCGGTG‑T‑T‑TGC‑GGGCACAC‑CCACTCCCACAGACTGGTTTCCCTTCAGCGGAACTGTTGCTGCAGTGAAGGgcc                                                                         >  1:18470448/1‑137 (MQ=9)
                                                                    aaaCGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGTAGTGCATGGC‑ACCAGCTTTGCCGGTG‑T‑TCG‑C‑GGGCAC‑GCCCGCTTCCACAGAATGATTGCCATTAAGCGGGACCACCGCTACAGTGACGGCCCGtttt                                                                    >  2:11839400/1‑137 (MQ=17)
                                                                    aaaCGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGC‑ACCGGCTTTGCCGGTG‑T‑TCG‑C‑GGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCACGtttt                                                                    >  2:9333810/1‑137 (MQ=21)
                                                                    aaaCGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGC‑ACCGGCTTTGCCGGTG‑T‑T‑CGC‑GGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTCAGCGGGCCTCTCGCCACAGTGACGGCCCGtttt                                                                    >  2:13888989/1‑137 (MQ=18)
                                                                     aaCGGTGGGAGCGGGCGTGCCCGCGAAACAAGCGCCGCGGTGCACGGC‑ACCGGCCTTTCCGGTG‑T‑TCG‑C‑GGGCAC‑GCCCGCTCCCACAAAACGATTGACCTTAAGCGGGCCTATCGCTCCAGTGACGGCCCGttttt                                                                   >  2:6925553/1‑137 (MQ=17)
                                                                        tttGGGAGCGGGCGTCCCCGCGAAATAGGCGTCGTGGGGCAGGGC‑ACTGGTTTTGCTGGTG‑T‑TCG‑C‑GGGCAC‑TCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTttgt                                                                <  1:5784409/135‑1 (MQ=255)
                                                                          gggggGCGGGCTTGCCCGCGAAACAGGCGCCGTGGTGCATGGC‑ACTGTCTTTGCTGGTG‑T‑TCG‑T‑GGGCAC‑CCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTttgttg                                                              <  1:13768348/136‑1 (MQ=255)
                                                                            ggAGCGGGTGTGCCCGCGAAACGTTTGCCGCGGTGCAGGGC‑ACCGGTTTTGCCGTTG‑T‑TCG‑C‑GGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGcc                                                            <  1:12398725/137‑1 (MQ=25)
                                                                             gAGCGGGCGTGCCCGCGAAACAAGCGTCCCAGTGCCTGGC‑CCCAGCTTTGTCCGTG‑T‑TCG‑C‑CGGCACAC‑CCGCTCCCACCAACTGATTGTCATTAAGCGGGCCTCCCTCTCCAGTGGCGGACCGTTTTTTGTTGCCt                                                           >  2:400945/1‑137 (MQ=9)
                                                                              agtgggCGTACCCGAGAAACGGGCGCCGTGGTGAATGGC‑ACTGTTTTTGCTGGTT‑T‑TCG‑G‑GGGAAC‑TCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTc                                                          <  2:1195313/134‑1 (MQ=255)
                                                                               gCGGGCGTGCCCACGAAACAGGCGCCGCGGTGCCTGGC‑ACCGGCTTTGCCGGTG‑T‑TCG‑C‑GGGCAC‑GCCCGACCCCACAAAATGGTTGCCATTAAGCGGGGCTATCGCTACAGTGAAGGCCCGGTTTTTGTTGCCTCg                                                         >  2:10089338/1‑137 (MQ=18)
                                                                                 gggTGTGCCTGAGAAACAGGCGCCGCGGTGCGTGGC‑ACCGGCTTTGCCGGTG‑T‑T‑CGC‑GGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCt                                                       <  2:16479936/137‑1 (MQ=25)
                                                                                   gCGTGCCCGCGAAACGTGCGCTGCGGTGCATGGC‑ACTGTTTTTGCTGGTG‑T‑T‑TGT‑GGGCAC‑TCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTaa                                                     <  2:10368299/137‑1 (MQ=255)
                                                                                          tgcgAAACAGGTTCCGTGGGGCAGGGC‑ACTGTTTTTGTGGGTG‑T‑T‑TGG‑GGGCAC‑CCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                              <  2:12828264/136‑1 (MQ=255)
                                                                                          cgcgAAACAGGTGCCGCGGTGCATGGC‑ACCGTCTTTACCGTTG‑T‑TCG‑C‑GGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                              <  2:1977617/137‑1 (MQ=37)
                                                                                          cgcgAAACAGGCGCCGCGGTGCATGGC‑ACCGGCTTTCCCGGTG‑T‑T‑C‑CTGGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                              <  2:15905412/137‑1 (MQ=37)
                                                                                          cgcgAAACAGGCGCCGCGGTGCATGAC‑ACCGGCTTTTCCGTTG‑T‑TCG‑C‑GGGCAC‑CCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                              <  2:15655868/137‑1 (MQ=255)
                                                                                          cgcaAAACGGGCGCCGCGGTGCATGCC‑ACCGTTTTTTCCGTTG‑T‑TCT‑C‑GGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                              <  1:18116023/137‑1 (MQ=255)
                                                                                          cgaaaaaCAGGCGCCGTGGTGCGTGGC‑ACTGGTTTTGCTGGTG‑T‑TCG‑C‑GGACAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                              <  2:18785069/133‑1 (MQ=255)
                                                                                          agagAAACGGGCGCTGGGGTGCAGGGC‑ACTGTTTTTGCTGGTG‑T‑TCG‑C‑GGGCAC‑CCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                              <  1:14902472/134‑1 (MQ=255)
                                                                                           gcgAAACAGGCGCCGCGGTGCATGGC‑ACCGGCTTTGCCGGTG‑T‑TCG‑C‑GGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCAt                                             <  2:17732192/137‑1 (MQ=25)
                                                                                                aaaGGCGCCGCGGTGCATGGC‑ACCGGCTTTGCCGGTG‑T‑TCG‑C‑GGACAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTc                                        <  2:17704421/135‑1 (MQ=37)
                                                                                                     ctccgcGGTGCATGGC‑ACCGGCTTTTCCGGTG‑T‑T‑CGT‑GGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa                                   <  2:6490289/135‑1 (MQ=37)
                                                                                                     cgctgcGGTGCATGGC‑ACCGGCTTTGCTGGTG‑T‑TCG‑C‑GGGCAC‑GCCCGCTCCCACATACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa                                   <  2:17446809/137‑1 (MQ=25)
                                                                                                     cgccgtGGTGCGTGAC‑ACCGTCTTTGCCGGTG‑T‑T‑CGC‑GGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa                                   <  2:7882999/137‑1 (MQ=37)
                                                                                                     cgccgcGGTGCGTGGC‑ACGGGTTTTGCAGGTG‑T‑T‑CGT‑GGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa                                   <  1:10425024/137‑1 (MQ=37)
                                                                                                     cgccgcGGTGCATGGC‑ACCGGCTTTGCCGGTG‑T‑TCG‑C‑GGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa                                   >  2:10681687/1‑137 (MQ=37)
                                                                                                     cgccgcGGGGCATGGC‑ACCGGCTTGTCTGGTG‑T‑TCG‑C‑GGGAAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa                                   <  1:13527620/137‑1 (MQ=255)
                                                                                                               cATGGC‑ACCGGCTTTGCCGGTG‑T‑T‑CGC‑GGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATc                         <  2:6417235/137‑1 (MQ=37)
                                                                                                                      accGGCTTTGCCGGTG‑T‑TCG‑C‑GGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAg                   >  1:14210516/3‑137 (MQ=38)
                                                                                                                        cGGCTTTGCCGGTG‑T‑TCG‑C‑GGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGAt                 >  2:48268/1‑137 (MQ=38)
                                                                                                                        cGGCTTTGCCGGTG‑T‑TCG‑C‑GGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGAt                 >  1:16984923/1‑137 (MQ=38)
                                                                                                                        cGGCTTTGCCGGTG‑T‑TCG‑C‑GGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGAt                 >  1:8418263/1‑137 (MQ=38)
                                                                                                                         ggCTTTGCCGGTG‑T‑TCG‑C‑GGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                              <  1:438360/107‑1 (MQ=37)
                                                                                                                           cTTTGCCGGTG‑T‑TCG‑C‑GGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCGTGACTACCTTAGAACGCTATCAAGCAGATCTg              >  2:15637092/1‑137 (MQ=38)
                                                                                                                           cTTTGCCGGTG‑T‑TCG‑C‑GGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg              >  1:2152626/1‑137 (MQ=39)
                                                                                                                           cTTTGCCGGTG‑T‑TCG‑C‑GGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg              >  2:3657076/1‑137 (MQ=39)
                                                                                                                           cTTTGCCGGTG‑T‑TCG‑C‑GGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg              >  1:15668606/1‑137 (MQ=39)
                                                                                                                           cTTTGCCGGTG‑T‑TCG‑C‑GGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg              >  1:14339335/1‑137 (MQ=39)
                                                                                                                           cTTTGCCGGTG‑T‑TCG‑C‑GGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg              >  1:1260398/1‑137 (MQ=39)
                                                                                                                           cTTTGCCGGTG‑T‑TCG‑C‑GGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg              >  1:12067149/1‑137 (MQ=39)
                                                                                                                           cTTTGCCGGTG‑T‑TCG‑C‑GGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg              >  1:11478245/1‑137 (MQ=39)
                                                                                                                           cTTTGCCGGTG‑T‑TCG‑C‑GGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg              >  1:10052060/1‑137 (MQ=39)
                                                                                                                           cTTTGCCGGTG‑T‑TCG‑C‑GGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg              >  2:12616548/1‑137 (MQ=39)
                                                                                                                           cTTTGCCGGTG‑T‑TCG‑C‑GGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg              >  2:15471196/1‑137 (MQ=39)
                                                                                                                           cTTTGCCGGTG‑T‑TCG‑C‑GGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg              >  2:14727701/1‑137 (MQ=39)
                                                                                                                           cTTTGCCGGTG‑T‑TCG‑C‑GGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg              >  2:13319116/1‑137 (MQ=39)
                                                                                                                           cTTTGCCGGTG‑T‑TCG‑C‑GGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg              >  1:2714525/1‑137 (MQ=39)
                                                                                                                           cTTTGCCGGTG‑T‑TCG‑C‑GGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg              >  1:7480151/1‑137 (MQ=39)
                                                                                                                           cTTTGCCGGTG‑T‑TCG‑C‑GGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg              >  1:5531673/1‑137 (MQ=39)
                                                                                                                           cTTTGCCGGTG‑T‑TCG‑C‑GGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg              >  1:404536/1‑137 (MQ=39)
                                                                                                                           cTTTGCCGGTG‑T‑TCG‑C‑GGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg              >  1:290695/1‑137 (MQ=39)
                                                                                                                             ttGCCGGTG‑T‑TCG‑C‑GGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCACGCAGATCTGaa            >  2:13675086/1‑137 (MQ=39)
                                                                                                                             ttGCCGGTG‑T‑TCG‑C‑GGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGaa            >  2:18726152/1‑137 (MQ=39)
                                                                                                                             ttGCCGGTG‑T‑TCG‑C‑GGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGaa            >  1:1866242/1‑137 (MQ=39)
                                                                                                                             ttGCCGGTG‑T‑TCG‑C‑GGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGaa            >  2:13717972/1‑137 (MQ=39)
                                                                                                                             ttGCCGGTG‑T‑TCG‑C‑GGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGaa            >  1:1923103/1‑137 (MQ=39)
                                                                                                                             ttGCCGGTG‑T‑TCG‑C‑GGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGaa            >  1:10546976/1‑137 (MQ=39)
                                                                                                                             ttGCCGGTG‑T‑TCG‑C‑GGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACG‑ATCCTGCCTCCCTTAGAACGCTATCAAGCAGATCTGaaa           >  1:11209056/1‑137 (MQ=38)
                                                                                                                             ttGCCGGTG‑T‑T‑CGC‑GGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGaa            >  1:10897406/1‑137 (MQ=39)
                                                                                                                                       tgttcgc‑GGGCAC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTccc     <  2:13460277/131‑1 (MQ=255)
                                                                                                                                         tct‑C‑GGGCCC‑GCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGAc  <  1:15480738/135‑1 (MQ=255)
                                                                                                                                           |||                                                                                                                                       
TGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGC‑CACGGCTTTGCCGGTG‑T‑T‑T‑C‑GGGCAC‑G‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGAC  >  AE015451_phaM1/1499363‑1499631

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