Predicted mutation
evidence seq id position mutation annotation gene description
RA AE015451_phaM1 5,659,798 G→A R281H (CGC→CAC)  metK CDS (methionine adenosyltransferase) →

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451_phaM15,659,7980GA99.7% 2283.6 / NA 631R281H (CGC→CAC) metK CDS (methionine adenosyltransferase)
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (334/295);  total (335/296)
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 654 total aligned reads displayed.

CCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCGCTACGTGGCCAAGAACATCGTCGCCGCCGGTTTGGCCGAGCGTTGCGAGATCCAGGTGTCGTACGCCATTGGCGTGGCCCAGCCGACCTCCATCTCGATCAACACCTTCGGTACCGGCAAGGTCTCCGACGAC  >  AE015451_phaM1/5659664‑5659931
                                                                                                                                      |                                                                                                                                     
  cGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCACTAc                                                                                                                                   >  1:13133508/1‑137 (MQ=255)
   ggtgggtgACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCACTACg                                                                                                                                  >  1:16124230/1‑137 (MQ=255)
      gggtgACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCACTACGTgg                                                                                                                               >  2:1195867/1‑137 (MQ=255)
         tgACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTCCCGCTGCCTACGCCGGCCACTACGTGGcaa                                                                                                                            >  2:16803861/1‑135 (MQ=255)
          gACTGCGGCCTGACTGGCCGCAAGATCATTGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACACGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCACTACGTGGCCaa                                                                                                                           >  2:14678179/1‑137 (MQ=255)
                  ccTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAACAtcg                                                                                                                   >  1:7718721/1‑137 (MQ=255)
                  ccTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAACAtcg                                                                                                                   >  2:14993804/1‑137 (MQ=255)
                  ccTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAACAtcg                                                                                                                   >  1:11709977/1‑137 (MQ=255)
                  ccTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAACAtcg                                                                                                                   >  2:17461686/1‑137 (MQ=255)
                  ccTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAACAtcg                                                                                                                   >  2:13624603/1‑137 (MQ=255)
                    tGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGATCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCGCTACGTGGCCAAGAACAtcgtc                                                                                                                 >  2:14218059/1‑137 (MQ=255)
                    tGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAACAtcgtc                                                                                                                 >  1:12253575/1‑137 (MQ=255)
                     gACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAACAtcgtcg                                                                                                                >  1:2019215/1‑137 (MQ=255)
                      aCTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAACATCGTcgc                                                                                                               >  2:16498539/1‑137 (MQ=255)
                         ggCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAACATCGTcgccgc                                                                                                            >  1:15531870/1‑137 (MQ=255)
                         ggCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAACATCGTcgccgc                                                                                                            >  2:859832/1‑137 (MQ=255)
                         ggCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCAATACGTGGCCAAGAACATCGTcgccgc                                                                                                            >  2:8220681/1‑137 (MQ=255)
                           ccGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAACATCGTCGccgccg                                                                                                          <  2:15873537/137‑1 (MQ=255)
                            cGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAACATCGTCGCCGCCgg                                                                                                         >  1:10675512/1‑137 (MQ=255)
                             gCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAACATCGTCGCCGCCGGt                                                                                                        >  1:3836519/1‑137 (MQ=255)
                              cAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAACATCGTCGCCGCCGGtt                                                                                                       >  2:9638910/1‑137 (MQ=255)
                               aaGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGTTCGACCGTTCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAACATCGTCGCCGCCGGttt                                                                                                      <  1:9124131/137‑1 (MQ=255)
                                aGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAACATCGTCGCCGCCGGTTTg                                                                                                     <  2:4570436/137‑1 (MQ=255)
                                   tcatcGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGNTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAACATCGTCGCCGCCGGTTTGGcc                                                                                                  <  1:14248208/137‑1 (MQ=255)
                                     atcGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAACATCGTCGCCGCCGGTTTGGCCGa                                                                                                <  2:9827990/137‑1 (MQ=255)
                                     atcGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAACATCGTCGCCGCCGGTTTGGCCGa                                                                                                <  2:752394/137‑1 (MQ=255)
                                     atcGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAACATCGTCGCCGCCGGTTTGGCCGa                                                                                                <  1:12961785/137‑1 (MQ=255)
                                        gtcgACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAACATCGTCGCCGCCGGTTTGGCCGAGCg                                                                                             <  2:13903598/137‑1 (MQ=255)
                                        gtcgACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAACATCGTCGCCGCCGGTTTGGCCGAGCg                                                                                             <  2:6098808/137‑1 (MQ=255)
                                        gtcgACCCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAACATCGTCGCCGCCGGTTTGGCCGAGCg                                                                                             <  1:6085920/137‑1 (MQ=255)
                                         tcgACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAACATCGTCGCCGCCGGTTTGGCCGAGCGt                                                                                            <  2:15383774/137‑1 (MQ=255)
                                          cgACTCCTACGGTTGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCTGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAACATCGTCGCCGTCGGTTTGGCCGAGCGtt                                                                                           <  2:1158798/137‑1 (MQ=255)
                                          cgACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAACATCGTCGCCGCCGGTTTGGCCGAGCGtt                                                                                           <  1:8786492/137‑1 (MQ=255)
                                          cgACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAACATCGTCGCCGCCGGTTTGGCCGAGCGtt                                                                                           <  1:4876618/137‑1 (MQ=255)
                                           gACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAACATCGTCGCCGCCGGTTTGGCCGAGCGTTg                                                                                          <  1:4306422/137‑1 (MQ=255)
                                              tCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAACATCGTCGCCGCCGGTTTGGCCGAGCGTTGCga                                                                                       >  2:6787288/1‑137 (MQ=255)
                                                cTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAACATCGTCGCCGCCGGTTTGGCCGAGCGTTGCgaga                                                                                     >  1:6423280/1‑137 (MQ=255)
                                                  aCGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAACATCGTCGCCGCCGGTTTGGCCGAGCGTTGCGAGATc                                                                                   <  2:13141136/137‑1 (MQ=255)
                                                  aCGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAACATCGTCGCCGCCGGTTTGGCCGAGCGTTGCGAGATc                                                                                   <  1:10276313/137‑1 (MQ=255)
                                                   cGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAACATCGTCGCCGCCGGTTTGGCCGAGCGTTGCGAGATcc                                                                                  <  2:4110547/137‑1 (MQ=255)
                                                          aTGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAACATCGTCGCCGCCGGTTTGGCCGAGCGTTGCGAGATCCAGGTGTc                                                                           <  1:7264460/137‑1 (MQ=255)
                                                          aTGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAACATCGTCGCCGCCGGTTTGGCCGAGCGTTGCGAGATCCAGGTGTc                                                                           <  1:1114345/137‑1 (MQ=255)
                                                          aTGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAACATCGTCGCCGCCGGTTTGGCCGAGCGTTGCGAGATCCAGGTGTc                                                                           <  2:5418859/137‑1 (MQ=255)
                                                          aTGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAACATCGTCGCCGCCGGTTTGGCCGAGCGTTGCGAGATCCAGGTGTc                                                                           <  2:3402467/137‑1 (MQ=255)
                                                          aTGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAACATCGTCGCCGCCGGTTTGGCCGAGCGTTGCGAGATCCAGGTGTc                                                                           <  1:17665800/137‑1 (MQ=255)
                                                           tGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAACATCGTCGCCGCCGGTTTGGCCGCGCGTTGCGAGATCCAGGTGTCg                                                                          >  2:3009376/1‑137 (MQ=255)
                                                           tGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAACATCGTCGCCGCCGGTTTGGCCGAGCGTTGCGAGATCCAGGTGTCg                                                                          >  2:1003855/1‑137 (MQ=255)
                                                           tGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTAGACCGTTCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAACATCGTCGCCGCCGGTTTGGCCGAGCGTTGCGAGATCCAGGTGTCg                                                                          >  1:8501913/1‑137 (MQ=255)
                                                            ggCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAACATCGTCGCCGCCGGTTTGGCCGAGCGTTGCGAGATCCAGGTGTCGt                                                                         >  2:10882201/1‑137 (MQ=255)
                                                               ccGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAACATCTTCGCCGCCGGTTTGGCCGAGCGTTGCGAGATCCAGGTGTCGTACg                                                                      <  1:5432225/137‑1 (MQ=255)
                                                               ccGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAACATCGTCGCCGCCGGTTTGGCCGAGCGTTGCGAGATCCAGGTGTCGTACg                                                                      <  1:14805261/137‑1 (MQ=255)
                                                               ccGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAACATCGTCGCCGCCGGTTTGGCCGAGCGTTGCGAGATCCAGGTGTCGTACg                                                                      <  2:215674/137‑1 (MQ=255)
                                                               ccGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAACATCGTCGCCGCCGGTTTGGCCGAGCGTTGCGAGATCCAGGTGTCGTACg                                                                      <  2:14519407/137‑1 (MQ=255)
                                                               ccGTCACGGGGGCGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAACATCGTCGCCGCCGGTTTGGCCGAGCGTTGCGAGATCCAGGTGTCGTACg                                                                      <  2:9299185/137‑1 (MQ=255)
                                                                cGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAACATCGTCGCCGCCGGTTTGGCCGAGCGTTGCGAGATCCAGGTGTCGTACGc                                                                     <  1:10022919/137‑1 (MQ=255)
                                                                     cGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAACATCGTCGCCGCCGGTTTGGCCGAGCGTTGCGAGATCCAGGTGTCGTACGCCATTg                                                                >  2:12541320/1‑137 (MQ=255)
                                                                     cGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCATACGCCGGCCACTACGTGGCCAAGAACATCGTCGCCGCCGGTTTGGCCGAGCGTTGCGAGATCCAGGTGTCGTACGCCATTg                                                                >  2:2621776/1‑137 (MQ=255)
                                                                        gggTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAACATCGTCGCCGCCGGTTTGGCCGAGCGTTGCGAGATCCAGGTGTCGTACGCCATTGGCg                                                             <  1:16826320/137‑1 (MQ=255)
                                                                           tGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAACATCGTCGCCGCCGGTTTGGCCGAGCGTTGCGAGATCCAGGTGTCGTACGCCATTGGCGTgg                                                          >  2:5230808/1‑137 (MQ=255)
                                                                               gcgTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAACATCGTCGCCGCCGGTTTGGCCGAGCGTTGCGAGATCCAGGTGTCGTACGCCATTGGCGTGGCCCa                                                      <  1:5868811/137‑1 (MQ=255)
                                                                               gcgTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAACATCGTCGCCGCCGGTTTGGCCGAGCGTTGCGAGATCCAGGTGTCGTACGCCATTGGCGTGGCCCa                                                      <  1:8089362/137‑1 (MQ=255)
                                                                               gcgTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAACATCGTCGCCGCCGGTTTGGCCGAGCGTTGCGAGATCCAGGTGTCGTACGCCATTGGCGTGGCCCa                                                      <  2:10561183/137‑1 (MQ=255)
                                                                               gcgTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAACATCGTCGCCGCCGGTTTGGCCGAGCGTTGCGAGATCCAGGTGTCGTACGCCATTGGCGTGGCCCa                                                      <  1:172499/137‑1 (MQ=255)
                                                                                cgTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAACATCGTCGCCGCCGGTTTGGCCGAGCGTTGCGAGATCCAGGTGTCGTACGCCATTGGCGTGGCCCAg                                                     >  1:12803065/1‑137 (MQ=255)
                                                                                cgTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAACATCGTCGCCGCCGGTTTGGCCGAGCGTTGCGAGATCCAGGTGTCGTACGCCATTGGCGTGGCCCAg                                                     >  2:707738/1‑137 (MQ=255)
                                                                                 gTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAACATCGTCGCCGCCGGTTTGGCCGAGCGTTGCGAGATCCAGGTGTCGTACGCCATTGGCGTGGCCCAGc                                                    >  1:3584180/1‑137 (MQ=255)
                                                                                   tctcCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAACATCGTCGCCGCCGGTTTGGCCGAGCGTTGCGAGATCCAGGTGTCGTACGCCATTGGCGTGGCCCAGCCg                                                  >  2:11314230/1‑137 (MQ=255)
                                                                                    ctcCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAACATCGTCGCCGCCGGTTTGGCCGAGCGTTGCGAGATCCAGGTGTCGTACGCCATTGGCGTGGCCCAGCCGa                                                 <  1:16892726/137‑1 (MQ=255)
                                                                                     tcCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACTCCGGCCACTAAGTGGCCAAGAACATCGTCGCCGCCGGTTTGGCCGCGCGTTGCGAGATCCAGGCGTCGTACGCCATTGGCGTGGCCCAGCCGAc                                                >  2:12324001/1‑137 (MQ=255)
                                                                                     tcCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAACATCGTCGCCGCCGGTTTGGCCGAGCGTTGCGAGATCCAGGTGTCGTACGCCATTGGCGTGGCCCAGCCGAc                                                >  2:1755134/1‑137 (MQ=255)
                                                                                     tcCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAACATCGTCGCCGCCGGTTTGGCCGAGCGTTGCGAGATCCAGGTGTCGTACGCCATTGGCGTGGCCCAGCCGAc                                                >  1:6942276/1‑137 (MQ=255)
                                                                                       cGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAACATCGTCGCCGCCGGTTTGGCCGAGCGTTGCGAGATCCAGGTGTCGTACGCCATTGGCGTGGCCCAGCCGACCt                                              >  1:2439089/1‑137 (MQ=255)
                                                                                         gCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAACATCGTCGCCGCCGGTTTGGCCGAGCGTTGCGAGATCCAGGTGTCGTACGCCATTGGCGTGGCCCAGCCGACCTcc                                            <  1:9519169/137‑1 (MQ=255)
                                                                                          cAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAACATCGTCGCCGCCGGTTTGGCCGAGCGTTGCGAGATCCAGGTGTCGTACGCCATTGGCGTGGCCCAGCCGACCTCCa                                           <  2:11827416/137‑1 (MQ=255)
                                                                                          cAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAACATCGTCGCCGCCGGTTTGGCCGAGCGTTGCGAGATCCAGGTGTCGTACGCCATTGGCGTGGCCCAGCCGACCTCCa                                           <  2:5827198/137‑1 (MQ=255)
                                                                                          cAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAACATCGTCGCCGCCGGTTTGGCCGAGCGTTGCGAGATCCAGGTGTCGTACGCCATTGGCGTGGCCCAGCCGACCTCCa                                           <  1:14013827/137‑1 (MQ=255)
                                                                                          cAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAACATCGGCGCCGCCGGTTTGGCCGAGCGTTGCGAGATCCAGGTGTCGTACGCCATTGGCGTGGCCCAGCCGACCTCCa                                           <  1:14604732/137‑1 (MQ=255)
                                                                                                    tCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAACATCGTCGCCGCCGGTTTGGCCGAGCGTTGCGAGATCCAGGTGTCGTACGCCATTGGCGTGGCCCAGCCGACCTCCATCTCGATCaa                                 >  1:15120827/1‑137 (MQ=255)
                                                                                                       aaaatCGACCGTTCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAAGATCGTCGCCGCCGGTTTGGCCGAGCGCTGCGAGATCCAGGTGTCGTGCGCCATTGGCGTGGCCCAGCCGACCGCCATCTCGATCAacac                              <  2:16995156/133‑1 (MQ=255)
                                                                                                       aaGGTCGACCGTTCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAACATCGTCGCCGCCGGTTTGGCCGAGCGTTGCGAGATCCAGGTGTCGTACGCCATTGGCGTGGCCCAGCCGACCTCCATCTCGATCAacac                              <  2:6609005/137‑1 (MQ=255)
                                                                                                       aaGGTCGACCGTTCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAACATCGTCGCCGCCGGTTTGGCCGAGCGTTGCGAGATCCAGGTGTCGTACGCCATTGGCGTGGCCCAGCCGACCTCCATCTCGATCAacac                              >  1:13713559/1‑137 (MQ=255)
                                                                                                           tCGACCGTTCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAACATCGTCGCCGCCGGTTTGGCCGAGCGTTGCGAGATCCAGGTGTCGTACGCCATTGGCGTGGCCCAGCCGACCTCCATCTCGATCAACACCTTc                          <  1:2571135/137‑1 (MQ=255)
                                                                                                           tCGACCGTTCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAACATCGTCGCCGCCGGTTTGGCCGAGCGTTGCGAGATCCAGGTGTCGTACGCCATTGGCGTGGCCCAGCCGACCTCCATCTCGATCAACACCTTc                          <  1:5137615/137‑1 (MQ=255)
                                                                                                           tCGACCGTTCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAACATCGTCGCCGCCGGTTTGGCCGAGCGTTGCGAGATCCAGGTGTCGTACGCCATTGGCGTGGCCCAGCCGACCTCCATCTCGATCAACACCTTc                          <  1:10910509/137‑1 (MQ=255)
                                                                                                                 gTTCCGCTGTCTACGCCGGCCACTACGTGGCCAAGAAAATCGTCGCCGCCGGTTTGGCCGAGCGTTGCGAGATCCAGGTGTCGTACGCCATTGGCGTGGCCCAGCCGACCTCCATCTCGATCAACACCTTCGGTAcc                    >  2:15570934/1‑137 (MQ=255)
                                                                                                                 gTTCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAACATCGTCGCCGCCGGTTTGGCCGAGCGTTGCg                                                                                        >  1:326100/1‑69 (MQ=255)
                                                                                                                 gTTCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAACATCGTCGCCGCCGGTTTGGCCGAGCGTTGCGAGATCCAGGTGTCGTACGCCATTGGCGTGGCCCAGCCGACCTCCATCTCGATCAACACCTTCGGTAcc                    >  1:2857039/1‑137 (MQ=255)
                                                                                                                 gTTCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAACATCGTCGCCGCCGGTTTGGCCGAGCGTTGCGAGATCCAGGTGTCGTACGCCATTGGCGTGGCCCAGCCGACCTCCATCTCGATCAACACCTTCGGTAcc                    >  2:487957/1‑137 (MQ=255)
                                                                                                                 gTTCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAACATCGTCGCCGCCGGTTTGGCCGAGCGTTGCGAGATCCAGGTGTCGTACGCCATTGGCGTGGCCCAGCCGACCTCCATCTCGATCAACACCTTCGGTAcc                    >  2:16248310/1‑137 (MQ=255)
                                                                                                                 gTTCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAACATCGTCGCCGCCGGTTTGGCCGAGCGTTGCGAGATCCAGGTGTCGTACGCCATTGGCGTGGCCCAGCCGACCTCCATCTCGATCAACACCTTCGGTAcc                    >  1:4489336/1‑137 (MQ=255)
                                                                                                                 gTTCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAACATCGTCGCCGCCGGTTTGGCCGAGCGTTGCGAGATCCAGGTGTCGTACGCCATTGGCGTGGCCCAGCCGACCTCCATCTCGATCAACACCTTCGGTAcc                    >  1:9000250/1‑137 (MQ=255)
                                                                                                                  ttCCGCTGCCTACGCCGGCCACTACGTGGCCAAGAACATCGTCGCCGCCGGTTTGGCCGAGCGTTGCGAGATCCAGGTGTCGTACGCCATTGGCGTGGCCCAGCCGACCTCCATCTCGATCAACACCTTCGGTACCg                   <  1:16590884/137‑1 (MQ=255)
                                                                                                                     cGCTGCCTACGCCGGCCACTACGTGGCCAAGAACATCGTCGCCGCCGGTTTGGCCGAGCGTTGCGAGATCCAGGTGTCGTACGCCATTGGCGTGGCCCAGCCGACCTCCATCTCGATCAACACCTTCGGTACCGGCa                >  1:16693827/1‑137 (MQ=255)
                                                                                                                          ccTACGCCGGCCACTACGTGGCCAAGAACATCGTCGCCGCCGGTTTGGCCGAGCGTTGCGAGATCCAGGTGTCGTACGCCATTGGCGTGGCCCAGCCGACCTCCATCTCGATCAACACCTTCGGTACCGGCAAGGtc           <  2:7833956/137‑1 (MQ=255)
                                                                                                                          ccTACGCCGGCCACTACGTGGCCAAGAACATCGTCGCCGCCGGTTTGGCCGAGCGTTGCGAGATCCAGGTGTCGTACGCCATTGGCGTGGCCCAGCCGACCTCCATCTCGATCAACACCTTCGGTACCGGCAAGGtc           <  1:149658/137‑1 (MQ=255)
                                                                                                                          ccTACGCCGGCCACTACGTGGCCAAGAACATCGTCGCCGCCGGTTTGGCCGAGCGTTGCGAGATCCAGGTGTCGTACGCCATTGGCGTGGCCCAGCCGACCTCCATCTCGATCAACACCTTCGGTACCGGCAAGGtc           <  2:17298035/137‑1 (MQ=255)
                                                                                                                               gccggccaCTACGTGGCCAAGAACATCGTCGCCGCCGGTTTGGCCGAGCGTTGCGAGATCCAGGTGTCGTACGCCATTGGCGTGGCCCAGCCGACCTCCATCTCGATCAACACCTTCGGTACCGGCAAGGTCTCcga      <  2:3752886/137‑1 (MQ=255)
                                                                                                                                ccggccaCTACGTGGCCAAGAACATCGTCGCCGCCGGTTTGGCCGAGCGTTGCGAGATCCAGGTGTCGTACGCCATTGGCGTGGCCCAGCCGACCTCCATCTCGATCAACACCTTCGGTACCGGCAAGGTCTCcgac     >  1:8325592/1‑137 (MQ=255)
                                                                                                                                ccggccaCTACGTGGCCAAGAACATCGTCGCCGCCGGTTTGGCCGAGCGTTGCGAGATCCAGGTGTCGTACGCCATTGGCGTGGCCCAGCCGACCTCCATCTCGATCAACACCTTCGGTACCGGCAAGGTCTCcgac     >  2:12949463/1‑137 (MQ=255)
                                                                                                                                  ggccaCTACGTGGCCAAGAACATCGTCGCCGCCGGTTTGGCCGAGCGTTGCGAGATCCAGGTGTCGTACGCCATTGGCGTGGCCCAGCCGACCTCCATCTCGATCAACACCTTCGGTACCGGCAAGGTCTCcgacga   <  1:4037328/137‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                     
CCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCGCTACGTGGCCAAGAACATCGTCGCCGCCGGTTTGGCCGAGCGTTGCGAGATCCAGGTGTCGTACGCCATTGGCGTGGCCCAGCCGACCTCCATCTCGATCAACACCTTCGGTACCGGCAAGGTCTCCGACGAC  >  AE015451_phaM1/5659664‑5659931

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