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.C91.0% 783.9 / 84.8 258intergenic (+83/‑81)trpS CDS (Tryptophan‑‑tRNA ligase)/zapE CDS (putative enzyme with nucleoside triphosphate hydrolase domain)–/–
Reads supporting (aligned to +/- strand):  ref base . (17/6);  new base C (167/67);  total (185/73)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 270 total aligned reads displayed.

GATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTGTTTCGGGCACGCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTC  >  AE015451_phaM1/1499373‑1499637
                                                                                                                                              |                                                                                                                                    
gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGTGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCAcgcc                                                                                                                                    >  2:8526317/1‑137 (MQ=25)
gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAAGCGCCACGGTGCACGGCACCGGCTTTTCCGGTGTTTTCGGGCAcgcc                                                                                                                                    >  2:2232358/1‑137 (MQ=37)
gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAAAAGGCGCCACGGTGCCCGGCCCCCGCTTTTCCGGTGTTCGCGGGCActcc                                                                                                                                    >  2:9334172/1‑137 (MQ=37)
gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCACGAAAACGGCGCCAAGGTTCACTGCCCCAGATTTGCCGGTGTTCTCAGGCAagcc                                                                                                                                    >  1:16820284/1‑137 (MQ=255)
  tGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCACGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCAcgcccg                                                                                                                                  >  2:16520937/1‑135 (MQ=25)
         aTTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCACGAAACAAGCGCCACGGTGCATGGCAACGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCcac                                                                                                                           >  2:2335353/1‑137 (MQ=21)
              gagCCTGAACGAATCCTTGGAACCACAGCGGGACTGAGGCTGCCAAATGGTGGGAGTGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGACACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACt                                                                                                                      <  2:2156795/137‑1 (MQ=18)
                  cTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTTCCCGCGAAACAAGCGCCGCGGTGTATGGCACCGGATTTTCCGGTGTTCTCGGGCACGCCCGCTCCCACAGAATGAtt                                                                                                                  >  1:6997782/1‑137 (MQ=21)
                    gAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTTCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGc                                                                                                                >  2:9123087/1‑137 (MQ=21)
                    gAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCACGAAACAAGCGCCACGGTGCACGGCACCAGCTTTTCCGGTTTTCGCGGGCACGCCCACCCCCACAAACTGATTGc                                                                                                                >  1:8043030/1‑137 (MQ=255)
                       cGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGGGCACGCCCGCTCCCACAGACTGATTGCCATTa                                                                                                           >  1:16329588/1‑137 (MQ=25)
                       cGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCACCGCGGTGCATAGCACCGGCTTTGCCGGTGTTCGCGGGCAAACCAGTTCCCAAAGAATGATTGCCAt                                                                                                             >  2:14535565/1‑137 (MQ=18)
                        gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGGGTGTCCGCGAAAAAAACGCCCCGGTGCATGGCACCAACTTTTCCGGTGTTTTTGGGCACGCCCACCCCCCCAGAATGATTGTCAtt                                                                                                            >  1:11041004/1‑137 (MQ=255)
                        gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGTGCCACGGTGCCTGGCACCAGCTTTGCCGGTGTTCTCGGGCACACCCGCTTCCACAGAATGATTGACAtt                                                                                                            >  1:13053385/1‑137 (MQ=21)
                        gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCAAAAAAAGAGCCACGGGGCATGGCACCGGCTTTTCCAGTGTTCGCGGGAACACCCGCTCCCCCAAACTTATTGCCact                                                                                                            >  2:15674160/1‑135 (MQ=12)
                        gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGATTGCCCTCCAAACCAGCGCCACGGTGTCTGGCACCGGCTCTGTCCGTGTTCTCGGGAAAACCCACTCCCCCAAACCGATTGCCact                                                                                                            >  1:4714848/1‑135 (MQ=255)
                        gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCAAAAAGGTGGGAGCGGGCGCGCCCACGAAACAGGCACCGCGGTGCATGGCACCGGATTTGCCGGTGTTCTCGGGCACGCCCGCTCCCACGGACGGATTACCAtt                                                                                                            >  2:13146225/1‑137 (MQ=18)
                          aGCCTTGTCACCACAGCAGGACTGGGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACGGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTaa                                                                                                          <  2:16948129/137‑1 (MQ=18)
                            ccTTGTCACCACAACAGAACGGAGGCTTACAAACGGTGGAAGCGGGCGTGCCCGCGAAACAGGCACCACGGTGCATGGCACCGACTTTTCCGTTGTTCCCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGc                                                                                                        <  2:12043863/137‑1 (MQ=12)
                              ttGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGGGCCGAGGTGTATTGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGAATGAGTGACATTCAGCgg                                                                                                      >  1:4206868/1‑137 (MQ=17)
                              ttGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCAAAACAAGCACAGCAGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAGAATGATTGCCCTTCAGTgg                                                                                                      >  2:13255893/1‑137 (MQ=12)
                                 gCACCACAGCAGGACTGAGTTTGCCAAACGTTGGGAGCGGGCGTGCCCGCGAAAAGGGCGCCGTGGTGCATGGCACCGGTTTTGCCGTTTTTCGCGGGCACTCCCGCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                                                   <  2:11945443/137‑1 (MQ=18)
                                 ccaccacAGAGGACCTGCGTCTGCAAAACGGTGGGAGCGGGTGTGCCCGTAAAAAAGGTCCCATGGTCCCTGACCCTGGCTTTGCCGGTGTTCGTGGGCACACCCGCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                                                   <  2:5791/136‑1 (MQ=255)
                                  caccacAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGATCCCACAGACTGATTGCCATTAAGCGGGACt                                                                                                  >  2:10402044/1‑137 (MQ=18)
                                  caccacAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGAGTGTCCACGAAAAAAGAACCACGGTGTCTGGCACCCGATTTTCCGGTGTTTTCGGGGACACCCCCCCCCCCAAAATGATTGTCACTCAGCGGGGCt                                                                                                  >  1:14649028/1‑137 (MQ=255)
                                         gcagGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCCTGGCACCGGCCTTGCCGGTGTTCGCGGGCACGCCCGCTTCCACAGAATGATTGCCATTCAGCGGGCCTATCGCTa                                                                                           >  2:15912533/1‑137 (MQ=18)
                                          cagGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGGAACAGGCGCCGCGGTGCATGGCACCAGATTTGCCAGTGTTCGCGGGCACGCCCGCTCCCCCAGACTGATTGCCATTAAGCGGGTCTATCGCtcc                                                                                          >  2:6118829/1‑135 (MQ=18)
                                                cTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCCCAGAATGATTGCCATTAAACCGGCCCCTCGCTGCAGTgg                                                                                     >  2:13436845/1‑136 (MQ=14)
                                                 tGAGGCTGCCAAACGGTGGGAGCGGGGGTGTCCGCGAAAAAAGCGCCGCGGTGCACGGAACCGGCTTTGTCGGTGTTCGCGGGCACGCCCGCCCCCAAAGAATGACTGCCATTCAGCGGGTCTATCGCTCCAATGAc                                                                                    >  1:12049481/1‑137 (MQ=12)
                                                   gggCTGCCAAATGGGGGGAGTGGGTTTGCCTGGGAAAAGGGTGCCGCGGTGCGGGGCACTGGTTTTGCCGGTGTTTGTGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACgg                                                                                  <  1:10344492/136‑1 (MQ=255)
                                                   aGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCCTGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCACTCCCACAGAACGGGTGCCATTAAGCGGGCCTCTCATTACAATGACgg                                                                                  >  2:4534108/1‑137 (MQ=14)
                                                      cTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGAATGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGccc                                                                               >  2:1747795/1‑137 (MQ=18)
                                                      cTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCAGCTTTGCCGGTGTTCGCGGGCCCGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGccc                                                                               >  2:8418754/1‑137 (MQ=21)
                                                           aaaCGGTGGGAGCGGGCGTGCCCGCGAAACAGGTGCCGCGGTGCATGGCCCCGGATTTGCCGGTGTTCGCGGGCACGCCCGCTCCCCCAAACTGATTGCCATTCAGCGGGCCTATCGCTACATTGACGGCCCTtttt                                                                          >  2:11132556/1‑137 (MQ=18)
                                                           aaaCGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGAGGTGCAGGGCACCGGCTTTGCCGGTGTTCGAGGGCAAGCCCGCCCCCACAGAATGATTGCCATTAAGCGGGACTCTCGCTACAATGACGGCCCGtttt                                                                          >  2:16489255/1‑137 (MQ=17)
                                                           aaaCGGTGGGAGCGGGCGTGCCCGCGAAACAAGCGCCACGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAAACTGATTGCCATTAAGCGGGACTATCGCTACAGTGACGGCCCGtttt                                                                          >  2:6236668/1‑137 (MQ=17)
                                                             aCGTTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCTGGGTATGGCACCGGTTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGtttttt                                                                        <  2:2605183/137‑1 (MQ=21)
                                                               ggggggAGTGGGTGTGCCCGAGAAAAGGGCGCTGGGGGGCGTGCCACTGGTTTTGCTGGTGTTCGGGGGCACCCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTttgt                                                                      <  1:6448523/134‑1 (MQ=255)
                                                               ggTGGGAGAGGGCGTGCCTGCGAAACGGGCGCTGCGGTGAATGGCACCGGCTTTGCGGGTGTTCGCGGGAACCCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTttgt                                                                      <  2:631284/137‑1 (MQ=255)
                                                                   ggAGCGGGCGTGCCCGCGAAACAGGCGCCGTGGTGCATGGCACCGTCTTTGCCGGTGTTAGCGGACACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGcc                                                                  <  2:6774309/137‑1 (MQ=21)
                                                                      gCGGGCGTGCCCGCGAAACAAGCGCCACGGTGCACGGCACCAGCTTTTCCAGTGTTCGCGGGCACACCCACTCCCCAAAACTGACTGCCATTAAGCGGGCCTCTCGCTCCAATGACGGCCCGTTTTTTTTTTCCTCg                                                               >  2:5656582/1‑137 (MQ=12)
                                                                           cGTGCCCGCGAAACAGGCGCCACGGTGCATTGCACCGGCTTTTCCCATGTTCTCGGGCACACCCACTCCCACAGACTGATTGTCATTCAGCGGGCCTCCTGCTACAATGACAGCCCGTTTTTTGTTTCCCCACTAAc                                                          >  2:12646643/1‑137 (MQ=255)
                                                                            gTGCCCGAGAAACGGGCGCCGTGGTGCAGGGCACTGGTTTTGCCGGTGTTTGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACg                                                         <  1:12740911/137‑1 (MQ=37)
                                                                               cccGCGAAACAGGCGCCGCGGTGCATGGCACCGTCTTTGCAGGTTTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAAt                                                      <  2:17178478/137‑1 (MQ=25)
                                                                                 cgcgAAACGTGCGCCGCGGGGCATGGCACTGGTTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                    <  2:1851752/137‑1 (MQ=25)
                                                                                 cgcgAAACGGGCGCCGTGGTGCGTGGCACCGTTTTTGCTGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                    <  2:484541/137‑1 (MQ=25)
                                                                                 cgcgAAACGGGCGCCGCGGTGAGTGGCACCGGCTTTGCCGGTGTTCGGGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                    <  2:2627045/137‑1 (MQ=25)
                                                                                 cgcgAAACAGGCGCCGCGTTGCATGGCACCGGTTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                    <  2:2016868/137‑1 (MQ=37)
                                                                                 cgcgAAACAGGCGCCGCGGTGAATGGCACTGGTTTTGCCGGTGTTCGCGGACACCCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                    <  1:3080486/137‑1 (MQ=255)
                                                                                 cgcgAAAAAGTCGCTGTGGTGAGTGGCACTGGTTTTGCGGGTGTTTGCGGACACCCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                    <  1:558772/137‑1 (MQ=255)
                                                                                 cacgAAAAAGGCGCCGTGGGGCGTGGCACTGTCTTTTCCGGTGTTTGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                    <  1:6906198/135‑1 (MQ=25)
                                                                                            cgccgtGGTGCATGGCACCGGTTTTGTTGGTGTTCGCGGACACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa                                         <  1:9522607/137‑1 (MQ=255)
                                                                                            cgccggGGGGCATGGCACTGGTTTTGCTGGTGTTCGCGGGCACCCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa                                         <  1:14915996/137‑1 (MQ=255)
                                                                                            caccgcGGCGCATGGCACCGGTTTTGCTGGTGTTTGGGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa                                         <  2:10198965/135‑1 (MQ=37)
                                                                                                      cATGGCACCGGCTTTGCCGTTGTTCTCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATc                               <  1:17280017/137‑1 (MQ=38)
                                                                                                            accGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAg                         >  2:8961495/3‑137 (MQ=38)
                                                                                                            accGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAg                         >  1:2375555/3‑137 (MQ=38)
                                                                                                              cGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGAt                       >  1:16772680/1‑137 (MQ=38)
                                                                                                              cGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCAACAGACTGATTGCCCTTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCATCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGTAGAt                       >  2:2758059/1‑137 (MQ=38)
                                                                                                                  cTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACCGACTGCTTGCCCTTAAGCGGGCCTCTCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGAATCCCTTAGAACGCTCTCAAGCAGATCTg                    >  2:4245061/1‑137 (MQ=37)
                                                                                                                  cTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTATAACGCTATCAAGCAGATCTg                    >  2:4080456/1‑137 (MQ=38)
                                                                                                                  cTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                    >  2:3008480/1‑137 (MQ=39)
                                                                                                                  cTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                    >  2:5083728/1‑137 (MQ=39)
                                                                                                                  cTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                    >  2:6969436/1‑137 (MQ=39)
                                                                                                                  cTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                    >  2:7205293/1‑137 (MQ=39)
                                                                                                                  cTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                    >  2:9808887/1‑137 (MQ=39)
                                                                                                                  cTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                    >  1:2841450/1‑137 (MQ=39)
                                                                                                                  cTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                    >  1:10093460/1‑137 (MQ=39)
                                                                                                                  cTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                    >  1:10669867/1‑137 (MQ=39)
                                                                                                                  cTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                    >  1:3519442/1‑137 (MQ=39)
                                                                                                                  cTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                    >  1:4445103/1‑137 (MQ=39)
                                                                                                                  cTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                    >  1:15601554/1‑137 (MQ=39)
                                                                                                                  cTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                    >  1:13705614/1‑137 (MQ=39)
                                                                                                                  cTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                    >  1:14392878/1‑137 (MQ=39)
                                                                                                                  cTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                    >  1:6843968/1‑137 (MQ=39)
                                                                                                                  cTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                    >  1:9759177/1‑137 (MQ=39)
                                                                                                                  cTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                    >  1:15093673/1‑137 (MQ=39)
                                                                                                                  cTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                    >  2:12197098/1‑137 (MQ=39)
                                                                                                                  cTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCACGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                    >  2:17556756/1‑137 (MQ=38)
                                                                                                                  cTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGACCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCACTTAGAACGCTATCAAGCAGATCTg                    >  2:8762547/1‑137 (MQ=38)
                                                                                                                  cTTTGCCGGTGTTCGTGGGCACGCCAGCTCCCACAGACTGATTGCCATTAAGCGGGCATATCGCTACAGTGACGGCCCGTTTTTTATTGCCTCGCTAACGAATCATGACTCACTTAGAACGCTATCGAGCAGATCTg                    >  2:10855645/1‑137 (MQ=255)
                                                                                                                  cTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg                    >  2:14100283/1‑137 (MQ=39)
                                                                                                                    ttGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATc                               <  1:9042898/124‑1 (MQ=38)
                                                                                                                    ttGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGaa                  >  2:7768816/1‑137 (MQ=39)
                                                                                                                    ttGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGaa                  >  2:7814361/1‑137 (MQ=39)
                                                                                                                    ttGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGaa                  >  1:13157796/1‑137 (MQ=39)
                                                                                                                    ttGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGaa                  >  1:5559036/1‑137 (MQ=39)
                                                                                                                    ttGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGaa                  >  2:1256290/1‑137 (MQ=39)
                                                                                                                    ttGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGaa                  >  2:12859773/1‑137 (MQ=39)
                                                                                                                    ttGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGaa                  >  2:14633632/1‑137 (MQ=39)
                                                                                                                    ttGCCGGTGTTCGCGGGCACGCCCGCTCCCACAAACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGGCTCCCTTAGAACGCTATAAAGCAGGTCTGaa                  >  2:7114176/1‑137 (MQ=38)
                                                                                                                                tcgcGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                    >  2:15102024/4‑93 (MQ=255)
                                                                                                                                tcgcGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGAc        >  2:3335374/4‑137 (MQ=255)
                                                                                                                                tcgcGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGACCTGAAACGTCCCGAc        >  2:1789650/4‑137 (MQ=255)
                                                                                                                                 cgcGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGAct       >  2:3444218/3‑137 (MQ=255)
                                                                                                                                 cgcGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGAct       >  2:1613518/3‑137 (MQ=255)
                                                                                                                                 cgcGGGCACGCCCGCTCCCACAGACTGACTGCCATTAAGCGGGCCTATCGCTACCGTGACGGCCCGTTTTTTGTTGCCTCGCTACCGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCAGAct       >  2:2931630/3‑137 (MQ=255)
                                                                                                                                    gggCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGa                                                        <  1:10015718/85‑1 (MQ=255)
                                                                                                                                    gggCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGa                                       >  2:3752797/1‑102 (MQ=255)
                                                                                                                                      gCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACttcttc  >  2:8153757/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: