Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451441,3260G.25.0% 136.6 / 47.4 57coding (207/1173 nt)bioF8‑amino‑7‑oxononanoate synthase
Reads supporting (aligned to +/- strand):  ref base G (20/22);  new base . (6/8);  total (26/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.51e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AGCACGCCGACGAAGTGCACCAGGCAGCCATCGAGCAAGCGCTGGTCGAGCAGCGCAGCAGCGAGATGTTCTACGACGCCGCGACAGCCTGACATGGCTTTCGACCTGGCAGCGCGCCTGGCCGAACGGCGCGCCGCAGACCTGTATCGGCAGCGGCCATTGCTGGAAAGCCCGCAAGGGCCGGAAGTGGTGGTCGACGGCCAGCGGTTGCTGGCCTTCTGCAGCAATGACTACCTGGGCCTGGCCAACCACCCCGAGGTGATCGCCGCCTGGCAGGCCGGCGCCGAGCGCTGGGGTGTGGGGGGGGGGGCCTCGCACCTGGTAGTTGGCCACAGCACCCCGCACCATCAGGTTGAAGAGGCACTGGCTGAACTCACCGGCCGCCCGCGCGCCCTGTTGTTCTCCACCGGCTATATGGCCAATCTGGGCGCCATCACCGCGCTGGTGGGGCAGGGCGACACGGTATTGCAAGACCGCCTCAACCATGCCTCGTTGCTGGATGGCGGGTTGCTCAGTGGTGCCCGCTTCAACCGCTACCTGCACAACGACCCGGCCAGCCTGGCCAGCCGCCTGGACAAGGCGGTCGGCAATACCCTGGTG  >  AE015451/441027‑441626
                                                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                                                              
agCACGCCGACGAAGTGCACCAGGCAGCCATCGAGCAAGCGCTGGTCGAGCAGCGCAGCAGCGAGATGTTCTACGACGCCGCGACAGCCTGACATGGCTTTCGACCTGGCAGCGCGCCTGGCCGAACGGCGCGCCGCAGACCTGTATCGGCAGCGGCCATTGCTGGAAAGCCCGCAAGGGCCGGAAGTGGTGGTCGACGGCCAGCGGTTGCTGGCCTTCTGCAGCAATGACTACCTGGGCCTGGCCAACCACCCCGAGGTGATCGCCGCCTGGCAGGCCGGCGCCGAGCGCTGGGGTGTgg                                                                                                                                                                                                                                                                                                              >  1:690609/1‑301 (MQ=255)
            aaGTGCACCAGGCAGCCATCGAGCAAGCGCTGGTCGAGCAGCGCAGCAGCGAGATGTTCTACGACGCCGCGACAGCCTGACATGGCTTTCGACCTGGCAGCGCGCCTGGCCGAACGGCGCGCCGCAGACCTGTATCGGCAGCGGCCATTGCTGGAAAGCCCGCAAGGGCCGGAAGTGGTGGTCGACGGCCAGCGGTTGCTGGCCTTCTGCAGCAATGACTACCTGGGCCTGGCCAACCACCCCGAGGTGATCGCCGCCTGGCAGGCCGGCGCCGAGCGCTGGGGTGGTgggggggggcccc                                                                                                                                                                                                                                                                                                   >  2:126164/1‑297 (MQ=255)
                    cAGGCAGCCATCGAGCAAGCGCTGGTCGAGCAGCGCAGCAGCGAGATGTTCTACGACGCCGCGACAGCCTGACATGGCTTTCGACCTGGCAGCGCGCCTGGCCGAACGGCGCGCCGCAGACCTGTATCGGCAGCGGCCATTGCTGGAAAGCCCGCAAGGGCCGGAAGTGGTGGTCGACGGCCAGCGGTTGCTGGCCTTCTGCAGCAATGACTACCTGGGCCTGGCCAACCACCCCGAGGTGATCGCCTCCTGGCAGGCCGGCGCCGAGCGCTGGGGGGTGGGGGGGGGGCCTCGCACCTgg                                                                                                                                                                                                                                                                                         >  2:309592/1‑301 (MQ=255)
                         aGCCATCGAGCAAGCGCTGGTCGAGCAGCGCAGCAGCGAGATGTTCTACGACGCCGCGACAGCCTGACATGGCTTTCGACCTGGCAGCGCGCCTGGCCGAACGGCGCGCCGCAGACCTGTATCGGCAGCGGCCATTGCTGGAAAGCCCGCAAGGGCCGGAAGTGGTGGTCGACGGCCAGCGGTTGCTGGCCTTCTGCAGCAATGACTACCTTGGCCTGGCCAACCACCCCGCGGTGATCGCCGCCTGGCCGGCCCGCGCCGAGCGCTGGGGTGTGGGGGGGGGGGCCTCGCACCTGTTAGt                                                                                                                                                                                                                                                                                     >  2:352247/1‑301 (MQ=255)
                         aGCCATCGAGCAAGCGCTGGTCGAGCAGCGCAGCAGCGAGATGTTCTACGACGCCGCGACAGCCTGACATGGCTTTCGACCTGGCAGCGCGCCTGGCCGAACGGCGCGCCGCAGACCTGTATCGGCAGCGGCCATTGCTGGAAAGCCCGCAAGGGCCGGAAGTGGTGGTCGACGGCCAGCGGTTGCTGGCCTTCTGCAGCAATGACTACCTGGGCCTGGCCAACCACCCCGAGGTGATCGCCGCCTGGCAGGCCGGCGCCGAGCGCTGGGGTGTGGGGGGGGGGCCTCGCACCTGGTAGtt                                                                                                                                                                                                                                                                                    >  1:377677/1‑301 (MQ=255)
                            cATCGAGCAAGCGCTGGTCGAGCAGCGCAGCAGCGAGATGTTCTACGACGCCGCGACAGCCTGACATGGCTTTCGACCTGGCAGCGCGCCTGGCCGAACGGCGCGCCGCAGACCTGTATCGGCAGCGGCCCTTTCGGGAAAGCCCGCAAGGGCCGGAAGGTGTGGTCGACGGCCAGCGGTTGGTGGCCTTCTGCCGCCATTACCACCTGGGCCGTGCCAAACACCCCGCTGTGCTCGCCGCCTGGTAGGCCGGCGCTCAGCGCGTGGGTGTGGGGGGGGGGGACTCGcgcgtgtgtgggtg                                                                                                                                                                                                                                                                                  >  2:13099/1‑288 (MQ=255)
                             aTCGAGCAAGCGCGGGGCGAGCAGCGCATGAGCGAGATGTTCTACGACGCCGCGACAGCCTGACATGGCGTTAGACCTGGCAGCGCGCCGGGCCGGACGGCGCGCCGCAGACCTGTATCGGCAGCGGCCATTGCTGGAAAGCCCGCAAGGGCCGGAAGTGGTGGTCGACGGCCAGCGGTTGCTGGCCTTCTGCAGCAATGACTACCTGGGCCTGGCCAACCACCCCGAGGTGATCGCCGCCTGGCAGGCCGGCGCCGAGCGCTGGGGTGTGGGGGGGGGGGCCTCGCACCTGGTAGTTGGc                                                                                                                                                                                                                                                                                 <  2:311121/301‑1 (MQ=255)
                                                 ggagcgcagcAGCGGGATGTTCTACGACGCCGCGACAGCCTGACATGGCTTTCGACCTGGCAGCGCGCCTGGCCGAACGGCGCGCCGCAGACCTGTATCGGCAGCGGCCATTGCTGGAAAGCCCGCAAGGGCCGGAAGTGGTGGTCGACGGCCAGCGGTTGCTGGCCTTCTGCAGCAATGACTACCTGGGCCTGGCCAACCACCCCGAGGTGATCGCCGCCTGGCAGGCCGGCGCCGGGCGCTGGGGTGTGGGGGGGGGGGCCTCGCACCTGGTAGTTGGCCACAGCACCCCGCACCATCa                                                                                                                                                                                                                                                             <  1:719616/299‑1 (MQ=255)
                                                   agcgcagcAGCGAGATGTTCTACGACGCCGCGACAGCCTGACATGGCTTTCGACCTGGCAGCGCGCCTGGCCGAACGGCGCGCCGCAGACCTGTATCGGCAGCGGCCATTGCTGGAAAGCCCGCAAGGGCCGGAAGTGGTGGTCGACGGCCAGCGGTTGCTGGCCTTCTGCAGCAATGACTACCTGGGCCTGGCCAACCACCCCGAGGTGATCACCGCCTGGCAGGCCGGCGCCGCGCGCTGGGGTGTGGGGGGGGGGGCCTCGCCCCTGGTAGATGGcctaggtcccccgcccgacacgg                                                                                                                                                                                                                                                           >  2:161868/1‑280 (MQ=255)
                                                          cagcGAGATGTTCTACGACGCCGCGACAGCCTGACATGGCTTTCGACCTGGCAGCGCGCCTGGCCGAACGGCGCGCCGCAGACCTGTATCGGCAGCGGCCATTGCTGGAAAGCCCGCAAGGGCCGGAAGTGGTGGTCGACGGCCAGCGGTTGCTGGCCTTCTGCAGCAATGACTACCTGGGCCTGGCCAACCACCCCGAGGTGATCGCCGCCTGGCAGGCCGGCGCCGAGGGCTGGGGTGTGGGGGGGGGGGCCTCGCACCGGGTGGTTGGCcaacggcccccgccccaaccggttgaagg                                                                                                                                                                                                                                                    >  2:6783/1‑274 (MQ=255)
                                                                 ttgtcgtccacgcccgcGACGTCCTGACATGGTTTCCGCCCTGCCAGCGCGCCTGGCCGAACGGCGCGCCGAGAACCGTTATCGGCAGCGGCCATTGCGGGAAAGCCCGCAAGGGCCGGAGGTGGTGGTCGACGGCCAGCGGTTGCTGGCCTTCTGCAGCAATGACTACCTGGGCCTGGCCAACCACCCCGAGGTGATCGCCGCCTGGCAGGCCGGCGCCGAGCGCTGGGGTGTGGGGGGGGGGGCCTCGCACCTGGTAGTTGGCCACAGCACCCCGCACCATCAGGTTGAAGAGGCACTg                                                                                                                                                                                                                                             <  1:283186/288‑1 (MQ=255)
                                                                     tCTACGACGCCGCGACAGCCTGACATGGCTTTCGACCTGGCAGCGCGCCTGGCCGAACGGCGCGCCGCAGACCTGTATCGGCAGCGGCCATTGCTGGAAAGCCCGCAAGGGCCGGAAGTGGTGGTCGACGGCCAGCGGTTGCTGGCCTTCTGCAGCAATGACTACCTGGGCCTGGCCAACCACCCCGAGGTGATCGCCGCCTGGCAGGCCGGCGCCGAGCGCTGGGGTGTGGGGGGGGGGCCTCGCACCTGGTAGTTGGCCACAGCACCCCGCACCATCAGGTTGAAGAGGCACTGGCTGa                                                                                                                                                                                                                                        <  1:132853/301‑1 (MQ=255)
                                                                     ggcgcggcgccgcGACGGCCTGACCTGGCTTTCGACCTGGCAGCGCGCCTGGCCGAACGGCGCGCCGCAGACCTGTATCGGGAGCGGCCATTGCGGGAAAGCCCGCAAGGGCCGGAAGTGGTGGTCGACGGCCAGCGGTTGCTGGCCTTCTGCAGCAATGACTACCTGGGCCTGGCCAACCACCCCGGGGTGATCGCCGCCTGGCAGGCCGGCGCCGAGCGCTGGGGCGTGGGGGGGGGGGCCTCGCACCTGGTAGTTGGCCACAGCACCCCGCACCATCAGGTTGAAGAGGCACTGGCTg                                                                                                                                                                                                                                         <  1:4342/294‑1 (MQ=255)
                                                                         cgacgCCGCGACAGCCTGACATGGCTTTCGACCTGGCAGCGCGCCTGGCCGAACGGCGCGCCGCAGACCTGTATCGGCAGCGGCCATTGCTGGAAAGCCCGCAAGGGCCGGAAGTGGTGGTCGACGGCCAGCGGTTGCTGGCCTTCTGCAGCAATGACTACCTGGGCCTGGCCAACCACCCCGAGGTGATCGCCGCCTGGCAGGCCGGCGCCGAGCGCTGGGGTGTGGGGGGGGGCCTCGCACCTGGTAGTTGGCCACAGCACCCCGCACCATCAGGTTGAAGAGGCACTGGCTGAACTcc                                                                                                                                                                                                                                   >  1:650499/1‑300 (MQ=255)
                                                                            cgccgcGACAGCCTGACATGGCTTTCGACCTGGCAGCGCGCCTGGCCGAACGGCGCGCCGCAGACCTGTATCGGCAGCGGCCATTGCTGGAAAGCCCGCAAGGGCCGGAAGTGGTGGTCGACGGCCAGCGGTTGCTGGCCTTCTGCAGCAATGACTACCTGGGCCTGGCCAACCACCCCGAGGTGATCGCCGCCTGGCAGGCCGGCGCCGAGCGCTGGGGTGTGGTGGGGGGGGCCTCGCACCTTGAAGTTGGCAACAGCACCCCGCACCTTCAGTTTCATGAGCCACTggtggacctccc                                                                                                                                                                                                                                  >  2:375838/1‑291 (MQ=255)
                                                                             ggcgcgACAGCCTGACATGGGCTTGCGGCTGGGAGGGGCGCTGGCCGGACGGGGGGCCGCAGACCTGTATCGGCAGGGGCCATTGCGGGAAAGCCCGCAAGGGCGGGAAGGGGTGGTCGACGGCCAGCGGTTGCTGGCCTTCTGCAGCAATGACTACCTGGGCCTGGCCAACCACCCCGAGGTGATCGCCGCCTGGCAGGCCGGCGCCGAGCGCTGGGGTGTGGGGGGGGGGGCCTCGCACCTGGTAGTTGGCCACAGCACCCCGCACCATCAGGTTGAAGAGGCACTGGCTGAACTCAcc                                                                                                                                                                                                                                 <  2:248436/299‑1 (MQ=255)
                                                                                             aTGGCTTTCGACCTGGCAGCGCGCCTGGCCGAACGGCGCGCCGCAGACCTGTATCGGCAGCGGCCATTGCTGGAAAGCCCGCAAGGGCCGGAAGTGGTGGTCGACGGCCAGCGGTTGCTGGCCTTCTGCAGCAATGACTACCTGGGCCTGGCCAACCACCCCGAGGTGATCGCCGCCTGGCAGGCCGGCGCCGAGCGCTGGGGTGTGGGGGGGGGGGCCTCGCACCTGGTAGTTGGGCACAGCACCCCGCACCATCAGGTTGAAGAGGCACTGGCTGAACTCACCGGCCGCCCGCGCGccc                                                                                                                                                                                                                 >  1:381499/1‑301 (MQ=255)
                                                                                                ggCTTTCGACCTGGCAGCGGGCCTGGCCGAACGGCGGGCGGCCGACCTGTATCGGCAGCGGCCATTGCTGGAAAGCCCGCAAGGGCCGGAAGTGGTGGTCGACGGCCAGCGGTTGCTGGCCTTCTGCAGCAATGACTACCTGGGCCTGGCCAACCACCCCGAGGTGATCGCCGCCTGGCAGGCCGGCGCCGAGCGCTGGGGTGTGGGGGGGGGCCTCGCACCTGGTAGTTGGCCACAGCACCCCGCACCATCAGGTTGAAGAGGCACTGGCTGAACTCACCGGCCGCCCGCGCGCCCtgtt                                                                                                                                                                                                             <  2:650227/301‑1 (MQ=255)
                                                                                                      cGACCTGGCAGCGCGCCTGGCCGAACGGCGCGCCGCAGACCTGTATCGGCAGCGGCCATTGCTGGAAAGCCCGCAAGGGCCGGAAGTGGTGGTCGACGGCCAGCGGTTGCTGGCCTTCTGCAGCAATGACTACCTGGGCCTGGCCAACCACCCCGAGGTGATCGCCGCCTGGCAGGCCGGCGCCGAGCGCTGGGGTGTGGGGGGGGGGGCCTCGCACCTGGTAGTTGGCCACAGCCCCCCGCCCCATCAGGTTGAAGAGGCACTGGCTGAACTCACCGGCCGCCCCCGCGCCCTGTTGTtc                                                                                                                                                                                                         >  2:193329/1‑301 (MQ=255)
                                                                                                       gACCTGGCAGCGCGCCTGGCCGAACGGCGCGCCGCAGACCTGTATCCGCAGCGGCCATTGCTGGAAAGCCCGCAAGGGCCGGCAGTGGTGGTCGACGGCCAGCGGTTGCTGGCCTTCTGCAGCAATGACTACCTGGGCCTGGCCACCCACCCCGAGGTGATCGCCGCCTGGCAGGCCGGCCCCGCGCGCTCGGGTGTGGGGGGGGGGGCCTCCCACCTGTTAGTTGGCCCCAGCCCCCCCCACCATCCGGTTGATGAGCCACGGGCTGAACCCACCGGCCGCCCGCCCGCCCTGTTGTtcc                                                                                                                                                                                                        >  2:582051/1‑300 (MQ=255)
                                                                                                               cacgcggcCTGGCAGAACGGGGCCCCCCAGACCTGGATCGGCAGGGGCCATTGCTGGAAAGCCCGCAAGGGCCGGAAGTGGGGGTCGACGGCCAGCGGTTGCTGGCCTTCTGCAGCAATGACTACCTGGGCCTGGCCAACCACCCCGAGGTGATCGCCGCCTGGCAGGCCGGCGCCGAGCGCTGGGGTGTGGGGGGGGGGGCCTCGCACCTGGTAGTTGGCCACAGCACCCCGCACCATCAGGTTGAAGAGGCACTGGCTGAACTCACCGGCCGCCCGCGCGCCCTGTTGTTCTCCACCgg                                                                                                                                                                                                 <  1:323029/295‑1 (MQ=255)
                                                                                                                      ccgggCCGAACGGCGCGCCGCAGACCTGTATCGGCAGCGGCCATTGCTGGAAAGCCCGCAAGGGCCGGAAGTGGTGGTCGACGGCCAGCGGTTGCTGGCCTTCTGCAGCAATGACTACCTGGGCCTGGCCAACCACCCCGAGGTGATCGCCGCCTGGCAGGCCGGCGCCGAGCGCTGGGGTGTGGGGGGGGGGGCCTCGCACCTGGTAGTTGGCCACAGCACCCCGCACCATCAGGTTGAAGAGGCACTGGCTGAACTCACCGGCCGCCCGCGCGCCCTGTTGTTCTCCACCGGCTATATg                                                                                                                                                                                          <  1:6783/298‑1 (MQ=255)
                                                                                                                      ccTGGCCGAACGGCGCGCCGCAGACCTGTATCGGCAGCGGCCATTGCTGGAAAGCCCGCAAGGGCCGGAAGTGGTGGTCGACGGCCAGCGGTTGCTGGCCTTCTGCAGCAATGACTACCTGGGCCTGGCCAACCACCCCGAGGTGATCGCCGCCTGGCAGGCCGGCGCCGAGCGCTGGGGTGTGGGGGGGGGGCCTCGCACCTGGTAGTTGGCCACAGCACCCCGCACCATCCGGTTGAAGAGGCACTGGCTGAACTCACCGGCCGCCCGCGCGCCCTGTTGTTCTCCACCGGCTatattg                                                                                                                                                                                         >  1:338543/1‑299 (MQ=255)
                                                                                                                                 ggCGCGCCGCAGACCTGTATCGGCAGCGGCCATTGCTGGAAAGCCCGCAAGGGCCGGAAGTGGTGGTCGACGGCCAGCGGTTGCTGGCCTTCTGCAGCAATGACTACCTGGGCCTGGCCAACCACCCCGAGGTGATCGCCGCCTGGCAGGCCGGCGCCGAGCGCTGGGGTGTGGGGGGGGGGGCCTCGCACCTGGTAGTTGGCCACAGCACCCCGCCCCCTCAGGTTGAAGAGGCCCTGGCTGAACTCCCCGCCCGCCCGCGCCCCCTTTTGTTCCCCACCGGCtgttgggccaatgcggg                                                                                                                                                                               >  2:389085/1‑285 (MQ=255)
                                                                                                                                 ggCGCGCCGCAGACCTGTATCGGCAGCGGCCATTGCTGGAAAGCCCGCAAGGGCCGGAAGTGGTGGTCGACGGCCAGCGGTTGCTGGCCTTCTGCAGCAATGACTACCTGGGCCTGGCCAACCACCCCGAGGTGATCGCCGCCTGGCAGGCCGGCGCCGAGCGCTGGGGTGTGGGGGGGGGGGCCTCGCACCGGGTAGTTGGCCACAGCACCCCGCACCATCAGGTTGAAGAGGCACTGGCTGAACTCACCGGCCGCCCGCGCGCCCTGTTGTTCTCCACCGGCTATATGGCCAAACTggg                                                                                                                                                                               >  1:638184/1‑301 (MQ=255)
                                                                                                                                         cccagccctgtgtgggccgCGGCCATTGCTGGGAAGCCCGCAAGGGCCGGAAGTGGTGGTCGACGGCCAGCGGTTGCTGGCCTTCTGCAGCAATGACTACCTGGGCCTGGCCAACCACCCCGAGGTGATCGCCGCCTGGCAGGCCGGCGCCGAGCGCTGGGGTGTGGGGGGGGGGCCTCGCACCTGGTAGTTGGCCACAGCACCCCGCACCATCAGGTTGAAGAGGCACTGGCTGAACTCACCGGCCGCCCGCGCGCCCTGTTGTTCTCCACCGGCTATATGGCCAATCTGGGCGCCATCa                                                                                                                                                                       <  2:167100/283‑1 (MQ=255)
                                                                                                                                          cAGACCTGTATCGGCAGCGGCCATTGCTGGAAAGCCCGCAAGGGCCGGAAGTGGTGGTCGACGGCCAGCGGTTGCTGGCCTTCTGCAGCAATGACTACCTGGGCCTGGCCAACCACCCCGAGGTGATCGCCGCCTGGCAGGCCGGCGCCGAGCGCTGGGGTGTGGGGGGGGGGCCTCGCACCTGGTAGTTGGCCACAGCACCCCGCACCATCCGGTTGAAGAGGCACTGGCTGCACTCCCCGGCCGCCCGCGCGCCCTGTTTTTCCCCCCCGCCGTTATGGCCAAttcggcccacaccccc                                                                                                                                                                     >  2:571583/1‑286 (MQ=255)
                                                                                                                                            aggcccgaagcgccagcagacatggatGGACAGCCCGCAATGGCCGGAACTGGTGCTGGACGGCAAGCGGATGCTGTCCTTATCCCGCAATTCCTACCTGCGCCTAGCCAACCACCCCGAGGTGATCGCCGCCTGGCAGGCCGGCGCCGAGCGCTGGGGTGTGGGGGGGGGGCCTCGCACCTGGTAGTTGGCCACAGCACCCCGCACCATCAGGTTGAAGAGGCACTGGCTGAACTCACCGGCCGCCCGCGCGCCCTGTTGTTCTCCACCGGCTATATGGCCAATCTGGGCGCCATCACcg                                                                                                                                                                    <  2:480913/275‑1 (MQ=255)
                                                                                                                                                 ggTATCGGCAGCGGCCATTGCGGGAAAGCCCGCAAGGGCCGGAAGTGGTGGTCGACGGCCAGCGGTTGCTGGCCTTCTGCAGCAATGACTACCTGGGCCTGGCCAACCACCCCGAGGTGATCGCCGCCTGGCAGGCCGGCGCCGAGCGCTGGGGTGTGGGGGGGGGGCCTCGCACCTGGTAGTTGGCCACAGCACCCCGCACCATCAGGTTGAAGAGGCACTGGCTGAACTCACCGGCCGCCCGCGCGCCCTGTTGTTCTCCACCGGCTATATGGCCAATCTGGGCGCCATCACCGCGCtg                                                                                                                                                               <  1:174394/300‑1 (MQ=255)
                                                                                                                                                       ggCAGCGGTCATTGCGCGAAAGCCCGCAAGGGGCGGAAGTGGTGGGAGACGGCCAGCGGTGGCTGGCCTTTTGCAGCAATGACCACCTGGGCCTGGCCAACCCCCCCGAGGTGATCGCCGCCTGGCAGGCCGGCGCCGAGCGCTGGGGTGTGGGGGGGGGGGCCTCGCACCTGGTAGTTGGCCACAGCACCCCGCACCATCAGGTTGAAGAGGCACTGGCTGAACTCACCGGCCGCCCGCGCGCCCTGTTGTTCTCCACCGGCTATATGGCCAATCTGGGCGCCATCACCGCGCTGGTggg                                                                                                                                                          <  1:126221/301‑1 (MQ=255)
                                                                                                                                                         cAGCGGCCATTGCTGGAAAGCCCGCAAGGGCCGGAAGTGGTGGTCGACGGCCAGCGGTTGCTGGCCTTCTGCAGCAATGACTACCTGGGCCTGGCCAACCACCCCGAGGTGATCGCCGCCTGGCAGGCCGGCGCCGAGCGCTGGGGTGTGGGGGGGGGGGCCTCGCACCTGGTAGTTGGCCACAGCACCCCGCACCATCAGGTTGAAGAGGCACTGGCTGAACTCACCGGCCGCCCGCGCGCCCTGTTGTTCTCCACCGGCTATATGGCCAATCTGGGCGCCATCACCGCGCTGGTGGGGc                                                                                                                                                        >  1:701841/1‑301 (MQ=255)
                                                                                                                                                                 aTTGCTGGAAAGCCCGCAAGGGCCGGAAGTGGTGGTCGACGGCCAGCGGTTGCTGGCCTTCTGCAGCAATGACTACCTGGGCCTGGCCAACCACCCCGAGGTGATCGCCGCCTGGCAGGCCGGCGCCGAGCGCTGGGGTGTGGGGGGGGGGGCCTCGCACCTGGTAGTTGGCCACAGCACCCCGCACCATCAGGTTGAAGAGGCACTGGCTGAACTCACCGGCCGCCCGCGCCCCCTGTTGTTCGCCACCGGCTATATGGCCAACCTGGGCGCCATCACCGCCCCGGTGGGGAGGGccgac                                                                                                                                               >  2:496320/1‑297 (MQ=255)
                                                                                                                                                                    gcgtgAAAGCCCGCAAGGGCCGGAAGGGGTGGTCGACGGCCATCGGTTGCTGGCCTTCTGCAGCAATGACTACCTGGGCCTGGCCAACCACCCCGAGGTGATCGCCGCCTGGCAGGCCGGCGCCGAGCGCTGGGGTGTGTGGGGGGGGGGCCTCGCACCTGGTAGTTGGCCACAGCACCCCGCACCATCAGGTTGAAGAGGCACTGGCTGAACTCACCGGCCGCCCGCGCGCCCTGTTGTTCTCCACCGGCTATATGGCCAATCTGGGCGCCATCACCGCGCTGGTGGGGCAGGGCGacac                                                                                                                                              <  1:730283/297‑1 (MQ=255)
                                                                                                                                                                          aaGCCCGCAAGGGCCGGAAGTGGTGGTCGACGGCCAGCGGTTGCTGGCCTTCTGCAGCAATGACTACCTGGGCCTGGCCAACCACCCCGAGGTGATCGCCGCCTGGCAGGCCGGCGCCGAGCGCTGGGGTGTGGGGGGGGGGGCCTCGCACCTGGTAGTTGGCCACAGCACCCCGCACCATCAGGTTGAAGAGGCACTGGCTGAACTCACCGGCCGCCCGCGCGCCCTGTTGTTCTCCACCGGCTATATGGCCAATCTGGGCGCCATCACCGCGCTGGTGGGGCAGGGCGACACGGTATTg                                                                                                                                       >  2:452432/1‑301 (MQ=255)
                                                                                                                                                                                gCAAGGGCCGGAAGTGGTGGTCGACGGCCAGCGGTTGCTGGCCTTCTGCAGCAATGACTACCTGGGCCTGGCCAACCACCCCGAGGTGATCGCCGCCTGGCAGGCCGGCGCCGAGCGCTGGGGTGTGGGGGGGGGGCCTCGCACCTGGTAGTTGGCCACAGCACCCCGCACCATCAGGTTGAAGAGGCACTGGCTGAACTCACCGGCCGCCCGCGCGCCCTGTTGTTCTCCACCGGCTATATGGCCAATCTGGGCGCCATCACCGCGCTGGTGGGGCAGGGCGACACGGTATTGCAAGAcc                                                                                                                                >  1:417512/1‑301 (MQ=255)
                                                                                                                                                                                  aaGGGCCGGAAGTGGTGGTCGACGGCCAGCGGTTGCTGGCCTTCTGCAGCAATGACTACCTGGGCCTGGCCAACCACCCCGAGGTGATCGCCGCCTGGCAGGCCGGCGCCGAGCGCTGGGGTGTGGGGGGGGGGGCCTCGCACCTGGTAGTTGGCCACAGCACCCCGCACCATCAGGTTGAAGAGGCACTGGCTGAACTCACCGGCCGCCCGCGCGCCCTGTTGTTCTCCACCGGCTATATGGCCAATCTGGGCGCCATCACCGCGCTGGTGGGGCAGGGCGACACGGTATTGCAAGACCg                                                                                                                               <  1:751430/301‑1 (MQ=255)
                                                                                                                                                                                   ggggccggaggggtgggtcgcgggCCGGCGGTGGCTGCCCTTCTTCAGAAATGACTACCTGGGCCGGGCAAACCACCCCGGGGTGATGGCCGCCTGGCGGGCGGGGGCCGGGCGCTGGGGTGGGGGGGGGGGGCCTCGCACCTGGTAGTTGGCCACAGCACCCCGCACCATCAGGTTGAAGAGGCACTGGCTGAACTCACCGGCCGCCCGCGCGCCCTGTTGTTCTCCACCGGCTATATGGCCAATCTGGGCGCCATCACCGCGCTGGTGGGGCAGGGCGACACGGTATTGCAAGACCGcc                                                                                                                             <  1:309711/279‑1 (MQ=255)
                                                                                                                                                                                         ggacgtggtgCTGGACGGCCAGCGGCTGCTGTCCTTCTCCAGCACTGACTACCTGAGCCTGTCCAACCACCCCGAGGTGATCGCCGCCTGGCAGGGCGGCGCCGAGCGCTGGGGTGTGGGGGGGGGGGCCTCGCACCTGGTAGTTGGCCACAGCACCCCGCACCATCAGGTTGAAGAGGCACTGGCTGAACTCACCGGCCGCCCGCGCGCCCTGTTGTTCTCCACCGGCTATATGGCCAATCTGGGCGCCATCACCGCGCTGGTGGGGCAGGGCGACACGGTATTGCAAGACCGCCTCAAc                                                                                                                        <  1:161941/297‑1 (MQ=255)
                                                                                                                                                                                                        gtcccggCGGTTGCTGGCCTCTGCACACACTGACGTCCAGGGACTGGCCAACCCCCCCGAGGTGATCGCCGCCTGGCAGGCCGGCGCCGGGCGCTGGGGTGTGGGGGGGGGGGCCTCGCACCTGGTAGTTGGCCACAGCACCCCGCACCATCAGGTTGAAGAGGCACTGGCTGAACTCACCGGCCGCCCGCGCGCCCTGTTGTTCTCCACCGGCTATATGGCCAATCTGGGCGCCATCACCGCGCTGGTGGGGCAGGGCGACACGGTATTGCAAGACCGCCTCAACCATGCCTCGTTGCTg                                                                                                         <  2:65426/295‑1 (MQ=255)
                                                                                                                                                                                                        gggCCAGCGGTTGCTGGCCTGTTGCAGCAATGACGACCTGGGCCTGGCCAACCACCCCGAGGTGATCGCCGCCTGGCAGGCCGGCGCCGAGCGCTGGGGTGTGGGGGGGGGCCTCGCACCTGGTAGTTGGCCACAGCACCCCGCACCATCAGGTTGAAGAGGCACTGGCTGAACTCACCGGCCGCCCGCGCGCCCTGTTGTTCTCCACCGGCTATATGGCCAATCTGGGCGCCATCACCGCGCTGGTGGGGCAGGGCGACACGGTATTGCAAGACCGCCTCAACCATGCCTCGTTGCTGGa                                                                                                       <  2:133879/300‑1 (MQ=255)
                                                                                                                                                                                                          gCCAGCGGTTGCGGGCCTTCTGCAGCAATGACTACCTGGGCCTGGCCAACCACCCCGAGGTGATCGCCGCCTGGCAGGCCGGGGCCGAGCGCTGGGGTGTGGGGGGGGGGCCTCGCACCTGGTAGTTGGCCACAGCACCCCGCACCATCAGGTTGAAGAGGCACTGGCTGAACTCACCGGCCGCCCGCGCGCCCTGTTGTTCTCCACCGGCTATATGGCCAATCTGGGCGCCATCACCGCGCTGGTGGGGCAGGGCGACACGGTATTGCAAGACCGCCTCAACCATGCCTCGTTGCTGGAt                                                                                                      <  2:266030/301‑1 (MQ=255)
                                                                                                                                                                                                               cGGTTGCTGGCCTTCTGCAGCAATGACTACCTGGGCCTGGCCAACCACCCCGAGGTGATCGCCGCCTGGCAGGCCGGCGCCGAGCGCTGGGGTGTGGGGGGGGGGGCCTCGCACCTGGTAGTTGGCCACAGCACCCCGCACCATCAGGTTGAAGAGGCACTGGCTGAACTCACCGGCCGCCCGCGCGCCCTGTTGTTCTCCACCGGCTATATGGCCAATCTGGGCGCCATCACCGCGCTGGTGGGGCAGGGCGCCACGGTATTGCAAGACAGCATCACCCATGCCtgtgtgcttggtgggg                                                                                                  >  1:403740/1‑286 (MQ=255)
                                                                                                                                                                                                                         ccttctcgagcaCTGACTCCCCGGGCCCGGCCAACCACACCGAGGTGATCGCCGCCTGGCCGGCCGGCGCCGAGCGCTGGGGTGGTGGGGGGCGGGCCTCGCACCTGGTAGTTGGCCACAGCACCCCGCACCATCAGGTTGAAGAGGCACTGGCTGAACTCCCCGGCCGCCCGCGCGCCCTGTTGTTCTCCACCGGCTATATGGCCAATCTGGGCGCCATCACCGCGCTGGTGGGGCAGGGCGACACGGTATTGCAAGACCGCCTCAACCATGCCTCGTTGCTGGATGGCGGGTTGCTCAg                                                                                        <  2:381351/293‑1 (MQ=255)
                                                                                                                                                                                                                           ttCTGCAGCAATGACTACCTGGGCCTGGCCAACCACCCCGAGGTGATCGCCGCCTGGCAGGCCGGCGCCGAGCGCTGGGGTGTGGGGGGGGGGGCCTCGCACCTGGTAGTTGGCCACAGCACCCCGCACCATCAGGTTGAAGAGGCACTGGCTGAACTCACCGGCCGCCCGCGCGCCCTGTTGTTCTCCACCGGCTATATGGCCAATCTGGGCGCCATCACCGCGCTGGTGGGGCAGGGCGACACGGTATTGCAAGACCGCCTCAACCATGCGTCGTTGCTGGCTGGCGGGTTGCTCAgtg                                                                                      >  2:518474/1‑301 (MQ=255)
                                                                                                                                                                                                                                 gaggaatgcgtACCTGGGCCGGGACAAACACCCCGGGGGGATCGACGCCTGGCAGGCCGGCGCCGAGCGCTGGGGTGTGGGGGGGGGGGCCTCGCACCTGGTAGTTGGCCACAGCACCCCGCACCATCAGGTTGAAGAGGCACTGGCTGAACTCACCGGCCGCCCGCGCGCCCTGTTGTTCTCCACCGGCTATATGGCCAATCTGGGCGCCATCACCGCGCTGGTGGGGCAGGGCGACACGGTATTGCAAGACCGCCTCAACCATGCCTCGTTGCTGGATGGCGGGTTGCTCAGTGGTGcc                                                                                <  2:701542/291‑1 (MQ=255)
                                                                                                                                                                                                                                   ccaATGACTCCCTGGGCCTGGCCAACCACCCCGCGGGGAGCGCCGCCTGGCAGGCCGGCGCCGAGCGCTGGGGTGTGGGGGGGGGGGCCTCGCACCTGGTAGTTGGCCACAGCACCCCGCACCATCAGGTTGAAGAGGCACTGGCTGAACTCACCGGCCGCCCGCGCGCCCTGTTGTTCTCCACCGGCTATATGGCCAATCTGGGCGCCATCACCGCGCTGGTGGGGCAGGGCGACACGGTATTGCAAGACCGCCTCAACCATGCCTCGTTGCTGGATGGCGGGTTGCTCAGTGGTGCCCg                                                                              <  2:637912/300‑1 (MQ=255)
                                                                                                                                                                                                                                         ccTACCTGGGCCTGGCCAACCACCCCGAGGTGATCGCCGCCTGGCAGGCCGGCGCCGAGCGCTGGGGTGTGGGGGGGGGGCCTCGCACCTGGTAGTTGGCCACAGCACCCCGCACCATCAGGTTGAAGAGGCACTGGCTGAACTCACCGGCCGCCCGCGCGCCCTGTTGTTCTCCACCGGCTATATGGCCAATCTGGGCGCCATCACCGCGCTGGTGGGGCAGGGCGACACGGTATTGCAAGACCGCCTCAACCATGCCTCGTTGCTGGATGGCGGGTTGCTCAGTGGTGCCCGCTTCAAc                                                                       <  1:589904/300‑1 (MQ=255)
                                                                                                                                                                                                                                         aCTACCTGGGCCTGGCCAACCACCCCGAGGTGATCGCCGCCTGGCAGGCCGGCGCCGAGCGCTGGGGTGTGGGGGGGGGGGCCTCGCACCTGGTAGTTGGCCACAGCACCCCGCACCATCAGGTTGAAGAGGCACTGGCTGAACTCACCGGCCGCCCGCGCGCCCTGTTGTTCTCCACCGGCTATATGGCCAATCGGGGCGCCATCACCGCGCTGGTGGGGCAGGGCGACACGGTATTGCAAGACCGCCTCAACCATGCCTCGTTGCTGGATGGCTGGTTGCCCAGTGTGTCCCGCTTCaa                                                                        >  2:189661/1‑301 (MQ=255)
                                                                                                                                                                                                                                            cgcacggcctaccggtctcccccGAGGGGCTCGCCGCCGGGCAGGCGGGCGCCGGGGGCGGGGGGGGGGGGGGGGGGGCCTCGCACCTGGTAGTTGGCCACAGCACCCCGCACCATCAGGTTGAAGAGGCACTGGCTGAACTCACCGGCCGCCCGCGCGCCCTGTTGTTCTCCACCGGCTATATGGCCAATCTGGGCGCCATCACCGCGCTGGTGGGGCAGGGCGACACGGTATTGCAAGACCGCCTCAACCATGCCTCGTTGCTGGATGGCGGGTTGCTCAGTGGTGCCCGCTTCAACCg                                                                     <  2:640822/282‑1 (MQ=255)
                                                                                                                                                                                                                                                gggCCTGGCCAACCACCCCGGGGGGAGCGGCGGCTGGGCGGCCGGCGCCGGGCGCGGGGGTGGGGGGGGGGGGGCCTCGCACCTGGTAGTTGGCCACAGCACCCCGCACCATCAGGTTGAAGAGGCACTGGCTGAACTCACCGGCCGCCCGCGCGCCCTGTTGTTCTCCACCGGCTATATGGCCAATCTGGGCGCCATCACCGCGCTGGTGGGGCAGGGCGACACGGTATTGCAAGACCGCCTCAACCATGCCTCGTTGCTGGATGGCGGGTTGCTCAGTGGTGCCCGCTTCAACCGCTAc                                                                 <  1:496525/301‑1 (MQ=255)
                                                                                                                                                                                                                                                 ggCCTGGCCAACCACCCCGAGGTGATCGCCGCCTGGCAGGCCGGCGCCGAGCGCTGGGGTGTGGGGGGGGGGGCCTCGCACCTGGTAGTTGGCCACAGCACCCCGCACCATCAGGTTGAAGAGGCACTGGCTGAACTCACCGGCCGCCCGCGCGCCCTGTTGTTCTCCACCGGCTATATGGCCAATCTGGGCGCCATCACCGCGCTGGTGGGGCAGGGCGACACGGTATTGCAAGACCGCCTCAACCATGCCTCGTTGCTGGATGGCGGGTTGCTCAGTGGTGCCCGCTTCAACCGCTAcc                                                                <  1:452622/301‑1 (MQ=255)
                                                                                                                                                                                                                                                         cacaccccccGCGGGGATCGCCGCCTGGCAGGGCGGCGACGGGCGCTGGGGGGGGGGGGGGGGGGCCTCGCACCTGGTAGTGGGCCACAGCACCCCGCACCATCAGGTTGAAGAGGCACTGGCTGAACTCACCGGCCGCCCGCGCGCCCTGTTGTTCTCCACCGGCTATATGGCCAATCTGGGCGCCATCACCGCGCGGGTGGGGCAGGGCGACACGGTATTGCAAGACCGCCTCAACCATGCCTCGTTGCTGGATGGCGGGTTGCTCAGTGGTGCCCGCTTCAACCGCTACCTGCACAAc                                                        <  1:18089/295‑1 (MQ=255)
                                                                                                                                                                                                                                                          aaCCACCCCGGGGGGCTTGCCGCCTGGCAGGCCGGGGCCGGGGGCTGGGGTGGGGGGGGGGGGGCCTCGCACCTGGTAGTTGGCCACAGCACCCCGCACCATCAGGTTGAAGAGGCACTGGCTGAACTCACCGGCCGCCCGCGCGCCCTGTTGTTCTCCACCGGCTATATGGCCAATCTGGGCGCCATCACCGCGCTGGTGGGGCAGGGCGACACGGTATTGCAAGACCGCCTCAACCATGCCTCGTTGCTGGATGGCGGGTTGCTCAGTGGTGCCCGCTTCAACCGCTACCTGCACAACg                                                       <  1:389243/301‑1 (MQ=255)
                                                                                                                                                                                                                                                             caccCCGAGGTGATCGCCGCCTGGCAGGCCGGCGCCGAGCGCTGGGGTGTGGGGGGGGGGGCCTCGCACCTGGTAGTTGGCCACAGCACCCCGCACCATCAGGTTGAAGAGGCACTGGCTGAACTCACCGGCCGCCCGCGCGCCCTGTTGTTCTCCACCGGCTATATGGCCAATCTGGGCGCCATCACCGCGCTGGTGGGGCAGGGCGACACGGTATTGCAAGACCGCCTCAACCATGCCTCGTTGCTGGATGGCGGGTTGCTCAGTGGTGCCCGCTTCAACCGCTACCGGCACAACGAcc                                                    >  1:525056/1‑301 (MQ=255)
                                                                                                                                                                                                                                                              ccccgggggggatggcggccgggccggccggggccggGCGCTGGGGTGTGGGGGGGGGGGCCTCGCACCTGGTAGTTGGCCACAGCACCCCGCACCATCAGGTTGAAGAGGCACTGGCTGAACTCACCGGCCGCGCGCGCGCCCTGTTGTTCTCCACCGGCTATATGGCCAATCTGGGCGCCATCACCGCGCTGGTGGGGCAGGGCGACACGGTATTGCAAGACCGCCTCAACCATGCCTCGTTGCTGGATGGCGGGTTGCTCAGTGGTGCCCGCTTCAACCGCTACCTGCACAACGAccc                                                   <  1:582301/276‑1 (MQ=255)
                                                                                                                                                                                                                                                              accCCGAGGTGATCGCCGCCTGGCAGGCCGGCGCCGAGCGCTGGGGTGTGGGGGGGGGGGCCTCGCACCTGGTAGTTGGCCACAGCACCCCGCACCATCAGGTTGAAGAGGCACTGGCTGAACTCACCGGCCGCCCGCGCGCCCTGTTGTTCTCCACCGGCTATATGGCCAATCTGGGCGCCATCACCGCGCTGGTGGGGCAGGGCGACACGGTATTGCAAGACCGCCTCAACCATGCCTCGTTGCTGGATGGCGGGTTGCTCAGTGGTGCCCGCTTCAACCgagatcggaagagcacacg                                                   >  1:641094/1‑283 (MQ=255)
                                                                                                                                                                                                                                                                            gcggccTGGCAGGCGGGAGCGGAGGGCTGGGGGGGGGGGGGGGGGGCCGCGCACCTGGTAGTTGGCCACAGCACCCCGCACCATCAGGTTGAAGGGGCACTGGCTGAACTCACCGGCCGCCCGCGCGCCCTGTTGTTCTCCACCGGCTATATGGCCAATCTGGGCGCCATCACCGCGCTGGTGGGGCAGGGCGACACGGTATTGCAAGACCGCCTCAACCATGCCTCGTTGCTGGATGGCGGGTTGCTCAGTGGTGCCCGCTTCAACCGCTACCTGCACAACGACCCGGCCAGCCTGGCCa                                     <  2:483809/298‑1 (MQ=255)
                                                                                                                                                                                                                                                                                ccTGGCAGGCCGGCGCCGAGCGCTGGGGTGTGGGGGGGGGGGCCTCGCACCTGGTAGTTGGCCACAGCACCCCGCACCATCAGGTTGAAGAGGCACTGGCTGAACTCACCGGCCGCCCGCGCGCCCTGTTGTTCTCCACCGGCTATATGGCCAATCTGGGCGCCATCCCCGCGCTGGTGGGGCAGGGCGACACGGTATTGCAAGACCGCCTCAACCATGCCTCGTTGCTGGATGGCGGGTTGCTCAGTGGTGCCCGCTTCAACCGCTACCTGCCCAACGACCCGGCCCGccgtgccacccc                                 >  2:727793/1‑291 (MQ=255)
                                                                                                                                                                                                                                                                                                      cTGGGGTGTGGGGGGGGGGGCCTCGCACCTGGTAGTTGGCCACAGCACCCCGCACCATCAGGTTGAAGAGGCACTGGCTGAACTCACCGGCCGCCCGCGCGCCCTGTTGTTCTCCACCGGCTATATGGCCAATCTGGGCGCCATCACCGCGCTGGTGGGGCAGGGCGACACGTTATTGCAAGACCGCCTCAACCATCCCTCGTTGCTGGATGGCGGGTTGCTCAGTGGTGCCTGCTTCAACCGCTACCTGCACAACGACCCGGCCAGCCTGGCCAGCCGCCTGGCAAAGCCGTTCGGCAAt           >  2:346083/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                                                gggggggggggCCTGGCACCTGGTAGTTGGCCACAGCACCCCGCACCATAAGGTTGAAGAGGCACTGGCTGAACTCACCGGCCGCCCGCGCGCCCTGTTGTTCTCCACCGGCTATATGGCCAATCTGGGCGCCATCACCGCGCTGGTGGGGCAGGGCGACACGGTATTGCAAGACCGCCTCAACCATGCCTCGTTGCTGGATGGCGGGTTGCTCAGTGGTGCCCGCTTCAACCGCTACCTGCACAACGACCCGGCCAGCCTGGCCAGCCGCCTGGACAAGGCGGTCGGCAATACCCtggtg  <  2:690317/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                                                              
AGCACGCCGACGAAGTGCACCAGGCAGCCATCGAGCAAGCGCTGGTCGAGCAGCGCAGCAGCGAGATGTTCTACGACGCCGCGACAGCCTGACATGGCTTTCGACCTGGCAGCGCGCCTGGCCGAACGGCGCGCCGCAGACCTGTATCGGCAGCGGCCATTGCTGGAAAGCCCGCAAGGGCCGGAAGTGGTGGTCGACGGCCAGCGGTTGCTGGCCTTCTGCAGCAATGACTACCTGGGCCTGGCCAACCACCCCGAGGTGATCGCCGCCTGGCAGGCCGGCGCCGAGCGCTGGGGTGTGGGGGGGGGGGCCTCGCACCTGGTAGTTGGCCACAGCACCCCGCACCATCAGGTTGAAGAGGCACTGGCTGAACTCACCGGCCGCCCGCGCGCCCTGTTGTTCTCCACCGGCTATATGGCCAATCTGGGCGCCATCACCGCGCTGGTGGGGCAGGGCGACACGGTATTGCAAGACCGCCTCAACCATGCCTCGTTGCTGGATGGCGGGTTGCTCAGTGGTGCCCGCTTCAACCGCTACCTGCACAACGACCCGGCCAGCCTGGCCAGCCGCCTGGACAAGGCGGTCGGCAATACCCTGGTG  >  AE015451/441027‑441626

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: