Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451_phaM15,737,7560AG38.2% 209.3 / 86.5 213intergenic (+45/+26)hutF CDS (probable formiminoglutamate deiminase)/PP_5037 CDS (lipocalin family lipoprotein)–/–
Reads supporting (aligned to +/- strand):  ref base A (38/93);  new base G (0/81);  total (38/175)
Fisher's exact test for biased strand distribution p-value = 1.17e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 218 total aligned reads displayed.

TGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCC‑TTCGCGTGCAAGCCCGCGAA‑GGGTTGCAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAG  >  AE015451_phaM1/5737621‑5737892
                                                                                                                                        |                                                                                                                                          
        gggTGTTGCGCCAGGGTCTACATGCCCAGGAAGAGGAAAGTGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACTGGGGTTGTGTTACAGCCC‑TTCGCGTGCAAGCCCGCGAA‑GGGTTgcaa                                                                                                                                  <  2:16496537/137‑1 (MQ=255)
                        gcggcATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCTTTTCCGCCC‑CTCTCAGTCCAGACCTCGAAGGGGTTTCCACAC‑ACACCCaaaaaa                                                                                                                  >  2:2478289/1‑137 (MQ=255)
                         cggcGTCCCCGGGAGGGTGAAATTGGGCGGGCTTTTGCAAGGGTTTTGCTCCAGTTTTTGGGCTGAGCGTACTGGGGTGTCGTGGCATCCC‑TTTGGGTGCAAGCCCGCGAA‑GGGTTCCAACGC‑AACCCCAAAAAACt                                                                                                                <  2:12654766/137‑1 (MQ=255)
                         cggcATGCCCGGGAAGGGGAAAGTGGGCGGGCGTGTGCCAAGTTGTTGCACCACCTTTTGGACTAAGCGTACCGGGGTTGCGTTGCAGCCC‑TTCGCGTGCAACCCCGCGAG‑GGGTTGAAACGC‑AACCCCAAAAAACt                                                                                                                <  2:4390681/137‑1 (MQ=255)
                           gcATGCCAAGGAAGAGGAAAGCGGGCGGGCGTTTGCAAAGGTGTTACACCAGTTTTGGGGCGGAGCGGACTGGGGTTGCGTGGCAACCT‑TTCGCGGGCAACCCCGCGAA‑GGGTGGCAACGC‑AACCCCAAAAAACTCa                                                                                                              <  2:5372505/137‑1 (MQ=255)
                                ccagggAAGATGAAAGTGGGCGGGCGTTGGCCAGGGTGTTGCGCCAGTTTTTGGGCGGAGCGTACTGGGGTTGCTTGGCAGCCC‑TTCGTGTGCAAGCCCGCGAG‑GGGTGGCAACGC‑AACCCCAAAAAACTCAGTGCa                                                                                                         <  2:8047496/133‑1 (MQ=255)
                                     gAAGATAAAAGCGGGCGGGCGTTTGCCAAGGTGTTTCGCCATTTTTTGGGTTGAGTGTACCGGGTTTGCGTTGCAGCTT‑TTCGTGGGAAAGCCCGTGAG‑GGGTTGAAACGC‑AACCCCAAAAAACTCAGTGCAttttt                                                                                                    <  1:14380911/137‑1 (MQ=255)
                                                  ggcgggcgTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTTTGTTTCAACCC‑CTCTCGTTTAAACCCACAAA‑AAGGTGCCACACAAACCCCAAAAAACCCA‑TGGATTTTTTCTCtttttttat                                                                                       >  2:12760770/1‑131 (MQ=255)
                                                     ggtcgTTGGCCAAGGTGTTACGCCACTTTTTGTGCTGAGTGTACCGGGGTTGCGTTTCAGCCC‑TTCACGTGCAGCCCTGCGAA‑GGGTTGCAACGC‑AACCCCAAAACACTCAGTGCATTTTTACTATTTGCTGATCgg                                                                                    <  2:15756094/134‑1 (MQ=255)
                                                        cgTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACTGGGGTTGCTTTGCAGCCC‑TTCGCGTGCAAGCCCGCGAA‑GGGTTGAAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGAcg                                                                                 <  2:16288515/137‑1 (MQ=255)
                                                         ttttGCCAGGGTTTTCCACCAGTTTTTGGGTTAAGCGTACTGGGGTTGTTTTGAATCCT‑TTCGCGGAAAAGCCCGCGAG‑GGGTTGAAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGAcgc                                                                                <  2:9199014/136‑1 (MQ=255)
                                                          tttGCAAAGGTGTTGCGCCATCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCATCCC‑TTCGCGGGAAAGCCCGCGAA‑GGGTTGAAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGAcgcc                                                                               <  2:8133947/137‑1 (MQ=255)
                                                             gcaaaGGGGTTGCACAGGTTTTTGGGTGGAGTGGACGGGGGTTGTGTTGCGGCCT‑TTTGTGGGAAAGCCCGCGAG‑GGGTTGCAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGAcgcccgc                                                                            <  2:3678511/134‑1 (MQ=255)
                                                                 ggttgtTGCGCCAGCTTTTGGGTTGAGGGTATTGGGGTTGTTTTGCGCTCT‑TTTGCGGGCAGGCCCGCGAG‑GGGTTGCAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGa                                                                        <  1:17281805/134‑1 (MQ=255)
                                                                   gttttGCACCAGCTTTGGGGCGGAGAGTACTGGGTTTGCTTTGAAACCC‑TTCTAGTACAAGCCCGCGAA‑GGGTGACAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATc                                                                      <  1:15128439/134‑1 (MQ=255)
                                                                   gtgtTTCGCTGGCTTTGGGGCGGGGCGCGCTGGGGGTGCGTTACACCCT‑TTTGCGTGCAATCCCGCGAG‑GGGTTGCAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATc                                                                      <  1:3319366/137‑1 (MQ=255)
                                                                   gtgtTGCGCCAGCTTTTGGGCTGAGCGTACGGGGGTTGTGTGGCAGCCC‑TTTGCGTGCAAGCCTGCGAA‑GGGTGGCAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATc                                                                      <  2:6615793/137‑1 (MQ=255)
                                                                       tGCTCCAGCTTTTGGCCTGATTTCACTTGGGTTTCTTTTCGGCCT‑TTCTCGTGCAAGCCCGCGAA‑GGGTTGCAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGc                                                                  <  2:6243145/137‑1 (MQ=255)
                                                                       tGCGCCATCTTTTGGGCGGAGCGGA‑CGGGGTTGGTGTGAAACCCTTTTGCGGGAAACCCCGGGAG‑GGGTGGAAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGc                                                                  <  1:6057443/137‑1 (MQ=255)
                                                                       tGCGCCAGTTTTGGGGTGGGGCGGATGGGGGGGGTGTTGCGGCCT‑TTCGTGGGCAAGCCCGGGAG‑GGGTTGAAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGc                                                                  <  2:4253102/137‑1 (MQ=255)
                                                                       tGCGCAAGTTTTTGGGTGGGGTGTACTGGGGTTGTGTTGCATCCT‑TTCGCGGGCAAGCCCGCGAG‑GGGTGGCAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGc                                                                  <  1:12926681/137‑1 (MQ=255)
                                                                       tGCGCAAGCTTTGGGGCTGAGCGGACTGGGGTTGCGGTGCAGCCC‑TTCGTGTGAAAGCCCGCGAA‑GGGTGGCAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGc                                                                  <  2:5858115/137‑1 (MQ=255)
                                                                       ggcgcCAGCTTTTGGGCTGAGCGTACCGGGTTTGTTTTGCAGCCC‑TTCGTGGGCAATCCCGCGAG‑GGGTGGCAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGc                                                                  <  2:13359102/136‑1 (MQ=255)
                                                                        gcgcCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCC‑TTCGCGTGCAAGCCCGCGAA‑GGGTTGCAACGC‑AATCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCg                                                                 <  2:11526417/137‑1 (MQ=255)
                                                                         cgcCATCTTTTGGGCTGAGCGTACCGGGGTTGCTTTGCATCCT‑TTCGCGTGCAAGCCCACGAG‑GGGTTACAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCgg                                                                <  2:341140/137‑1 (MQ=255)
                                                                         cgcCAGTTTTTGGGTGGAGCGTACGGGGTTTGCGTTGCAGCCC‑TTCGCGGGAAAGCCGGAGAG‑GGGTTGCAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCgg                                                                <  1:7340228/137‑1 (MQ=255)
                                                                         cgcAAGTTTTTGGGTGGAGCGTACTGGGGTGGCGTTGCAGCCC‑TTCACGTAAAAGCCCGAGAA‑GGGTGGCAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCgg                                                                <  2:4482257/137‑1 (MQ=255)
                                                                         cgcAAGTTTTGGGGTGGAGCGTACGGGGGTGGGGGTGAACCCC‑TTCGTGGACAAACCAGAGAA‑GGGTGGAAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCgg                                                                <  1:13668673/137‑1 (MQ=255)
                                                                          tccAGTTTTTGGTCTGAGCGTACCGGGGTTGCGTTGCAACCC‑TTCGCGTACAAGCCCGCGAA‑GGGTGGCAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGc                                                               <  2:15144816/136‑1 (MQ=255)
                                                                          gccggttttgGGGCTGAGCGTACCGGGGTTGGGTTGCAGCCC‑TTCGCGTGAAAGCCCGCGAA‑GGGTTGCAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGc                                                               <  2:11343784/131‑1 (MQ=255)
                                                                          cccAGCTTTTGGGCTGAGCGTACCGGGGTGGCGTTACAGCCC‑TTCGCGTGCAATCCCGCGAA‑GGGTTGCAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGc                                                               <  2:3260605/136‑1 (MQ=255)
                                                                             ggttttgGGGTTGAGTGTACGGGGGTTTTGTTGCATCCT‑TTTGTGGGAAAGCCTGGGAG‑GGGTTGCAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTgg                                                            <  1:9822630/134‑1 (MQ=255)
                                                                             agtttttGGGCGGAGCGTACGGGGGTTGTGTTGCATCCC‑TTCGCGGGAAAGCCCGCGAA‑GGGTGACAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTgg                                                            <  2:12256157/134‑1 (MQ=255)
                                                                             aGCTTTTGGGCTGAGCGGACTGGGGTTGCGTGGCAGCCC‑TTCGCGTGAAAGCCCGCGAA‑GGGTTAAAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTgg                                                            <  2:12751328/137‑1 (MQ=255)
                                                                                          agtgtcccGGGTTTGCTTTGAATCCC‑TTTGCGTGCAGGCCCGCGAG‑GGGTTGCAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCctg                                               <  2:6340902/131‑1 (MQ=255)
                                                                                            tttACGGGGGTTGTGTTGCGGCCT‑TTCGTGGGAAATCCTGAGAG‑GGGTTAAAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCctgct                                             <  1:15648114/135‑1 (MQ=255)
                                                                                            tggactggggTTGTGTTGAATCCC‑TTCGTGTGAAAGCCCGCGAG‑GGGTTGAAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCctgct                                             <  1:3640678/131‑1 (MQ=255)
                                                                                            tggacgggggTTGTGTTGCAGCCT‑TTCGTGTGAAAGACCGGGAG‑GGGTGGCAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCctgct                                             <  1:4360123/131‑1 (MQ=255)
                                                                                            tgTACCTGGTTTGCGTTGAAGCCC‑TTCGCGTGCAAGCCCGCGAA‑GGGTTGCAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTGGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCctgct                                             <  2:12040269/136‑1 (MQ=255)
                                                                                            gttACTGGGGTTGTTTTGCAGCTC‑TTCACGTACAAGCCCGCGAA‑GGGTTACAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCctgct                                             <  2:87106/135‑1 (MQ=255)
                                                                                            cttACCGGGGTTTCTTTGAACCCC‑TTCTCGGACAACCCCACGAG‑GGGTGAAAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCctgct                                             <  1:14953298/135‑1 (MQ=255)
                                                                                            cggattggggTGGTGTTGCATCCT‑TTTGCGGGAAATCCTGGGAG‑GGGTTAAACCGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCctgct                                             <  1:12236334/131‑1 (MQ=255)
                                                                                            cGTACTGGGGTTGCGTTGCAGCCC‑TTCGCGTTCAAGCCCGCGAA‑GGGTTGCAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCctgct                                             <  2:2233203/137‑1 (MQ=255)
                                                                                            cGTACTGGGGTTGCGTGGCGTCCC‑TTCGATAGCAAGCCCCCGAA‑GGGTTGAAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCctgct                                             <  1:8426980/137‑1 (MQ=255)
                                                                                            cGTACCGGGTTTGCGTGGAATCCC‑TTCGCGTGCAGCCCCGCGAA‑GGGTTACAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCctgct                                             <  1:15461542/137‑1 (MQ=255)
                                                                                            cGTACCGGGGTTGTGTTGCAGCCC‑TTCGCGGGAAAGCCTGCGAG‑GGGTGGCAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCctgct                                             <  1:8790885/137‑1 (MQ=255)
                                                                                            cGTACCGGGGTTGCGTTGCAGCCC‑TTCGCGTACAAGCCCGCGAG‑GGGTTGCAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCctgct                                             <  2:13497954/137‑1 (MQ=255)
                                                                                            cGTACCGGGGTGGCGTGACAGCTC‑TTCGCGTGAAAGCTCGCGAG‑GGGTTGCAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCctgct                                             <  1:16363971/137‑1 (MQ=255)
                                                                                               aCCGGGGTTGCGTTACAGCCC‑TTCGTGTACAAGCCCGCGAA‑GGGTGGCAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGcc                                          <  2:2928715/137‑1 (MQ=255)
                                                                                                 tggggTTGCGTTGCAGCCC‑TTCGCGTGCAATCCCGCGAG‑GGGTTGCAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCcgc                                        <  2:3628234/136‑1 (MQ=255)
                                                                                                 tggggTGGCGTTGAAGCCC‑TTCGTGGGCAAGCCCGCGAG‑GGGTTGCAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCcgc                                        <  2:10446286/136‑1 (MQ=255)
                                                                                                 cGGGGTTGCGTTGCAGCCC‑TTCGCTTGCACGCCCGCGAA‑GGTTTGCAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCcgc                                        <  2:17293311/137‑1 (MQ=255)
                                                                                                    ggTTGCGTTGCAGCCC‑TTCGCGTGCAAGCCCGCGAA‑GGGTTGCAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGcc                                     <  2:3876484/137‑1 (MQ=255)
                                                                                                    ggTTGCGTTGCAGCCC‑TTCGCGTACAAGCCCGCGAG‑GGGTTGCAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGcc                                     <  2:8328769/137‑1 (MQ=255)
                                                                                                         cgttgcAGCCC‑TTCGCTTACAAACCCCCGAA‑GGGTTGCAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGa                                <  2:6000838/137‑1 (MQ=255)
                                                                                                            tgcACCCT‑TTTGCGTGCAAGCCTGTGAG‑GGGTGGCAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGAcag                             <  1:16861588/137‑1 (MQ=255)
                                                                                                                 ccc‑TTCGCGTGCAAGCCCGCGAA‑GGGTTGCAACGC‑AACCCCCAAAAACTCAGTGCATTTTTACTGTTTGCTGATCGGACGGACACCAGATCAAGCGGATGGTGTCATGCCCCTGCTGCCGCGCCCTGGACAGAAGGt                        >  2:13573750/1‑137 (MQ=255)
                                                                                                                 ccc‑TTCGCATGCCAGCCCACGAA‑GGGTTTCAACGC‑AACCCCAAAAAAATCCGTGTATTTTTACCATTTGCTGATCGGACGCCCGTCAGATCAAGCGGCTGGTGTCCTCCCCCTGCTGCCGCGCCCTGGACAGCAGGt                        >  2:1311156/1‑137 (MQ=255)
                                                                                                                   c‑TTCGCGTGCAAGCCTGCGAG‑GGGTTGCAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTtc                      <  2:2098348/137‑1 (MQ=255)
                                                                                                                   c‑TTCGCGTGCAAGCCCGCGAA‑GGTTTGCAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTtc                      <  2:3534365/137‑1 (MQ=255)
                                                                                                                   c‑TTCGCGTGCAAGCCCGCGAA‑GGTTTACAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTtc                      <  1:11357004/137‑1 (MQ=255)
                                                                                                                   c‑TTCGCGTGCAAGCCCGCGAA‑GGGTTGCAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTtc                      <  2:6640824/137‑1 (MQ=255)
                                                                                                                   c‑TTCGCGTGCAAGCCCGCGAA‑GGGTTGCAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTtc                      <  2:913238/137‑1 (MQ=255)
                                                                                                                   c‑TTCGCGTGCAAGCCCACGAA‑GGGTTGCAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTtc                      <  1:15356696/137‑1 (MQ=255)
                                                                                                                    ctttgtgTGAAAGTCCGTGAA‑GGGTTGCAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTtc                      <  1:10578362/131‑1 (MQ=255)
                                                                                                                       cgcgTGCAAGCCCGCGAA‑GGGTTGCAACGC‑AACCCCAAAAAACTCAGTGCCTTTTTACCATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCCGCTGCCGCGCCCTGGACAGCAGGTTCTCa                   >  1:1445693/1‑137 (MQ=255)
                                                                                                                        acatGCCAGCCCGCGAA‑GGATTGCAACGT‑ACCCCCAAAAAACCCAGTGCATTTTTTCTCTTTGCTGACCTGACGCCCGCCAAATCAAACTGTTGGTTTCATACCCCTGCTGCCGCGCCCTGGGGAGCAGGTCCTCAc                  >  2:360824/4‑137 (MQ=255)
                                                                                                                          ttGCAATCCCGCAAA‑GGGTTTAAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGc                <  1:14138356/136‑1 (MQ=255)
                                                                                                                          gTGCAAGCCGGCGAG‑GGGTTGCAACGC‑AACCCCAAAAAACTCATTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGc                <  2:2878476/137‑1 (MQ=255)
                                                                                                                          gTGCAAGCCCGCGAA‑GGGTTGCAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGc                >  1:13314916/1‑137 (MQ=255)
                                                                                                                          gTGCAAGCCCGCGAA‑GGGTTGCAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGc                <  2:7525204/137‑1 (MQ=255)
                                                                                                                          gTGCAAGCCCGCGAA‑GGGTTGCAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGc                >  2:7443663/1‑137 (MQ=255)
                                                                                                                          gTGCAAGCCCGCGAA‑GGGTTGCAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGc                >  2:5154933/1‑137 (MQ=255)
                                                                                                                          gTGCAAGCCCGCGAA‑GGGTTGCAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGc                >  2:14050331/1‑137 (MQ=255)
                                                                                                                          gTGCAAGCCCGCGAA‑GGGTTGCAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGc                >  2:14840792/1‑137 (MQ=255)
                                                                                                                          gTGCAAGCCCGCGAA‑GGGTTGCAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTAGACAGCAGGTTCTCACGc                >  2:1714147/1‑137 (MQ=255)
                                                                                                                          gTGCAAGCCCGCGAA‑GGGTTGCAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCAAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACgt                >  2:14160697/1‑136 (MQ=255)
                                                                                                                          gTGCAAGCCCGCGAA‑GGGTTGCAACGC‑AACCCCAAAAAACTCAGTGCATCTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGc                >  2:3061622/1‑137 (MQ=255)
                                                                                                                                gCCCGCGAA‑GGGTTGCAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCgg                                                                                    <  2:15216984/63‑1 (MQ=255)
                                                                                                                                gCCCGCGAA‑GGGTTGCAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAg          >  2:5734193/1‑137 (MQ=255)
                                                                                                                                gCCCGCGAA‑GGGTTGCAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAg          >  2:10029927/1‑137 (MQ=255)
                                                                                                                                 cccGCGAA‑GGGTTGCAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAgg         >  2:15831140/1‑137 (MQ=255)
                                                                                                                                  cctcgAC‑GGGTTGCAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCACCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGc        <  2:6681707/134‑1 (MQ=255)
                                                                                                                                  cctcgAA‑GGGTTGCAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGc        <  2:12478781/134‑1 (MQ=255)
                                                                                                                                  ccGCGAG‑GGGTTGCAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGc        <  1:11847681/137‑1 (MQ=255)
                                                                                                                                  ccGCGAA‑GGGTTGCAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGc        >  2:11933211/1‑137 (MQ=255)
                                                                                                                                  ccGCGAA‑GGGTTGCAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGc        <  2:7459171/137‑1 (MQ=255)
                                                                                                                                  ccGCGAA‑GGGTTGCAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGc        <  2:14698620/137‑1 (MQ=255)
                                                                                                                                  ccGCGAA‑GGGTTGCAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGc        >  2:10300817/1‑137 (MQ=255)
                                                                                                                                  ccGCGAA‑GGGTTGCAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGc        <  2:4058033/137‑1 (MQ=255)
                                                                                                                                  ccGCGAA‑GGGTTGCAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGc        <  2:10895183/137‑1 (MQ=255)
                                                                                                                                  ccGCGAA‑GGGTTGAAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGc        <  2:5430060/137‑1 (MQ=255)
                                                                                                                                  ccGCGAA‑GGGTTACAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGc        <  2:4776651/137‑1 (MQ=255)
                                                                                                                                  ccGCGAA‑GGGTTACAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTAATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGc        <  2:7290299/137‑1 (MQ=255)
                                                                                                                                    gcgAG‑GGGTTGCAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCg      <  2:98475/137‑1 (MQ=255)
                                                                                                                                    gcgAA‑GGGTTGCAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCg      >  1:11535933/1‑137 (MQ=255)
                                                                                                                                    gcgAA‑GGGTTGCAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCg      >  1:8458835/1‑137 (MQ=255)
                                                                                                                                    gcgAA‑GGGTGGCAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCg      <  2:15226051/137‑1 (MQ=255)
                                                                                                                                     cgAA‑GGTTTACACCGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCgg     <  1:1862513/137‑1 (MQ=255)
                                                                                                                                        a‑GGGTTGCAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCCTACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAg  >  2:5545105/1‑137 (MQ=255)
                                                                                                                                        |                                                                                                                                          
TGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCC‑TTCGCGTGCAAGCCCGCGAA‑GGGTTGCAACGC‑AACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAG  >  AE015451_phaM1/5737621‑5737892

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