Predicted mutation
evidence seq id position mutation annotation gene description
RA AE015451_phaM1 5,659,665 C→A P237T (CCG→ACG)  metK CDS (methionine adenosyltransferase) →

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451_phaM15,659,6650CA99.9% 3209.8 / NA 879P237T (CCG→ACG) metK CDS (methionine adenosyltransferase)
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (431/447);  total (432/447)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 879 total aligned reads displayed.

CCGGAGGTTTCGCAGAAAGACCTGCAAGAAGCCGTGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTA‑‑‑CTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCG  >  AE015451_phaM1/5659533‑5659798
                                                                                                                                       |                                                                                                                                     
  ggAGGTTTCGCAGAAAGACCTGCAAGAAGCCGTGATGGAGCTGATCGTCAAGCACGCGCTGCCGGCCGAACTA‑‑‑CTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCACggtgg                                                                                                                                 >  1:13726008/1‑137 (MQ=255)
  ggAGGTTTCGCAGAAAGACCTGCAAGAAGCCGTGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTA‑‑‑CTGCACCCAGGCACCCAGTCCCACATCCACCCGACCGTCCCCTTCATCATCGGTGGCACgttgg                                                                                                                                 >  1:11455545/1‑137 (MQ=255)
  ggAGGTTTCGCAGAAAGACCTGCAAGAAGCCGTGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTA‑‑‑CTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCACggtgg                                                                                                                                 >  2:2124001/1‑137 (MQ=255)
  ggAGGTTTCGCAGAAAGACCTGCAAGAAGCCGTGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTA‑‑‑CTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCACggtgg                                                                                                                                 >  2:4773803/1‑137 (MQ=255)
  ggAGGTTTCGCAGAAAGACCTGCAAGAAGCCGTGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTA‑‑‑CTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCACggtgg                                                                                                                                 >  2:16856875/1‑137 (MQ=255)
  ggAGGTTTCGCAGAAAGACCTGCAAGAAGCCGTGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTA‑‑‑CTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCACggtgg                                                                                                                                 >  2:1561606/1‑137 (MQ=255)
  ggAGGTTTCGCAGAAAGACCTGCAAGAAGCCGTGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTA‑‑‑CTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCACggtgg                                                                                                                                 >  2:416854/1‑137 (MQ=255)
  ggAGGTTTCGCAGAAAGACCTGCAAGAAGCCGTGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTA‑‑‑CTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCACggtgg                                                                                                                                 >  1:5158324/1‑137 (MQ=255)
  ggAGGTTTCGCAGAAAGACCTGCAAGAAGCCGTGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTA‑‑‑CTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCACggtgg                                                                                                                                 >  2:3346698/1‑137 (MQ=255)
  ggAGGTTTCGCAGAAAGACCTGCAAGAAGCCGTGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTA‑‑‑CTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCACggtgg                                                                                                                                 >  1:13016861/1‑137 (MQ=255)
  ggAGGTTTCGCAGAAAGACCTGCAAGAAGCCGTGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTA‑‑‑CTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCACggtgg                                                                                                                                 >  1:3251032/1‑137 (MQ=255)
  ggAGGTTTCGCAGAAAGACCTGCAAGAAGCCGTGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTA‑‑‑CTGCACAAAGGCACCCAGTACCACATCAAACCGACCGGCAACTTCATCATCGGTGGCACggtgg                                                                                                                                 >  1:9752074/1‑137 (MQ=255)
   gAGGTTTCGCAGAAAGACCTGCAAGAAGCCGTGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTA‑‑‑CTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCACggtggg                                                                                                                                >  1:12327195/1‑137 (MQ=255)
      gTTTCGCAGAAAGACCTGCAAGAAGCCGTGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTA‑‑‑CTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCACGGTGGGTGa                                                                                                                             <  2:3812745/137‑1 (MQ=255)
       tttCGCAGAAAGACCTGCAAGAAGCCGTGATGGAGCTGATCGTCAAGCACACGCTGCCGACCAAACTA‑‑‑CTGCACCAAAGCACCCAGTACCCACGCAACCCGCCCGGAAAATTCATCATGGGTGGCACGGTGGGTGAc                                                                                                                            >  2:19484154/1‑137 (MQ=255)
            cAGAAAGACCTGCAAGAAGCCGTGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTA‑‑‑CTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCATGGTGGGTGACTGCgg                                                                                                                       >  1:15331308/1‑137 (MQ=255)
            cAGAAAGACCTGCAAGAAGCCGTGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTA‑‑‑CTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCACGGTGGGTGACTGCgg                                                                                                                       >  2:1415303/1‑137 (MQ=255)
             aGAAAGACCTGCAAGAAGCCGTGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTA‑‑‑CTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCACGGTGGGTGACTGCGGc                                                                                                                      >  2:9620588/1‑137 (MQ=255)
             aGAAAGACCTGCAAGAAGCCGTGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTA‑‑‑CTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCACGGTGGGTGACTGCGGc                                                                                                                      <  1:6517882/137‑1 (MQ=255)
              gAAAGACCTGCAAGAAGCCGTGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTA‑‑‑CTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCACGGTGGGTGACTGCGGcc                                                                                                                     <  1:14465097/137‑1 (MQ=255)
               aaaGACCTACAAGAAGCTGTGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTA‑‑‑CTGCACAAATGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCACGGTGGGTGACTGCGGCCt                                                                                                                    <  2:9838633/137‑1 (MQ=255)
                   aCCTGCAAGAAGCCGTGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTA‑‑‑CTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCACGGTGGGTGACTGCGGCCTGACt                                                                                                                <  2:2372911/137‑1 (MQ=255)
                   aCCTGCAAGAAGCCGTGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTA‑‑‑CTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCACGGTGGGTGACTGCGGCCTGACt                                                                                                                <  2:5111494/137‑1 (MQ=255)
                    ccTGCAAGAAGCCGTGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTA‑‑‑CTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCACGGTGGGTGACTGCGGCCTGACTg                                                                                                               >  2:17984105/1‑137 (MQ=255)
                                cGTGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTA‑‑‑CTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCACGGTGGGTGACTGCGGCCTGACTGGCCGCAAGatca                                                                                                   <  1:9280809/137‑1 (MQ=255)
                                  tGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTA‑‑‑CTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCACGGTGGGTGACTGCGGCCTGACTGGCCGCAAGatcatc                                                                                                 <  2:7765761/137‑1 (MQ=255)
                                  tGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTA‑‑‑CTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCACGGTGGGTGACTGCGGCCTGACTGGCCGCAAGatcatc                                                                                                 <  2:12752351/137‑1 (MQ=255)
                                  tGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTA‑‑‑CTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCACGGTGGGTGACTGCGGCCTGACTGGCCGCAAGatcatc                                                                                                 <  1:15009573/137‑1 (MQ=255)
                                    aTGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTA‑‑‑CTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCACGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCAtcgt                                                                                               <  1:6360088/137‑1 (MQ=255)
                                            gATCGTCAAGCACACGCTGCCGGCCGAACTA‑‑‑CTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCACGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCt                                                                                       <  2:7193203/137‑1 (MQ=255)
                                            gATCGTCAAGCACACGCTGCCGGCCGAACTA‑‑‑CTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCACGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCt                                                                                       <  2:11003149/137‑1 (MQ=255)
                                               cGTCAAGCACACGCTGCCGGCCGAACTACTGCTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCACGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCt                                                                                       <  1:5741434/137‑1 (MQ=255)
                                               cGTCAAGCACACGCTGCCGGCCGAACTA‑‑‑CTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCACGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACg                                                                                    >  2:6100288/1‑137 (MQ=255)
                                                gTCAAGCACACGCTGCCGGCCGAACTA‑‑‑CTGCACAACGGCACCCAGTACCACACCAACCCGACCGGCAACCTCATCATCGGTGGCACGGTGGGTGACTGCGGACTGACTGGCCGCGAGATCATCGTCGACTCATACgg                                                                                   >  2:9214877/1‑137 (MQ=255)
                                                gTCAAGCACACGCTGCCGGCCGAACTA‑‑‑CTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCACGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACgg                                                                                   >  1:12840483/1‑137 (MQ=255)
                                                 tCAAGCACACGCTGCCGGCCGAACTA‑‑‑CTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCACGGTGGGTGACTGCGGCCTGACTGGCCGCAAGatcatc                                                                                                 <  2:18442666/122‑1 (MQ=255)
                                                 tCAAGCACACGCTGCCGGCCGAACTA‑‑‑CTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCACGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGt                                                                                  >  1:10011999/1‑137 (MQ=255)
                                                  cAAGCACACGCTGCCGGCCGAACTA‑‑‑CTGCACAAAGGCACCCAGTACCACATCAACCCGACCAGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTg                                                                                 >  2:6247657/1‑137 (MQ=255)
                                                      cacaCGCTGCCGGCCGAACTA‑‑‑CTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCACGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCAt                                                                             <  1:6927298/137‑1 (MQ=255)
                                                      cacaCGCTGCCGGCCGAACTA‑‑‑CTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCACGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCAt                                                                             <  2:12376884/137‑1 (MQ=255)
                                                       acacGCTGCCGGCCGAACTA‑‑‑CTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCACGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATg                                                                            <  1:550469/137‑1 (MQ=255)
                                                       acacGCTGCCGGCCGAACTA‑‑‑CTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCACGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATg                                                                            <  1:8171733/137‑1 (MQ=255)
                                                       acacGCTGCCGGCCGAACTA‑‑‑CTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCACGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATg                                                                            <  2:8666008/137‑1 (MQ=255)
                                                       acacGCTGCCGGCCGAACTA‑‑‑CTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCACGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATg                                                                            <  1:16100522/137‑1 (MQ=255)
                                                            cTGCCGGCCGAACTA‑‑‑CTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCACGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCg                                                                       <  2:6528779/137‑1 (MQ=255)
                                                            cTGCCGGCCGAACTA‑‑‑CTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCACGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCg                                                                       <  1:1194642/137‑1 (MQ=255)
                                                            cTGCCGGCCGAACTA‑‑‑CTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCACGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCg                                                                       <  2:12005917/137‑1 (MQ=255)
                                                              gccggccgAACTA‑‑‑CTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCACGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTc                                                                     >  2:18741415/1‑137 (MQ=255)
                                                                  gccgAACTA‑‑‑CTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCACGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACgg                                                                 >  2:738061/1‑137 (MQ=255)
                                                                      aaCTA‑‑‑CTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCACGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGt                                                             >  2:1024959/1‑137 (MQ=255)
                                                                       acta‑‑‑ctGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCACGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTg                                                            <  1:7820534/137‑1 (MQ=255)
                                                                         ta‑‑‑ctGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCACGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGgc                                                          <  1:8500612/137‑1 (MQ=255)
                                                                                 cacaAATGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCACGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGtt                                                     <  1:1055198/137‑1 (MQ=255)
                                                                                 cacaAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCACGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGtt                                                     <  1:4686234/137‑1 (MQ=255)
                                                                                   caAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCACGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCAAGGGGGTGGCGCGTtct                                                   >  1:11667100/1‑137 (MQ=255)
                                                                                      aGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCACGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCg                                                >  1:17786641/1‑137 (MQ=255)
                                                                                      aGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCACGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCg                                                <  1:18111410/137‑1 (MQ=255)
                                                                                         cACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCACGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGTATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCa                                             >  1:16938994/1‑137 (MQ=255)
                                                                                         cACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCACGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCa                                             >  1:18629360/1‑137 (MQ=255)
                                                                                          aCCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCACGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCGACGGTGACATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCaa                                            <  1:17556431/137‑1 (MQ=255)
                                                                                           cccAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCACGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAg                                           <  1:2927637/137‑1 (MQ=255)
                                                                                           cccAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCACGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAg                                           <  1:16497404/137‑1 (MQ=255)
                                                                                           cccAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCACGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAg                                           <  2:9074700/137‑1 (MQ=255)
                                                                                           cccAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCACGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAg                                           <  1:2004905/137‑1 (MQ=255)
                                                                                           cccAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCACGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAg                                           <  1:2424567/137‑1 (MQ=255)
                                                                                           cccAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCACGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAg                                           <  2:2756978/137‑1 (MQ=255)
                                                                                           cccAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCACGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAg                                           <  2:13239211/137‑1 (MQ=255)
                                                                                           cccAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCACGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAg                                           <  2:13711149/137‑1 (MQ=255)
                                                                                           cccAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCACGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAg                                           >  1:9562664/1‑137 (MQ=255)
                                                                                           cccAGTACCACATCAACCCGACCGGCAACTTCATCATCGGCGGCACGGTGTGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACTGGGGTGGCGCGGGCTCCGGCAAg                                           <  2:15195223/137‑1 (MQ=255)
                                                                                           cccAGTACCACATAAATCCGACCGGCAACTGCATTATCTGTGGCACGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAg                                           <  1:1674612/137‑1 (MQ=255)
                                                                                             cAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCACGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCAGCAAGGa                                         >  2:17145248/1‑137 (MQ=255)
                                                                                              aGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCACGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGtt                                                     <  2:17374525/124‑1 (MQ=255)
                                                                                              aGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCACGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGAc                                        >  2:9472670/1‑137 (MQ=255)
                                                                                              aGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCACGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGAc                                        >  2:8361050/1‑137 (MQ=255)
                                                                                              aGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCACGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGAc                                        >  1:6070678/1‑137 (MQ=255)
                                                                                              aGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCACGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGAc                                        >  2:18887271/1‑137 (MQ=255)
                                                                                                 aCCACATCAACCCGACCGGCAACTTCATCATCGGTGGCACGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCg                                     >  2:14857086/1‑137 (MQ=255)
                                                                                                        cAACCCGACCGGCAACTTCATCATCGGTGGCACGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGt                                                      <  2:13140709/113‑1 (MQ=255)
                                                                                                        cAACCCGACCGGCAACTTCATCATCGGTGGCACGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAgg                              >  2:10707127/1‑137 (MQ=255)
                                                                                                         aaCCCGACCGGCAACTTCATCATCGGTGGCACGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATg                                                                            >  1:3598436/1‑90 (MQ=255)
                                                                                                         aaCCCGACCGGCAACTTCATCATCGGTGGCACGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGt                             >  2:14448843/1‑137 (MQ=255)
                                                                                                         aaCCCGACCGGCAACTTCATCATCGGTGGCACGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGt                             >  2:6857577/1‑137 (MQ=255)
                                                                                                         aaCCCGACCGGCAACTTCATCATCGGTGGCACGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCAAGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGt                             >  2:17616363/1‑137 (MQ=255)
                                                                                                           cccGACCGGCAACTTCATCATCGGTGGCACGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCg                           >  1:4286774/1‑137 (MQ=255)
                                                                                                           cccGACCGGCAACTTCATCATCGGTGGCACGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCg                           >  1:3964956/1‑137 (MQ=255)
                                                                                                             cGACCGGCAACTTCATCATCGGTGGCACGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGAc                         >  2:11685459/1‑137 (MQ=255)
                                                                                                             cGACCGGCAACTTCATCATCGGTGGCACGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGAc                         >  1:17339044/1‑137 (MQ=255)
                                                                                                             cGACCGGCAACTTCATCATCGGTGGCACGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCACGTCAAGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGAc                         >  1:14076929/1‑137 (MQ=255)
                                                                                                               aCCGGCAACTTCATCATCGGTGGCACGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCg                       >  2:16372617/1‑137 (MQ=255)
                                                                                                                ccGGCAACTTCATCATCGGTGGCACGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGt                      >  2:4641402/1‑137 (MQ=255)
                                                                                                                ccGGCAACTTCATCATCGGTGGCACGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGt                      <  1:7409142/137‑1 (MQ=255)
                                                                                                                 cGGCAACTTCATCATCGGTGGCACGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGtt                     <  2:11261123/137‑1 (MQ=255)
                                                                                                                 cGGCAACTTCATCATCGGTGGCACGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGtt                     <  1:6703011/137‑1 (MQ=255)
                                                                                                                 cGGCAACTTCATCATCGGTGGCACGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGtt                     <  2:3691878/137‑1 (MQ=255)
                                                                                                                 cGGCAACTTCATCATCGGTGGCACGGNGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGtt                     <  1:19283356/137‑1 (MQ=255)
                                                                                                                   gCAACTTCATCATCGGTGGCACGGTGGTTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTcc                   >  1:2599467/1‑137 (MQ=255)
                                                                                                                   gCAACTTCATCATCGGTGGCACGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCACGGACCCGTCCAAGGTCGACCGTTcc                   >  2:5564441/1‑137 (MQ=255)
                                                                                                                               tcGTTGGCACGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACgcc       <  1:9027210/137‑1 (MQ=255)
                                                                                                                               tcGTTGGCACGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACgcc       <  1:18945218/137‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                     
CCGGAGGTTTCGCAGAAAGACCTGCAAGAAGCCGTGATGGAGCTGATCGTCAAGCACACGCTGCCGGCCGAACTA‑‑‑CTGCACAAAGGCACCCAGTACCACATCAACCCGACCGGCAACTTCATCATCGGTGGCCCGGTGGGTGACTGCGGCCTGACTGGCCGCAAGATCATCGTCGACTCCTACGGTGGCATGGCCCGTCACGGGGGTGGCGCGTTCTCCGGCAAGGACCCGTCCAAGGTCGACCGTTCCGCTGCCTACGCCGGCCG  >  AE015451_phaM1/5659533‑5659798

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: