Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451_phaM15,199,6550CT29.2% 84.6 / 28.0 110intergenic (+30/+43)PP_4579 CDS (Transcriptional regulator_ LysR family)/PP_4580 CDS (conserved protein of unknown function)–/–
Reads supporting (aligned to +/- strand):  ref base C (20/55);  new base T (21/10);  total (41/69)
Fisher's exact test for biased strand distribution p-value = 1.43e-04
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 123 total aligned reads displayed.

TGGCAGACCGGGGTGGCGTTGGTCGATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTCTCCGCCTGTAGCGGCCCTTTCGCGGGCATGCCCGCGAAAGGGCCGGTACAGGTACAGCAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGCGCCGCCCGCCCAGGATCAGCCAGCCAATCCCCAGCAACAGGCTTTCGACCACAAAGGCCAGC  >  AE015451_phaM1/5199521‑5199788
                                                                                                                                        |                                                                                                                                      
tggcAGACCGGGGTGGCGTTGGTCGATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGTGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTCTCCGCCTGTAGCGGCCCCTTcac                                                                                                                                      >  2:12337635/1‑135 (MQ=255)
 ggcAGACCGGGGTGGCGTTGGTCGATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTCTCCGCCTGTAGCGGCCCTTTcttt                                                                                                                                     >  1:8390013/1‑134 (MQ=255)
 ggcAGACCGGGGTGGCGTTGGTCGATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTCTCCGCCTGTAGCGGCCCTTTcgcg                                                                                                                                     <  1:1127343/137‑1 (MQ=255)
 ggcAGACCGGGGTGGCGTTGGTCGATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTCTCCGCCTGTAGCGGCCCTTTcgcg                                                                                                                                     <  1:7407377/137‑1 (MQ=255)
 ggcAGACCGGGGTGGCGTTGGTCGATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTCTCCGCCTGTAGCGGCCCTTTcgcg                                                                                                                                     <  2:1384142/137‑1 (MQ=255)
 ggcAGACCGGGGTGGCGTTGGTCGATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTCTCCGCCTGTAGCGGCCCTTTcgcg                                                                                                                                     <  2:3634835/137‑1 (MQ=255)
 ggcAGACCGGGGTGGCGTTGGGCGATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTCTCCGCCTGTAGCGGCCCTTTcgcg                                                                                                                                     <  2:2100955/137‑1 (MQ=255)
  gcAGACCGGGGTGGCGTTGGTCGATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTCTCCGCCTGTAGCGGCCCTTTCGCgg                                                                                                                                    <  2:1611112/137‑1 (MQ=255)
  gcAGACCGGGGTGGCGTTGGTCGATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTCTCCGCCTGTAGCGGCCCTTTCGCgg                                                                                                                                    <  2:5777486/137‑1 (MQ=255)
    aGACCGGGGTGGCGTTGTTCGATGAGGTGGTTGGGGTTGGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTCTCCGCCTGTAGCGGCCCTTTCGCGGGc                                                                                                                                  <  1:17177489/137‑1 (MQ=255)
    aGACCGGGGTGGCGTTGGTCGATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTCTCCGCCTGTAGCGGGCCTTTcttgggg                                                                                                                                  >  1:3733985/1‑131 (MQ=255)
    aGACCGGGGTGGCGTTGGTCGATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTCTCCGCCTGTAGCGGGCCTTTcttggga                                                                                                                                  >  2:5868516/1‑131 (MQ=255)
    aGACCGGGGTGGCGTTGGTCGATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTCTCCGCCTGTAGCGGCCCTTTCTCgggg                                                                                                                                  >  2:8154983/1‑136 (MQ=255)
    aGACCGGGGTGGCGTTGGTCGATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTCTCCGCCTGTAGCGGCCCTTTCGCGGGc                                                                                                                                  <  2:3242572/137‑1 (MQ=255)
    aGACCGGGGTGGCGTTGGTCGATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTCTCCGCCTGTAGCGGCCCTTTCGCGGGc                                                                                                                                  <  2:10271180/137‑1 (MQ=255)
    aGACCGGGGTGGCGTTGGTCGATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTCTCCGCCTGTAGCGGCCCCTTCCCgggg                                                                                                                                  >  2:3090859/1‑136 (MQ=255)
    aGACCGGGGTGGCGTTGGTCGATGAGGTGGTTGGGGGTTGCCGGCGGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTCTCCGCCTGTAGCGGCCCTTTCGCGGGc                                                                                                                                  <  1:17177488/137‑1 (MQ=255)
     gACCGGGGTGGCGTTGGTCGATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTCTCCGCCTGTAGCGGGCCCTTTctggggg                                                                                                                                  >  1:8025449/1‑131 (MQ=255)
     gACCGGGGTGGCGTTGGTCGATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTCTCCGCCTGTAGCGGCCCTTTTGCGGGCa                                                                                                                                 >  2:9014778/1‑137 (MQ=255)
     gACCGGGGTGGCGTTGGTCGATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTCTCCGCCTGTAGCGGCCCTTTCGCgggaa                                                                                                                                 >  2:16755862/1‑135 (MQ=255)
      aCCGGGGTGGCGTTGGTCGATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTCTCCGCCTGTAGCGGCCCTTTCGAGGGCAt                                                                                                                                >  2:1520/1‑137 (MQ=255)
       ccGGGGTGGCGTTGGTAGAGGAGGGGGTTGGGGTGTGCCGGCAGGTGTTGGAGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGGGCTGGTCTGACAAACAGTCTCCGCCTGTAGCGGCCCTTTCGCGGGCATg                                                                                                                               <  2:5192937/137‑1 (MQ=255)
                cGTTGGTCGATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTCTCCGCCTGTAGCGGGCCCTTTGTGGGAATGGCCACGaaa                                                                                                                      >  1:16008643/1‑137 (MQ=255)
                   tGGTCGATGAGGTGGTTGGGGTTTGTCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTCTCCGCCTGTAGCGGACCTTTTGTGGGAATGTCCTCGAAAggg                                                                                                                   >  2:7129448/1‑137 (MQ=255)
                   tGGTCGATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTCTCCGCCTGTAGCGGCCCTTTTGCGGGCCTTTCCGCGAAAggg                                                                                                                   >  2:17370819/1‑137 (MQ=255)
                   tGGTCGATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTCTCCGCCTGTAGCGGCCCCTTCGGAGGGATGGCCGCGAAAggg                                                                                                                   >  2:3538040/1‑137 (MQ=255)
                   tGGTCGATGAGGTGGTTGGGGTTTGCAGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTCTCCGCATGTAGCGGACCTTTCGTAGGCATGCCCACAAAAagg                                                                                                                   >  2:11985796/1‑137 (MQ=255)
                     gTCGATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTCTCCGCCTGTAGCGGCCCTTTTTCGGGCATGCCCACAAAAaggac                                                                                                                 >  1:12156354/1‑135 (MQ=255)
                       cGATGAGGTGGTTGGGGTTTTCCGGCAGGTGTTGGCGCGGTATGCGTGGGGTGTGGGTGGGAAGCGGGTTGGTCTGGTTTGACAAAAGTTTTCCGCTTGTAGTGGCCTTTTCGTGGACGTACCCGAGAAGGGGCCgg                                                                                                               <  2:13453105/137‑1 (MQ=255)
                       cGATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTCTCCGCCTGTAGCGGCCCTTTCTCGGGCCTTCCCACAAAAAGGGCgg                                                                                                               >  1:13533522/1‑137 (MQ=255)
                        gATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTCTCCGCCTGTAGCGGACCCCTTTTTGGGATGCCCACGAAAgggggggg                                                                                                              >  1:14729207/1‑132 (MQ=255)
                        gATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTCTCCGCCTGTAGCGGCCCTTTTTCGGGGATGCCCACGAAAAGGCCggt                                                                                                              >  2:14950574/1‑137 (MQ=255)
                        gATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTCTCCGCCTGTAGCGGCCCTTTCGTGGGGATGTCCGCAAAAgggatggg                                                                                                              >  2:678598/1‑132 (MQ=255)
                        gATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTCTCCGCCTGTAGCGGCCCCTTTTCCGGCATTCCCCCCAAAagggcagg                                                                                                              >  2:2339097/1‑132 (MQ=255)
                        gATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAACGGATCGTGCTGGTCTGACCAACAGTCTCCGCCTGTAGCGGCCCCTTCCCGTGCACGCCCACGAAAAGGccagg                                                                                                              >  1:13968430/1‑134 (MQ=255)
                          tggggtgttTGGGGTTTGCTTGCGGGTGTTGGCGCGGTATGCGCGGGCTGGGGGTGGTCAGGGAGTCGGGCTGTTCGGCCCAAAAGTCTCCGCCTGTAGCGGCCCTTTCGCGGGCATTCCCGCGAAAGGGCCggtac                                                                                                            <  2:3638937/134‑1 (MQ=255)
                          ggggggggtTGGGTTTTGCTGGGAGGTTTGGGTGTGGGATGCGCGGGGTGTGGGTGGGCGGGGGATTGTGTGGGTCTGACAAAAGTTCTCCGCTTGGAGCGCCCTTTTCGTGGGCGTTCCCGCGAAAGGGCCggtac                                                                                                            <  2:7558157/131‑1 (MQ=255)
                            agttggtggGGGTTTATTGGGGGGTGTTGGCGCGGGGGACGTGGGATGTGGGTGGGCAGCGGGTTTTGCTGGCCGGACCAACGTTCTCCGCCTGGAGCGGCCCTTTCGTGGGCATGCCCGAAAAAGGGCCggtacag                                                                                                          <  2:6334761/134‑1 (MQ=255)
                            aGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTCTCCGCCTGTAGCGGACCTTTCGCGGGCATGCCCGCGAAAGGGACagtacaa                                                                                                          >  2:11073742/1‑136 (MQ=255)
                                 tgtGGGGTTGGCGGGAGGGTTTTGGCGAGGGGTGCGGGGGATGTGGGTGTGCAAGGGATCGTGTGGGTCGAACAAACGTTCCCCGCCTGCAGTGGCCTTTTGGCGGGCATACCAGGGAAAGGGCCggtacaggtaca                                                                                                     <  2:2889968/135‑1 (MQ=255)
                                 gtTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTCTCCGCCTGTAGCGGCCCCTTTTCTGGCACTCCCCCGAAAAGGACagtgcaggtaca                                                                                                     >  2:14774455/1‑137 (MQ=255)
                                  ttGGGGTTTGCCGGCAGTTGTTGTTGTGTTATGCGCGGGATGTGGGTGGGCGGCGGGTCGTGCTGTTCTGACCAACATTCTCCGCTTGCAGCGTCCCTTTTGCGGGCGTGCCCGCGAAAGGGCCGgtacaggtacag                                                                                                    <  2:16113394/137‑1 (MQ=255)
                                  tgggggTTTGCTGCTAGGGTTGGGCGCGGTATGCGCGGGATGTGGGTGGGCAGTGGAGCGTGTTGTTTTGACCAACATTCTCTGCCTTTAGCGGCCTTTTCGCGGGCATGCCCGCGAAAGGGCCGgtacaggtacag                                                                                                    <  2:10829613/135‑1 (MQ=25)
                                          tGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTCTCCGCCTGTAGCGGCCCTTTTACGGGCACTCCCACGAAAAGGCCAGTACAAGTACCACAGAACtt                                                                                            >  2:7219350/1‑137 (MQ=255)
                                           gCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTCTCCGCCTGTAGCGGACCTTTTGCGGGGCTTTCCGCGAAAAGGCCAGTACAAATACAGCAAAAcctc                                                                                           >  2:3968328/1‑134 (MQ=255)
                                            ccGGCAGGTGTTGGCGCGGTATGCGCGGAATGTGGGTGGGCAGCGGATCGTGCGGGCCTAACCAACAGCCTCCGCCTGGAGCGACCCTTTCGCGGGCATTCCCGCGAAAGGGCCGGTACAGGTACAGCAGAACTTAg                                                                                          <  2:2003253/137‑1 (MQ=255)
                                             cGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTCTCCGCCTGTAGCGGCCCCTTCTCAGGGACGCCCACGAAAAGGGCAGTACAAGGAAAGAAGAAcctcac                                                                                         >  2:557835/1‑132 (MQ=255)
                                              ggCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTCTCCGCCTGTAGCGGCCCTTTTCCGGGGACGCCCGCGGAAAGGACGGTACAAGTACAAAAAAACttcacg                                                                                        >  2:37943/1‑133 (MQ=255)
                                               gcggttgttgtCGCGGGATGTGTGGGATGTGGGTGGGCAGAGGATCGTTCTGGTTGGACCAACAGTCTCCGCCTGGAGCGGCCTTTTCGGGGGCAGGCCCGCGAAAGGGCCGGTACAGGTACAGCAGAACTTAgcgc                                                                                       <  2:2951204/132‑1 (MQ=37)
                                               gCAGGTTTTGGCCCGGGAGGCGTGGGGTGGGGGTGGGCGGCGGATTGTGTTGGTCTGACCAACGGTCTCCGCCTGTAGCGGCCCTTTCTCGGGCAGGCCCGCGAAAGGGCCGGTACAGGTACAGCAGAACTTAgcgc                                                                                       <  2:16649553/137‑1 (MQ=255)
                                                catttgtTGGCGTGGTATGCGCGGAGTGTGGGTGGAAAGTGGATCGTGTTGTTCGGACCAACATTCCCCTCTTGTAGTGTCCCTTTCGCGGTCGTCCCCGCGAAAGGGCCGGTACAGGTACAGCAGAACTTAGCGcc                                                                                      <  2:9373048/133‑1 (MQ=255)
                                                cAGGTTTTGGCGCGGTGTGCGTGGGGGGGGGGTGGGCAGCGGATCGTGTGGGTCGGACCAACAGTCTCCGCTTGTGGCGTCCCTTTCGCGGGCATGCCCGCGAAAGGTCCGGTACAGGTACAGCAGAACTTAGCGcc                                                                                      <  2:13113106/137‑1 (MQ=255)
                                                aggttgttgtCGTGGTGTGTGTGGGGTGTGGGTGGGCGGGGGGTTGTGTTGTTCGGACCAACGTTCTCCTCTTTTGGTGTCCTTTTCGTGGGCATTCCTGCGAAAGGGCCGGTACAGGTACAGCAGAACTTAGCGcc                                                                                      <  2:16066212/133‑1 (MQ=255)
                                                    tgtTGGCGTGTTATGCGCGGGATGTGGTTGGCCATTGGATCGTGCTGGTCTGACCAACGTTCTCCGCCTGTAGCGGCCCTTTTGTGGGCAGCCCCAAGAAAGGGCCGGTACAGGTACAGCAGAACTTAGCGCCAGct                                                                                  <  2:15442116/137‑1 (MQ=255)
                                                         tcgtggTATGCGCGGGGTGGGGTTGGAAGGCGGATGGTGCTGGTCTGACCAATGGTTTCTGCTTGTGGTGGCCTTTTCGTGGGCATGCCCGCGAAGGGGCCGGTACAGGTACAGCAGAACTTAGCGCCAGCTCTCtt                                                                             <  2:3336329/133‑1 (MQ=255)
                                                         tcgcgGCATGCGTGGGGTGTGGTTTGGCAGTGGATCGTGTTGGTCTAACCAACATTCTCCTCTTGGAGTGTCCCTTTCGCGGGCATACCCGCGAAAGGGCCGGTACAGGTACAGCAGAACTTAGCGCCAGCTCTCtt                                                                             <  2:17396818/136‑1 (MQ=255)
                                                         gcgcGTTATGCGCGGGGTGGGGGTGGGCAGCGGATCGTGCTGGTCTAACCAACGTTCTCCGCTTGGAGCGGCCCTTTCGCGGGCATGCCCGAGAAAGGGCCGGTACAGGTACAGCAGAACTTAGCGCCAGCTCTCtt                                                                             <  2:16945525/137‑1 (MQ=38)
                                                         gcgcGGGGTGTGTGGGATGTGGGTGGGCAGGGGATCGTGCTGTTTTGACCAACGTTCTCCGCTTGGAGCGGCCCTTCTGTGGGAATTCCCGCGAAAGGGCCGGTACAGGTACAGCAGAACTTAGCGCCAGCTCTCtt                                                                             <  2:5350754/137‑1 (MQ=255)
                                                                   cggggggTGTGGGTGGGGGGAGGATTGTGTGGGTCGGACAAAAATTCTCCGCTTGGAGCGTCCTTTTTGTGGGAATACCCGCGAAAGGGCCGGTACAGGTACAGCAGAACTTAGCGCCAGCTCTCTTTCACCACCCt                                                                   <  2:8990936/134‑1 (MQ=255)
                                                                        gATGTGGTTGGTCATCGGATTTTGCTGTGCGGACAAACAGTTTCCTCCTGTAGTTGCCCTTTCGCTGGCATGCTTGGAAAGGGGCCGGTACAGGTACAGCAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGcgcc                                                              <  2:14200110/137‑1 (MQ=255)
                                                                            ttggtgggCGGAGGATCGTGTTGGTCTGACAAACGTCCTCCGCCTGTAGCGGCCCTTTCTTGGGCGTGCCCGCGAAAGGGCCGGTACAGGTACAGCAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGCGccgccc                                                         <  2:16755046/135‑1 (MQ=255)
                                                                            tggttgggCAGAGGGTCGTGTTGGTCCGACCAACAGTCTCCTCCTGTAGCGGCCCTTTCGCGGGCAGGCCCACGAAAGGGCCGGTACAGGTACAGCAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGCGccgccc                                                         <  2:10093760/137‑1 (MQ=255)
                                                                            tgggtgggCGGCGGATCGTGTTGGTCGGACCAAAAGTCTCCGCCTGTAGCGGCCTTTTCGCGGGCATGCCCGCGAAAGGGCCGGTACAGGTACAGCAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGCGccgccc                                                         <  2:7678523/137‑1 (MQ=39)
                                                                            tgggtgggCGATGGATCGGATTGGTCTGACCAACAGACTCCGCCTGGATCGGCCCTTTAGCGGGCAGGCCCGCGAAAGGGCCGGTACAGGTACAGCAGAGCTTAGCGCCAGCTCTCTTTTACCACCCTGCGccgccc                                                         <  2:11359015/137‑1 (MQ=255)
                                                                            tgggtgggCAGCGGATCGTGCTGGTCTGACCAACAGTCTCCGCCTGTAGCGGCCCTTTTGCCGGAATGCCCGCGAAAGGGGCGGTACCGGTCCCGGAAAACTTAGCGCCAGCCCTCTTTCACCAACCTTCGccgccc                                                         >  2:12272415/1‑137 (MQ=255)
                                                                            tgggggggAAGCGGATTGTGTGGGCCGGACAAAAATTCTCCGCCGGGAGCGCCCTTTTCGCGGGCATGCCCGAGAAAGGGCCGGTACAGGTACAGCAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGCGccgccc                                                         <  2:1316928/137‑1 (MQ=255)
                                                                            gggttgggCAGTGGGTTGTGTGGTTCTGACAAACATTCCCCTCTTGGATCTCCCCTTTTGTGGGCGTCCCCAAGAAAGGGCCGGTACAGGTACAGCAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGCGccgccc                                                         <  1:11999949/133‑1 (MQ=255)
                                                                                tgggAAGGGGATTGTGTTGGTCGGACAAACGTTTTCCTCTTGGAGCGTCCTTTTTGTGGGCGTTCCCGCGAAGGGTCCGGTACAGGTACAGCAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGCGCCgcccgccc                                                     <  1:8579341/137‑1 (MQ=255)
                                                                                     aGCGGATCGTGCTGTTCTGACAAACAGTCTCCGCTTGTAGTGTCCCTTTCGCGGGCGTGCCCGTGAAAGGGCCGGTACAGGTACAGCAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGCGCCGCCCGCCCAGGAt                                                <  2:8173235/137‑1 (MQ=255)
                                                                                     aGCGGATCGTGCTGGTCTGACCAACGGTCTCCGCTTGTAGCGCCCCTTTCTCGGGCATGCCTGCGAAAGGGCCGGTACAGGTACAGCAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGCGCCGCCCGCCCAGGAt                                                <  2:3240188/137‑1 (MQ=255)
                                                                                     aGCGGATCGTGCTGGTCTGACCAACAGTCTCCGCCTGTAGCGGCCCTTTTTGGGGGATTTCCCCCAAAAGGACGATACAAATACCACAAAACTTCACGCCCACTCTCCTTCCCCCCCCCACACCCCCCCCCCacaaa                                                >  1:7258114/1‑133 (MQ=255)
                                                                                     aGCGGATCGTGCTGGTCTGACCAACAGTCTCCGCCTGTAGCGGCCCCTTCGCAGGGATGCCCGCAAAAAGGGCGGGGCAAGGAAAGAAGGAAGTCGAGCCAGACCTCTTTCCACAACACGCGCCGCCCGACCaaggg                                                >  2:17025240/1‑133 (MQ=255)
                                                                                       cGGATCGTGTTGGCCTGACCAACAGTCTCCGCCTTTAGCGTCCCTTTCTCGGGCATACCCGCGAAGGGGCCGGTACAGGTACAGCAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGCGCCGCCCGCCCAGGATca                                              <  2:11277498/137‑1 (MQ=255)
                                                                                       cGGATCGTGCTGGTCTGACCAACAGTCTCCGCCTGTAGCGGGCCTTTCGCGGGCATTCCCGCGAAAAGGAAAGTACAAGTGCCACAAAACCTCGCTTCCGATTTCTTTCCCCACCCCTCGCCGGCCCCCCCAgggca                                               >  2:3765032/1‑134 (MQ=255)
                                                                                          gttgttcTGTTCTGACAAACATCCTCCGCTGGTGGTGCCCTTTTCGGGGGCATCCCTGCGAAAGGGCCGGTACAGGTACAGCAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGCGCCGCCCGCCCAGGATcagcc                                           <  2:4939646/131‑1 (MQ=255)
                                                                                          attgTGCTGTTCTGACCAACGGTCTCCTCTGGTGGCGTCCTTTTTGCGGGCATTCCCGCGAAAGGGCCGGTACAGGTACAGCAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGCGCCGCCCGCCCAGGATcagcc                                           <  1:17339343/134‑1 (MQ=255)
                                                                                            cttGCTGGTCTGACCAACATTCTCCGCCTGTGGCGGCCCTTCCGCGGGCATTCCCGCGAAAGGTCCGGTACAGGTACAGCAGAACTTAGCGCCAGTTCTCTTTCACCACCCTGCGCCGCCCGCCCAGGATcagccag                                         <  2:16237063/135‑1 (MQ=255)
                                                                                             gTGCGGTTCTAACCAACATTCTCTTCCTGTGGTGGCCCTTTCGCGGACATGCCCGCGAAAGGGCCGGTACAGGTACAGCAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGCGCCGCCCGCCCAGGATcagccagc                                        <  2:16301206/137‑1 (MQ=38)
                                                                                                        aCCAACATTCTCCGCCTGTAGCGGCCCTTTCGCGGGCATGCCCACGAAAGGGCCGGTACAGGTACAGCAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGCGCCGCCCGCCCAGGATCAGCCAGCCAATCCCCAGc                             <  2:4832772/137‑1 (MQ=255)
                                                                                                         ccagcggTCTCAGCTGGTAGTGGCCCTTTCGGGGGCGGACCCCCAAAAGGTCCGGTACAGGTACAGCAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGCGCCGCCCGCCCAGGATCAGCCAGCCAATCCCCAGCa                            <  2:11248308/131‑1 (MQ=255)
                                                                                                         ccaccaGCCTCCACATGTAGCGTCCCTTTAGCGGGCATTCCCGCGAAGGGGCCGGTACAGGTACAGCAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGCGCCGCCCGCCCAGGATCAGCCAGCCAATCCCCAGCa                            <  2:3927727/137‑1 (MQ=255)
                                                                                                         ccAACGTTCTCCTCTTGGAGCGTCCTTTTTGGGGGCGTTCCCGCAAAAGGGCCGGTACAGGTACAGCAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGCGCCGCCCGCCCAGGATCAGCCAGCCAATCCCCAGCa                            <  1:11258457/137‑1 (MQ=255)
                                                                                                         ccAACGTTCTCCGCCTGTAGTGCCCCTTTGGCGGGCATCCCCACGAAAGGGCCGGTACAGGTACAGCAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGCGCCGCCCGCCCAGGATCAGCCAGCCAATCCCCAGCa                            <  1:11258427/137‑1 (MQ=255)
                                                                                                         ccAACATTCTCCTCTTGGAGTGTCCCTTTCGGGGACGTGCCAGAGAAAGGGCCGGTACAGGTACAGCAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGCGCCGCCCGCCCAGGATCAGCCAGCCAATCCCCAGCa                            <  2:11530936/137‑1 (MQ=255)
                                                                                                         ccAACATTCTCCTCCTGTAGTGCCCCTTTCGCGGGCATGCCAGAAAAAGGGCCGGTACAGGTACAGCAGAACTTAGCGCCAGCTCTCTTTCACCAACCTGCGCCGCCCGCCCAGGATCAGCCAGCCAATCCCCAGCa                            <  2:4479463/137‑1 (MQ=255)
                                                                                                           aaCAGTCTCCGCCTGTAGCGGCCCTTTTGCGGGTACGTCCCAAAAAAGGGCGGTACCGGGACAGAAAAAATTCGGAGCAAATCTCTTTTCCCCCCCCGTGGCACCCGCCCAGGAGCCGCCCGACAAGcccccaaaaa                          >  2:17631707/1‑131 (MQ=255)
                                                                                                           aaCAGTCTCCGCCTGTAGCGGACCCTTCGAGGGGAGGGCCACGAAAAGGGCAGTAAAAGGACAGGAGAACCTCGTTCCAGATCTCTTTCCCCAACCTGTGCCATCCGAACCAGATCAGCCAGGCAACcccccacaaa                          >  2:4190532/1‑131 (MQ=255)
                                                                                                               ttctcCTCTTGGAGTGCCCTTTTTGCGGGCATCCCAGAGAGAGGGCCGGTACAGGTACAGCAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGCGCCGCCCGCCCAGGATCAGCCAGCCAATCCCCAGCAACAGGc                      <  1:2906754/136‑1 (MQ=255)
                                                                                                                tctcCGCCTGTAGCGGCCCTTTCGCGGGCATGCCCGCGACGAGGCCGGTACAGGTCCAGCAGAACTTAGCGCCAGCTCTCTTTCACCACCACGCGCCGCCCGCCCAGGATCAGCCAGCCAATCCCCAGCAACAGGCt                     >  2:17636850/1‑137 (MQ=255)
                                                                                                                  tcCGCCTGTAGTGGCCCTTTCGTGGGCATTCCCGCGAAAGGGCCGGTACAGGTACAGCAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGCGCCGCCCGCCCAGGATCAGCCAGCCAATCCCCAGCAACAGGCttt                   <  2:4261544/137‑1 (MQ=255)
                                                                                                                  tcCGCCTGGAGCGGCCCTTTCGCGGGCATGCCCGCGAAAGGGCCGGTACAGGTACAGCAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGCGCCGCCCGCCCAGGATCAGCCAGCCAATCCCCAGCAACAGGCttt                   <  2:2790443/137‑1 (MQ=255)
                                                                                                                       ttGTAGTGTCCTTTTCGCGGTCGGCCCCACGAAAGGGCCGGTACAGGTACAGCAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGCGCCGCCCGCCCAGGATCAGCCAGCCAATCCCCAGCAACAGGCTTTCGAcc              <  2:5455954/136‑1 (MQ=255)
                                                                                                                       ttGTAGCGCCCCTTTCTCGGGCATACCCGCGAAAGGGCCGGTACAGGTACAGCAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGCGCCGCCCGCCCAGGATCAGCCAGCCAATCCCCAGCAACAGGCTTTCGAcc              <  2:11625379/136‑1 (MQ=255)
                                                                                                                       ttGTAGCGCCCCTTCCGCGGGCATTCCCGCGAAAGGGCCGGTACAGGTACAGCAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGCGCCGCCCGCCCAGGATCAGCCAGCCAATCCCCAGCAACAGGCTTTCGAcc              <  2:9204303/136‑1 (MQ=255)
                                                                                                                       cggTATCGGCCCTTCCGGGGGCATTCCCGTGAAAGGGCCGGTACAGGTACAGCAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGCGCCGCCCGCCCAGGATCAGCCAGCCAATCCCCAGCAACAGGCTTTCGAcc              <  2:10120390/135‑1 (MQ=255)
                                                                                                                       cTGTAGTGCCCTTTTCGTGGACGTCCCCCAGAAAGGTCCGGTACAGGTACAGCAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGCGCCGCCCGCCCAGGATCAGCCAGCCAATCCCCAGCAACAGGCTTTCGAcc              <  1:6335334/137‑1 (MQ=255)
                                                                                                                       cTGTAGCGGCCCTTTCGCGGGCATACCCGCGAAAGGGCCGGTACAGGTACAGCAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGCGCCGCCCGCCCAGGATCAGCCAGCCAATCCCCAGCAACAGGCTTTCGAcc              <  2:15621682/137‑1 (MQ=255)
                                                                                                                            gCGGCCCTTTCGCGGGCATGCCCGCGAAAGGGCCGGTACAGGTACAGCAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGCGCCGCCCGCCCAGGATCAGCCAGCCAATCCCCAGCAACGGGCTTTCGACCACaaa         <  2:6876422/137‑1 (MQ=255)
                                                                                                                               cgtcttttttGCGGGCATGCCCGCGAAAGGGCCGGTACAGGTACAGCAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGCGCCGCCCGCCCAGGATCAGCCAGCCAATCCCCAGCAACAGGCTTTCGACCACAAAg        <  2:13391907/131‑1 (MQ=255)
                                                                                                                                     tttCGCGGACATGCCCGCGAAAGGGCCGGTACAGGTACAGCAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGCGCCGCCCGCCCAGGATCAGCCAGCCAATCCCCAGCAACAGGCTTTCGACCACAAAGGCCAGc  <  2:13674430/137‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                      
TGGCAGACCGGGGTGGCGTTGGTCGATGAGGTGGTTGGGGTTTGCCGGCAGGTGTTGGCGCGGTATGCGCGGGATGTGGGTGGGCAGCGGATCGTGCTGGTCTGACCAACAGTCTCCGCCTGTAGCGGCCCTTTCGCGGGCATGCCCGCGAAAGGGCCGGTACAGGTACAGCAGAACTTAGCGCCAGCTCTCTTTCACCACCCTGCGCCGCCCGCCCAGGATCAGCCAGCCAATCCCCAGCAACAGGCTTTCGACCACAAAGGCCAGC  >  AE015451_phaM1/5199521‑5199788

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: