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.C77.7% 573.3 / 277.5 313intergenic (+74/‑90)trpS CDS (Tryptophan‑‑tRNA ligase)/zapE CDS (putative enzyme with nucleoside triphosphate hydrolase domain)–/–
Reads supporting (aligned to +/- strand):  ref base . (31/38);  new base C (172/68);  total (203/110)
Fisher's exact test for biased strand distribution p-value = 8.19e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*AE015451_phaM11,499,4980TG84.1% 532.9 / 36.3 313intergenic (+75/‑89)trpS CDS (Tryptophan‑‑tRNA ligase)/zapE CDS (putative enzyme with nucleoside triphosphate hydrolase domain)–/–
Reads supporting (aligned to +/- strand):  ref base T (33/11);  new base G (159/74);  total (203/110)
Fisher's exact test for biased strand distribution p-value = 4.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 324 total aligned reads displayed.

TGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTGTTTCGGGCACGCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACT  >  AE015451_phaM1/1499363‑1499632
                                                                                                                                           |||                                                                                                                                        
  aatcccTTGATGCGGCCGTTAAGAGCCTGAACGAATCCTTGGAACCACAGCAGGACTGAGGCTGCAAACCGGTGGGAGCGGGCGTGCCTGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGtttgcg                                                                                                                                         <  2:358662/134‑4 (MQ=21)
       cGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCAc                                                                                                                                   >  2:6058457/1‑137 (MQ=37)
          gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGGGTGTCCGCGAAACAAGAGCCGCGGTGCATGGCACCGGTTTTTCCGGTGTTTACGGGCAcgcc                                                                                                                              >  1:6044929/1‑137 (MQ=255)
          gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGTCGGTGTTCACGGGCAcgcc                                                                                                                              >  2:9254722/1‑137 (MQ=25)
          gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCACGAAACAGGCGCCGCGGTGCACGGCACCGGCCTTGCCAGTGTTCGCGGGCAcgcc                                                                                                                              >  2:16777977/1‑137 (MQ=25)
          gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGAGTGGCCGCGAAAAAAGAGCCACAGTGCCCGGCACCGAATTTTCCAGTGTTTTCGGGCAcacc                                                                                                                              >  1:7548722/1‑137 (MQ=255)
            tGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGTCGCGGTGCATGGCACCGGGCTTGCCGGTGTTCGCGGGCAcgtccg                                                                                                                             >  2:11577767/1‑133 (MQ=25)
                 ctaTTCAGAGCTTGAACGAAGCCTGGGCACCACAGCAGGACTGAGGCTGCCAAACGGGGGGAGTGGGTGTGCCCGCGAAACGGGCGCTGCTGGGCCTGGCACTGGTTTTGCTGGTGTTAGCGGGCACGCCCGCTccc                                                                                                                        <  2:12141575/135‑1 (MQ=255)
                   aTTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCAGTGCACGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCcac                                                                                                                      >  2:14394989/1‑137 (MQ=21)
                   aTTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCAAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGTCCGCTCCcac                                                                                                                      >  2:6933881/1‑137 (MQ=21)
                              gAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCAAAAAAAGCGCCACGGTGCACGGCCCCTGATCTGCCGGTTTTCGCGGGCACGCCCGCTCACCCAGAATGATTGc                                                                                                           >  1:14574739/1‑137 (MQ=19)
                                aCGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCATGCCCACGAAACAAGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGTTCCCACAAACTGAGTGCCa                                                                                                         >  2:16935062/1‑137 (MQ=18)
                                 cGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCAt                                                                                                        >  2:4475608/1‑137 (MQ=21)
                                 cGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAAAAAGCGCCACGGTGCCTGGCACCGGCTTTGTCGGTGTTCTCGGGCCCGCCCGCTCCCACAGACTGATTGCCAt                                                                                                        >  2:10207892/1‑137 (MQ=21)
                                 cGAAGCCTTGGCACCAAAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGCAACAAGCGCCTCGGTGCATGGCAACGGTTTTGCCAGTATTCGCGGGCACAACCGCTCCCACAGACCGGTTTCCAt                                                                                                        >  2:12339615/1‑137 (MQ=18)
                                  ggcaccTTGTCACCACAATAATACTGAGTTGACCAAATGTTGGAGGCGGGAGTGCCCGCGAAAAAGGCGCCGTGGTGCAGGGCCCTGGTTGTTTAGGTGTTCGGGGACAAGCCCGCTCCCACAGACTGATTGCCAtt                                                                                                       <  2:5112977/133‑1 (MQ=255)
                                  gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGTGTGTCCGCGAAAAAGGCGCCACGGTGCCTGGCACCAGCTTTTTCGGTGTTCTCAGGCACACCCACCCCCACAAAATGATTGCCAtt                                                                                                       >  1:14486795/1‑137 (MQ=255)
                                  gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGGGTGTCCGCGAAACAAGAGCCGCAGTGCCTGGCACCGGCTTTTCCGGTGTTTTTGGGCACGCCCGCTCCCCCAAAACGAGTGCCact                                                                                                       >  1:13031249/1‑135 (MQ=255)
                                  gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGGGTGTCCGCGAAAAAAGCGCCGTGGTGCCCGGCACCTGCTTTGCCGGTGTTTGCGGGCACGCCCGCTCCCACCAACCGAGTGCCAtt                                                                                                       >  1:17826374/1‑137 (MQ=19)
                                  gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCTCGAAAAAGGAGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCAtt                                                                                                       >  1:13518502/1‑137 (MQ=21)
                                  gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGTATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCAtt                                                                                                       >  2:17841217/1‑137 (MQ=21)
                                  gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCACGAAACAGGCGCCGCGGTGCCTGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGAATGATTGCCAtt                                                                                                       >  2:17929630/1‑137 (MQ=21)
                                      ccTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAAAAGGCGCCGCGGTGCATGGCAACGGATTTGCCAGTGTTCGCGGGCACGCCCGCTCCCACAGACCGATTGCCATTCAGc                                                                                                   >  2:17165786/1‑137 (MQ=21)
                                      ccTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCACGAAAAAAGCGCCACGGTGCACGGCACCGGATTTTTCGGTGTTCTCGGGCACGCCCACTCCCACAAACTGATTGCCCTTAAga                                                                                                   >  1:9116470/1‑136 (MQ=255)
                                        ttGGCACCACAGAAGGATTGGGGTGGCAAAATGGTGGGAGTGGGTGTTCCCGAGAAACAGGCGCTGGGGGGAGGGGCACGGGTTTTGCTGGTGTTTGGGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCgg                                                                                                 <  1:1105792/137‑1 (MQ=255)
                                           gCACCACAGCAGGACTGAGGCTGCCAAACGGGGGGAGCGGGCGGGCCCGCGAAAAAGGCGCCGCGGTGCATGGCACTGGCTTTCCCGGTGTTCTCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                                              <  2:2605680/137‑1 (MQ=25)
                                           gCACCAAAGCAGGACTGAGGTTGCCAAACGGTGGGAGCGGGCGTGCCTGCGAAAAGGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                                              <  2:4859200/137‑1 (MQ=18)
                                           gCACCAAAGAAGGATTGAGGCGGCAAAATGGTGGGAGTGGGTGTTCCCGCGAAAAGGGCCCGGTGGTGAGGGGAACTGGTTTTGCTGGTGTTTGTGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                                              <  1:18697115/137‑1 (MQ=255)
                                           gCACAACAGCAGGATTGGGGTTGCCAAATGGTGGGACTGGGCGGGCCCGGGAAACGGGCTCTGGGGGGCGGGGCACTGGTTTTGCTGGTGTTTGTGGGCACTCCCGCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                                              <  1:13896797/137‑1 (MQ=255)
                                            caccacAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGGCCGTGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCAGGCACACCCGCCCCCACAGACTGATTGCCATTCAGCGGGCCt                                                                                             >  2:3166273/1‑137 (MQ=21)
                                            caccacAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGTCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCCTTCAGCGGGCCt                                                                                             >  2:13300600/1‑137 (MQ=18)
                                            caccacAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGACt                                                                                             >  2:19088371/1‑137 (MQ=21)
                                            caccacAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAAAAGGCGCCGCGGTGCCTGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGATCCCACAGACTGATTGCCCTTAAGCGGGctt                                                                                             >  2:8626533/1‑135 (MQ=18)
                                            aaccacAGAAGGACGGAGGCTACAAAATGGGGGGAGCGGGTGTGCCTGCGAAACGGGCGCCGCGGTGCATGGAACCGGTTTTGCTGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCt                                                                                             <  2:7794666/136‑1 (MQ=16)
                                             accacAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAAGCGCCACGGTGCACGGCACCAAATTTTCCAGTGTTCGCGGGCACGCCCACCCCCCCAGACTGATTGCCATTCAGCGGGCCTa                                                                                            >  1:7062510/1‑137 (MQ=16)
                                                acaGCAGGACTAAGGCTGCCAAACGGTGGGAATGGGCGTTCCCGCAAAACATGCGCCGCGGTACATGACACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCg                                                                                         <  2:7616976/137‑1 (MQ=21)
                                                acaGCAGAACGGAGGTGGCAAAAGGGGGGGGGTGGGTGTGCCCGAAAAACAGGCACCGGGGGGCATGGCACCGGTTTTGCTGGTGTCCGAGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCg                                                                                         <  1:18358297/137‑1 (MQ=12)
                                                   gcagGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTa                                                                                      >  2:9120134/1‑137 (MQ=21)
                                                   gcagGACTGAGGCTGCCAAACGGTGGGAGCGGGAGTGCCCGCGAAACAAGCGCCGCGGTGTACGGCACCAGCTTTTCCGGTGTTCGCGGGCACACCCGCTCCCCCAAACCGACTTCCACTCAGCGGGCCTATCGCtc                                                                                      >  2:1613244/1‑136 (MQ=255)
                                                       gACTGAGGCTGCCAAACGGTGGGAGCGGGAGTGCCCACGAAAAAAGCGCCGCGGTGCATGGCCCCGGCCTTTCCGGTGTTCGCGGGCACACCCACCCCCACAAACCGATTTCCCTTCAACAGGCCTCTCGCTAccat                                                                                  >  1:9433310/1‑134 (MQ=255)
                                                          tGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAAGAGCCACGGTTTATTGCCCCCGCTTTTCCGGTGTTCGCTGGCACGCCCGCTCCCACAGACCGATTGCCATTCAGCGGGCCTGTTGCTACCGTGAc                                                                               >  2:14682684/1‑137 (MQ=12)
                                                               ctacaaaaCGGGGGGGGGGGGGGTGCCCGAGAAACAGGCGCTGCGGGGAATGGCACTGTTTTTTCCGGTGTTCGTGGGCACCCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGccc                                                                          <  1:17479890/132‑1 (MQ=255)
                                                               cTGCCAAACGGTGGGAGCGGGCGTGCCCTCGAAACAGGCGCCACGGTGTATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAAACTGATTGCCATTCAGCGGGCCTCTCGCTCCAATGACGGccc                                                                          >  2:14412365/1‑137 (MQ=17)
                                                               cTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGccc                                                                          >  2:19273542/1‑137 (MQ=21)
                                                               cTGCCAAACGGTGGGAGCGGGCGTGCCCACGAAACAAGCGCCGCGGTGCATGGCACCGGCTTTTCCGGTGTTCTCGGGCACACCCGCTCCCACAGACTGATTGCCACTCAACGGGCCTCTCGCTCCAGTGACAGccc                                                                          >  2:13153187/1‑137 (MQ=17)
                                                                 gCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGATTTGTCGGTGTTCGCGGGAACACCCACTCCCAAAAACTGAGTGCCATTAAACGGGCCTCTCGCTACAGTGGCGGcccat                                                                        >  2:13582731/1‑135 (MQ=12)
                                                                 gCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCACCGCGGTGCATGGCACCGGATTTGCCGGTGTTCGCGGGCACGCCCGCTCCAAAAGAATGATTGCCATTAAGCGGGCCTATCGTTACAGTGGAGGCCCGt                                                                        >  2:18591763/1‑137 (MQ=17)
                                                                    aaaCGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCAGTGCATGGCCCCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGACTGCCATTAAGAAGGACTATTGCTAAAGTGACGGTCCGtttt                                                                     >  2:9729314/1‑137 (MQ=17)
                                                                     aaCGGTGGGAGCGGGCGTGTCCACGAAACAGGCGCCGCGGTGCATGGCACCGGCCTTGCCGGTGTTCGCAGGCACACCCACTCCCACAAACCGATTGCCACTCAACAGGCCCATCACTCCAGTGACAGCCCAttttt                                                                    >  2:6589226/1‑137 (MQ=9)
                                                                            ggAGCGGGCGGGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTTCCGGTGTTCGCGGGCACTCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGcc                                                             <  2:9485577/137‑1 (MQ=21)
                                                                                   gtgTGCCTTCGAAACGTGCGCCGTGGTGCGTGGCACTGTCTTTGCTGGTGTTCGTGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTaa                                                      <  1:2380521/135‑1 (MQ=19)
                                                                                        cccGCGAAACAGGCGCCGTGGCGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAAt                                                 <  2:12603422/137‑1 (MQ=25)
                                                                                          cgtgAAACAGGCGCCATGGTGCATGGCACCGGTTTTTCTGGTGTTCGTGGGCATTCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                               <  2:8974696/134‑1 (MQ=255)
                                                                                          cgcgAAACAGGCGCCGCGGTGCATGCCACCGGTTTTGCCGTTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCTTCGCTAACGAATCa                                               <  2:16155409/137‑1 (MQ=25)
                                                                                          cgcgAAACAGCCTCCGCGGTGCATGCCACCGGTTTTGCTGTTGTTCGCGGCCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                               <  1:4241438/137‑1 (MQ=255)
                                                                                          cgcgAAAAAGGCGCCGGGTGGCATGGCACCGGTTTTGCCGTTGTTCGCGGTCATGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGTCCGTTTTTTGTTGCCTCGCTAACGAATCa                                               <  2:16563317/137‑1 (MQ=255)
                                                                                          cgcaAAACGGGCGCTTCGGTGCATGGCACCGGTTTTGTTGGGGTTCGTGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                               <  2:18421284/137‑1 (MQ=25)
                                                                                          cgagAAAAAGGCGTCGTGGTGCAGGGCACTGGTTTTGCTGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                               <  1:16714405/134‑1 (MQ=25)
                                                                                          agagAAACAGTCGCCGCGGTGCATGGCACTGGCTTTGCAGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                               <  2:7872206/134‑1 (MQ=25)
                                                                                           gcgAAACGGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCAt                                              <  2:13636591/137‑1 (MQ=25)
                                                                                           gcgAAACAGGCGCCGCGGGGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCAt                                              <  2:19466132/137‑1 (MQ=37)
                                                                                           gagAAACAGGCGCCGCGGGGCATGGCACCGGCTTTGCCGGTGTTCGAGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCAt                                              <  2:9210811/135‑1 (MQ=37)
                                                                                               aaCAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                               >  1:540212/1‑132 (MQ=37)
                                                                                               aaCAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                               <  2:540212/132‑1 (MQ=25)
                                                                                                acggtcgccgcGGTGCATTGCACCGGTTTTGCTGGGGTTCGTGGGCACCCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTc                                         <  1:19440212/132‑1 (MQ=255)
                                                                                                aCAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACta                                         >  2:10908353/1‑136 (MQ=37)
                                                                                                aCAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTc                                         >  2:2124070/1‑137 (MQ=37)
                                                                                                aCAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTc                                         >  1:7424000/1‑137 (MQ=37)
                                                                                                     cgccgtGGGGCATGGCACTGGTTTTGCTGGTGTTCGGGGACACCCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa                                    <  1:1905471/137‑1 (MQ=255)
                                                                                                     cgccgcGGTGCATGGCACCGGCTTTGCTGGTTTCCACGGGCACACCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa                                    <  2:13934561/137‑1 (MQ=37)
                                                                                                       ctgcgGGGCAGGGCACCGGTTTTGCCGGTGTTCGAGGGCACTCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGa                                  <  1:16125544/135‑1 (MQ=255)
                                                                                                       ccgcGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGa                                  >  1:10129997/1‑137 (MQ=37)
                                                                                                                  ggcaccGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAg                       >  1:3699995/6‑137 (MQ=37)
                                                                                                                     accGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAg                    >  2:5839534/3‑137 (MQ=38)
                                                                                                                     accGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAg                    >  2:5691744/3‑137 (MQ=38)
                                                                                                                       cGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGAt                  >  2:11730762/1‑137 (MQ=38)
                                                                                                                       cGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGAt                  >  1:7617038/1‑137 (MQ=38)
                                                                                                                         ggCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATc                 >  1:12129885/1‑137 (MQ=38)
                                                                                                                           cTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAATGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg               >  2:8436614/1‑137 (MQ=38)
                                                                                                                           cTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCtt                                     >  2:18932193/1‑115 (MQ=38)
                                                                                                                           cTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGCACGCTATCAAGCAGATCTg               >  1:1231647/1‑137 (MQ=38)
                                                                                                                           cTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg               >  1:14072479/1‑137 (MQ=39)
                                                                                                                           cTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg               >  1:14747819/1‑137 (MQ=39)
                                                                                                                           cTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg               >  2:1515140/1‑137 (MQ=39)
                                                                                                                           cTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg               >  1:13438155/1‑137 (MQ=39)
                                                                                                                           cTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg               >  2:15959612/1‑137 (MQ=39)
                                                                                                                           cTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg               >  2:14936048/1‑137 (MQ=39)
                                                                                                                           cTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg               >  1:8721701/1‑137 (MQ=39)
                                                                                                                           cTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg               >  1:2849596/1‑137 (MQ=39)
                                                                                                                           cTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg               >  1:15344151/1‑137 (MQ=39)
                                                                                                                           cTTTGCAGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGTCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCACTTAGAACGCTATCAAGCAGATCTg               >  2:17095145/1‑137 (MQ=38)
                                                                                                                            tttGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGa              >  1:5238019/1‑137 (MQ=39)
                                                                                                                             ttGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGaa             >  2:10649380/1‑137 (MQ=39)
                                                                                                                             ttGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGaa             >  1:18545355/1‑137 (MQ=39)
                                                                                                                             ttGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGaa             >  2:11982096/1‑137 (MQ=39)
                                                                                                                             ttGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGaa             >  2:12869465/1‑137 (MQ=39)
                                                                                                                             ttGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACAAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGaa             >  2:16670651/1‑137 (MQ=39)
                                                                                                                                       tgttcgcGGGCACTCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTccc      <  2:11222890/131‑1 (MQ=255)
                                                                                                                                          tcgcGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGAc   >  1:10605428/4‑137 (MQ=255)
                                                                                                                                           cgcGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGAct  >  2:15337319/3‑137 (MQ=255)
                                                                                                                                           |||                                                                                                                                        
TGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTGTTTCGGGCACGCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACT  >  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 < 24 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: