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.9% 2143.0 / 101.9 676intergenic (‑237/+77)PP_0224 CDS (Monooxygenase_ DszC family)/sctC CDS (sulfur compound ABC transporter ‑ ATP‑binding subunit)–/–
Reads supporting (aligned to +/- strand):  ref base . (4/29);  new base A (313/303);  total (344/332)
Fisher's exact test for biased strand distribution p-value = 7.91e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 789 total aligned reads displayed.

GCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGC‑TCTTCG‑C‑GGGT‑A‑ACCCGCTCCTACACACATGCCTCAAGATCTGTG‑CACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTC  >  AE015451_phaM1/278606‑278875
                                                                                                                                            |                                                                                                                                       
gCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGC‑TCTTCG‑C‑GGGt                                                                                                                                          <  1:4413133/137‑1 (MQ=255)
gCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGC‑TCTTCG‑C‑GGGt                                                                                                                                          >  2:13777300/1‑137 (MQ=255)
  tGATCGCGACAGACGTGTGCTGTTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGC‑TCTTCG‑C‑GGGT‑a‑a                                                                                                                                      <  1:13801011/137‑1 (MQ=255)
  tGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGC‑TCTTCG‑C‑GGGT‑a‑a                                                                                                                                      <  2:6723070/137‑1 (MQ=255)
  tGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGC‑TCTTCG‑C‑GGGT‑a‑a                                                                                                                                      <  1:7051562/137‑1 (MQ=255)
  tGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGC‑TCTTCG‑C‑GGGT‑a‑a                                                                                                                                      <  2:1095714/137‑1 (MQ=255)
  tGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGC‑TCTTCG‑C‑GGGT‑a‑a                                                                                                                                      <  2:13211749/137‑1 (MQ=255)
  tGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGC‑TCTTCG‑C‑GGGT‑a‑a                                                                                                                                      <  1:1310547/137‑1 (MQ=255)
  tGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGC‑TCTTCG‑C‑GGGT‑a‑a                                                                                                                                      <  2:290196/137‑1 (MQ=255)
    aTCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGC‑TCTTCG‑C‑GGGT‑a‑aac                                                                                                                                    <  2:8555490/137‑3 (MQ=255)
    aTCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGC‑TCTTCG‑C‑GGGT‑a‑aac                                                                                                                                    <  1:6734277/137‑3 (MQ=255)
       gcgACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGC‑TCTTCG‑C‑GGGT‑a‑aacccg                                                                                                                                 >  1:10008109/1‑132 (MQ=255)
       gcgACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGC‑TCTTCG‑C‑GGGT‑a‑aacccg                                                                                                                                 >  2:6236100/1‑132 (MQ=255)
       gcgACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGC‑TCTTCG‑C‑GGGT‑a‑aac                                                                                                                                    <  1:16621012/134‑3 (MQ=255)
         gACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGC‑TCTTCG‑C‑GGGTAA‑ACCCGCt                                                                                                                                <  2:14011247/137‑1 (MQ=255)
          aCAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGC‑TCTTCG‑C‑GGGTAA‑ACCCGCTc                                                                                                                               <  1:1517679/137‑1 (MQ=255)
           cAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGC‑TCTTCG‑C‑GGGTAA‑ACCCGCTcc                                                                                                                              >  2:4968602/1‑137 (MQ=255)
               cGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGC‑TCTTCG‑C‑GGGTAA‑ACCCGCTCCTaca                                                                                                                          <  1:12685231/137‑1 (MQ=255)
               cGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGC‑TCTTCG‑C‑GGGTAA‑ACCCGCTCCTaca                                                                                                                          <  1:1822191/137‑1 (MQ=255)
                gtgtGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGC‑TCTTCG‑C‑GGGTAA‑ACCCGCTCCTacac                                                                                                                         <  2:1132894/137‑1 (MQ=255)
                 tgtgCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGATGTGCCGGGC‑TCTTCG‑C‑GGGTAA‑ACCCGCTCCTacaca                                                                                                                        <  1:7249503/137‑1 (MQ=255)
                 tgtgCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGC‑TCTTCG‑C‑GGGTAA‑ACCCGCTCCTacaca                                                                                                                        <  1:5217234/137‑1 (MQ=255)
                 tgtgCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGC‑TCTTCG‑C‑GGGTAA‑ACCCGCTCCTacaca                                                                                                                        <  1:11792662/137‑1 (MQ=255)
                 tgtgCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGC‑TCTTCG‑C‑GGGTAA‑ACCCGCTCCTacaca                                                                                                                        <  2:5502708/137‑1 (MQ=255)
                      tggtccggtggtccggTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGC‑TCTTCG‑C‑GGGTAA‑ACCCGCTCCTACACACATGc                                                                                                                   >  1:9281646/1‑137 (MQ=255)
                           cggtggtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGC‑TCTTCG‑C‑GGGTAA‑ACCCGCTCCTACACACATGCCTCaa                                                                                                              >  1:10290843/1‑137 (MQ=255)
                           cggtggtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGC‑TCTTCG‑C‑GGGTAA‑ACCCGCTCCTACACACATGCCTCaa                                                                                                              <  2:11949558/137‑1 (MQ=255)
                           cggtggtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGC‑TCTTCG‑C‑GGGTAA‑ACCCGCTCCTACACACATGCCTCaa                                                                                                              <  2:4380559/137‑1 (MQ=255)
                           cggtggtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGC‑TCTTCG‑C‑GGGTAA‑ACCCGCTCCTACACACATGCCTCaa                                                                                                              <  2:2332192/137‑1 (MQ=255)
                           cggtggtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGC‑TCTTCG‑C‑GGGTAA‑ACCCGCTCCTACACACATGCCTCaa                                                                                                              >  1:4162881/1‑137 (MQ=255)
                           cggtggtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGC‑TCTTCG‑C‑GGGTAA‑ACCCGCTCCTACACACATGCCTCaa                                                                                                              >  1:14121737/1‑137 (MQ=255)
                           cggtggtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGC‑TCTTCG‑C‑GGGTAA‑ACCCGCTCCTACACACATGCCTCaa                                                                                                              <  1:10779042/137‑1 (MQ=255)
                                gtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGC‑TCTTCG‑C‑GGGTAA‑ACCCGCTCCTACACACATGCCTCAAGATCt                                                                                                         <  2:17061300/137‑1 (MQ=255)
                                gtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGC‑TCTTCG‑C‑GGGTAA‑ACCCGCTCCTACACACATGCCTCAAGATCt                                                                                                         <  1:3719343/137‑1 (MQ=255)
                                gtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGC‑TCTTCG‑C‑GGGTAA‑ACCCGCTCCTACACACATGCCTCAAGATCt                                                                                                         <  2:12845976/137‑1 (MQ=255)
                                gtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGC‑TCTTCG‑C‑GGGTAA‑ACCCGCTCCTACACACATGCCTCAAGATCt                                                                                                         <  1:15848061/137‑1 (MQ=255)
                                gtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGC‑TCTTCG‑C‑GGGTAA‑ACCCGCTCCTACACACATGCCTCAAGATCt                                                                                                         <  1:9658433/137‑1 (MQ=255)
                                gtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGC‑TCTTCG‑C‑GGGTAA‑ACCCGCTCCTACACACATGCCTCAAGATCt                                                                                                         <  2:5684906/137‑1 (MQ=255)
                                gtccggtATATATACCCTAATGATTTATAAAAGTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGC‑TCTTCG‑C‑GGGTAA‑ACCCGCTCCTACACACATGCCTCAAGATCt                                                                                                         <  2:15311042/137‑1 (MQ=255)
                                     gtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGC‑TCTTCG‑C‑GGGTAA‑ACCCGCTCCTACACACACGCCTCAAGATCTGTG‑Ca                                                                                                   >  1:2648261/1‑137 (MQ=255)
                                      tatatataCCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGC‑TCTTCG‑C‑GGGTAA‑ACCCGCTCCTACACACATGCCTCAAGGTCTGTG‑CAc                                                                                                  >  2:3205416/1‑137 (MQ=255)
                                         atataCCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGC‑TCTTCG‑C‑GGGTCA‑ACCTGCTCCTCCACCCATACCTCAAGATCTGTG‑CACTGt                                                                                               >  1:14607584/1‑137 (MQ=255)
                                          tataCCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGC‑TCTTCG‑C‑GGGTAA‑ACCCGCTCCTACACACATGCCTCAAGATCTGTG‑CACTGTa                                                                                              <  2:16235021/137‑1 (MQ=255)
                                          tataCCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGC‑TCTTCG‑C‑GGGTAA‑ACCCGCTCCTACACACATGCCTCAAGATCTGTG‑CACTGTa                                                                                              <  2:17669999/137‑1 (MQ=255)
                                          tataCCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGC‑TCTTCG‑C‑GGGTAA‑ACCCGCTCCTACACACATGCCTCAAGATCTGTG‑CACTGTa                                                                                              <  1:14945572/137‑1 (MQ=255)
                                          tataCCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGC‑TCTTCG‑C‑GGGTAA‑ACCCGCTCCTACACACATGCCTCAAGATCTGTG‑CACTGTa                                                                                              <  1:497892/137‑1 (MQ=255)
                                          tataCCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGC‑TCTTCG‑C‑GGGTAA‑ACCCGCTCCTACACACATGCCTCAAGATCTGTG‑CACTGTa                                                                                              <  1:111131/137‑1 (MQ=255)
                                          tataCCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGC‑TCTTCG‑C‑GGGTAA‑ACCCGCTCCTACACACATGCCTCAAGATCTGTG‑CACTGTa                                                                                              <  2:9323675/137‑1 (MQ=255)
                                                  aaTGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGCAGGGGGGTTTGTTACAGAGGTGCCGAGC‑TCTTCG‑C‑AGGTAA‑AACCGCTCCTACCCACATGCCTCAAGATCTGTG‑CACTGTAGGAGCggg                                                                                      >  2:17419935/1‑137 (MQ=255)
                                                  aaTGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGC‑TCTTCG‑C‑GGGTAA‑ACCCGCTCCTACACACATGCCTCAAGATCTGTG‑CACTGTAGGAGCggg                                                                                      >  2:15630421/1‑137 (MQ=255)
                                                  aaTGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGC‑TCTTCG‑C‑GGGTAA‑ACCCGCTCCTACACACATGCCTCAAGATCTGTG‑CACTGTAGGAGCggg                                                                                      >  2:15052302/1‑137 (MQ=255)
                                                  aaTGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGC‑TCTTCG‑C‑GGGTAA‑ACCCGCTCCTACACACATGCCTCAAGATCTGTG‑CACTGTAGGAGCggg                                                                                      >  2:892599/1‑137 (MQ=255)
                                                   aTGATTTATAACCTTTTTATGAACGAAATTTTTGGAGTATTTATAGAAAGGGGACTTGTTACCGAGGTGCCGGGC‑TCTTCC‑C‑GGGC‑ACACCCGCTCCTACACACATGCTGTAAGATCTGTG‑CACTGTAGTAGCgggg                                                                                     >  2:16603111/1‑136 (MQ=255)
                                                   aTGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGC‑TCTTCG‑C‑GGGTCA‑ACCCGCTCCTACACACATGCCTCAAGATTTGTG‑CACTGTAGGAGCGGGt                                                                                     >  1:6088024/1‑137 (MQ=255)
                                                     gATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTCTAGAAGGGGGATTTGTTACAGAGGTGCCGGGC‑TCTTCG‑C‑AGGTCA‑ACCCGCTCCTACACACATGCCTCCAGATCTGTG‑CACTGTAGGAGCGGGctt                                                                                   >  2:4003554/1‑137 (MQ=255)
                                                     gATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGC‑TCTTCG‑C‑GGGTAA‑ACCCGCTCCTacaca                                                                                                                        >  2:13516442/1‑101 (MQ=255)
                                                     gATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGC‑TCTTCG‑C‑GGGTAA‑ACCCGCTCCTACACACATGCCTCAAGATCTGTG‑TACTGTAGGAGCGGGttt                                                                                   >  1:5925038/1‑137 (MQ=255)
                                                     gATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGC‑TCTTCG‑C‑GGGTAA‑ACCCGCTCCTACACACATGCCTCAAGATCTGTG‑CACTGTAGGAGCGGGttt                                                                                   >  2:4622679/1‑137 (MQ=255)
                                                     gATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGC‑TCTTCG‑C‑GGGTAA‑ACCCGCTCCTACACACATGCCTCAAGATCTGTG‑CACTGTAGGAGCGGGttt                                                                                   >  2:7543521/1‑137 (MQ=255)
                                                     gATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGC‑TCTTCG‑C‑GGGTAA‑ACCCGCTCCTACACACATGCCTCAAGATCTGTG‑CACTGTAGGAGCGGGttt                                                                                   >  2:6907022/1‑137 (MQ=255)
                                                     gATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGC‑TCTTCG‑C‑GGGTAA‑ACCCGCTCCTACACACATGCCTCAAGATCTGTG‑CACTGTAGGAGCGGGttt                                                                                   >  2:10586326/1‑137 (MQ=255)
                                                     gATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGC‑TCTTCG‑C‑GGGTAA‑ACCCGCTCCTACACACATGCCTCAAGATCTGTG‑CACTGTAGGAGCGGGttt                                                                                   >  2:102886/1‑137 (MQ=255)
                                                     gATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGC‑TCTTCG‑C‑GGGTAA‑ACCCGCTCCTACACACATGCCTCAAGATCTGTG‑CACTGTAGGAGCGGGttt                                                                                   >  1:8797505/1‑137 (MQ=255)
                                                     gATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGC‑TCTTCG‑C‑GGGTAA‑ACCCGCCCCTACACACACGCCTCAAGATCTGTG‑CACTGTAGGAGCGGGttt                                                                                   >  1:7859940/1‑137 (MQ=255)
                                                     gATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGC‑TCTTCG‑C‑GGGT‑a‑aac                                                                                                                                    <  2:3525727/88‑3 (MQ=255)
                                                      aTTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGC‑TCTTCG‑C‑GGGTAA‑ACCCGCTCCTCCACACATGCCTCAAGATCTGTG‑CACTGTAGGAGCGGGtttc                                                                                  >  1:12349721/1‑136 (MQ=255)
                                                       tttATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGC‑TCTTCG‑C‑GGGTAA‑ACCCGCTCCTACACACATGCCTCAAGATCTGTG‑CACTGTAGGAGCGGGTTTAc                                                                                 >  2:12136621/1‑137 (MQ=255)
                                                           taAAATTTTTGTGAACTAACTTTTTGGAAGATTTATAGAAGGGGGACTTGTCACAGAGTTCCTGGGC‑TTTTCG‑C‑GGGGAA‑ACCAGCTCCTACACACATGCCTCAAGACCTGTG‑CACTGTAGGAGTGGGTTTACCcgc                                                                             <  2:9542636/137‑1 (MQ=255)
                                                           taAAATTTTTATGAACTAACTTTCTGGACTATTGATAGAAGGGGGATTTGTTACAGAGGTGCCAGGC‑TCTTCG‑C‑GGGTAA‑ACCCGCTCCTACACACATGCCTCAAGATCTTTG‑TACTGTAGGAGCGTGCTTACCcga                                                                             >  2:8242828/1‑136 (MQ=255)
                                                           taAAATTTTGGGGAATTAATTTTTTGGAAGATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGT‑CCTTCGTG‑GGGA‑A‑ACCCGCTCCTACACACGGGCCTCAAGATCGGGG‑CACTGTAGGAGCGGGTTTACCcgc                                                                             <  2:3763229/137‑1 (MQ=255)
                                                            aaaaTTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGC‑TCTTCG‑C‑GGGTAA‑ACCCGCTCCTACACACATGCCTCAAGATCt                                                                                                         >  1:17329728/1‑109 (MQ=255)
                                                            aaaaTTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGC‑TCTTCG‑C‑GGGTAA‑ACCCACTCCTCCACACCTGCCTCAAAATCTGTG‑CACTGTAGGAGCGGGTTTCCCcacg                                                                            >  1:6551174/1‑137 (MQ=255)
                                                                     gtGAACTAATTTTTTGGAGTATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGC‑TTTTCG‑C‑GGGGAA‑ACCCGCTCCTACACACATGCCTCAAGATCTGTG‑CACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCg                                                                   <  2:2419057/136‑1 (MQ=255)
                                                                     aTGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGC‑TCTTCG‑C‑GGGGAA‑ACCCGCTCCTACACACATGCCTCAAGATCTGTG‑CACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCg                                                                   <  2:6356217/137‑1 (MQ=255)
                                                                     aTGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGC‑TCTTCG‑C‑GGGGAA‑ACCCGCTCCTACACACATGCCTCAAGATCTGTG‑CACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCg                                                                   <  2:8101168/137‑1 (MQ=255)
                                                                     aTGAACTAACTTTTTGGAATATTGAGAGAAGGGGGGTTTGTTACAGAGGTGTTGGGT‑TTTTTG‑G‑GGGTAA‑ACCTTCTCCTACACACGTGCCTCAAGATTTTTG‑CACTGGAGGAGTGGGTTTACCCGCGAAAGGGCCg                                                                   <  1:11674188/137‑1 (MQ=255)
                                                                                ttttGGAATATTTATAGAGGGGGGATTTGTTACAGAGGTGCCTGGC‑TCTTCG‑C‑GGGTAA‑ACCCGTTCCTACACACATGCCTCAAGATCTGTG‑CACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTa                                                        <  2:11717514/137‑1 (MQ=39)
                                                                                   tGGAATATTGATAGAGGGGGGATTTTTTACAGAGGTGCCGGGC‑TCTTCG‑CGGGGA‑A‑ACCCGCTCCTACACACATGCCTCAAGATCTGTG‑CATTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTc                                                     <  1:8102508/137‑1 (MQ=38)
                                                                                     gAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGC‑TCTTCG‑C‑GGGTGA‑AACCACTCCCCCCCACACGCCCCCAAAGCTGTG‑CACCGTAGGAGAGGGTTTCCCCCCAAAAAGGACAGAACAAGTTTCCTCgg                                                   >  1:6348950/1‑137 (MQ=255)
                                                                                     gAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGC‑TCTTCG‑C‑GGGTAA‑ACCCGCTCCTACACACATGCCTCAAGATCTGTG‑CACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCgg                                                   >  1:16316407/1‑137 (MQ=255)
                                                                                        gatTTAGAGAAGGGGGATTTGTTACAGAGGTGCCGGGC‑TCTTCG‑C‑GGGTAA‑ACCCGCTCCTACACACATGCCTCAAGATCTGTG‑CACTGTGGGAGCGGTTTCACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAgg                                                <  2:14240342/136‑1 (MQ=39)
                                                                                          tttATAGAAGGGGGATTTGTTACAGAGGTGCCGGGT‑TCTTCG‑C‑GGGTAA‑ACCCGCTCCTACACACATGCCTCAAGATCTGTG‑CACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGcc                                              <  2:10002410/137‑1 (MQ=255)
                                                                                           ttAGAGAGGGGGGATTTGTGAAGGAGGTGCTGGGCTTTTTGG‑G‑GGGA‑A‑ACCCTCTCCGACACACATTCCTAAAGATTTTTG‑CACTGTGGGAGTGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCg                                             <  1:3477422/137‑1 (MQ=255)
                                                                                            tgtaGAAGGGGGATTTGTTACGGAGGTGCCGGGC‑TTTTCG‑C‑GGGTAA‑ACCCACTCCTACACACATGCCTCAAGATCTGTG‑CATTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGc                                            <  2:5294034/135‑1 (MQ=38)
                                                                                                       gATTTGTTACAGAGGTGCCGGGC‑TCTTCG‑C‑GGGTAA‑ACCCGCTCCTCCACACATGCCTCAAGATCTGTG‑CACTGTAGGAGCGGGTTTCCCCGCGAAAAGGCCGGAACAGGTTTCCTCGGAGGCCGCAGTCCCAAGGa                                 >  1:3784148/1‑137 (MQ=38)
                                                                                                       gATTTGTTACAGAGGTGCCGGGC‑TCTTCG‑C‑GGGTAA‑ACCCGCTCCTCCACACATGCCTCAAGATCTGTG‑CACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGa                                 >  2:1466499/1‑137 (MQ=255)
                                                                                                       gATTTGTTACAGAGGTGCCGGGC‑TCTTCG‑C‑GGGTAA‑ACCCGCTCCTACACACATGCCTCAAGATCTGTG‑CACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGa                                 >  2:15873970/1‑137 (MQ=255)
                                                                                                        aTTTGTTAAAGAGGTGCCGGGC‑TCTTCG‑C‑GGGTAA‑ACCCGCTCCTCCACACATGCCTCAAGATCTGTG‑CACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAAAAGGTTTACTCGGAGGCCGCAGTCCCGAGGaa                                >  2:7103556/1‑137 (MQ=39)
                                                                                                              tACAGAGGTGCCGGGC‑TCTTCG‑C‑GGGTGA‑ACCCACTCCCCCACCCCCACCCCCAAACCTTTGTCAC‑GTCGGAGCGGGTTTTCCCCCGAAAAGGCCGGAACAGGGTTCCTCTGAGGCCACAATCCCAAAGAAccttcg                          >  1:1352789/1‑132 (MQ=255)
                                                                                                              tACAGAGGTGCCGGGC‑TCTTCG‑C‑GGGTAA‑ACCCGCTCCTACACACATGCCTCAAGATCTGTG‑CACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTgcg                          >  2:7661718/1‑137 (MQ=255)
                                                                                                               aCAGAGGTGCCGGGC‑TCTTCG‑C‑GGGTAA‑ACCCGCTCCTACACACATGCCTCAAGATCTGTG‑CACTGGAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTgcgc                         <  2:6558080/137‑1 (MQ=255)
                                                                                                                        ccGGGC‑TCTTCG‑C‑GGGTAA‑ACCCGCTCCTACACACATGCATCAAGGTCTGTG‑CCTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTCCTCGGAGGCCGCAGTCCCGAGGAACTTGCACAAGAACCGc                >  2:2020677/1‑137 (MQ=38)
                                                                                                                        ccGGGC‑TCTTCG‑C‑GGGTAA‑ACCCACCCCTCCACACATGCCTCCAAATCTTTG‑CACCTTAGGAGCGGGTTTCCCCGCGAAAGGGCCGGAACAGGTTTCCTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAAAACTGc                >  1:8308567/1‑137 (MQ=37)
                                                                                                                         cGGGC‑TCTTCGTG‑GGGA‑A‑ACCTGCTCCTACACACATGCCTCAAGGTCTGTG‑CATTGGAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGcc               <  1:240769/137‑1 (MQ=38)
                                                                                                                           ggC‑TCTTCG‑C‑GGGGAA‑ACCCGCTCCTACACACGTGCCTCAAGATCTGGG‑CACTGGAGGAGTGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCcgc             <  1:3130340/137‑1 (MQ=255)
                                                                                                                                tttCG‑C‑GGGTAA‑ACCCGCTCCTACACACATGCCTAAAGATCTGTG‑CACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTAc         <  2:1386275/136‑1 (MQ=255)
                                                                                                                                ctTCG‑C‑GGGTAA‑ACCCGCTCCTACACACATGCCTCAAGATCTGTT‑CACTGGAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTAc         <  2:16112022/137‑1 (MQ=255)
                                                                                                                                   cg‑c‑gGGTAA‑ACCCGCTCCTCCACACCTGCCTCAAAATCTGTG‑CACTGTAGGAGCGGGTTTCCCCGCGAAAGGGCCGGAACAAGTTTCCTCGGAGGCCGCAGTCCCAAAGAACTTGCGCAAGGACTGCCGCGTACAct      >  1:15683189/1‑137 (MQ=255)
                                                                                                                                    g‑c‑gGGTAA‑ACCCGCTCCTCCACACATGTCTCAAGATCTGTG‑TAATGTAGGAGCGGGTTTACCCGAGAAAGGGCCGGAACAGGGTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGctt     >  1:11341745/1‑136 (MQ=38)
                                                                                                                                        gggA‑A‑ACCCGCTCCTACACACATGCCTCAAGATCTGTG‑CACTTGAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCAGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGttcttc  <  1:11965475/137‑1 (MQ=39)
                                                                                                                                            |                                                                                                                                       
GCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGC‑TCTTCG‑C‑GGGT‑A‑ACCCGCTCCTACACACATGCCTCAAGATCTGTG‑CACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTC  >  AE015451_phaM1/278606‑278875

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: ‑