Predicted mutation
evidence seq id position mutation annotation gene description
RA AE015451_phaM1 1,070,248 +GA intergenic (+48/+94) aroP‑I CDS (aromatic amino acid transport protein) → / ← PP_0928 CDS (K+‑dependent Na+_Ca+ exchanger related‑protein) –/–

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451_phaM11,070,2461.G97.4% 2437.0 / 53.9 648intergenic (+46/+96)aroP‑I CDS (aromatic amino acid transport protein)/PP_0928 CDS (K+‑dependent Na+_Ca+ exchanger related‑protein)–/–
Reads supporting (aligned to +/- strand):  ref base . (4/13);  new base G (289/340);  total (293/355)
Fisher's exact test for biased strand distribution p-value = 8.43e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*AE015451_phaM11,070,2462.A97.4% 2287.8 / 47.8 648intergenic (+46/+96)aroP‑I CDS (aromatic amino acid transport protein)/PP_0928 CDS (K+‑dependent Na+_Ca+ exchanger related‑protein)–/–
Reads supporting (aligned to +/- strand):  ref base . (4/13);  new base A (288/342);  total (293/355)
Fisher's exact test for biased strand distribution p-value = 8.51e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 655 total aligned reads displayed.

ATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGT‑‑GA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGG  >  AE015451_phaM1/1070111‑1070373
                                                                                                                                          ||                                                                                                                                 
     cATCTCGGTAATGCTGATCCCGGCATGGATCGCGGTGCTGTGGGTGGCGTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGAGCGGGTTTGTTCGGAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGG‑‑Gagatc                                                                                                                             <  2:15090954/137‑5 (MQ=255)
             tAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGcgccg                                                                                                                     >  2:6471207/1‑137 (MQ=255)
                tGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCcgcgc                                                                                                                  <  1:16589081/137‑1 (MQ=255)
                 gCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCGGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCcgcgcg                                                                                                                 <  2:304245/137‑1 (MQ=255)
                     aTCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCgcgggcg                                                                                                             >  2:11145539/1‑137 (MQ=255)
                        ccGGTATGGATCGCTGTGCTGTGGTTGGCTTTCAAGCTGAAGAAGGCTCGCCGGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGgcggcg                                                                                                          <  2:4747614/137‑1 (MQ=38)
                        ccGGCATGGATCGCTGTGCTGTGGTTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGgcggcg                                                                                                          <  1:9514164/137‑1 (MQ=255)
                        ccGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGgcggcg                                                                                                          >  1:873174/1‑137 (MQ=255)
                        ccGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGgcggcg                                                                                                          >  2:6241789/1‑137 (MQ=255)
                        ccGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCAGCGCGGgcggcg                                                                                                          >  1:13496864/1‑137 (MQ=255)
                        acGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGgcggcg                                                                                                          >  2:976242/2‑137 (MQ=255)
                             aTGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGa                                                                                                     >  2:10598878/1‑137 (MQ=255)
                               ggATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGAtc                                                                                                   <  2:1991338/137‑1 (MQ=255)
                               ggATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGGGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGAtc                                                                                                   <  1:5116902/137‑1 (MQ=255)
                               ggATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGGGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGAtc                                                                                                   <  1:11686234/137‑1 (MQ=255)
                                    gctgtgctgtGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTcaca                                                                                              >  1:595827/1‑137 (MQ=255)
                                    gctgtgctgtGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTcaca                                                                                              >  1:16242546/1‑137 (MQ=255)
                                    gctgtgctgtGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTcaca                                                                                              >  2:3538515/1‑137 (MQ=255)
                                        tgctgtgGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCg                                                                                          <  2:8589797/137‑1 (MQ=255)
                                        tgctgtgGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTAATAGGGTTTGTCCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCA‑‑TGTGAGA‑‑TCGAATGCCGCGGGGGCGGCGCCCGATCTAACAGGCg                                                                                          >  2:13034255/1‑137 (MQ=255)
                                         gctgtgGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGa                                                                                         <  2:1721493/137‑1 (MQ=255)
                                         gctgtgGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGa                                                                                         <  2:13717762/137‑1 (MQ=255)
                                         gctgtgGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGa                                                                                         <  1:6093202/137‑1 (MQ=255)
                                         gctgtgGGTGGCTTTCAAGCTAAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGa                                                                                         <  1:3138209/137‑1 (MQ=255)
                                         gctgtgGGTGGCTTTCAAGCGGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGa                                                                                         <  1:152500/137‑1 (MQ=255)
                                         gatgtgGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCTGGCGGCGCTCGATCTCACAGGCGa                                                                                         <  2:16783717/135‑1 (MQ=38)
                                          ctgtgGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGAt                                                                                        >  1:6643147/1‑137 (MQ=255)
                                                ttggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGc                                                                                  <  2:9385892/136‑1 (MQ=255)
                                                gtggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGc                                                                                  <  2:19073115/137‑1 (MQ=255)
                                                 tggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGcc                                                                                 <  1:12726659/137‑1 (MQ=255)
                                                  ggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCt                                                                                >  2:14915018/1‑137 (MQ=255)
                                                     tttCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTaaa                                                                             >  1:18140052/1‑137 (MQ=255)
                                                     tttCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTaaa                                                                             >  2:1527639/1‑137 (MQ=255)
                                                     tttCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTaaa                                                                             >  2:15893641/1‑137 (MQ=255)
                                                     tttCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTaaa                                                                             >  2:18764273/1‑137 (MQ=255)
                                                      ttCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAg                                                                            >  1:13165297/1‑137 (MQ=255)
                                                       tCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCCCAGGCGATGAAAGCCTAAAgg                                                                           >  1:11358219/1‑137 (MQ=255)
                                                       tCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAgg                                                                           >  1:2242791/1‑137 (MQ=255)
                                                       tCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAgg                                                                           >  1:14404420/1‑137 (MQ=255)
                                                       tCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAgg                                                                           >  2:11675815/1‑137 (MQ=255)
                                                       tCAAGATGACGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCAGCACGGGCGGAGCTCGCTCTCACAGGCGATGAAAGCCTAAAgg                                                                           >  1:8384825/1‑137 (MQ=38)
                                                          aGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCAt                                                                        >  1:15554219/1‑137 (MQ=255)
                                                          aGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCAt                                                                        >  1:7905594/1‑137 (MQ=255)
                                                            cTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACCGGCGATGAAAGCCTAAAGGCATGc                                                                      >  1:14544202/1‑137 (MQ=255)
                                                            cTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGc                                                                      >  1:14907633/1‑137 (MQ=255)
                                                               aagaagGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTc                                                                   <  1:9710198/137‑1 (MQ=255)
                                                               aagaagGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTc                                                                   <  2:3897879/137‑1 (MQ=255)
                                                                 gaagGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTg                                                                 >  1:4855448/1‑137 (MQ=255)
                                                                  aagGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTgg                                                                >  1:9917551/1‑137 (MQ=255)
                                                                      cTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTggcggc                                                            <  2:12341796/137‑1 (MQ=255)
                                                                          ccAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTg                                                        <  1:2871803/137‑1 (MQ=255)
                                                                            aGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGAtc                                                                                                   >  2:8180343/1‑92 (MQ=37)
                                                                            aGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGAc                                                      >  1:9252930/1‑137 (MQ=255)
                                                                            aGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGAc                                                      >  1:6218405/1‑137 (MQ=255)
                                                                            aGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGAc                                                      <  2:19373630/137‑1 (MQ=255)
                                                                               ccAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGAcgcg                                                   <  1:10077594/137‑1 (MQ=255)
                                                                                cAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGt                                                  <  2:1402089/137‑1 (MQ=255)
                                                                                 aaGGGCTGATCGGGTTTGTTCGTAAGGGGTGTGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGtt                                                 <  1:3516136/137‑1 (MQ=255)
                                                                                 aaGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGtt                                                 >  1:10277640/1‑137 (MQ=255)
                                                                                 aaGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACACGtt                                                 >  1:4653549/1‑137 (MQ=255)
                                                                                      cTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCt                                            >  2:11932370/1‑137 (MQ=255)
                                                                                      cTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCt                                            >  2:17534850/1‑137 (MQ=255)
                                                                                       tGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCtt                                           <  1:12405126/137‑1 (MQ=255)
                                                                                         aTCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGTCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTcc                                         <  2:18334938/137‑1 (MQ=255)
                                                                                         aTCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTcc                                         <  1:15942077/137‑1 (MQ=255)
                                                                                         aTCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTcc                                         <  1:7655087/137‑1 (MQ=255)
                                                                                         aTCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTcc                                         <  2:5413249/137‑1 (MQ=255)
                                                                                           cGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCAt                                       >  2:6734298/1‑137 (MQ=255)
                                                                                           cGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCAt                                       >  2:4187253/1‑137 (MQ=255)
                                                                                           cGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCAt                                       >  1:17447879/1‑137 (MQ=255)
                                                                                           cGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCAt                                       >  1:12235919/1‑137 (MQ=255)
                                                                                            gggTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCCAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATg                                      >  1:19369205/1‑137 (MQ=255)
                                                                                                   ttCGTAATGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCTCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTTGCGGCCCTGACGCGTTCCCCCTCCCTTccccccc                               <  2:449839/137‑4 (MQ=25)
                                                                                                             gTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTc                     <  1:3731659/137‑1 (MQ=255)
                                                                                                             gTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTc                     <  1:18665922/137‑1 (MQ=255)
                                                                                                             gTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTc                     <  2:12808879/137‑1 (MQ=255)
                                                                                                               aTGCCGTCAGGTTTGCGGCGCCTTTGA‑‑GA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCg                   <  1:19073609/137‑1 (MQ=37)
                                                                                                                   cGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCaa               >  2:16601224/1‑137 (MQ=255)
                                                                                                                   cGTTAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCaa               >  1:9000543/1‑137 (MQ=255)
                                                                                                                     ttAGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgc             <  2:4995268/137‑1 (MQ=255)
                                                                                                                       aGGTTTTCGGCGTC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcac           <  2:2707354/137‑1 (MQ=255)
                                                                                                                       aGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGATCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcac           <  2:16277229/137‑1 (MQ=255)
                                                                                                                       aGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCTGGCTGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcac           <  2:14445409/137‑1 (MQ=38)
                                                                                                                       aGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGTCGCTCGATCTCACAGGCGATGAAAGTCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcac           <  2:7083159/137‑1 (MQ=38)
                                                                                                                       aGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcac           <  1:10737580/137‑1 (MQ=255)
                                                                                                                       aGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcac           <  1:17013439/137‑1 (MQ=255)
                                                                                                                       aGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcac           <  1:7398116/137‑1 (MQ=255)
                                                                                                                       aGGTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcac           <  1:8192697/137‑1 (MQ=255)
                                                                                                                        ggTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcacg          >  2:9089553/1‑137 (MQ=255)
                                                                                                                         gTTTTCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcacgc         >  1:5588953/1‑137 (MQ=255)
                                                                                                                          ttttCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCTGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcacgca        <  2:18033110/137‑1 (MQ=255)
                                                                                                                          ttttCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATNCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcacgca        <  1:17976549/137‑1 (MQ=255)
                                                                                                                          ttttCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcacgca        <  1:15086744/137‑1 (MQ=255)
                                                                                                                          ttttCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcacgca        <  2:5800771/137‑1 (MQ=255)
                                                                                                                          ttttCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcacgca        <  1:2007331/137‑1 (MQ=255)
                                                                                                                          ttttCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcacgca        <  2:7466839/137‑1 (MQ=255)
                                                                                                                          ttttCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcacgca        <  1:10924580/137‑1 (MQ=255)
                                                                                                                          ttttCGGCGCC‑‑TGT‑‑GAGGTCGAGCGCCGCGCGGGCGGCGCTCGATCCCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcacgca        <  1:1401395/137‑1 (MQ=255)
                                                                                                                            ttCGGCGCC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGcaca      <  1:4123655/137‑1 (MQ=255)
                                                                                                                                gcgcC‑‑TGTGAGA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAgtgg  >  1:13796708/1‑137 (MQ=255)
                                                                                                                                          ||                                                                                                                                 
ATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGT‑‑GA‑‑TCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGG  >  AE015451_phaM1/1070111‑1070373

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