Predicted mutation
evidence seq id position mutation annotation gene description
RA AE015451_phaM1 4,586,030 2 bp→TC intergenic (+113/+101) PP_4061 CDS (conserved protein of unknown function) → / ← PP_4063 CDS (putative Long‑chain‑fatty‑acid‑CoA ligase) –/–

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451_phaM14,586,0300CT99.6% 800.5 / NA 229intergenic (+113/+102)PP_4061 CDS (conserved protein of unknown function)/PP_4063 CDS (putative Long‑chain‑fatty‑acid‑CoA ligase)–/–
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base T (9/218);  minor base G (0/1);  total (10/219)
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
*AE015451_phaM14,586,0310TC98.7% 872.4 / ‑3.6 229intergenic (+114/+101)PP_4061 CDS (conserved protein of unknown function)/PP_4063 CDS (putative Long‑chain‑fatty‑acid‑CoA ligase)–/–
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base C (8/218);  total (10/219)
Fisher's exact test for biased strand distribution p-value = 5.05e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 229 total aligned reads displayed.

CCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTCTGC‑GGTAAACCCGCTCCTACC‑AGGGA‑CGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCC  >  AE015451_phaM1/4585900‑4586116
                                                                                                                                  ||                                                                                        
        gAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGTTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCg                                                                             <  1:5147863/137‑1 (MQ=255)
        gAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCg                                                                             <  2:11674892/137‑1 (MQ=255)
        gAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCg                                                                             <  1:16363617/137‑1 (MQ=255)
        gAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCg                                                                             <  1:3578133/137‑1 (MQ=255)
        gAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCg                                                                             <  1:4473869/137‑1 (MQ=255)
        gAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCg                                                                             <  1:6980443/137‑1 (MQ=255)
        gAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCg                                                                             <  1:9354394/137‑1 (MQ=255)
        gAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCg                                                                             <  2:10427638/137‑1 (MQ=255)
        gAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCg                                                                             <  2:11177551/137‑1 (MQ=255)
        gAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCg                                                                             <  2:11283124/137‑1 (MQ=255)
        gAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCg                                                                             <  1:11030933/137‑1 (MQ=255)
        gAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCg                                                                             <  2:11816361/137‑1 (MQ=255)
        gAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCg                                                                             <  2:14693338/137‑1 (MQ=255)
        gAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCg                                                                             <  2:15239826/137‑1 (MQ=255)
        gAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCg                                                                             <  2:1666824/137‑1 (MQ=255)
        gAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCg                                                                             <  2:5105306/137‑1 (MQ=255)
        gAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCg                                                                             <  2:5153891/137‑1 (MQ=255)
        gAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCg                                                                             <  2:7901211/137‑1 (MQ=255)
        gAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCg                                                                             <  2:8712919/137‑1 (MQ=255)
        gAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCg                                                                             <  1:11989265/137‑1 (MQ=255)
        gAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCg                                                                             <  1:14046632/137‑1 (MQ=255)
        gAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCg                                                                             <  1:11560753/137‑1 (MQ=255)
        gAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCACTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCg                                                                             <  1:1726846/137‑1 (MQ=255)
        gAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCGGATCTGCGCTGCTCCGTTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCg                                                                             <  1:481524/137‑1 (MQ=255)
             cTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCt                                                                        <  1:12878374/137‑1 (MQ=255)
                  gggTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACC‑Ag                                                                  <  2:12432108/137‑1 (MQ=39)
                  gggTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACC‑Ag                                                                  <  2:13996754/137‑1 (MQ=255)
                       gACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACC‑AGGGA‑cc                                                            <  1:12393249/137‑2 (MQ=39)
                       gACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACC‑AGGGA‑cc                                                            <  1:150845/137‑2 (MQ=39)
                       gACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACC‑AGGGA‑cc                                                            <  1:5930025/137‑2 (MQ=39)
                        aCCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACC‑AGGGA‑ccg                                                           <  2:15524731/137‑3 (MQ=38)
                        aCCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACC‑AGGGA‑ccg                                                           <  1:367101/137‑3 (MQ=39)
                         ccACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGNTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACC‑AGGGA‑ccgt                                                          <  1:15766464/137‑4 (MQ=38)
                          cACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACC‑AGGGA‑ccgtg                                                         <  1:2439389/137‑5 (MQ=38)
                            cGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACC‑AGGGACCGTGcc                                                        <  2:2178149/137‑1 (MQ=38)
                             ggCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCg                                                       <  2:3245452/137‑1 (MQ=38)
                               cACCGGCGTGAACAGCTGGTCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGc                                                     <  2:1987213/137‑1 (MQ=38)
                                aCCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGcc                                                    >  2:12758416/1‑137 (MQ=38)
                                aCCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGcc                                                    >  1:18350365/1‑137 (MQ=38)
                                    gCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGa                                                <  2:16182943/137‑1 (MQ=38)
                                    gCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGa                                                <  2:16029880/137‑1 (MQ=38)
                                     cGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGaa                                               <  2:4223198/137‑1 (MQ=38)
                                        gAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATgc                                            <  1:4651443/137‑1 (MQ=38)
                                         aaCAGCTGATCTGCGCTGCTCCATTGAGAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACTCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATgcg                                           <  2:18584990/137‑1 (MQ=38)
                                         aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACC‑ATGGACCGTGCCGGCCTTGAATgcg                                           <  2:9282942/137‑1 (MQ=38)
                                         aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATgcg                                           <  1:1249979/137‑1 (MQ=38)
                                         aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATgcg                                           <  1:18052642/137‑1 (MQ=38)
                                         aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATgcg                                           <  2:13852447/137‑1 (MQ=38)
                                         aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATgcg                                           <  2:17480120/137‑1 (MQ=38)
                                         aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATgcg                                           <  2:8221450/137‑1 (MQ=38)
                                         aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATgcg                                           <  1:17939341/137‑1 (MQ=38)
                                         aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATgcg                                           <  2:246620/137‑1 (MQ=38)
                                         aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATgcg                                           <  1:13527451/137‑1 (MQ=38)
                                         aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATgcg                                           <  2:2570642/137‑1 (MQ=38)
                                         aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATgcg                                           <  2:2753710/137‑1 (MQ=38)
                                         aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATgcg                                           <  2:16958767/137‑1 (MQ=38)
                                         aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATgcg                                           <  1:8440283/137‑1 (MQ=38)
                                         aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATgcg                                           <  1:11691004/137‑1 (MQ=38)
                                         aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATgcg                                           <  1:5572718/137‑1 (MQ=38)
                                         aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATgcg                                           <  1:5550407/137‑1 (MQ=38)
                                         aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATgcg                                           <  1:15912548/137‑1 (MQ=38)
                                         aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATgcg                                           <  1:16040855/137‑1 (MQ=38)
                                         aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATgcg                                           <  1:17009431/137‑1 (MQ=38)
                                         aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGACGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATgtg                                           <  1:16201017/137‑3 (MQ=38)
                                         aaCAGCTGATCTGCGCTGCTCCATCGAAAGTTTGAGGAGGACCGGGTGCTTGGCAAGTTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACC‑ATGGACCGTACCGGCCTTGAATgcg                                           <  2:11782704/137‑1 (MQ=37)
                                          aCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATgcgg                                          <  1:13361205/137‑1 (MQ=38)
                                          aCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATgcgg                                          <  1:16001660/137‑1 (MQ=38)
                                           cAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATgcggc                                         <  1:19345172/137‑1 (MQ=38)
                                            aGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATgcggcg                                        <  1:1876407/137‑1 (MQ=38)
                                            aGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATgcggcg                                        <  2:6590620/137‑1 (MQ=38)
                                            aGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATgcggcg                                        <  2:16496945/137‑1 (MQ=38)
                                              cTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATGCGGCGCt                                      <  2:1361752/137‑1 (MQ=38)
                                              cTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATGCGGCGCt                                      <  2:5853843/137‑1 (MQ=38)
                                              cTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATGCGGCGCt                                      <  1:14868769/137‑1 (MQ=38)
                                              cTGAGCTGCGCTACTCCATTGAAAGGTTGAGGATCACCGGGTACGTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTTGCGGGTAAACCCGCTCCTACC‑ATGGACCGTGCCGGCCTTGAATGCGGCGCt                                      <  2:3941882/137‑1 (MQ=37)
                                               ttaTCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTGGAAGGCGgcgcgg                                     <  2:4597407/135‑3 (MQ=37)
                                               tGATCTGCGCTGCTCCATTGAAAGGTTTAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACC‑ATGGACCGTGCCGGCCTTGAATGCGGCGCTg                                     <  2:9671274/137‑1 (MQ=38)
                                               tGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATGCGGCGCTg                                     <  1:3824808/137‑1 (MQ=38)
                                                   cTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATGCGGCGCTGcccc                                 <  2:18424292/137‑1 (MQ=38)
                                                    tGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACC‑ATGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCt                                <  2:18278688/137‑1 (MQ=38)
                                                    tGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCt                                <  2:9750241/137‑1 (MQ=38)
                                                    tGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCt                                <  2:15645242/137‑1 (MQ=38)
                                                    tGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCt                                <  1:18960738/137‑1 (MQ=38)
                                                    tGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCt                                <  1:3504712/137‑1 (MQ=38)
                                                    tGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCt                                <  1:8677588/137‑1 (MQ=38)
                                                    tGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACC‑AGGGACCATGCCGGCCTCGAATGCGGCGCTGCCCCt                                <  2:14604619/137‑1 (MQ=38)
                                                        ctgctCCATTGAAAGTTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTa                            <  2:14786664/137‑1 (MQ=38)
                                                        ctgctCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTACACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTa                            <  1:19128751/137‑1 (MQ=38)
                                                        ctgctCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTa                            <  1:9605382/137‑1 (MQ=38)
                                                        ctgctCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTa                            <  2:1225246/137‑1 (MQ=38)
                                                        ctgctCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTa                            <  2:6281913/137‑1 (MQ=38)
                                                        ctgctCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTa                            <  1:10608784/137‑1 (MQ=38)
                                                        ctgctCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTa                            <  1:17400254/137‑1 (MQ=38)
                                                         tgctCCATTGAAAGGTTGGGGGTGGCAGGGTGCTTGGCGAGGGGAAGAGCTGGAATGGATCGGGTCGGCCTCTGCGCGGGTAAACCCGCTCCTACC‑GTGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGGAg                           <  1:2291774/137‑1 (MQ=255)
                                                         tgctCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAg                           <  1:3454311/137‑1 (MQ=38)
                                                              cATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATgcggcg                                        <  2:5014335/119‑1 (MQ=37)
                                                              cATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAAACCGCTCCCCCCAAGGGA‑ACTGCCGGACCTGGACGCTGTGCTGTCCCCGGTAGGAgag                      >  2:5259717/1‑135 (MQ=255)
                                                               agtaaaaGGTGGGGGAGGACCGGGTGCTGGACAAGGTGAAGATCTTGCATTGATTGGGCTGTCCTCTTCGCGGGGAAACCCGCCCCTACC‑AGGGACCGTGCCGTCCTGGAATGCGACGCTGCCCCGGGTAGGAGCgg                     <  2:14357948/133‑1 (MQ=255)
                                                                        ttGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAg                           >  2:11569844/1‑122 (MQ=37)
                                                                        ttGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACC‑AGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAg                           <  1:11569844/122‑1 (MQ=37)
                                                                                                                                  ||                                                                                        
CCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTCTGC‑GGTAAACCCGCTCCTACC‑AGGGA‑CGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCC  >  AE015451_phaM1/4585900‑4586116

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