Predicted mutation
evidence seq id position mutation annotation gene description
RA AE015451_phaM1 1,499,508 +C intergenic (+85/‑79) 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,5061.C92.8% 903.6 / 77.2 308intergenic (+83/‑81)trpS CDS (Tryptophan‑‑tRNA ligase)/zapE CDS (putative enzyme with nucleoside triphosphate hydrolase domain)–/–
Reads supporting (aligned to +/- strand):  ref base . (18/4);  new base C (171/111);  total (193/115)
Fisher's exact test for biased strand distribution p-value = 6.62e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 325 total aligned reads displayed.

GATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTGTTTCGGGCACGCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTC  >  AE015451_phaM1/1499373‑1499637
                                                                                                                                            |                                                                                                                                   
gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGGGTGTCCGCGAAACAAGAGCCGCGGTGCATGGCACCGGTTTTTCCGGTGTTTACGGGCAcgcc                                                                                                                                   >  1:6044929/1‑137 (MQ=255)
gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTTCCCACTAAACAAGCACCGCATTGCATAGCACCGACCTTGCCGGTGTTCGCGGGCAcgcc                                                                                                                                   >  1:15383535/1‑137 (MQ=25)
gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCAcacc                                                                                                                                   >  2:19005410/1‑137 (MQ=25)
gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTACCGGTGTTCGCGGGCAcgcc                                                                                                                                   >  2:19349402/1‑137 (MQ=25)
gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCACGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCAcacc                                                                                                                                   >  2:13168356/1‑137 (MQ=25)
gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCACGAAACAGGCGCCACGGTGCACGGCCCCGGCTTTTCCAGTGTTCTCGGGTAcgcc                                                                                                                                   >  2:2472108/1‑137 (MQ=38)
gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGAGTGGCCGCGAAAAAAGAGCCACAGTGCCCGGCACCGAATTTTCCAGTGTTTTCGGGCAcacc                                                                                                                                   >  1:7548722/1‑137 (MQ=255)
 aTGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAAGCGCCGCGGTGCATGGCCCCGGCTTTGCCGGTGTTCGCGGGCAcgccc                                                                                                                                  >  2:3069346/1‑136 (MQ=25)
       ccATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCAAAAAAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTccc                                                                                                                             >  2:13985222/1‑137 (MQ=21)
         aTTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGTCTGGCACCGGGTTTGCCGGTGTTCGCGGGCACGCCCACTCCcac                                                                                                                           >  2:1792298/1‑137 (MQ=21)
         aTTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGAGGTGCATGGCACCGGCTTTGCCAGTGTTCGCGGGCACGCCCGCTCCcac                                                                                                                           >  2:5956804/1‑137 (MQ=21)
         aTTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCAAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGTCCGCTCCcac                                                                                                                           >  2:6933881/1‑137 (MQ=21)
             agagCCGAAACGAAGCTTTGGCACCACAGCAGGACGAAGGCTGCCAAACGGGGGGAGCGGTTGTACCCCCGAAAAAGGCGCCGCGGGGCATGACACCGTCTTTGCTGGGGTTCGAGGACACACCCGCTCCCACAGAc                                                                                                                       <  2:18096541/137‑1 (MQ=255)
                   tGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGACTGATTg                                                                                                                 >  2:16317676/1‑137 (MQ=18)
                    gAACGACGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCACGAAACAGGCGCCACGGTGCACGGCACCGGATTTGCCGGTGTTCGCGGGAACGCCCGCTCCCACAGAATGATTGc                                                                                                                >  2:12003050/1‑137 (MQ=21)
                    gAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCTCTCCCACAGACTGATTGc                                                                                                                >  2:12420983/1‑137 (MQ=21)
                       cGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCCCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGccct                                                                                                             >  2:15877932/1‑135 (MQ=21)
                       cGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGGCTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCac                                                                                                             >  2:1661530/1‑136 (MQ=21)
                       cGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAAACTGATTGCCAt                                                                                                             >  2:13814346/1‑137 (MQ=18)
                       cGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAAAAAGCGCCACGGTGCCTGGCACCGGCTTTGTCGGTGTTCTCGGGCCCGCCCGCTCCCACAGACTGATTGCCAt                                                                                                             >  2:10207892/1‑137 (MQ=21)
                        gAAGCTTTGGCACCACAGCAGGATTGGGGCGGCAAAACGGTGGGAACGGGCGTTCCCGCGAAGCAGTCGCCGTGGTGCGTGGCACCGGCTTTGCTGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCAtt                                                                                                            <  2:9231633/137‑1 (MQ=12)
                        gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGGGTGTCCGCGAAACAAGAGCCGCAGTGCCTGGCACCGGCTTTTCCGGTGTTTTTGGGCACGCCCGCTCCCCCAAAACGAGTGCCact                                                                                                            >  1:13031249/1‑135 (MQ=255)
                        gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGGGTGTCCGCGAAAAAAGCGCCGTGGTGCCCGGCACCTGCTTTGCCGGTGTTTGCGGGCACGCCCGCTCCCACCAACCGAGTGCCAtt                                                                                                            >  1:17826374/1‑137 (MQ=19)
                        gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGTCCGCGAAAAAGGAGCCGCGGTGCATGGTACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGAATGGTTGCCact                                                                                                            >  1:8670328/1‑135 (MQ=21)
                        gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCAtt                                                                                                            >  2:9291196/1‑137 (MQ=21)
                        gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCACGGTGCATGGCCCCGGCTTTGCCAGTGTTCGCAGGCACGCCCTCCCCCACAGACCGATTGCCact                                                                                                            >  2:5698429/1‑135 (MQ=21)
                        gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAAGCGCCACGGTGCACGGCACCGGCTTTTCCGGTGTTCGCGGGCACGCCCACTCCCACAAACTGATTGCCCtt                                                                                                            >  2:16740160/1‑137 (MQ=18)
                        gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAAAGGTGGGAGCGGGCGTGTCCGCGAAACAGACGCTGAGATGTCTGTCACCGGCTTTGCCGGTGTTCACGGGTACGCCCGCTCCCAAAGACCGATTGCCAtt                                                                                                            >  2:10588537/1‑137 (MQ=19)
                              ttGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGTCGCGGTGCATGGCACCGGCTTTGCCAGTGTTCGCGGGCACGCCCGCTCCCACAGAATGATTGCCACTAAGCgg                                                                                                      >  2:15292133/1‑137 (MQ=18)
                              ttGGCACCACAGAAGGATTGGGGTGGCAAAATGGTGGGAGTGGGTGTTCCCGAGAAACAGGCGCTGGGGGGAGGGGCACGGGTTTTGCTGGTGTTTGGGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCgg                                                                                                      <  1:1105792/137‑1 (MQ=255)
                                 gCACCACAGCAGGACGGAGGCTGCAAAACGGGGGGAGCGGGTGTGCCCGCGAAACGGGCGCCGTGTTGCGTGGCACCGGTTTTTCCGGTGTTTGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                                                   <  2:8721340/137‑1 (MQ=21)
                                 gCACCACAGAAGGATGGAGGCTGCAAAATGGGGGGAGCGGGCGTGCCCGAGAAAAAGGCACAGTGGGGCGGGAAACCGGTTTTGCCGGTGTTCGGGGACACCCCCGCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                                                   <  1:8793533/137‑1 (MQ=255)
                                 gCACCAAAGCAGGACTGAGGCTGCCAAACGGTGGGGACGGGCGTGCCCGCGAAACATGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTATTCGTGGACACGCCCGCTCCCACAGACTGATTACCATTAAGCGAGcc                                                                                                   <  2:18855286/137‑1 (MQ=17)
                                  caccacAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGGGTGTCCGCGAAACAGGAGACACGGTGCATGGCACCAGCTTTGCCGGTGTTTGCGGGTACACCCTCTCCCCCAGACTGACTGCCCTTCCGCGGGCCt                                                                                                  >  1:7326042/1‑137 (MQ=255)
                                  caccacAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGTATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGAACGATTGCCATTAAGCGGGCCt                                                                                                  >  2:16908586/1‑137 (MQ=18)
                                  aaccacAGAAGGACGGAGGCTACAAAATGGGGGGAGCGGGTGTGCCTGCGAAACGGGCGCCGCGGTGCATGGAACCGGTTTTGCTGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCt                                                                                                  <  2:7794666/136‑1 (MQ=16)
                                      acaGCAGGACTGAGGCTACTAAATGGTGGGAGCGGGCGTTCCCTCGAAACAGGCCCCGCGGTACGTGGCACCGGCTTTGCCGGTGTTCGTGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCg                                                                                              <  2:4428047/137‑1 (MQ=18)
                                      acaGCAGGACTGAGGCTACCAAACGGTGGGAGTGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGTTTGGCCGGGTTTAGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCg                                                                                              <  2:17594403/137‑1 (MQ=18)
                                      acaGCAGGACTAAGGCTGCCAAACGGTGGGAATGGGCGTTCCCGCAAAACATGCGCCGCGGTACATGACACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCg                                                                                              <  2:7616976/137‑1 (MQ=21)
                                      acaGCAGAACGGAGGTGGCAAAAGGGGGGGGGTGGGTGTGCCCGAAAAACAGGCACCGGGGGGCATGGCACCGGTTTTGCTGGTGTCCGAGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCg                                                                                              <  1:18358297/137‑1 (MQ=12)
                                        agcagGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCACGGTGCATGGCACCGGCCTTGCCGGTGTTCGCGGGCACGCCCACTCCCACAGACTGAGTGCCATTCAAAGGGACTCTCGCt                                                                                            >  2:3263343/1‑137 (MQ=17)
                                             gACTGAGGCTGCCAAACGGTGGGAGCGGGAGTGCCCACGAAAAAAGCGCCGCGGTGCATGGCCCCGGCCTTTCCGGTGTTCGCGGGCACACCCACCCCCACAAACCGATTTCCCTTCAACAGGCCTCTCGCTAccat                                                                                       >  1:9433310/1‑134 (MQ=255)
                                                ggAGGCTGCCAAACGGGGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGACACCGGCTTTGCCGTTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGAc                                                                                    <  2:14410429/136‑1 (MQ=18)
                                                     cTGCCAAACGGTGGGAGCGGGCGTTCCCGCGAAACAGGGGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACACCCGATCCCACAGGATGATTGTCATTAAGTGGGCCTATCGCTGAAGTGGCGGtcc                                                                               >  2:5326698/1‑137 (MQ=9)
                                                     cTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACCAGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTCAGCGGGCCTATCGCTCCAGTGACGGccc                                                                               >  2:1500370/1‑137 (MQ=17)
                                                     cTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCGTTGCCGGTGTTCACGGGCAAGCCCGCTCCCACAAACTGATTGCCATTCAGCGGGCCTATCGCTCAAGTGACGGccc                                                                               >  2:9024319/1‑137 (MQ=18)
                                                       gCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCCTGGCACCGGCTTTGCCGGTGTTCGCGGGCACGACCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTCCAGTGACGGCCCGt                                                                             >  2:18547803/1‑137 (MQ=21)
                                                          aaaCGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCAGTGCATGGCCCCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGACTGCCATTAAGAAGGACTATTGCTAAAGTGACGGTCCGtttt                                                                          >  2:9729314/1‑137 (MQ=17)
                                                           aaCGGTGGGAGCGGGCGTGTCCGCGGAACAGGCGCCGCGGTGCATGGCACCGGATTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGAGGGCCTATCGCTACAATGACGGCCCGctttt                                                                         >  2:11547479/1‑137 (MQ=18)
                                                           aaCGGTGGGAGCGGGCGTGTCCACGAAACAGGCGCCGCGGTGCATGGCACCGGCCTTGCCGGTGTTCGCAGGCACACCCACTCCCACAAACCGATTGCCACTCAACAGGCCCATCACTCCAGTGACAGCCCAttttt                                                                         >  2:6589226/1‑137 (MQ=9)
                                                               gTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTTCCATTAAGCGGGCCTCTCGCTACAATGACGGCCCGGTTTttgtt                                                                     >  2:16147671/1‑137 (MQ=18)
                                                                ggggAGCGGGCGTGCCCGCGAAACAGGCGCCGTGGGGAATGGAACCGTCCTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTttgttg                                                                    <  2:9127012/136‑1 (MQ=21)
                                                                  ggAGCGGGCGGGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTTCCGGTGTTCGCGGGCACTCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGcc                                                                  <  2:9485577/137‑1 (MQ=21)
                                                                     gCGGGTGTGTCCGCGAAATAGGCGCCGCGGTGCATGGAACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGAATGATTGCCATTAAGCGGGACTATCGCTACAGTGACGGCCCGTTTTTTGTTGCATCg                                                               >  2:19213887/1‑137 (MQ=21)
                                                                     gCGGGCGTGCCCACGAAAAAAGCGCCACGGTGCATGGCACCAGCTTTTCCGGTGTTCGCGGGCACACCCGCTCCCACAGAATGAGTGCCATTAAGCGGGCCTATCGCTACAATGACGGACCATTTTTTGTTGCCCCg                                                               >  2:16086950/1‑137 (MQ=255)
                                                                      tgtttttGCCCGTGAAACAGGCCCCGTGGGGCGCGACACCGGCTTTGCCGGTGTTCGAGGACACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTAGCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGc                                                              <  2:10692027/131‑1 (MQ=19)
                                                                      cGGGCTTGCCCGCGAAACAGGCGCCGCGGGGCATGACACCGGCTTGGCTGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGc                                                              <  1:14500790/137‑1 (MQ=21)
                                                                         gtgTGCCTTCGAAACGTGCGCCGTGGTGCGTGGCACTGTCTTTGCTGGTGTTCGTGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTaa                                                           <  1:2380521/135‑1 (MQ=19)
                                                                                cgtgAAACAGGCGCTGTGGTGCATGGCACCGTCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                    <  1:14066113/134‑1 (MQ=25)
                                                                                cgcgAAACAGGCGCTACGGTGCATGGCACCGGCTTTGCCGGTGTTCGTGGACACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACNGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                    <  1:18785082/137‑1 (MQ=25)
                                                                                cgcgAAACAGGCGCCGCGGTGCATGGCACTGGTTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                    <  2:14559781/137‑1 (MQ=25)
                                                                                cgcaAAACAGGCGCCACGGTGCATGGCACCGGCTTTGCCGGTGTTCGTGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                    <  1:13625716/137‑1 (MQ=25)
                                                                                cgagAAAAAGGCGTCGTGGTGCAGGGCACTGGTTTTGCTGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                    <  1:16714405/134‑1 (MQ=25)
                                                                                agagAAACAGTCGCCGCGGTGCATGGCACTGGCTTTGCAGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                    <  2:7872206/134‑1 (MQ=25)
                                                                                 tcgAAAAAGGCTCCGCGGTGCGGGGAACCGGTTTTTCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCAt                                                   <  2:4996916/136‑1 (MQ=25)
                                                                                 gggAAACGGTCGCCGTGGTGCATGGCACTGTCTTTGCTGTTTTTCGCGGACACTCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCAt                                                   <  1:13475424/135‑1 (MQ=255)
                                                                                 gcgAAACAGGCGCCGCGGTGCAGGGCACCGGCTTTTCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCAt                                                   <  2:11731394/137‑1 (MQ=37)
                                                                                 gcaaaaCAGGCCCCGCGGTACAGGGCACCGGCTTTCCCGGTTTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCAt                                                   <  2:5654137/134‑1 (MQ=25)
                                                                                 gagAAACGGGCGCCGTGGGGCGTGGCACCGGTTTGGCTGTTGTTTGGGGGCACACCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCAt                                                   <  1:5203727/135‑1 (MQ=255)
                                                                                   gAAACAGGCGCCGCGGGGCATGGCACCGGTTTTGCTGGTGTTCGTGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGa                                                 <  1:17650694/137‑1 (MQ=37)
                                                                                     aaCAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                    >  1:540212/1‑132 (MQ=37)
                                                                                     aaCAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACt                                               >  2:13618278/1‑137 (MQ=37)
                                                                                      acggtcgccgcGGTGCATTGCACCGGTTTTGCTGGGGTTCGTGGGCACCCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTc                                              <  1:19440212/132‑1 (MQ=255)
                                                                                           cgccgtGGTACGTGGCACCGGCTTTTCCGGTTTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa                                         <  2:13136217/137‑1 (MQ=25)
                                                                                           cgccgtGGGGCATGGCACTGGTTTTGCTGGTGTTCGGGGACACCCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa                                         <  1:1905471/137‑1 (MQ=255)
                                                                                             ctgtggTGCGTGGCCCTGGTTTTTCCGGTGTTTGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGa                                       <  1:3389339/133‑1 (MQ=255)
                                                                                             ctgcgGGGCAGGGCACCGGTTTTGCCGGTGTTCGAGGGCACTCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGa                                       <  1:16125544/135‑1 (MQ=255)
                                                                                             ccgcGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGa                                       >  1:10129997/1‑137 (MQ=37)
                                                                                                     cATGGCACCGGCTTTGCCGGTGTTCGTGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATc                               <  2:12333394/137‑1 (MQ=37)
                                                                                                     cATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATc                               <  2:12693782/137‑1 (MQ=37)
                                                                                                     cATGGCACCGGCTTTGCCGGTGTTCGCGGGCACTCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATc                               <  2:11666021/137‑1 (MQ=37)
                                                                                                        ggcaccGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAg                            >  1:3699995/6‑137 (MQ=37)
                                                                                                        ggcaccGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAg                            >  2:11200355/6‑137 (MQ=37)
                                                                                                           accGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAg                         >  2:4722300/3‑137 (MQ=38)
                                                                                                           accGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAg                         >  2:17053346/3‑137 (MQ=38)
                                                                                                             cGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGAt                       >  1:7617038/1‑137 (MQ=38)
                                                                                                               ggCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATc                      >  1:12129885/1‑137 (MQ=38)
                                                                                                                 cTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAATGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                    >  2:8436614/1‑137 (MQ=38)
                                                                                                                 cTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCGTGACTCCCTTAGGACGCTATCAAGCAGATCTg                    >  2:14887565/1‑137 (MQ=38)
                                                                                                                 cTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCCAGCAGATCTg                    >  1:14072510/1‑137 (MQ=38)
                                                                                                                 cTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                    >  1:1463526/1‑137 (MQ=39)
                                                                                                                 cTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                    >  1:14975006/1‑137 (MQ=39)
                                                                                                                 cTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                    >  1:12739307/1‑137 (MQ=39)
                                                                                                                   ttGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGaa                  >  2:1342042/1‑137 (MQ=39)
                                                                                                                   ttGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGaa                  >  1:6894022/1‑137 (MQ=39)
                                                                                                                   ttGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGaa                  >  1:18545355/1‑137 (MQ=39)
                                                                                                                               tcgcGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGAc        >  2:1711713/4‑137 (MQ=255)
                                                                                                                               tcgcGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGAc        >  1:10605428/4‑137 (MQ=255)
                                                                                                                                cgcGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGAct       >  2:2107025/3‑137 (MQ=255)
                                                                                                                                   gggCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGActtct    >  1:3929480/1‑137 (MQ=255)
                                                                                                                                            |                                                                                                                                   
GATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTGTTTCGGGCACGCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTC  >  AE015451_phaM1/1499373‑1499637

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: