Predicted mutation
evidence seq id position mutation annotation gene description
RA AE015451_phaM1 278,744 +A intergenic (‑239/+75) PP_0224 CDS (Monooxygenase_ DszC family) ← / ← sctC CDS (sulfur compound ABC transporter ‑ ATP‑binding subunit) –/–

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451_phaM1278,7421.A94.8% 1846.4 / 85.7 550intergenic (‑237/+77)PP_0224 CDS (Monooxygenase_ DszC family)/sctC CDS (sulfur compound ABC transporter ‑ ATP‑binding subunit)–/–
Reads supporting (aligned to +/- strand):  ref base . (6/22);  new base A (251/259);  total (269/281)
Fisher's exact test for biased strand distribution p-value = 5.62e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 662 total aligned reads displayed.

GCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTC  >  AE015451_phaM1/278606‑278872
                                                                                                                                           |                                                                                                                                   
  tGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTaa                                                                                                                                   <  1:8398831/137‑1 (MQ=255)
  tGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTaa                                                                                                                                   <  1:9766472/137‑1 (MQ=255)
  tGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTaa                                                                                                                                   <  1:10954085/137‑1 (MQ=255)
  tGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTaa                                                                                                                                   <  1:5825708/137‑1 (MQ=255)
  tGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTaa                                                                                                                                   <  2:7778655/137‑1 (MQ=255)
  tGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTaa                                                                                                                                   <  1:4793847/137‑1 (MQ=255)
  tGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTaa                                                                                                                                   <  2:13974292/137‑1 (MQ=255)
  tGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTaa                                                                                                                                   <  2:14421409/137‑1 (MQ=255)
  tGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTaa                                                                                                                                   <  2:223413/137‑1 (MQ=255)
  tGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTaa                                                                                                                                   <  2:3318460/137‑1 (MQ=255)
   gATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTaaa                                                                                                                                  <  2:3625360/137‑2 (MQ=255)
    aTCGCGACAGACGTGTGCTGGTCCGGTGGTTCGGTACATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTACCGGGCTCTTCGCGGGTaaac                                                                                                                                 <  2:11406458/137‑3 (MQ=255)
    aTCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTaaac                                                                                                                                 <  1:10205277/137‑3 (MQ=255)
    aTCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTaaac                                                                                                                                 <  2:10784078/137‑3 (MQ=255)
    aTCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTaaac                                                                                                                                 <  2:11625805/137‑3 (MQ=255)
    aTCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTaaac                                                                                                                                 <  2:16348108/137‑3 (MQ=255)
    aTCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTaaac                                                                                                                                 <  2:5205673/137‑3 (MQ=255)
    aTCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTaaac                                                                                                                                 <  1:11323664/137‑3 (MQ=255)
    aTCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTaaac                                                                                                                                 <  1:13174138/137‑3 (MQ=255)
    aTCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTaaac                                                                                                                                 <  1:13465929/137‑3 (MQ=255)
    aTCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTAGAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTaaac                                                                                                                                 <  1:389804/137‑3 (MQ=255)
       gcgACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTaaacccg                                                                                                                              >  2:1513403/1‑132 (MQ=255)
       gcgACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAacccg                                                                                                                             >  1:878354/1‑135 (MQ=255)
       gcgACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCACGGGTaaaccca                                                                                                                              >  2:9548711/1‑132 (MQ=255)
              aCGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTac                                                                                                                       <  1:10045251/137‑1 (MQ=255)
                gtgtGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTacac                                                                                                                     <  1:6897276/137‑1 (MQ=255)
                gtgtGCTGGTCCGGTGGTCCGGTATANATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTacac                                                                                                                     <  1:16671447/137‑1 (MQ=255)
                 tgtgCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTacaca                                                                                                                    <  2:7587523/137‑1 (MQ=255)
                 tgtgCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTacaca                                                                                                                    <  1:5179885/137‑1 (MQ=255)
                      tggtccggtggtccggTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGc                                                                                                               <  2:7288613/137‑1 (MQ=255)
                      tggtccggtggtccggTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGc                                                                                                               <  1:1629685/137‑1 (MQ=255)
                      tggtccggtggtccggTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGc                                                                                                               <  1:5588927/137‑1 (MQ=255)
                      tggtccggtggtccggTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGc                                                                                                               >  1:1464984/1‑137 (MQ=255)
                      tggtccggtggtccggTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGc                                                                                                               <  2:6683326/137‑1 (MQ=255)
                        gtccggtggtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCt                                                                                                             <  2:9151557/137‑1 (MQ=255)
                           cggtggtccggtATATATACCCTAATGATTTATAAAATTTTTCTGAACTAACTTTTTGGAATATTTATAGAAGGGGGGTTTGTTACAGAGGTGCCGGGCTCTTCGCGGGCAAACCCGCTCCTACACACATGCCTCaa                                                                                                          >  2:2785838/1‑137 (MQ=255)
                           cggtggtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCaa                                                                                                          >  1:9026887/1‑137 (MQ=255)
                           cggtggtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCaa                                                                                                          <  2:9687885/137‑1 (MQ=255)
                           cggtggtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCaa                                                                                                          <  2:4523982/137‑1 (MQ=255)
                           cggtggtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCaa                                                                                                          <  2:12370284/137‑1 (MQ=255)
                           cggtggtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCaa                                                                                                          >  2:6902793/1‑137 (MQ=255)
                           cggtggtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCaa                                                                                                          >  2:15661209/1‑137 (MQ=255)
                           cggtggtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCaa                                                                                                          <  1:8537448/137‑1 (MQ=255)
                           cggtggtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCaa                                                                                                          >  1:7555406/1‑137 (MQ=255)
                                gtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCt                                                                                                     <  1:7824427/137‑1 (MQ=255)
                                gtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCt                                                                                                     <  1:10418538/137‑1 (MQ=255)
                                gtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCt                                                                                                     <  1:10699488/137‑1 (MQ=255)
                                gtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCt                                                                                                     <  2:6341884/137‑1 (MQ=255)
                                gtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCt                                                                                                     <  1:15661209/137‑1 (MQ=255)
                                gtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCt                                                                                                     <  1:2706589/137‑1 (MQ=255)
                                gtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCt                                                                                                     <  2:14981287/137‑1 (MQ=255)
                                gtccggtATATATACCCTAATGATTTATAAAATTTTTAGGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCt                                                                                                     <  2:15815176/137‑1 (MQ=255)
                                    ggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGc                                                                                                 <  1:6477953/137‑1 (MQ=255)
                                          tataCCCTAATGATTTATAAAATTTTTATGGACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGATGTGCTGGGCTCTTCGCGGGTAAACCCGCTCCTACACACAGGCCTCAAGATCTGTGCACTGTa                                                                                           <  2:3094174/137‑1 (MQ=39)
                                          tataCCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGATGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTa                                                                                           <  2:1021844/137‑1 (MQ=255)
                                          tataCCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTa                                                                                           <  1:7697725/137‑1 (MQ=255)
                                          tataCCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTa                                                                                           <  1:4583754/137‑1 (MQ=255)
                                          tataCCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATAGTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTGCGCGGGTAAACCCGCGCCTACACACATGCCTCAAGATCTGTGCACTGTa                                                                                           <  2:14060009/137‑1 (MQ=255)
                                            taCCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAgg                                                                                         <  1:12731414/137‑1 (MQ=255)
                                                  aaTGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTacac                                                                                                                     >  2:13425069/1‑103 (MQ=255)
                                                  aaTGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACACGCCTCAAGATCTGTGCACTGTAGGAGCggg                                                                                   >  2:46844/1‑137 (MQ=255)
                                                  aaTGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCACTCCTACACACATGCCTCAAGATTTGTGCACTGTAGGAGAggg                                                                                   >  1:3041319/1‑137 (MQ=255)
                                                   aTGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCAGGt                                                                                  >  1:16469348/1‑137 (MQ=255)
                                                     gATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGttt                                                                                >  2:5710138/1‑137 (MQ=255)
                                                     gATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGttt                                                                                >  1:1980468/1‑137 (MQ=255)
                                                     gATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCCTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGAGGGttt                                                                                >  2:16729967/1‑137 (MQ=255)
                                                      aTTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTa                                                                               >  2:10979690/1‑137 (MQ=255)
                                                           taAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCcgc                                                                          >  2:13202261/1‑137 (MQ=255)
                                                            aaaaTTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTGAACCCACTCCTCCACACACGCCTCAAAAGCTGTGCACCGTAAGAGAGGGTTTCCCcacg                                                                         >  1:4223570/1‑137 (MQ=255)
                                                            aaaaTTTTTAGGAGTTAGTTTTTGGGAATCTTTGCAGAGGGGGGGTTTGTGACGGAGGTGCTGGGCTCTTTGGGGGAAAACCCTCTCCCAAACCCAGCCCGCAAGCTCTGTGCCTTGTAGGAGCGGGTTTACCcgcg                                                                         <  2:11541790/137‑1 (MQ=255)
                                                             aaaTTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGTACTGTAGGAGCGGGTTTCCCCGCGa                                                                        >  1:14068758/1‑137 (MQ=255)
                                                                    tATGAACTAATTTTTTGGAAGATTTGTAGAGGGGGGATTTGTTACGGGGGTGCCGGGCTTTTTGCGGGGAAACCCGCTCCGACACACATGCCTCAAGATCTGTGCACTTTGGGAGTGGGTTTACCCGCGAAAGGGcc                                                                 <  1:15075668/137‑1 (MQ=255)
                                                                    tATGAACTAACTTTTTGGAATATTTAGAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGcc                                                                 <  2:15982230/137‑1 (MQ=255)
                                                                     agaaactaactTTTTGGAATATTGAGAGAAGGGGGATTTGTTACAGAGTTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACGTGCCTCAAGCTCTGTGCATTGGAGGAGTGGGTTTACCCGCGAAAGGGCCg                                                                <  1:11643638/134‑1 (MQ=255)
                                                                                   tGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTCCACACATGCCTCAAGATTTGTGCACTGTAGGGGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTc                                                  >  1:11404335/1‑137 (MQ=39)
                                                                                        tatTTATAGAAGGGGAATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCACTCCTACACACATTCCTCAAGACCTGTGCATTGGAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAgg                                             <  2:8480919/137‑1 (MQ=38)
                                                                                          tttATAGAGGGGGGGTTTGTTAAAGAGGTGCCGGGCTCTTTGCGGGTAAACCCTCTCCTACACACATCCCTCAAGATCTGTGCACTGGAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGcc                                           <  1:11987453/137‑1 (MQ=38)
                                                                                          ttgataGAAGGGGGATTTGTTACAGAGGTTCTGGGCTCTTTGCGGGTAACCCCGCTCCTCCACACATGCCTCAAGATCTGTACACTGGAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGcc                                           <  1:2111622/134‑1 (MQ=38)
                                                                                            tataGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGc                                         >  2:14263283/1‑137 (MQ=255)
                                                                                                   gggggATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTATAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAATCCCg                                  >  1:3319240/1‑137 (MQ=255)
                                                                                                        aTTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTCCCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGaa                             >  2:868966/1‑137 (MQ=255)
                                                                                                          ttGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTCCACACATGCCTCAAGGTCTGTGCACTGTAGGAGCGGGTTTACCCGCAAAAGGGCCGGGACAGGTTTCCTCGGAGGCCGCAGTCCCGAGGAACt                           >  1:12368898/1‑137 (MQ=38)
                                                                                                                    ggTGCCGGGCTCTTCGCGGGTAAACCCACTCCTCCACACATGCCTCAAAATCTGTGCACCGTAGGAGCGGGTTTACCCACGAAAGGGCCGGAACAGGTTTCCTCGGAGGCCGCAGTCCCGAGGAACTTTCGCAAGaa                 >  1:713616/1‑137 (MQ=38)
                                                                                                                        ccGGGCTCTTCGCGGGTAAACCCGCTCCTACACACCTGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGAGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGc             >  2:241056/1‑137 (MQ=255)
                                                                                                                        ccGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTCCTCGGAGGCCGCAGTCCCGAGGAACTTGTGCAAGAACTGc             >  2:12944885/1‑137 (MQ=255)
                                                                                                                        ccGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGc             >  2:10501860/1‑137 (MQ=255)
                                                                                                                        ccGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGc             >  2:11913672/1‑137 (MQ=255)
                                                                                                                        ccGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGc             >  2:824798/1‑137 (MQ=255)
                                                                                                                         cGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATTTGTGCACTGTAGGAGTGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGcc            <  1:16175949/137‑1 (MQ=255)
                                                                                                                          gggCTTTTCGTGGGGAAACCCGCCCCTACACACATGCCTCAAGATCTTTGCATTGTAGGAGTGGTTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCcg           <  1:2540044/137‑1 (MQ=255)
                                                                                                                          gggCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCcg           <  2:5365750/137‑1 (MQ=255)
                                                                                                                          gggCCCTTCGCGGGTAAACCCGCTTCTACACCCATGCCTCAAGATATGTGCACTGTAGGAGCGGTTTTCCCCGCGAAAGGGCCGGAACAGTTTACTCGGAGGCCGCAGTCCCGAGGAACTGGCGCAAGGACTGCcgc          <  2:14629683/137‑1 (MQ=255)
                                                                                                                           ggCTCTTCGCGGGTAAACCCGCTCCTCCACACATGCCTCAAGATCTGTGCAATGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTCCTCTGAGGCCGCAGTCCCAAGGAACTTTCGCAAAAACTGccac          >  1:527624/1‑135 (MQ=37)
                                                                                                                           ggCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAGCCGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGACCTGCcgc          >  2:15287053/1‑137 (MQ=39)
                                                                                                                               ctTCGCGGGGAAACCTGCTCCTACACACATCCCTCAAGATCTGTGCATTGTGGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTAc      <  2:4259882/137‑1 (MQ=38)
                                                                                                                               ctTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATTTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTAc      <  2:1248778/137‑1 (MQ=255)
                                                                                                                               ctTCGCGGGGAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTAc      <  2:5060283/137‑1 (MQ=255)
                                                                                                                                  cgcgGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGTCCTTTAGGGGCGGGTTTCCCCGCGAAAGGGCCGGAACAGGTTTCCTCGGAGGCCGCAGTCCCAAGGAACTTTCGCCAGAACCGCCCCGTCCGct   >  1:9285280/1‑137 (MQ=255)
                                                                                                                                  cgcgGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGGACTTGCGCAAGAACTGCCGCGTACGct   >  2:4805555/1‑137 (MQ=255)
                                                                                                                                   gcgGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTCCTCGGAGGCCGCAGTCCCGAGGAACTTGCGCACGAACTGCCGCGTACGctc  >  2:1945794/1‑137 (MQ=255)
                                                                                                                                           |                                                                                                                                   
GCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTC  >  AE015451_phaM1/278606‑278872

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