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.3% 595.0 / 224.2 272intergenic (+74/‑90)trpS CDS (Tryptophan‑‑tRNA ligase)/zapE CDS (putative enzyme with nucleoside triphosphate hydrolase domain)–/–
Reads supporting (aligned to +/- strand):  ref base . (27/29);  new base C (177/38);  total (204/68)
Fisher's exact test for biased strand distribution p-value = 9.66e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*AE015451_phaM11,499,4980TG84.1% 558.3 / 38.2 272intergenic (+75/‑89)trpS CDS (Tryptophan‑‑tRNA ligase)/zapE CDS (putative enzyme with nucleoside triphosphate hydrolase domain)–/–
Reads supporting (aligned to +/- strand):  ref base T (26/12);  new base G (163/38);  total (204/68)
Fisher's exact test for biased strand distribution p-value = 8.52e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 281 total aligned reads displayed.

TGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTG‑TT‑T‑C‑G‑G‑GCAC‑G‑CCGCTCCCACAGACTGA‑‑‑‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACT  >  AE015451_phaM1/1499363‑1499632
                                                                                                                                         |||                                                                                                                                              
tGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGAGTGCCCGCGAAACAAGCGCCACGGTGCACGGCACCGGATTTTTCGGTG‑tt‑tt                                                                                                                                                 >  1:2853386/1‑136 (MQ=255)
          gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGGGTGCCCGCGAAACAGGCGCCACGGTGCATGGAACCGGCTTTTTCGGTGTTT‑T‑C‑G‑G‑GCAc‑a‑cc                                                                                                                                  >  1:10467994/1‑137 (MQ=255)
          gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCCTGGCACCGGATTTGCCGGTG‑TTCG‑C‑G‑G‑GCAc‑g‑cc                                                                                                                                  >  2:16218334/1‑137 (MQ=25)
          gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCCCTGCTTTGTCGGTG‑TTCG‑C‑G‑G‑GCAc‑g‑cc                                                                                                                                  >  2:10817124/1‑137 (MQ=37)
          gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGAGTGCCCGCGAAACAGGCGCCACGGGGTATGGCACCCGCTTTTTCCGTG‑TTCT‑C‑G‑G‑GCAcaac                                                                                                                                    >  1:5258733/1‑134 (MQ=255)
          gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGAGTGCCCGCGAAAAAAGAGCCGCGGTGCATGGGACCAGATTTTCCGGTG‑TTCT‑T‑G‑G‑GCAc‑a‑cc                                                                                                                                  >  1:10175655/1‑137 (MQ=255)
            tGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAAGCGCCGCGGTGCATGGCAACAGCTTTTCCCGTG‑TTCG‑C‑G‑G‑GCAc‑g‑cccg                                                                                                                                >  1:10938027/1‑135 (MQ=25)
                 ccATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCCTGGCACCGGCTTTGCCGGTG‑TT‑CGC‑G‑G‑GCAC‑GCCCGCTccc                                                                                                                            >  2:6751384/1‑137 (MQ=21)
                 ccATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGAGTGCCCGCAAAACAAGCGCCATGGTGCACGGCACCAGCTTTTTCGGTGTTT‑T‑C‑G‑G‑GCAC‑GCCCGCTccc                                                                                                                            >  1:13872817/1‑137 (MQ=255)
                  cATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCACGGTGCATGGCACCGGATTTGCCGGTG‑TTCG‑C‑G‑G‑GAAC‑GCCCGCTcccc                                                                                                                           >  1:4175936/1‑136 (MQ=21)
                   aTTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGGGTGTCCGCGAAAAAGGAGCCGCGGTGTATGGAACCGGCTTTGTCGGTG‑TT‑TGC‑G‑G‑GCAC‑GACCGCTCCcac                                                                                                                          >  1:12722700/1‑137 (MQ=255)
                        gatccTGAATGAATCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAACGGGCTTGCCCGCGAAACAGGCGCCGCGGTACATGGCACCGGCTTTGCCGGTG‑TTCG‑C‑G‑G‑GCAC‑GCCCGCTCCCACAGACt                                                                                                                     <  2:10373625/134‑1 (MQ=18)
                        gagCCTGAACGAAGCCTTGGCACCACAGCAGGACGGAGGCTGCAAAATGGGGGGAGCGGGTGTGCCTGCGAAACGGGTGCCGTGGTGCGTGGCACTGGCTTTGCCGGTG‑TT‑CGG‑G‑G‑GAAC‑GCCCGCTCCCACAGACt                                                                                                                     <  2:1892168/137‑1 (MQ=255)
                        gagCCGGAACGAGGCCTTGGCACCACAGCAGGACTAAGGCTGCCAAACGGTGAGAGCGGGCGTGCCATCGAAACAGGCGCCGCGGTGCATGGCACCGTCTTTGCCGGTG‑TT‑CGC‑G‑G‑GCAC‑GCCCGCTCCCACAGACt                                                                                                                     <  2:5469515/137‑1 (MQ=18)
                              gAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑TTCG‑C‑G‑G‑GCAC‑GCCCGCTCCCACAGACTGA‑‑‑‑TTGc                                                                                                           >  2:10242089/1‑137 (MQ=21)
                              gAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑TTCG‑C‑G‑G‑GCAC‑GCCCGCTCCCACAGAATGA‑‑‑‑TTGc                                                                                                           >  2:3069670/1‑137 (MQ=21)
                              gAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑TT‑CGC‑G‑G‑GCAC‑GCCCGCTCCCAAAGACTGA‑‑‑‑TTGc                                                                                                           >  2:12613622/1‑137 (MQ=21)
                              gAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCAAAAAAGGCGCCACGGTGCCTGGCACCGGCTTTGCCGGTG‑TTCG‑C‑G‑G‑GCAC‑GCCCGCTCCCACAAACTGA‑‑‑‑TTGc                                                                                                           >  1:14132766/1‑137 (MQ=18)
                                  gaggtcTTGGCACCACAGAAGGACGGAGGCTGCAAAATGGTGGGATTGGGCGTTCCCGCGAAACAGGTTCCGTGGGGCATGGAACTGGTTTGGTTGGTG‑TTCG‑T‑G‑G‑GCAC‑GCCCGCTCCCACAGACTGA‑‑‑‑TTGCCAtt                                                                                                       <  2:8813186/132‑1 (MQ=255)
                                  gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGTGGGCGTGCCCGCGGAACAGGCGCCGAGGTGCATGGCACCGGCTTTGCCGGTG‑TTCG‑C‑G‑G‑GCAC‑GCCCGATCCCACGGGCTGA‑‑‑‑GTGCCagt                                                                                                       >  2:11254793/1‑135 (MQ=18)
                                  gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGTCCCCGAAACAGGAGCCACGGTGTCTGGCCCCAGCTTTGCCGGTG‑TT‑TGC‑G‑G‑GCCCAC‑CCACTCCCCCAGACCGG‑‑‑‑TTGCCact                                                                                                       >  1:13507491/1‑135 (MQ=255)
                                  gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑TTCG‑C‑G‑G‑GCAC‑GCCCGCTCCCACAGACTGA‑‑‑‑TTGCCact                                                                                                       >  2:3671859/1‑135 (MQ=21)
                                  gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAAAAGGCGCCGCAGGGTATGGCAACGGATTTGCCGG‑G‑AT‑T‑C‑GTGAGCGC‑GCCCGATCCCACAGACCGA‑‑‑‑TTGCCAtt                                                                                                       >  2:279590/1‑137 (MQ=19)
                                  gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAAAAGGCGCCGAAGTGCATGGCACCGGGTTTGACGGTG‑TTCT‑C‑G‑G‑GACC‑GTCCGCTCCCAAAGAACGA‑‑‑‑TTGCCAtt                                                                                                       >  2:7117617/1‑137 (MQ=21)
                                  gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCAAAATGGTTGAAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCGGGGCACCGGCTTTGCCGGTG‑TT‑CGC‑G‑G‑GCAC‑GCCCGCTCCCACAGACTGA‑‑‑‑TTGCCAtt                                                                                                       <  2:15150973/137‑1 (MQ=18)
                                      ccTTGGCACCACAACAGGACTAAGACTGCAAAATGGGGGGAGTGGGTGTGCCCGAAAAACAGGCGCCGTGGTGAGTGGAACTGGTTTTGCTGGTG‑TT‑CGT‑G‑G‑GCAC‑GCCCGCTCCCACAGACTGA‑‑‑‑TTGCCATTAAGc                                                                                                   <  1:7558198/137‑1 (MQ=255)
                                           tcaccacAGCAGGACTGGGGTTGCCAAACTGTGGGGGCGGGCGTACCCGCGAAATAGGCGCCGTGGTGCGGGGCACCGGTTTTGCTGGTG‑TTCG‑C‑G‑G‑GCAC‑GCCCGCTCCCACAGACTGA‑‑‑‑TTGCCATTAAGCGGGcc                                                                                              <  2:8464180/136‑1 (MQ=16)
                                                   gcagGACTGAGGCTGCCAAACGGTGGGAGCGGGGGTGCCCACGAAACAAAAGCCACGGTGTATGGCACCAACTTTGTCGGTT‑TTCT‑CAG‑G‑CCAA‑C‑CCACACCCACAAAATGA‑‑‑‑CTGCCACTAAGCTGGACTACCGCtc                                                                                      >  1:13032435/1‑136 (MQ=255)
                                                   gcagGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCCTGGCACCGGCTTTGCCGGTG‑TTCA‑C‑G‑G‑GAAC‑GCCCGCTCCCAAAGACTGA‑‑‑‑TTGCCATTCAGCGGGCCTATCGCTa                                                                                      >  2:3018150/1‑137 (MQ=21)
                                                        aCTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCCCGAAACAAGCGCCACGGTGCCTGGCACCAGCCTTTTCAGTG‑TTCG‑C‑A‑G‑GCAC‑GCCCGCTCCCACAAAATGA‑‑‑‑TTGCCACTCAGCGGGACCATCGCTCCAATg                                                                                 >  2:3006359/1‑137 (MQ=12)
                                                        aCTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCACGAAACAGGCGCCGCGGTGCATGGCACCGGATTTGCCGGTG‑TT‑TGC‑G‑G‑GCAC‑GCCCGCTCCCACAGACTGA‑‑‑‑TTGCCATTCAGCGGGACTATCGCTACAGTg                                                                                 >  2:7547334/1‑137 (MQ=21)
                                                            aGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCCTGGCCCCTGCTTTGTCGGTG‑TTCG‑C‑G‑G‑GCCC‑GCCCGCTCCCAAAGACTGA‑‑‑‑TTGCCATTAAGCGGGTCTACCGCTACAGTGACgg                                                                             >  2:10411353/1‑137 (MQ=21)
                                                               cTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGTGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑TT‑TGT‑G‑G‑GGAC‑GCCCGCTCCCACAGGCTGA‑‑‑‑TTGCCATTAAGCGGGCCTATCGCTACAGTGAAGGccc                                                                          >  2:8942531/1‑137 (MQ=21)
                                                               cTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑TTCG‑C‑G‑G‑GCAC‑GCCCGCTCCCACAGAATGA‑‑‑‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGccc                                                                          >  2:14058325/1‑137 (MQ=21)
                                                               cTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCACGGTGCACGGCACCAGCTTTGCCGGTG‑TTCG‑C‑G‑G‑GCAC‑GCCCGCTCCCACAGACTGA‑‑‑‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGccc                                                                          >  2:15911107/1‑137 (MQ=21)
                                                               cTGCCAAACGGTGGGAGCGGGCGTGCCCACGAAACAAGCACCGCGGTGCATGGCACCGGATTTGCCGGTG‑TTCG‑C‑G‑G‑GCAC‑GCCCGCTCCCAAAGAATGA‑‑‑‑TTGCAATTAAGCGGGCCGACCACTCAAGTGACGGacc                                                                          >  2:6453055/1‑137 (MQ=9)
                                                                    aaaCGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑TTCG‑C‑G‑G‑GCAC‑GCCCGCTCCCACAGAATGATTGGTTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCg                                                                         >  2:7210534/1‑137 (MQ=17)
                                                                    aaaCGGTGGGAGCGGGCGTGCCCGCAAAACAAGAGCCGCGGTGCACGGCACCGGCTTTGCCGGTG‑TTCG‑C‑G‑G‑GCAC‑GCCCGCTCCCCCCGACTGA‑‑‑‑TTGCCATTCAGGGGGACTATCGCCCCCGTGACGGCCCGtttt                                                                     >  2:5872676/1‑137 (MQ=9)
                                                                     aaCGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGTATGGCACCGGCTTTGCCGGTG‑TTCG‑C‑G‑G‑GCAC‑GCCCGCTCCCACAGACTGA‑‑‑‑TTGCCATTCAGCGGGCCTCTCGCTACAGTGACGGCCCGttttt                                                                    >  2:15598995/1‑137 (MQ=18)
                                                                               gCGGGCGTGCCCGCGAAACAGGGGCCACGGTGCCTGGCACCAGCTTTGCCCGTG‑TT‑TGC‑G‑G‑GAAC‑GCCCCCTCCCCCAAACTGA‑‑‑‑TTGCCATTAAGCAGGCCTACCTCTCCAACGAAAGCCCATTTTTTTTTTCCTcc                                                          >  2:9746565/1‑136 (MQ=9)
                                                                               gCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑TTCG‑C‑G‑G‑GCAC‑GCCCGCTCCCACAGACTGA‑‑‑‑TTGCCATTAAGCGGGCCTCTTGCTAACGTGAAGGACCGTTTTTTGTTGCCTCg                                                          >  2:12043368/1‑137 (MQ=18)
                                                                                cGGGCTTGCCCGCGAAACGGGCGCCGCGGTACATGGCACCGGTTTTTCTGGTG‑TT‑TGC‑G‑G‑GCAC‑GCCCGCTCCCACAGACTGA‑‑‑‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGc                                                         <  2:2506023/137‑1 (MQ=25)
                                                                                agggCGTACCAGAGAAACAGGCGCCGTGGTGCAGGGCACCGGTTTTGCTGGTG‑TTCG‑C‑G‑G‑ACAC‑TCCCGCTCCCACAGACTGA‑‑‑‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGc                                                         <  1:9071084/136‑1 (MQ=255)
                                                                                        cccGAGAAACAGTCGCCGCGGTGCATGGCACTGGCTTTGCTGGTG‑TTCG‑C‑G‑G‑GCAC‑GCCCGCTCCCACAGACTGA‑‑‑‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAAt                                                 <  2:4914145/137‑1 (MQ=25)
                                                                                          ggcgAAAAAGGCTCCGCGGGGCATGGCACCGTTTTTGCTGGTG‑TT‑CGC‑G‑G‑GAAC‑GCCCGCTCCCACAGACTGA‑‑‑‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                               <  2:4527606/136‑1 (MQ=255)
                                                                                          cgcgAAACAGGCGCTGCGGTGCATGGCACCGGCTTTGCCGGGG‑TTCG‑C‑G‑G‑GAAC‑GCCCGCTCCCACAGACTGA‑‑‑‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCTTTTTTTGTTGCCTCGCTAACGAATCa                                               <  2:10674676/137‑1 (MQ=25)
                                                                                          cgcgAAAAAGTCGCCGCGGGGCGTGGCACCGTTTTGGCCGGTG‑TTCG‑C‑G‑G‑GCAC‑GCCCGCTCCCACAGACTGA‑‑‑‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                               <  1:16120664/137‑1 (MQ=25)
                                                                                          cgaaaaaCAGGCGCCGTGGTGCGGGGCACTGGTTTTGCCGGTG‑TT‑TGC‑G‑G‑GCAC‑GCCCGCTCCCACAGACTGA‑‑‑‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                               <  1:200421/133‑1 (MQ=37)
                                                                                          cgaaaaaCAGGCGCCGTGGTGCATGGCACCGTCTTTGCTGGTG‑TTCG‑C‑G‑G‑ACAC‑CCCCGCTCCCACAGACTGA‑‑‑‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGATAACGAATCa                                               <  1:14170624/133‑1 (MQ=255)
                                                                                               aaCAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑TTCG‑C‑G‑G‑GCAC‑GCCCGCTCCCACAGACTGA‑‑‑‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAACCATGACt                                          >  1:4669988/1‑137 (MQ=25)
                                                                                                     cgccggGTTGCGTGGCACCGGTTTTTCGGGTG‑TT‑CGC‑G‑G‑GCAT‑GCCCGCTCCCACAGACTGA‑‑‑‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa                                    <  2:6252241/137‑1 (MQ=25)
                                                                                                       ccgcGGTGCATGGCACCGGCTTTGCCGGTG‑TTCG‑C‑G‑G‑GCAC‑GCCCGCTCCCACAGACTGA‑‑‑‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGa                                  >  2:5109519/1‑137 (MQ=37)
                                                                                                          cgGTGCATGGCACCGGCTTTGCCGGTG‑TTCG‑C‑G‑G‑GCAC‑GCCCGCTCCCACAGACTGA‑‑‑‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACg                               >  2:13207512/1‑137 (MQ=37)
                                                                                                          cgGTGCATGGCACCGGCTTTGCCGGTG‑TTCG‑C‑G‑G‑GCAC‑GCCCGCTCCCACAGACTGA‑‑‑‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACg                               >  1:3694153/1‑137 (MQ=37)
                                                                                                               cATGGCACCGGTTTTGCCGGTG‑TT‑CGC‑G‑G‑GCAC‑GCCCGCTCCCACAGACTGA‑‑‑‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATc                          <  2:6068448/137‑1 (MQ=37)
                                                                                                               cATGGCACCGGCTTTTCTGGGG‑TTCT‑C‑G‑G‑GCAC‑TCCCGCTCCCACAGACTGA‑‑‑‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATc                          <  2:2669011/137‑1 (MQ=255)
                                                                                                               cATGGCACCGGCTTTGCTGGTG‑TTCG‑C‑G‑G‑GCAC‑GCCCGCTCCCACAGACTGA‑‑‑‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATc                          <  2:4006194/137‑1 (MQ=37)
                                                                                                               aaTGGCACCGGCTTTGCTGGTG‑TTCG‑C‑G‑G‑CCAC‑GCCCGCTCCCACAGACTGA‑‑‑‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATc                          <  1:625757/136‑1 (MQ=255)
                                                                                                                  ggcaccGGCTTTGCCGGTG‑TTCG‑C‑G‑G‑GCAC‑GCCCGCTCCCACAGACTGA‑‑‑‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAg                       >  1:15781094/6‑137 (MQ=38)
                                                                                                                  ggcaccGGCTTTGCCGGTG‑TTCG‑C‑G‑G‑GCAC‑GCCCGCTCCCACAGACTGA‑‑‑‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAg                       >  2:13032248/6‑137 (MQ=37)
                                                                                                                  ggcaccGGCTTTGCCGGTG‑TT‑CGC‑G‑G‑GCAC‑GCCCGCTCCCACAGACTGG‑‑‑‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATTATGACTCCCTTAGAGCGCTATCAAg                       >  2:11642421/6‑137 (MQ=37)
                                                                                                                     accGGCTTTGCCGGTG‑TTCG‑C‑G‑G‑GCAC‑GCCCGCTCCCACAGACTGA‑‑‑‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                               <  2:9574558/108‑1 (MQ=37)
                                                                                                                     accGGCTTTGCCGGTG‑TTCG‑C‑G‑G‑GCAC‑GCCCGCTCCCACAGACTGA‑‑‑‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                               >  1:9574558/3‑110 (MQ=37)
                                                                                                                     accGGCTTTGCCGGTG‑TTCG‑C‑G‑G‑GCAC‑GCCCGCTCCCACAGACTGA‑‑‑‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAg                    >  2:1441595/3‑137 (MQ=38)
                                                                                                                     accGGCTTTGCCGGTG‑TT‑CGC‑G‑G‑GCAC‑GCCCGCTCCCACAGACTGA‑‑‑‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAg                    >  2:14656286/3‑137 (MQ=38)
                                                                                                                       cGGCTTTGCCGGTG‑TTCG‑C‑G‑G‑GCAC‑GCCCGCTCCCACAGACTGA‑‑‑‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGAt                  >  1:15548959/1‑137 (MQ=38)
                                                                                                                       cGGCTTTGCCGGTG‑TTCG‑C‑G‑G‑GCAC‑GCCCGCTCCCACAGACTGA‑‑‑‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGAt                  >  1:5827627/1‑137 (MQ=38)
                                                                                                                        ggttttGCTGGTG‑TT‑CGC‑G‑G‑GCAC‑GCCCGCTCCCACAGACTGA‑‑‑‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATc                 <  1:12254970/134‑1 (MQ=255)
                                                                                                                          cTTTGCCGGTG‑TTCG‑C‑G‑G‑GCAC‑GCCCGCTCCCACAGACTGA‑‑‑‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                               >  2:11788693/1‑105 (MQ=37)
                                                                                                                          cTTTGCCGGTG‑TTCG‑C‑G‑G‑GCAC‑GCCCGCTCCCACAGACTGA‑‑‑‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                               <  2:912627/105‑1 (MQ=37)
                                                                                                                          cTTTGCCGGTG‑TTCG‑C‑G‑G‑GCAC‑GCCCGCTCCCACAGACTGA‑‑‑‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                               <  1:11788693/105‑1 (MQ=37)
                                                                                                                          cTTTGCCGGTG‑TTCG‑C‑G‑G‑GCAC‑GCCCGCTCCCACAGACTGA‑‑‑‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg               >  2:4435104/1‑137 (MQ=39)
                                                                                                                          cTTTGCCGGTG‑TTCG‑C‑G‑G‑GCAC‑GCCCGCTCCCACAGACTGA‑‑‑‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg               >  1:15300761/1‑137 (MQ=39)
                                                                                                                          cTTTGCCGGTG‑TTCG‑C‑G‑G‑GCAC‑GCCCGCTCCCACAGACTGA‑‑‑‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg               >  2:10565938/1‑137 (MQ=39)
                                                                                                                          cTTTGCCGGTG‑TTCG‑C‑G‑G‑GCAC‑GCCCGCTCCCACAGACTGA‑‑‑‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg               >  2:12808071/1‑137 (MQ=39)
                                                                                                                          cTTTGCCGGTG‑TTCG‑C‑G‑G‑GCAC‑GCCCGCTCCCACAGACTGA‑‑‑‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg               >  2:5548392/1‑137 (MQ=39)
                                                                                                                          cTTTGCCGGTG‑TTCG‑C‑G‑G‑GCAC‑GCCCGCTCCCACAGACTGA‑‑‑‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg               >  2:13143096/1‑137 (MQ=39)
                                                                                                                          cTTTGCCGGTG‑TTCG‑C‑G‑G‑GCAC‑GCCCGCTCCCACAGACTGA‑‑‑‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg               >  1:12955628/1‑137 (MQ=39)
                                                                                                                          cTTTGCCGGTG‑TTCG‑C‑G‑G‑GCAC‑GCCCGCTCCCACAGACTGA‑‑‑‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg               >  2:14893039/1‑137 (MQ=39)
                                                                                                                          cTTTGCCGGTG‑TTCG‑C‑G‑G‑GCAC‑GCCCGCTCCCACAGACTGA‑‑‑‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg               >  2:4198291/1‑137 (MQ=39)
                                                                                                                          cTTTGCCGGTG‑TTCG‑C‑G‑G‑GCAC‑GCCCGCTCCCACAGACTGA‑‑‑‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg               >  2:3497560/1‑137 (MQ=39)
                                                                                                                          cTTTGCCGGTG‑TTCG‑C‑G‑G‑GCAC‑GCCCGCTCCCACAGACTGA‑‑‑‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg               >  2:15454048/1‑137 (MQ=39)
                                                                                                                          cTTTGCCGGTG‑TTCG‑C‑G‑G‑GCAC‑GCCCGCTCCCACAGACTGA‑‑‑‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg               >  1:322314/1‑137 (MQ=39)
                                                                                                                          cTTTGCCGGTG‑TTCG‑C‑G‑G‑GCAC‑GCCCGCTCCCACAGACTGA‑‑‑‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg               >  2:16350050/1‑137 (MQ=39)
                                                                                                                          cTTTGCCGGTG‑TTCG‑C‑G‑G‑GCAC‑GCCCGCTCCCACAGACTGA‑‑‑‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg               >  2:10004998/1‑137 (MQ=39)
                                                                                                                          cTTTGCCGGTG‑TTCG‑C‑G‑G‑GCAC‑GCCCGCTCCCACAGACTGA‑‑‑‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg               >  2:2347117/1‑137 (MQ=39)
                                                                                                                          cTTTGCCGGTG‑TT‑CGC‑G‑G‑GCAC‑GCCCGCTCCCACAGACTGA‑‑‑‑TTTCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGACCCATGACTCCCTTAGAACGCTATCAAGCAGATCTg               >  1:6705231/1‑137 (MQ=38)
                                                                                                                          cTTTGCCGGTG‑TT‑CGC‑G‑G‑GCAC‑GCCCGCTCCCACAGACTGA‑‑‑‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg               >  1:8389082/1‑137 (MQ=39)
                                                                                                                            ttGCCGGTG‑TTCG‑C‑G‑G‑GCAC‑GCCCGCTCCCACAGACTGA‑‑‑‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGaa             >  2:2451937/1‑137 (MQ=39)
                                                                                                                            ttGCCGGTG‑TTCG‑C‑G‑G‑GCAC‑GCCCGCTCCCACAGACTGA‑‑‑‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGaa             >  2:7877416/1‑137 (MQ=39)
                                                                                                                            ttGCCGGTG‑TTCG‑C‑G‑G‑GCAC‑GCCCGCTCCCACAGACTGA‑‑‑‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGaa             >  2:6664623/1‑137 (MQ=39)
                                                                                                                            ttGCCGGTG‑TTCG‑C‑G‑G‑GCAC‑GCCCGCTCCCACAGACTGA‑‑‑‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGaa             >  2:5330255/1‑137 (MQ=39)
                                                                                                                            ttGCCGGTG‑TTCG‑C‑G‑G‑GCAC‑GCCCGCTCCCACAGACTGA‑‑‑‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGaa             >  2:4091781/1‑137 (MQ=39)
                                                                                                                            ttGCCGGTG‑TTCG‑C‑G‑G‑GCAC‑GCCCGCTCCCACAGACTGA‑‑‑‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGaa             >  2:16615415/1‑137 (MQ=39)
                                                                                                                            ttGCCGGTG‑TTCG‑C‑G‑G‑GCAC‑GCCCGCTCCCACAGACTGA‑‑‑‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGaa             >  1:4015640/1‑137 (MQ=39)
                                                                                                                            gtGCCGGTG‑TTCG‑C‑G‑G‑GCAC‑GCCCGCTCCCACAGACTGA‑‑‑‑TTGCCATTACGCGGACCTATCGCTACAGTGACGGCCTGGTTTTTGTTGCCTAGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTgca             >  2:15303068/2‑135 (MQ=37)
                                                                                                                             tGCCGGTG‑TTCG‑C‑G‑G‑GCAC‑GCCCGCTCCCACAGACTGA‑‑‑‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGaaa            >  1:663542/1‑137 (MQ=39)
                                                                                                                                     ttttcgc‑G‑G‑GCAC‑GCCCGCTCCCACAGACTGA‑‑‑‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTccc      <  1:2533061/131‑1 (MQ=255)
                                                                                                                                        tcgc‑G‑G‑GCAC‑GCCCGCTCCCACAGACTGA‑‑‑‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGAc   >  2:9817597/4‑137 (MQ=255)
                                                                                                                                         cgc‑G‑G‑GCAC‑GCCCGCTCCCACAGACTGA‑‑‑‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGAct  >  1:7257831/3‑137 (MQ=255)
                                                                                                                                         |||                                                                                                                                              
TGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTG‑TT‑T‑C‑G‑G‑GCAC‑G‑CCGCTCCCACAGACTGA‑‑‑‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACT  >  AE015451_phaM1/1499363‑1499632

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