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.5% 705.3 / NA 201intergenic (+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 (0/1);  new base T (5/195);  total (5/196)
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,0310TC99.5% 776.2 / NA 201intergenic (+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 (0/1);  new base C (5/195);  total (5/196)
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 201 total aligned reads displayed.

GGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTCTGC‑GG‑TAAACCCGCTCCTACCAGGGA‑CGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGG  >  AE015451_phaM1/4585905‑4586092
                                                                                                                             ||                                                                
ggTGAATGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGTTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAAc                                                       <  2:14007307/137‑1 (MQ=255)
 gTGAAGGCTTAATGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCGTTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAAcc                                                      <  2:753845/137‑1 (MQ=255)
   gAAGGCTTAATGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCg                                                    <  2:8254500/137‑1 (MQ=255)
   gAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCg                                                    <  2:15595342/137‑1 (MQ=255)
   gAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCg                                                    <  1:2412092/137‑1 (MQ=255)
   gAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCg                                                    <  2:4145155/137‑1 (MQ=255)
   gAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCg                                                    <  1:3371137/137‑1 (MQ=255)
   gAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCg                                                    <  2:2933010/137‑1 (MQ=255)
   gAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCg                                                    <  2:18376681/137‑1 (MQ=255)
   gAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCg                                                    <  1:5492394/137‑1 (MQ=255)
   gAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCg                                                    <  1:580754/137‑1 (MQ=255)
   gAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCg                                                    <  1:1619155/137‑1 (MQ=255)
   gAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCg                                                    <  1:6018973/137‑1 (MQ=255)
   gAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCg                                                    <  1:15691620/137‑1 (MQ=255)
   gAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCg                                                    <  2:7302876/137‑1 (MQ=255)
   gAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCg                                                    <  2:8155589/137‑1 (MQ=255)
   gAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCg                                                    <  2:16960412/137‑1 (MQ=255)
   gAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCg                                                    <  1:13578099/137‑1 (MQ=255)
   gAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCg                                                    <  1:12973452/137‑1 (MQ=255)
   gAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCg                                                    <  2:15022931/137‑1 (MQ=255)
   gAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCg                                                    <  1:11317496/137‑1 (MQ=255)
   gAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCg                                                    <  1:11316579/137‑1 (MQ=255)
   gAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGACGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCg                                                    <  2:11678119/137‑1 (MQ=255)
   gAAGGCTGAAGGGTAGACCACGGCACNGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCg                                                    <  1:18623121/137‑1 (MQ=255)
   gAAGGCGTAAGGTTAGACCGCGGAACAGGCGTGAACAGCTGATCTGCGCTGCTCCATTTAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCg                                                    <  1:1020920/137‑1 (MQ=255)
        cTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCt                                               <  1:18704118/137‑1 (MQ=255)
             gggTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAg                                          <  1:2037773/137‑1 (MQ=39)
             gggTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAg                                          <  1:11211047/137‑1 (MQ=39)
             gggTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAg                                          <  2:5169704/137‑1 (MQ=39)
             gggTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAg                                          <  2:10608364/137‑1 (MQ=39)
             gggTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAg                                          <  2:8609342/137‑1 (MQ=39)
             gggTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAg                                          <  1:8792623/137‑1 (MQ=39)
             gggTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTACTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCat                                          <  2:13201384/137‑2 (MQ=39)
                  gACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGA‑cc                                    <  2:9516841/137‑2 (MQ=39)
                  gACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGA‑cc                                    <  1:11968916/137‑2 (MQ=39)
                   aCCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGA‑ccg                                   <  2:11191874/137‑3 (MQ=38)
                   aCCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGA‑ccg                                   <  2:7894061/137‑3 (MQ=38)
                   aCCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGA‑ccg                                   <  2:18266742/137‑3 (MQ=38)
                    ccACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGA‑ccgt                                  <  2:16238689/137‑4 (MQ=38)
                    ccACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGA‑ccgt                                  <  1:9780684/137‑4 (MQ=38)
                      aCGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGA‑cc                                    >  2:1055520/1‑132 (MQ=38)
                          cACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACCGTGCCGGc                             <  1:7931696/137‑1 (MQ=38)
                               gagTGAACAGCTTATCTGCACTGCTCCATTCAAAGTTTGAGCATTACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGC‑GGATAAACTCGCTCCTACCATGGACCGTGCCGGCCTTGa                        <  2:16707868/135‑1 (MQ=255)
                               gCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGa                        <  1:3334740/137‑1 (MQ=38)
                               gCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGa                        <  1:1441451/137‑1 (MQ=38)
                               gCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGa                        <  1:13949371/137‑1 (MQ=38)
                               gCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGa                        <  1:7762279/137‑1 (MQ=38)
                               gCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGa                        <  1:5564219/137‑1 (MQ=38)
                               gCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGa                        <  1:16532864/137‑1 (MQ=38)
                               gCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGa                        <  2:16363246/137‑1 (MQ=38)
                               gCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGa                        <  2:5701035/137‑1 (MQ=38)
                                   gAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTG‑GGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATgcg                   <  2:5960883/137‑1 (MQ=38)
                                    aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGATGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATgcg                   <  2:14101128/137‑1 (MQ=38)
                                    aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTTGCGGG‑TAAACCCGCTCCTACCATGGACCGTGCCGGCCTTGAATgcg                   <  1:3129174/137‑1 (MQ=255)
                                    aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATgcg                   <  2:1262533/137‑1 (MQ=38)
                                    aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATgcg                   <  1:12136196/137‑1 (MQ=38)
                                    aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATgcg                   <  1:15111406/137‑1 (MQ=38)
                                    aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATgcg                   <  1:15602379/137‑1 (MQ=38)
                                    aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATgcg                   <  1:4980314/137‑1 (MQ=38)
                                    aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATgcg                   <  2:4754608/137‑1 (MQ=38)
                                    aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACCGTGCCGGCCGTGAATgcg                   <  2:4007134/137‑1 (MQ=38)
                                     aCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCATGGACCGTGCCGGCCTTGAATgcgg                  <  1:4309318/137‑1 (MQ=38)
                                     aCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATgcgg                  <  2:12308180/137‑1 (MQ=38)
                                     aCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATgcgg                  <  1:17255788/137‑1 (MQ=38)
                                     aCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATgcgg                  <  2:15533458/137‑1 (MQ=38)
                                      cAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATgcggc                 <  2:18060759/137‑1 (MQ=38)
                                      cAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATgcggc                 <  2:7180687/137‑1 (MQ=38)
                                      cAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATgcggc                 <  2:8834838/137‑1 (MQ=38)
                                      cAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATgcggc                 <  1:848868/137‑1 (MQ=38)
                                      cAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATgcggc                 <  1:17179500/137‑1 (MQ=38)
                                      cAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATgcggc                 <  1:7431252/137‑1 (MQ=38)
                                       aGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATTGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATgcggcg                <  2:15805808/137‑1 (MQ=38)
                                       aGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATgcggcg                <  2:12466837/137‑1 (MQ=38)
                                       aGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATgcggcg                <  1:11519225/137‑1 (MQ=38)
                                       aGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATgcggcg                <  2:17225769/137‑1 (MQ=38)
                                       aGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATgcggcg                <  2:13101508/137‑1 (MQ=38)
                                       aGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATgcggcg                <  1:717299/137‑1 (MQ=38)
                                       aGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATgcggcg                <  1:14991365/137‑1 (MQ=38)
                                       aGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATgcggcg                <  1:15393197/137‑1 (MQ=38)
                                       aGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATgcggcg                <  1:15877148/137‑1 (MQ=38)
                                       aGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATgcggcg                <  1:16449261/137‑1 (MQ=38)
                                       aGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATgcggcg                <  2:4702355/137‑1 (MQ=38)
                                       aGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGACGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATgcggcg                <  1:4094817/137‑1 (MQ=38)
                                        gCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGacgc               <  1:13385654/137‑1 (MQ=38)
                                         cTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCt              <  1:14276065/137‑1 (MQ=38)
                                         cTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCt              <  1:10595647/137‑1 (MQ=38)
                                         cTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCt              <  1:17801543/137‑1 (MQ=38)
                                          tGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGTCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTg             <  1:18424699/137‑1 (MQ=38)
                                          tGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTg             <  1:13936499/137‑1 (MQ=38)
                                          tGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTg             <  2:14890810/137‑1 (MQ=38)
                                           gATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGc            <  1:14881481/137‑1 (MQ=38)
                                              cTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGcccc         <  2:6024623/137‑1 (MQ=38)
                                               tGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCt        <  1:4880650/137‑1 (MQ=38)
                                                   ctgctCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTa    <  2:5305406/137‑1 (MQ=38)
                                                   ctgctCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTa    <  1:3441629/137‑1 (MQ=38)
                                                   ctgctCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTa    <  2:3183505/137‑1 (MQ=38)
                                                   ctgcgCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTa    <  2:3694492/137‑1 (MQ=38)
                                                                                               gAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCt              >  1:2008840/1‑83 (MQ=21)
                                                                                               gAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCt              <  2:2008840/83‑1 (MQ=21)
                                                                                               gAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCt              <  2:11951655/83‑1 (MQ=21)
                                                                                                                             ||                                                                
GGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTCTGC‑GG‑TAAACCCGCTCCTACCAGGGA‑CGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGG  >  AE015451_phaM1/4585905‑4586092

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