Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451_phaM15,814,1780CT28.2% 409.2 / 47.2 223intergenic (+88/+41)PP_5092 CDS (putative exported NLP_P60 protein)/pilT CDS (twitching motility protein)–/–
Reads supporting (aligned to +/- strand):  ref base C (138/20);  new base T (0/62);  total (141/82)
Fisher's exact test for biased strand distribution p-value = 1.59e-37
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 227 total aligned reads displayed.

CAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGT‑AGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGG‑CCTTTGCAGGTACTACAAG‑CCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTCCGCCGCCCCGCCCTGCATCACTGAAG  >  AE015451_phaM1/5814042‑5814313
                                                                                                                                         |                                                                                                                                         
cAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCTTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGT‑AGCGCGTCGCTCTCGCGCGTAGCTCACATAcc                                                                                                                                           >  2:2652685/1‑137 (MQ=255)
cAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGT‑AGCGCGTCGCTCTCGCGCGCAGCTCACATAcc                                                                                                                                           <  2:7384639/137‑1 (MQ=255)
 aaGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGT‑AGCGCGTCGGTCTCGCGCGCAGCTCACATACCt                                                                                                                                          >  2:1530046/1‑137 (MQ=255)
  aGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGT‑AGCGCGTCGCTCTCGCGCGCAGCTCACATACCtg                                                                                                                                         <  1:11479076/137‑1 (MQ=255)
  aGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGT‑AGCGCGTCGCTCTCGCGCGCAGCTCACATACCtg                                                                                                                                         <  1:13394868/137‑1 (MQ=255)
           tgcCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGT‑AGTGCGTCGCTTTTGCGCGCAGCTTACATACCTGTGGGAGCgg                                                                                                                                >  2:10489748/1‑137 (MQ=255)
           tgcCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGT‑AGCGCGTCGCTTTCGCGCGCAGCTCAAATACCTGTGGGAGCgg                                                                                                                                >  2:16182608/1‑137 (MQ=255)
           tgcCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCACTCCCACAGGT‑AGCGCGGCGCTCCCACGCGCCGATCCCACACCTGTGGGAGCgg                                                                                                                                >  1:18883709/1‑137 (MQ=255)
                aaaGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCGCAGGT‑AGCGCGTCGTTCTCGCGCGCAGTTCACATAACTTCGGGGGCgggtaag                                                                                                                           >  2:17021130/1‑133 (MQ=255)
                    agcagcagcTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCAAAAGG‑AGCGCGTCGCTCTCGCGTACAGCTCACAGGCCTGTGGGAGCGGGAGAGCCcg                                                                                                                       >  2:7477215/1‑137 (MQ=255)
                     gcaccagcTGAGTCGGACCAAGGGGGAAGCGGATCGGCGCTATCCGGTACCGCCCCCTCGGCGGGCACGCCCCCCCCCACAGGT‑AAAGCGTCGCTTCCGCGCGCAGCCCACATCCTTGGGGGAGCGGCCGCGCCcgc                                                                                                                      <  2:5182167/137‑1 (MQ=255)
                          agcTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCCCCCACAGGT‑AGCGCGTCGCTCTTGCGCGCAGCTCCCACACCTGTGGGAGCGGGTGCGCCCCCGaaaa                                                                                                                 >  2:1950661/1‑136 (MQ=39)
                          agcTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCATTCGCGGGCACGCCCGCCCCCACAGGT‑AGCGCGTCGCTCTCGCGCGCCATTCACATACCTGTGGGAGTGGGCGCGACCGCaaaaa                                                                                                                 >  2:11314786/1‑136 (MQ=255)
                                   tGGCCAAGTGTGAATCTGGTTGGTGCTATCTGTGCCCTTCCCCTTCGCGGGCACCCCCGCTCCCAAGGGG‑ATCGCGTCGCTTCTGCAAGAAGCCAACATACCTGGGGGAGCGGCCGCGCCCGCGAAAGGG‑CCTTTGc                                                                                                       <  2:13499187/137‑1 (MQ=255)
                                   tGGCAAAGTGTGAAGCGGATTGGCGCTGTCTGGGACCGCCCCCTTCGCGGACACGCCCTCCCCCACGGGG‑GGCCTGTCCCTCTCACGCGCAGCCCACATATTTGTGGGAGCGGCCGCGCCCGCGAAAGGG‑CCTTTGc                                                                                                       <  2:14133044/137‑1 (MQ=255)
                                   tGGCAAAGCGTCAATCTAGTCGACGCTATCCGGCACCGTCCCTTCCGCGGACACCCCCTCCCCCACAGGT‑AGTTCTTCGCCCTCGCGCCCAGCCCACATACCTGGGGGGGCGGTCGCGCCCGCGAAAGGG‑CCTTTGc                                                                                                       <  2:18271780/137‑1 (MQ=255)
                                   gggCCAAGCGTGAAGGTGGTTGGCGCTCTCTGTGACTTTCCCTTTTGTGGGCACGCCCGTTCCCACGGGG‑ATGGTGTCGTTCTCGAGAGCGGCCAACATACTTGTGGGAGTGGGCGCGCCCGCGAAAGGG‑CCTTTGc                                                                                                       <  2:17503361/136‑1 (MQ=255)
                                   gggCAAAGCGGGAGGCTGATAGGCGCTATCTGGGACCGTCCCCGTCGTGGGCACTCCCTCTCCCACAGGT‑AGCGTGTCGCCCTCGCGCGCAGCCCACATACCTGGGGGGGTGGGCGCGCCCGCGAAAGGG‑CCTTTGc                                                                                                       <  2:17119326/136‑1 (MQ=255)
                                      ccAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCCCACCCTCTCCCACAAGT‑AGTGCGTTGCTCTCGAGGGGAAATTCAATACCTGTGGGGGAGGGGGCGGCCGCGAAAAGG‑ACTTTTCAgg                                                                                                    >  1:8161461/1‑137 (MQ=255)
                                      ccAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCACTCCCCCCGGT‑CGCACGTTGCTCCTTCTCGCAGGTTCCATTCATGTGGGAGCGGGGTTTCCCGCAAAAGGGACCttttagg                                                                                                     >  1:9643526/1‑133 (MQ=255)
                                              tAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGT‑AGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGG‑CCTTTGCAAGTACTACaa                                                                                            >  2:945596/1‑137 (MQ=255)
                                               aaGCTGGTCGGCGCTACCCGGTCCTGGCCCCTCCGCGGAAGCTCCCTCTCCCCCATGG‑ACCCCGTCTCTCTCGCGCACACCCCACCCACCTGTGGGAGCGTCCGCGCCCGCGAAAGGG‑CCTTTGCAGGTACTACAAg                                                                                           <  2:18752618/137‑1 (MQ=255)
                                                  cTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACACCCGCTCCCACCGGT‑GGCGTGTCGCTCTCACGCGCAACTTCAATACCTGGGGGAGTGGGCGAGCCCACTAAAGGG‑ACTTTGTCAGCACTACAAgact                                                                                        >  1:14706318/1‑134 (MQ=255)
                                                        ggCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCCCAGGT‑CGCGCGTCGCTCTCGCGCGCAAATCACATACATGTGGGAGCGGGTGCGCCCATGCAAAGG‑CCCTTGCCAGTACCCCAAA‑CCTTTTATc                                                                                 >  2:5385349/1‑137 (MQ=255)
                                                        ggCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCCCAGGT‑AGCGCGTTGCTCTCGCGCGCAGCTCAAATAACCGTGGGAGCGGGCGCGACCGCGAAAGGG‑CCTTTGCAGGTACTACAAA‑CCTTTGacc                                                                                 >  2:10785307/1‑135 (MQ=37)
                                                        ggCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCCCAGGT‑AGCGCGTCGCTCTCGAGCGTAACTCCAATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGG‑CCTTTGCAGGTACTACAAG‑TCTTTGATc                                                                                 >  2:13255315/1‑137 (MQ=37)
                                                        ggCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGT‑AGCGAGCCGCTCTCGCGGGCAGCTCACATAACCGTGGGAACGGGTGCGCACGCGACAAGG‑CCCCTACAGGGACTACAAG‑ACGTTGATc                                                                                 >  2:7683595/1‑137 (MQ=25)
                                                        ggCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACACGT‑AAAGAGTCACTCTCACGTGCAGCTTCCACACCTGTGGGGGCGGGGATGCCCACTAAAGGGACCTTTCCA‑GCCCTCCAAG‑CCTTTTATc                                                                                 >  1:12194828/1‑137 (MQ=12)
                                                                   gcaCCGCCCCCTTCGCGGAAACACCTTCTCACAAGGGG‑AGTGTTTTGCTCTCGCGCAAAGCTCACATCCTTGGGGGGGCGATCGCGCCCGCGAAAGGG‑CCTTTGCAGGTACTACAAG‑CCTTTGATCAACCCTCGCTc                                                                      <  1:13989645/135‑1 (MQ=255)
                                                                       cgttcccTTCGCGGGCCCTCCCGTCCCCACAGGT‑GGTGCTTTGCTTTCGCGCGCACCTCACGTCCTTGGGGGAGTGGGCGCGCCCGCGAAAGGG‑CCTTTGCAGGTACTACAAG‑CCTTTGATCAACCCTCGCTCCCTc                                                                  <  2:11366624/133‑1 (MQ=255)
                                                                       cgcccctTTCGCGGACACGCCCTCCCCCACAGGT‑GGCACGTTGTTCTCGCGCGCGGCTCCCATGTTGGGGGGGGTGGCCGCGCCCGCGAAAGGG‑CCTTTGCAGGTACTACAAG‑CCTTTGATCAACCCTCGCTCCCTc                                                                  <  2:2639573/134‑1 (MQ=255)
                                                                       cGGCCCCTTGGGGGGAACACCCGCTCCCAAGGGG‑AGCGTGTTGCTCTCGCGAGCGGCTCACATACCGGTGGGGTCGGGCGCGCCCGCGAAAGGG‑CCTTTGCAGGTACTACAAG‑CCTTTGATCAACCCTCGCTCCCTc                                                                  <  2:12130100/137‑1 (MQ=255)
                                                                       cGGCCCCTTCGCGGGCCCGCCCGCTCCCACAGGT‑AGCGCATTTCTCTCGCGCGCAGCTTACCTACCTGTGGGAGCGGGCCCGCCCACGAAAGGG‑CCTTTGTCAGTCCTACAAA‑CCTTTTGTCAACCCTCGCTccccc                                                                  >  2:7805247/1‑135 (MQ=255)
                                                                            ctttCGCGGGCACTCCCGCTCCAAAAGGG‑ATGGCTTTGCTCTCGAACGAATCTAACTAACTGGGGGGAGCGGCCGCGCCCGCGAAAGGG‑CCTTTGCAGGTACTACAAG‑CCTTTGATCAACCCTCGCTCCCTCAGCCt                                                             <  1:14845148/135‑1 (MQ=255)
                                                                            ccTTCGCGGGTACGCCCGCTCCCAAAGGT‑AGCGCGTCGCCCTCGCGAGAAGCTCACACGCCTGTGGGAGCGGGCGAACCCACGAAAGGG‑CCTTTTCCAGTACCAAAAGACCTTTGCT‑AACCCTCGCTCCCTCAGccc                                                             >  2:10141372/1‑136 (MQ=255)
                                                                            ccTTCGCGGGCACGCCCGCTCCCACAGGT‑AGCGCGTTGCTTTTGCGCGCAGATCACATACCTGGGGGAGAGGGGACGCCCGCGAAAGGG‑CCTTTGCAAGGGCCACAAG‑ACTCTTACCAACCCTCGCTCCCCCAgtcc                                                             >  2:15992062/1‑134 (MQ=19)
                                                                            ccTTCGCGGGCAAGCCCGCTCCCACAGGT‑AGCGCGTCACTCTCGGACGCAGATCAAATAACTGTGGGAGGGGGCGCGCCCATGAAAGGG‑CCCTTGCAGGTACTCAAAA‑CCTTTGATCAACCCCCGATCCATCAGCCt                                                             >  2:9109440/1‑137 (MQ=255)
                                                                              tttgggggCACTCCCTCTCCCA‑AGGTGAGCGCGTTCCTCTTGCGCGCACCCCACTCACTTGTGGGGGTGGCCGCGCCCGCGAAAGGG‑CCTTTGCAGGTACTACAAG‑CCTTTGATCAACCCTCGCTCCCTCAGCCTcc                                                           <  2:117611/132‑1 (MQ=255)
                                                                              tttgggggCACGCCCGCCCCCCCGGGG‑GGCGTTTCTCTTTCGCGCGCGTCTCACGTCCTTGTGGGGGTGGCCGCGCCCGCGAAAGGG‑CCTTTGCAGGTACTACAAG‑CCTTTGATCAACCCTCGCTCCCTCAGCCTcc                                                           <  2:7169928/132‑1 (MQ=255)
                                                                              ttCGCGGGAACGCCCTTTCCCAAAGGG‑GGCGCTTCGCTCTCGTATGAAGCTAACATACTGGGGGGAGAGGGCGCGCCCGCGAAAGGG‑CCTTTGCAGGTACTACAAG‑CCTTTGATCAACCCTCGCTCCCTCAGCCTcc                                                           <  2:18514216/137‑1 (MQ=255)
                                                                                       acacccTCTCCCACGGGT‑AGCGCTTCGTCCTCGCGCGCATCCCACCTCTTTGTGGGAGCGGCCGCGCCCGCGAAAGGG‑CCTTTGCAGGTACTACAAG‑CCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGccc                                                  <  1:4229264/134‑1 (MQ=255)
                                                                                            cGCTCCCACAGGT‑AGCAGGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCAAAAAGG‑CCTTTGCAGGTACTACAAG‑CCTTTGATCAACCCTCGCTCCCCCAGCCTCCGCAACGCCCCCTcc                                             >  1:1720586/1‑137 (MQ=255)
                                                                                             tcTCCCACATGG‑AGCGTTTCGTTCTTTTTCGCACCCCACATCTCTGGGGGGGCGGCCGCGCCCGCGAAAGGG‑CCTTTGCAGGTACTACAAG‑CCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCa                                            <  1:12077974/136‑1 (MQ=255)
                                                                                             ccccccACGGGT‑AGCGCTCCTCTCTCGCACGCACCCCACGTCCTTGGGGGGGTGTCCGCGCCCGCGAAAGGG‑CCTTTGCAGGTACTACAAG‑CCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCa                                            <  1:10433830/134‑1 (MQ=255)
                                                                                                  cacaGGT‑AGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGG‑CCTTTGCAGGTACTACAAG‑CCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTc                                       >  2:2397824/1‑137 (MQ=255)
                                                                                                      ggT‑AGTGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGG‑CCTTTGCAGGTACTACAAG‑CCTTTGACCACCCCTCACTCCCTCAGCCCCCGCAACGCCCCCTCCATCGTCAGCa                                   >  2:6486718/1‑137 (MQ=255)
                                                                                                             gcgTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGG‑CCTTTGCAGGTACTACAAG‑CCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCt                             >  2:5713146/1‑137 (MQ=255)
                                                                                                             gcgTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGG‑CCTTTGCAGGTACTACAAG‑CCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCt                             >  1:6211018/1‑137 (MQ=255)
                                                                                                             gcgTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGG‑CCTTTGCAGGTACTACAAG‑CCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCt                             >  2:1171353/1‑137 (MQ=255)
                                                                                                             gcgTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCACGAAAGGG‑CCTTTGCAAGTCCTACAAG‑CCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCCCCGCCCCCCCCATCGTCAGTCTCCCCt                             >  1:9784208/1‑137 (MQ=255)
                                                                                                             gcgTCGCTCTCGCGCGCAGCTCACATACCCGTGGGAGCGGGCGCGCCCGCGAAAGGG‑CCTTTGCAGGTACTACAAG‑CCTTTGATCAAACCTCGATCCCTCAGCCTCCGCAACGCCCCCTCCATCATCAGCATCCCCt                             >  2:13014613/1‑137 (MQ=255)
                                                                                                             gcgTCACTCTCTCGCGCAGCTCACAGACCTGTGGGAGTGGGCGCGCCCGCGAAAGGG‑CCTTTGCAGGTACTACAAG‑CCTTTGATCAACCCTCGCTCCCTCATCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCt                             <  2:13742954/137‑1 (MQ=255)
                                                                                                               gTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGG‑CCtt                                                                                                          <  1:8486529/59‑1 (MQ=255)
                                                                                                               gTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGG‑CCtt                                                                                                          >  2:8486529/1‑59 (MQ=255)
                                                                                                               gTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGG‑CCTTTGCAGGTACTACAAG‑CCTTTGATCAACCCTCGCTCCCTCCGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTcc                           >  1:5394136/1‑137 (MQ=255)
                                                                                                               gTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGG‑CCTTTGCAGGTACTACAAG‑CCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTcc                           >  1:7424652/1‑137 (MQ=255)
                                                                                                               gTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGG‑CCTTTGCAGGTACTACAAG‑CCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTcc                           >  2:18014198/1‑137 (MQ=255)
                                                                                                               gTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGG‑CCTTTGCAGGTACTACAAG‑CCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTcc                           >  1:18690063/1‑137 (MQ=255)
                                                                                                               gTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGG‑CCTTTGCAGGTACTACAAG‑CCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTcc                           >  1:201812/1‑137 (MQ=255)
                                                                                                               gTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGG‑CCTTTGCAGGTACTACAAG‑CCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTcc                           >  1:17570402/1‑137 (MQ=255)
                                                                                                               gTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGG‑CCTTTGCAGGTACTACAAG‑CCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTcc                           >  1:3844728/1‑137 (MQ=255)
                                                                                                               gTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGG‑CCTTTGCAGGTACTACAAG‑CCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTcc                           >  2:16383658/1‑137 (MQ=255)
                                                                                                               gTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGG‑CCTTTGCAGGTACTACAAG‑CCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTcc                           >  1:654047/1‑137 (MQ=255)
                                                                                                               gTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGG‑CCTTTGCAGGTACTACAAG‑CCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTcc                           >  1:11069378/1‑137 (MQ=255)
                                                                                                               gTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGG‑CCTTTGCAGGTACTACAAG‑CCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTcc                           >  2:1450502/1‑137 (MQ=255)
                                                                                                               gTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGG‑CCTTTGCAGGTACTACAAG‑CCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTcc                           >  1:858398/1‑137 (MQ=255)
                                                                                                               gTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGG‑CCTTTGCAGGTACTACAAG‑CCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTcc                           >  2:2996300/1‑137 (MQ=255)
                                                                                                               gTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGG‑CCTTTGCAGGTACTACAAG‑CCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTcc                           >  2:4761798/1‑137 (MQ=255)
                                                                                                               gTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGG‑CCCTTGCAGGTACTACAAG‑CCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTcc                           >  1:3405836/1‑137 (MQ=255)
                                                                                                               gTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAAGG‑CCTTTGCAGGTACTACAAG‑CCTTTGATCCACCCTCGCTCCCTCCACCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTcc                           >  2:3528374/1‑137 (MQ=255)
                                                                                                                 cGCTTTCGCGCGGAGCTCACGTCTTTGGGGGGGTGGTCGCGCCCGCGAAAGGG‑CCTTTGCAGGTACTACAAG‑CCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTccgc                         <  1:8965649/137‑1 (MQ=255)
                                                                                                                 cGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGT‑CCTTTGCAGGTACTACAAG‑CCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCCTCGTCAGCATCCCCTccgc                         >  2:6185247/1‑137 (MQ=255)
                                                                                                                 cGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGG‑CCTTTGCAGGTACTACAAG‑CCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTccgc                         >  1:12478284/1‑137 (MQ=255)
                                                                                                                 cGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGG‑CCTTTGCAGGTACTACAAG‑CCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAAAGCCCCCTCCATCGTCAGCATCCCCTccgc                         >  1:14019717/1‑137 (MQ=255)
                                                                                                                 cGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGG‑CCTTTGCAGGTACTAAAAG‑CCTTTGATCCACCCTCGCTCCCTCCGCCTCCGCAACGCCCCCTCCATCGCCAGCATCCCCTccgc                         >  2:8969265/1‑137 (MQ=255)
                                                                                                                  tctctCGAACGAAACTCACCCACCTGGGGAATCGGTCGCGCCCGCGAAAGGG‑CCTTTGCAGGTACTACAAG‑CCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTccgcc                        <  2:18112195/136‑1 (MQ=255)
                                                                                                                  gttctcGCGCGCAGCTCACATCCTTTGGGGGGGGGTCGCGCCCGCGAAAGGG‑CCTTTGCAGGTACTACAAG‑CCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTccgcc                        <  2:8530389/135‑1 (MQ=255)
                                                                                                                  gCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGG‑CCTTTGCAGGTACTACAAG‑CCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTccgcc                        >  1:451321/1‑137 (MQ=255)
                                                                                                                  gCTCTCGAGAGAAGCCCACATATTTGGGGGGGGGGGCGCGCCCGCGAAAGGG‑CCTTTGCAGGTACTACAAG‑CCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTccgcc                        <  2:10028927/137‑1 (MQ=255)
                                                                                                                   ctctCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGG‑CCTTTGCAGGTACTACAAG‑CCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTccgccg                       >  2:8334973/1‑137 (MQ=255)
                                                                                                                   ctctCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGG‑CCTTTGCAGGTACTACAAG‑CCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTccgccg                       >  2:9745728/1‑137 (MQ=255)
                                                                                                                   ctctCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGG‑CCTTTGCAGGTACTACAAG‑CCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTccgccg                       >  1:12927697/1‑137 (MQ=255)
                                                                                                                    tctcGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGG‑CCTTTGCAGGTACTACAAG‑CCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTCcgccgc                      >  1:15877152/1‑137 (MQ=255)
                                                                                                                     cttgcgcgaAGCCAACATACTTGGGGGAGTGGACGCGCCCGCGAAAGGG‑CCTTTGCAGGTACTACAAG‑CCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTCCgccgcc                     <  2:18511537/134‑1 (MQ=255)
                                                                                                                     ctcatgcgcAGCCCACATCCTGGGGGGAGCGTCCGCGCCCGCGAAAGGG‑CCTTTGCAGGTACTACAAG‑CCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTCCgccgcc                     <  1:5980879/132‑1 (MQ=255)
                                                                                                                      ccgcgcaCATCCCACATCCTTGGGGGGGCTTCCGCGCCCGCGAAAGGG‑CCTTTGCAGGTACTACAAG‑CCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTCCGccgccc                    <  1:10489748/136‑1 (MQ=255)
                                                                                                                       cgcgcgCAGCTCAAATACTTGTGGGAGCGGCCGCGCCCGCGAAAGGG‑CCTTTGCAGGTACTACAAG‑CCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTCCGccgcccc                   <  1:18330124/137‑1 (MQ=255)
                                                                                                                       cgcgcgCAGCCAACATATTTGGGGGAGCGTTCGCGCCCGCGAAAGGG‑CCTTTGCAGGTACTACAAG‑CCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTCCGccgcccc                   <  1:15174456/137‑1 (MQ=255)
                                                                                                                       cgcgagAAGCTCACATACTGGTGGGAGCGGTCGCGCCCGCGAAAGGG‑CCTTTGCAGGTACTACAAG‑CCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTCCGccgcccc                   <  2:15642104/137‑1 (MQ=255)
                                                                                                                        gcgcgaAGCTCAAGCACCGGGGGGAGCGGGCGCGCCCGCGAAAGGG‑CCTTTGCAGGTACTACAAG‑CCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTCCGccgccccg                  <  2:16642579/137‑1 (MQ=255)
                                                                                                                          gcgcAGCCCACATACCTGGGGGAGCGGGCGCGCCCGCGAAAGGG‑CCTTTGCAGGTACTACAAG‑CCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTCCGccgccccgcc                <  1:17361287/137‑1 (MQ=255)
                                                                                                                            gcAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGG‑CCTTTGCAGGTACTACAAG‑CCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTCCGCCGCCCCGCCCt              >  2:14881205/1‑137 (MQ=255)
                                                                                                                               tcTCACATCCTTGGGGGAGCGTCCGCGCCCGCGAAAGGG‑CCTTTGCAGGTACTACAAG‑CCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTCCGCCGCCCCGCCCTGCa           <  1:16762942/136‑1 (MQ=255)
                                                                                                                               gttaacaTCCTTGGGGGGGCGGCCGCGCCCGCGAAAGGG‑CCTTTGCAGGTACTACAAG‑CCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTCCGCCGCCCCGCCCTGCa           <  2:12294970/133‑1 (MQ=255)
                                                                                                                               gCTCACATATCGGGGGGGGCGGGCGCGCCCGCGAAAGGG‑CCTTTGCAGGTACTACAAG‑CCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTCCGCCGCCCCGCCCTGCa           <  2:12008349/137‑1 (MQ=255)
                                                                                                                               gCTCAAATACTGGGGGGGGGGGTCGCGCCCGCGAAAGGG‑CCTTTGCATGTACTACAAG‑CCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTCCGCCGCCCCGCCCTGCa           <  1:2964107/137‑1 (MQ=255)
                                                                                                                                cTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGG‑CCTTTGCAGGTACTACAAG‑CCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCCTCGTCAGCATCCCCTCCGCCGCCCCGCCCTGCAt          >  2:4227026/1‑137 (MQ=255)
                                                                                                                                cTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGG‑CCTTTGCAGGTACTACAAG‑CCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTCCGCCGCCCCGCCCTGCAt          >  1:5633017/1‑137 (MQ=255)
                                                                                                                                cTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGG‑CCCTTGCAGGTACTACAAG‑CCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTCCGCCGCCCCGCCCTGCAt          >  1:4931465/1‑137 (MQ=255)
                                                                                                                                   acaTACCTGTGGGAGCGGGCGCGCCCGCGAAAGGG‑CCTTTGCAGGTACTACAAG‑CCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTCCGCCGCCCCGCCCTGCATCAc       >  2:5639749/1‑137 (MQ=255)
                                                                                                                                         |                                                                                                                                         
CAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGT‑AGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGG‑CCTTTGCAGGTACTACAAG‑CCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTCCGCCGCCCCGCCCTGCATCACTGAAG  >  AE015451_phaM1/5814042‑5814313

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