Predicted mutation
evidence seq id position mutation annotation gene description
RA AE015451_phaM1 539,294 A→C T558P (ACC→CCC)  rpoB CDS (DNA‑directed RNA polymerase subunit beta) →

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451_phaM1539,2940AC99.8% 2165.3 / NA 536T558P (ACC→CCC) rpoB CDS (DNA‑directed RNA polymerase subunit beta)
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base C (247/288);  minor base T (1/0);  total (248/288)
Fisher's exact test for biased strand distribution p-value = 4.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 560 total aligned reads displayed.

CAGCCAGCTGTCCCAGTTCATGGACCAGAACAACCCTCTCTCGGAGATTACCCACAAGCGCCGCGTCTCCGCACTCGGCCCTGGCGGTCTGACCCGTGAGCGTGCCGGCTTCGAAGTCCGTGACGTACACCCGACCCACTATGGCCGTGTGTGCCCGATCGAGACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCCCTGGCGGCCTATGCCCGCACCAACCAGTACGGCTTCCTGGAAAGCCCGTACCGCGTGGTGAAGGAAG  >  AE015451_phaM1/539161‑539426
                                                                                                                                     |                                                                                                                                    
cagccagcTGTCCCAGTTCATGGACCAGAACAACCCTCTCTCGTAGATTGCCCACAAGCGCCGCGTCTCCGCACTCGGCCCTGGCGGTCTGACCCGTGAGCGTGCCGGCTTCGAAGTCCGTGACGTACACCCGcccc                                                                                                                                   <  2:2523919/137‑1 (MQ=255)
cagccagcTGTCCCAGTTCATGGACCAGAACAACCCTCTCTCGGAGATTACCCACAAGCGCCGCGTCTCCGCACTCGGCCCTGGCGGTCTGACCCGTGAGCGTGCCGGCTTCGAAGTCCGTGACGTACACCCGcccc                                                                                                                                   >  1:10326190/1‑137 (MQ=255)
cagccagcTGTCCCAGTTCATGGACCAGAACAACCCTCTCTCGGAGATTACCCACAAGCGCCGCGTCTCCGCACTCGGCCCTGGCGGTCTGACCCGTGAGCGTGCCGGCTTCGAAGTCCGTGACGTACACCCGcccc                                                                                                                                   <  1:9204318/137‑1 (MQ=255)
cagccagcTGTCCCAGTTCATGGACCAGAACAACCCTCTCTCGGAGATTACCCACAAGCGCCGCGTCTCCGCACTCGGCCCTGGCGGTCTGACCCGTGAGCGTGCCGGCTTCGAAGTCCGTGACGTACACCCGcccc                                                                                                                                   <  1:2088002/137‑1 (MQ=255)
cagccagcTGTCCCAGTTCATGGACCAGAACAACCCTCTCTCGGAGATTACCCACAAGCGCCGCGTCTCCGCACTCGGCCCTGGCGGTCTGACCCGTGAGCGTGCCGGCTTCGAAGTCCGTGACGTACACCCGcccc                                                                                                                                   <  1:14432606/137‑1 (MQ=255)
    cagcTGTCCCAGTTCATGGACCAGAACAACCCTCTCTCGGAGATTACCCACAAGCGCCGCGTCTCCGCACTCGGCCCTGGCGGTCTGACCCGTGAGCGTGCCGGCTTCGAAGTCCGTGACGTACACCCGCCCCACTa                                                                                                                               <  2:15044530/137‑1 (MQ=255)
    cagcTGTCCCAGTTCATGGACCAGAACAACCCTCTCTCGGAGATTACCCACAAGCGCCGCGTCTCCGCACTCGGCCCTGGCGGTCTGACCCGTGAGCGTGCCGGCTTCGAAGTCCGTGACGTACACCCGCCCCACTa                                                                                                                               <  2:13496211/137‑1 (MQ=255)
    cagcTGTCCCAGTTCATGGACCAGAACAACCCTCTCTCGGAGATTACCCACAAGCGCCGCGTCTCCGCACTCGGCCCTGGCGGTCTGACCCGTGAGCGTGCCGGCTTCGAAGTCCGTGACGTACACCCGCCCCACTa                                                                                                                               <  2:10451129/137‑1 (MQ=255)
    cagcTGTCCCAGTTCATGGACCAGAACAACCCTCTCTCGGAGATTACCCACAAGCGCCGCGTCTCCGCACTCGGCCCTGGCGGTCTGACCCGTGAGCGTGCCGGCTTCGAAGTCCGTGACGTACACCCGCCCCACTa                                                                                                                               <  1:1266720/137‑1 (MQ=255)
    cagcTGTCCCAGTTCATGGACCAGAACAACCCTCTCTCGGAGATTACCCACAAGCGCCGCGTCTCCGCACTCGGCCCTGGCGGTCTGACCCGTGAGCGTGCCGGCTTCGAAGTCCGTGACGTACACCCGCCCCACTa                                                                                                                               <  1:16403528/137‑1 (MQ=255)
     agcTGTCCCAGTTCATGGACCAGAACAACCCTCTCTCGGAGATTACCCACAAGCGCCGCGTCTCCGCACTCGGCCCTGGCGGTCTGACCCGTGAGCGTGCCGGCTTCGAAGTCCGTGACGTACACCCGCCCCACTAt                                                                                                                              >  2:3620982/1‑137 (MQ=255)
          tCCCAGTTCATGGACCAGAACAACCCTCTCTCGGAGATTACCCACAAGCGCCGCGTCTCCGCACTCGGCCCTGGCGGTCTGACCCGTGAGCGTGCCGGCTTCGAAGTCCGTGACGTACACCCGCCCCACTATGGCCg                                                                                                                         >  2:11161656/1‑137 (MQ=255)
          tCCCAGTTCATGGACCAGAACAACCCTCTCTCGGAGATTACCCACAAGCGCCGCGTCTCCGCACTCGGCCCTGGCGGTCTGACCCGTGAGCGTGCCGGCTTCGAAGTCCGTGACGTACACCCGCCCCACTATGGCCg                                                                                                                         >  1:14906125/1‑137 (MQ=255)
               gTTCATGGACCAGAACAACCCTCTCTCGGAGATTACCCACAAGCGCCGCGTCTCCGCACTCGGCCCTGGCGGTCTGACCCGTGAGCGTGCCGGCTTCGAAGTCCGTGACGTACACCCGCCCCACTATGGCCgtgtgt                                                                                                                    <  2:1690312/137‑1 (MQ=255)
                 tCATGGACCAGAACAACCCTCTCTCGGAGATTACCCACAAGCGCCGCGTCTCCGCACTCGGCCCTGGCGGTCTGACCCGTGAGCGTGCCGGCTTCGAAGTCCGTGACGTACACCCGCCCCCCTATGGCCGTGTGTGc                                                                                                                  >  2:7378258/1‑137 (MQ=255)
                   aTGGACCAGAACAACCCTCTCTCGGAGATTACCCACAAGCGCCGCGTCTCCGCACTCGGCCCTGGCGGTCTGACCCGTGAGCGTGCCGGCTTCGAAGTCCGTGACGTACACCCGCCCCACTATGGCCGTGTGTGccc                                                                                                                <  1:8681252/137‑1 (MQ=255)
                    tGGACCAGAACAACCCTCTCTCGGAGATTACCCACAAGCGCCGCGTCTCCGCACTCGGCCCTGGCGGTCTGACCCGTGAGCGTGCCGGCTTCGAAGTCCGTGACGTACACCCGCCCCACTATGGCCGTGTGTGCCCg                                                                                                               >  2:14642372/1‑137 (MQ=255)
                      gACCAGAACAACCCTCTCTCGGAGATTACCCACAAGCGCCGCGTCTCCGCACTCGGCCCTGGCGGTCTGACCCGTGAGCGTGCCGGCTTCGAAGTCCGTGACGTACACCCGCCCCACTATGGCCGTGTGTGCCCGAt                                                                                                             <  2:12996824/137‑1 (MQ=255)
                      gACCAGAACAACCCTCTCTCGGAGATTACCCACAAGCGCCGCGTCTCCGCACTCGGCCCTGGCGGTCTGACCCGTGAGCGTGCCGGCTTCGAAGTCCGTGACGTACACCCGCCCCACTATGGCCGTGTGTGCCCGAt                                                                                                             <  2:12606495/137‑1 (MQ=255)
                      gACCAGAACAACCCTCTCTCGGAGATTACCCACAAGCGCCGCGTCTCCGCACTCGGCCCTGGCGGTCTGACCCGTGAGCGTGCCGGCTTCGAAGTCCGTGACGTACACCCGCCCCACTATGGCCGTGTGTGCCCGAt                                                                                                             <  2:6733098/137‑1 (MQ=255)
                       aCCAGAACAACCCTCTCTCGGAGATTACCCACAAGCGCCGCGTCTCCGCACTCGGCCCTGGCGGTCTGACCCGTGAGCGTGCCGGCTTCGAAGTCCGTGACGTACACCCGCCCCACTa                                                                                                                               >  1:1130138/1‑118 (MQ=255)
                       aCCAGAACAACCCTCTCTCGGAGATTACCCACAAGCGCCGCGTCTCCGCACTCGGCCCTGGCGGTCTGACCCGTGAGCGTGCCGGCTTCGAAGTCCGTGACGTACACCCGCCCCACTATGGCCGTGTGTGCCCGATc                                                                                                            >  2:10657030/1‑137 (MQ=255)
                       aCCAGAACAACCCTCTCTCGGAGATTACCCACAAGCGCCGCGTCTCCGCACTCGGCCCTGGCGGTCTGACCCGTGAGCGTGCCGGCTTCGAAGTCCGTGACGTACACCCGCCCCACTATGGCCGTGTGTGCCCGATc                                                                                                            >  2:926745/1‑137 (MQ=255)
                       aCCAGAACAACCCTCTCTCGGAGATTACCCACAAGCGCCGCGTCTCCGCACTCGGCCCTGGCGGTCTGACCCGTGAGCGTGCCGGCTTCGAAGTCCGTGACGTACACCCGCCCCACTATGGCCGTGTGTGCCCGATc                                                                                                            >  2:8975974/1‑137 (MQ=255)
                       aCCAGAACAACCCTCTCTCGGAGATTACCCACAAGCGCCGCGTCTCCGCACTCGGCCCTGGCGGTCTGACCCGTGAGCGTGCCGGCTTCGAAGTCCGTGACGTACACCCGCCCCACTATGGCCGTGTGTGCCCGATc                                                                                                            <  1:11722910/137‑1 (MQ=255)
                         cAGAACAACCCTCTCTCGGAGATTACCCACAAGCGCCGCGTCTCCGCACTCGGCCCTGGCGGTCTGACCCGTGAGCGTGCCGGCTTCGAAGTCCGTGACGTACACCCGCCCCACTATGGCCGTGTGTGCCCGATCga                                                                                                          <  2:13760061/137‑1 (MQ=255)
                             acaacCCTCTCTCGGAGATTACCCACAAGCGCCGCGTCTCCGCACTCGGCCCTGGCGGTCTGACCCGTGAGCGTGCCGGCTTCGAAGTCCGTGACGTACACCCGCCCCACTATGGCCGTGTGTGCCCGATCGAGAcc                                                                                                      >  1:10525052/1‑137 (MQ=255)
                               aacCCTCTCTCGGAGATTACCCACAAGCGCCGCGTCTCCGCACTCGGCCCTGGCGGTCTGACCCGTGAGCGTGCCGGCTTCGAAGTCCGTGACGTACACCCGCCCCACTATGGCCGTGTGTGCCCGATCGAGAcccc                                                                                                    >  2:15509248/1‑137 (MQ=255)
                               aacCCTCTCTCGGAGATTACCCACAAGCGCCGCGTCTCCGCACTCGGCCCTGGCGGTCTGACCCGTGAGCGTGCCGGCTTCGAAGTCCGTGACGTACACCCGCCCCACTATGGCCGTGTGTGCCCGATCGAGAcccc                                                                                                    <  2:15369538/137‑1 (MQ=255)
                                    tctctcGGAGATTACCCACAAGCGCCGCGTCTCCGCACTCGGCCCTGGCGGTCTGACCCGTGAGCGTGCCGGCTTCGAAGTCCGTGACGTACACCCGCCCCACTATGGCCGTGTGTGCCCGATCGAGACCCCTGAAg                                                                                               >  2:7090213/1‑137 (MQ=255)
                                    tctctcGGAGATTACCCACAAGCGCCGCGTCTCCGCACTCGGCCCTGGCGGTCTGACCCGTGAGCGTGCCGGCTTCGAAGTCCGTGACGTACACCCGCCCCACTATGGCCGTGTGTGCCCGATCGAGACCCCTGAAg                                                                                               >  1:5745486/1‑137 (MQ=255)
                                    tctctcGGAGATTACCCACAAGCGCCGCGTCTCCGCACTCGGCCCTGGCGGTCTGACCCGTGAGCGTGCCGGCTTCGAAGTCCGTGACGTACACCCGCCCCACTATGGCCGTGTGTGCCCGATCGAGACCCCTGAAg                                                                                               >  1:4760803/1‑137 (MQ=255)
                                    tctctcGGAGATTACCCACAAGCGCCGCGTCTCCGCACTCGGCCCTGGCGGTCTGACCCGTGAGCGTGCCGGCTTCGAAGTCCGTGACGTACACCCGCCCCAATATGGCCGTGTGTGCCCGATCGAGACCCCTGAAg                                                                                               >  2:12198565/1‑137 (MQ=255)
                                     ctctcGGAGATTACCCACAAGCGCCGCGTCTCCGCACTCGGCCCTGGCGGTCTGACCCGTGAGCGTGCCGGCTTCGAAGTCCGTGACGTACACCCGCCCCACTATGGCCGTGTGTGCCCGATCGAGACCCCTGAAgg                                                                                              <  1:4545760/137‑1 (MQ=255)
                                         cGGAGATTACCCACAAGCGCCGCGTCTCCGCACTCGGCCCTGGCGGTCTGACCCGTGAGCGTGTCGGCTTCGAAGTCCGTGACGTACACCCGCCCCACTATGGCCGTGTGTGCCCGATCGAGACCCCTGAAGGTCCg                                                                                          >  1:15240117/1‑137 (MQ=255)
                                         cGGAGATTACCCACAAGCGCCGCGTCTCCGCACTCGGCCCTGGCGGTCTGACCCGTGAGCGTGCCGGCTTCGAAGTCCGTGACGTACACCCGCCCCACTATGGCCGTGTGTGCCCGATCGAGACCCCTGAAGGTCCg                                                                                          >  2:16346939/1‑137 (MQ=255)
                                         cGGAGATTACCCACAAGCGCCGCGTCTCCGCACTCGGCCCTGGCGGTCTGACCCGTGAGCGTGCCGGCTTCGAAGTCCGTGACGTACACCCGCCCCACTATGGCCGTGTGTGCCCGATCGAGACCCCTGAAGGTCCg                                                                                          >  1:3703530/1‑137 (MQ=255)
                                         cGGAGATTACCCACAAGCGCCGCGTCTCCGCACTCGGCCCTGGCGGTCTGACCCGTGAGCGTGCCGGCTTCGAAGTCCGTGACGTACACCCGCCCCACTATGGCCGTGTGTGCCCGATCGAGACCCCTGAAGGTCCg                                                                                          >  2:8832779/1‑137 (MQ=255)
                                         cGGAGATTACCCACAAGCGCCGCGTCTCCGCACTCGGCCCTGGCGGTCTGACCCGTGAGCGTGCCGGCTTCGAAGTCCGTGACGTACACCCGCCCCACTATGGCCGTGTGTGCCCGATCGAGACCCCTGAAGGTCCg                                                                                          >  1:11959922/1‑137 (MQ=255)
                                          ggAGATTACCCACAAGCGCCGCGTCTCCGCACTCGGCCCTGGCGGTCTGACCCGTGAGCGTGCCGGCTTCGAAGTCCGTGACGTACACCCGCCCCACTATGGCCGTGTGTGCCCGATCGAGACCCCTGAAGGTCCGa                                                                                         <  2:1445261/137‑1 (MQ=255)
                                          ggAGATTACCCACAAGCGCCGCGTCTCCGCACTCGGCCCTGGCGGTCTGACCCGTGAGCGTGCCGGCTTCGAAGTCCGTGACGTACACCCGCCCCACTATGGCCGTGTGTGCCCGATCGAGACCCCTGAAGGTCCGa                                                                                         <  1:2597766/137‑1 (MQ=255)
                                               ttACCCACAAGCGCCGCGTCTCCGCACTCGGCCCTGGCGGTCTGACCCGTGAGCGTGCCGGCTTCGAAGTCCGTGACGTACACCCGCCCCACTATGGCCGTGTGTGCCCGATCGAGACCCCTGAAGGTCCGAACATc                                                                                    <  1:3126370/137‑1 (MQ=255)
                                                tACCCACAAGCGCCGCGTCTCCGCACTCGGCCCTGGCGGTCTGACCCGTGAGCGTGCCGGCTTCGAAGTCCGTGACGTACACCCGCCCCACTATGGCCGTGTGTGCCCGATCGAGACCCCTGAAGGTCCGAACATCg                                                                                   <  2:6334896/137‑1 (MQ=255)
                                                 aCCCACAAGCGCCGCGTCTCCGCACTCGGCCCTGGCAGTCTGACCCGTGAGCGTGCCGGCTTCGAAGTCCGCGACGTACACCCGCCCCACTATGGCCGTGTGTGCCCGATCGAGACCCCTGAAGGTCCGAACATCgg                                                                                  >  2:14140424/1‑137 (MQ=255)
                                                   ccACAAGCGCCGCGTCTCCGCACTCGGCCCTGGCGGTCTGACCCGTGAGCGTGCCGGCTTCGAAGTCCGTGACGTACACCCGCCCCACTATGGCCGTGTGTGCCCGATc                                                                                                            >  1:7782458/1‑109 (MQ=255)
                                                   ccACAAGCGCCGCGTCTCCGCACTCGGCCCTGGCGGTCTGACCCGTGAGCGTGCCGGCTTCGAAGTCCGTGACGTACACCCGCCCCACTATGGCCGTGTGTGCCCGATCGAGACCCCTGAAGGTCCGAACATCGGTc                                                                                >  1:1213649/1‑137 (MQ=255)
                                                    cacaAGCGCCGCGTCTCCGCACTCGGCCCTGGCGGTCTGACCCGTGAGCGTGCCGGCTTCGAAGTCCGTGACGTACACCCGCCCCACTATGGCCGTGTGTGCCCGATCGAGACCCCTGAAGGTCCGAACATCGGTCt                                                                               >  1:15470001/1‑137 (MQ=255)
                                                    cacaAGCGCCGCGTCTCCGCACTCGGCCCTGGCGGTCTGACCCGTGAGCGTGCCGGCTTCGAAGTCCGTGACGTACACCCGCCCCACTATGGCCGTGTGTGCCCGATCGAGACCCCTGAAGGTCCGAACATCGGTCt                                                                               <  2:5423951/137‑1 (MQ=255)
                                                     acaAGCGCCGCGTCTCCGCACTCGGCCCTGGCGGTCTGACCCGTGAGCGTGCCGGCTTCGAAGTCCGTGACGTACACCCGCCCCACTATGGCCGTGTGTGCCCGATCGAGACCCCTGAAGGTCCGAACATCGGTCTg                                                                              >  2:280368/1‑137 (MQ=255)
                                                           gccgcGTCTCCGCACTCGGCCCTGGCGGTCTGACCCGTGAGCGTGCCGGCTTCGAAGTCCGTGACGTACACCCGCCCCACTATGGCCGTGTGTGCCCGATCGAGACCCCTGAAGGTCCGAACATCGGTCTGATCAAc                                                                        >  1:12486033/1‑137 (MQ=255)
                                                               cgTCTCCGCACTCGGCCCTGGCGGTCTGACCCGTGAGCGTGCCGGCTTCGAAGTCCGTGACGTACACCCGCCCCACTATGGCCGTGTGTGCCCGATCGAGACCCCTGAAGGTCCGAACATCGGTCTGATCAACTccc                                                                    <  2:1567768/137‑1 (MQ=255)
                                                               cgTCTCCGCACTCGGCCCTGGCGGTCTGACCCGTGAGCGTGCCGGCTTCGAAGTCCGTGACGTACACCCGCCCCACTATGGCCGTGTGTGCCCGATCGAGACCCCTGAAGGTCCGAACATCGGTCTGATCAACTccc                                                                    <  2:8487197/137‑1 (MQ=255)
                                                               cgTCTCCGCACTCGGCCCTGGCGGTCTGACCCGTGAGCGTGCCGGCTTCGAAGTCCGTGACGTACACCCGCCCCACTATGGCCGTGTGTGCCCGATCGAGACCCCTGAAGGTCCGAACATCGGTCTGATCAACTccc                                                                    <  2:16107962/137‑1 (MQ=255)
                                                               cgTCTCCGCACTCGGCCCTGGCGGTCTGACCCGTGAGCGTGCCGGCTTCGAAGTCCGTGACGTACACCCGCCCCACTATGGCCGTGTGTGCCCGATCGAGACCCCTGAAGGTCCGAACATCGGTCTGATCAACTccc                                                                    >  2:3869262/1‑137 (MQ=255)
                                                               cgTCTCCGCACTCGGCCCTGGCGGTCTGACCCGTGAGCGTGCCGGCTTCGAAGTCCGTGACGTACACCCGCCCCACTATGGCCGTGTGTGCCCGATCGAGACCCCTGAAGGTCCGAACATCGGTCTGATCAACTccc                                                                    <  1:7248399/137‑1 (MQ=255)
                                                               cgTCTCCGCACTCGGCCCTGGCGGTCTGACCCGTGAGCGTGCCGGCTTCGAAGTCCGTGACGTACACCCGCCCCACTATGGCCGTGTGTGCCCGATCGAGACCCCTGAAGGTCCGAACATCGGTCTGATCAACTccc                                                                    <  2:8196865/137‑1 (MQ=255)
                                                               cgTCTCCGCACTCGGCCCTGGCGGTCTGACCCGTGAGCGTGCCGGCTTCGAAGTCCGTGACGTACACCCGCCCCACTATGGCCGTGTGTGCCCGATCGAGACCCCTGAAGGTCCGAACATCGGTCTGATCAACTccc                                                                    <  1:6511155/137‑1 (MQ=255)
                                                               cgTCTCCGCACTCGGCCCTGGCGGTCTGACCCGTGAGCGTGCCGGCTTCGAAGTCCGTGACGTACACCCGCCCCACTATGGCCGTGTGTGCCCGATCGAGACCCCTGAAGGTCCGAACATCGGTCTGATCAACTccc                                                                    >  1:6032710/1‑137 (MQ=255)
                                                               cgTCTCCGCACTCGGCCCTGGCGGTCTGACCCGTGAGCGTGCCGGCTTCGAAGTCCGTGACGTACACCCGCCCCACTATGGCCGTGTGTGCCCGATCGAGACCCCTGAAGGTCCGAACATCGGTCTGATCAACTccc                                                                    >  2:4454896/1‑137 (MQ=255)
                                                               cgTCTCCGCACTCGGCCCTGGCGGTCTGACCCGTGAGCGTGCCGGCTTCGAAGTCCGTGACGTACACCCGCCCCACTATGGCCGTGTGTGCCCGATCGAGACCCATGAAGGTCCGAACATCGGTCTGATCAACTccc                                                                    >  2:12257481/1‑137 (MQ=255)
                                                                gTCTCCGCACTCGGCCCTGGCGGTCTGACCCGTGAGCGTGCCGGCTTCGAAGTCCGTGACGTACACCCGCCCCACTATGGCCGTGTGTGCCCGATCGAGACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCCCt                                                                   <  1:7936642/137‑1 (MQ=255)
                                                                    ccGCACTCGGCCCTGGCGGTCTGACCCGTGAGCGTGCCGGCTTCGAAGTCCGTGACGTACACCCGCCCCACTATGGCCGTGTGTGCCCGATCGAGACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCCCTggcg                                                               <  2:3144562/137‑1 (MQ=255)
                                                                              cccTGGCGTTCTGACCCGTGAGCGTGCCGGCTTCGAAGTCCGTGACGTACACCCGCCCCACTATGGCCGTGTGTGCCCGATCGAGACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCCCTGGCGGCCTATGccc                                                     <  2:4169733/137‑1 (MQ=255)
                                                                              cccTGGCGGTCTGACCCGTGAGCGTGCCGGCTTCGAAGTCCGTGACGTACACCCGCCCCACTATGGCCGTGTGTGCCCGATCGAGACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCCCTGGCGGCCTATGccc                                                     <  1:10011549/137‑1 (MQ=255)
                                                                              cccTGGCGGTCTGACCCGTGAGCGTGCCGGCTTCGAAGTCCGTGACGTACACCCGCCCCACTATGGCCGTGTGTGCCCGATCGAGACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCCCTGGCGGCCTATGccc                                                     <  2:5146968/137‑1 (MQ=255)
                                                                              cccTGGCGGTCTGACCCGTGAGCGTGCCGGCTTCGAAGTCCGTGACGTACACCCGCCCCACTATGGCCGTGTGTGCCCGATCGAGACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCCCTGGCGGCCTATGccc                                                     <  2:2018602/137‑1 (MQ=255)
                                                                               ccTGGCGGTCTGCCCCGTGAGCGTGCCGGCTTCGAAGTCCGTGACGTACACCCGCCCCACTATGGCCGTGTGTGCCCGATCGAGACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCCATGGCGGCCTATGCCCg                                                    >  2:8074597/1‑137 (MQ=255)
                                                                               ccTGGCGGTCTGACCCGTGAGCGTGCCGGCTTCGAAGTCCGTGACGTACACCCGCCCCACTATGGCCGTGTGTGCCCGATCGAGACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCCCTGGCGGCCTATGCCCg                                                    >  2:3306588/1‑137 (MQ=255)
                                                                               ccTGGCGGTCTGACCCGTGAGCGTGCCGGCTTCGAAGTCCGTGACGTACACCCGCCCCACTATGGCCGTGTGTGCCCGATCGAGACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCCCTGGCGGCCTATGCCCg                                                    >  1:1909803/1‑137 (MQ=255)
                                                                               ccTGGCGGTCTGACCCGTGAGCGTGCCGGCTTCGAAGTCCGTGACGTACACCCGCCCCACTATGGCCGTGTGTGCCCGATCGAGACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCCCTGGCGGCCTATGCCCg                                                    >  1:6815533/1‑137 (MQ=255)
                                                                                cTGGCGGTCTGACCCGTGAGCGTGCCGGCTTCGAAGTCCGTGACGTACACCCGCCCCACTATGGCCGTGTGTGCCCGATCGAGACCCCTGAAGGTCCGAACATCGGTCTGATCAACTccc                                                                    >  1:4447001/1‑120 (MQ=255)
                                                                                cTGGCGGTCTGACCCGTGAGCGTGCCGGCTTCGAAGTCCGTGACGTACACCCGCCCCACTATGGCCGTGTGTGCCCGATCGAGACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCCCTGGCGGCCTATGCCCGc                                                   <  1:15705975/137‑1 (MQ=255)
                                                                                   gCGGTCTGACCCGTGAGCGTGCCGGCTTCGAAGTCCGTGACGTACACCCGCCCCACTATGGCCGTGTGTGCCCGATCGAGACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCCCTGGCGGCCTATGCCCGCacc                                                <  2:9836387/137‑1 (MQ=255)
                                                                                        cTGACCCGGGCGCGTGCCGGCTTTGGAGTCCTTGACGTATACCCCCCCCACTATGGCCGTGTGTGCCCGATCGAGACCCCCGAAGGTCCGAACATCGGTCTGATCAACCCCCTGGCGGCCTATGCCCGCaccaacca                                           <  2:216435/137‑1 (MQ=255)
                                                                                          gACCCGTGAGCGTGCCGGCTTCGAAGTCCGTGACGTACACCCGCCCCACTATGGCCGTGTGTGCCCGATCGAGACCCCTGAAGGTCCGAACATCGGTCTNATCAACTCCCTGGCGGCCTATGCCCGCACCAACCAGt                                         <  1:16262456/137‑1 (MQ=255)
                                                                                           aCCCGTGAGCGTGCCGGCTTCGAAGTCCGTGACGTACACCCGCCCCACTATGGCCGTGTGTGCCCGATCGAGACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCCCTGGCGGCCTATGCCCGCACCAACCAGTa                                        >  2:13231061/1‑137 (MQ=255)
                                                                                                tGAGCGTGCCGGCTTCGAAGTCCGTGACGTACACCCGCCCCACTATGGCCGTGTGTGCCCGATCGAGACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCCCTGGCGGCCTATGCCCGCACCAACCAGTACGGCt                                   >  1:13945248/1‑137 (MQ=255)
                                                                                                   gCGTGCCGGCTTCGAAGTCCGTGACGTACACCCGCCCCACTATGGCCGTGTGTGCCCGATCGAGACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCCCTGGCGGCCTATGCCCGCACCAACCAGTACGGCTTcc                                <  2:5818746/137‑1 (MQ=255)
                                                                                                   gCGTGCCGGCTTCGAAGTCCGTGACGTACACCCGCCCCACTATGGCCGTGTGTGCCCGATCGAGACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCCCTGGCGGCCTATGCCCGCACCAACCAGTACGGCTTcc                                >  1:16017714/1‑137 (MQ=255)
                                                                                                   gCGTGCCGGCTTCGAAGTCCGTGACGTACACCCGCCCCACTATGGCCGTGTGTGCCCGATCGAGACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCCCTGGCGGCCTATGCCCGCACCAACCAGTACGGCTTcc                                >  1:3976393/1‑137 (MQ=255)
                                                                                                   gCGTGCCGGCTTCGAAGTCCGTGACGTACACCCGCCCCACTATGGCCGTGTGTGCCCGATCGAGACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCCCTGGCGGCCTATGCCCGCACCAACCAGTACGGCTTcc                                <  1:9466661/137‑1 (MQ=255)
                                                                                                   gCGTGCCGGCTTCGAAGTCCGTGACGTACACCCGCCCCACTATGGCCGTGTGTGCCCGATCGAGACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCCCTGGCGGCCTATGCCCGCACCAACCAGTACGGCTTcc                                >  2:11716082/1‑137 (MQ=255)
                                                                                                    cGTGCCGGCTTCGAAGTCCGTGACGTACACCCGCCCCACTATGGCCGTGTGTGCCCGATCGAGACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCCCTGGCGGCCTATGCCCGCACCAACCAGTACGGCTTCCt                               <  1:13363678/137‑1 (MQ=255)
                                                                                                    cGTGCCGGCTTCGAAGTCCGTGACGTACACCCGCCCCACTATGGCCGTGTGTGCCCGATCGAGACCCCTGAAGGTCCGAACAGCGGTCTGATCAACTCCCTGGCGGCCTATGCCCGCACCAACCAGTACTGCTTCCt                               <  2:7641597/137‑1 (MQ=255)
                                                                                                     gggCCGGCTTCGAAGTCCGTGACGTACACCCGCCCCACTATGGCCGTGTGTGCCCGATCGAGACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCCCTGGCGGCCTATGCCCGCACCAACCAGTACGGCTTCCTg                              <  1:8492806/135‑1 (MQ=255)
                                                                                                     gTGCCGGCTTCGAAGTCCGTGACGTACACCCGCCCCACTATGGCCGTGTGTGCCCGATCGAGACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCCCTGGCGGCCTATGCCCGCACCAACCAGTACGGCTTCCTg                              <  2:471021/137‑1 (MQ=255)
                                                                                                     gTGCCGGCTTCGAAGTCCGTGACGTACACCCGCCCCACTATGGCCGTGTGTGCCCGATCGAGACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCCCTGGCGGCCTATGCCCGCACCAACCAGTACGGCTTCCTg                              <  2:15534067/137‑1 (MQ=255)
                                                                                                     gTGCCGGCTTCGAAGTCCGTGACGTACACCCGCCCCACTATGGCCGTGTGTGCCCGATCGAGACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCCCTGGCGGCCTATGCCCGCACCAACCAGTACGGCTTCCTg                              <  1:9892048/137‑1 (MQ=255)
                                                                                                        ccGGCTTCGAAGTCCGTGACGTACACCCGCCCCACTATGGCCGTGTGTGCCCGATCGAGACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCCCTGGCGGCCTATGCCCGCACCAACCAGTACGGCTTCCTGGaa                           <  1:10957087/137‑1 (MQ=255)
                                                                                                                  aGTCCGTGACGTACACCCGCCCCACTATGGCCGTGTGTGCCCGATCGAGACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCCCTGGCGGCCTATGCCCGCACCAACCAGTACGGCTTCCTGGAAAGCCCGTAcc                 <  1:6305772/137‑1 (MQ=255)
                                                                                                                  aGTCCGTGACGTACACCCGCCCCACTATGGCCGTGTGTGCCCGATCGAGACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCCCTGGCGGCCTATGCCCGCACCAACCAGTACGGCTTCCTGGAAAGCCCGTAcc                 <  1:10657030/137‑1 (MQ=255)
                                                                                                                  aGTCCGTGACGTACACCCGCCCCACTATGGCCGTGTGTGCCCGATCGAGACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCCCTGGCGGCCTATGCCCGCACCAACCAGTACGGCTTCCTGGAAAGCCCGTAcc                 <  2:947738/137‑1 (MQ=255)
                                                                                                                   gTCCGTGACGTACACCCGCCCCACTATGGCCGTGTGTGCCCGATCGAGACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCCCTGGCGGCCTATGCCCGCACCAACCAGTACGGCTTCCTGGAAAGCCCGTActg                <  1:1725258/137‑3 (MQ=255)
                                                                                                                    tCCGTGACGTACACCCGCCCCACTATGGCCGTGTGTGCCCGATCGAGACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCCCTGGCGGCCTATGCCCGCACCAACCAGTACGGCTTCCTGGAAAGCCCGTACcgc               <  1:14613819/137‑1 (MQ=255)
                                                                                                                      ggTGACGTACACCCGCCCCACTATGGCCGTGTGTGCCCGATCGAGACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCCCTGGCGGCCTATGCCCGCACCAACCAGTACGGCTTCCTGGAAAGCCCGTACCGCgt             >  1:12971878/2‑137 (MQ=255)
                                                                                                                         gACGTACACCCGCCCCACTATGGCCGTGTGTGCCCGATCGAGACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCCCTGGCGGCCTATGCCCGCACAAACCAGTACGGCTTCCTGGAAAGCCCGTACCGCgtggt          >  1:1363262/1‑137 (MQ=255)
                                                                                                                           cGTACACCCGCCCCACTATGGCCGTGTGTGCCCGATCGAGACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCCCTGGCGGCCTATGCCCGCACCAACCAGTACGGCTTCCTGGAAAGCCCGTACCGCGTGGTga        <  1:3430037/137‑1 (MQ=255)
                                                                                                                           cGTACACCCGCCCCACTATGGCCGTGTGTGCCCGATCGAGACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCCCTGGCGGCCTATGCCCGCACCAACCAGTACGGCTTCCTGGAAAGCCCGTACCGCGTGGTga        >  1:4972052/1‑137 (MQ=255)
                                                                                                                           cGTACACCCGCCCCACTATGGCCGTGTGTGCCCGATCGAGACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCCCTGGCGGCCTATGCCCGCACCAACCAGTACGGCTTCCTGGAAAGCCCGTACCGCGTGGTga        <  1:5336401/137‑1 (MQ=255)
                                                                                                                                 cccGCCCCACTATGGCCGTGTGTGCCCGATCGAGACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCCCTGGCGGCCTATGCCCGCACCAACCAGTACGGCTTCCTGGAAAGCCCGTACCGCGTGGTgaaggaag  >  2:10274337/1‑137 (MQ=255)
                                                                                                                                     |                                                                                                                                    
CAGCCAGCTGTCCCAGTTCATGGACCAGAACAACCCTCTCTCGGAGATTACCCACAAGCGCCGCGTCTCCGCACTCGGCCCTGGCGGTCTGACCCGTGAGCGTGCCGGCTTCGAAGTCCGTGACGTACACCCGACCCACTATGGCCGTGTGTGCCCGATCGAGACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCCCTGGCGGCCTATGCCCGCACCAACCAGTACGGCTTCCTGGAAAGCCCGTACCGCGTGGTGAAGGAAG  >  AE015451_phaM1/539161‑539426

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: