Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154514,169,4720TC30.8% 45.3 / 20.2 40K141E (AAG→GAG) katGCatalase‑peroxidase
Reads supporting (aligned to +/- strand):  ref base T (14/13);  new base C (12/0);  total (26/14)
Fisher's exact test for biased strand distribution p-value = 2.97e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TGCACGGCGGCCAGCGGGTTTTCCAGGTTACGCTCGGCCGACAGGTCGCGATTCTGCTCTTCACCATGCTTGGCAGGCTCGGCCACCAGGTCACCCTGCCCCGGCGGCTGCGCCTTGACCTGATCCTTGCCGTAACGGGTATCACCCCCTAGCCAGACCTTTTCCGACCCCCAGTACACGTCTTCGTCCGGCTCCCACACATCGG‑CACGCCCGCCGGAGAAGCCGAAGGTCTTGAAGCCCATGGATTCGAGGGCAACATTACCCGTGAGTACGATCAGGTCGGCCCAGGAAATCTTG‑T‑TGCCGTACTTCTGCTTGATCGGCCACAGCAGGCGCCGGGCCTTGTCCAGGCTGACGTTGTCCGGCCAGCTGTTGAGCGGGGCGAAGCGTTGCTGGCCAGAGCCCGCACCACCGCGGCCATCACCAATGCGGTAGGTGCCGGCGCTATGCCAGGCCATGCGGATGAACAGCGGGCCGTAGTGGCCGAAGTCGGCCGGCCACCAATCCTGCGAGTCGGTCATCAGGGCGGTCAGGTCTTTCTTCAGGGCCTGAAAGTCGAGGCTCTTGAACGCCTTGGCATAGTCGAAGTCCGG  >  AE015451/4169177‑4169766
                                                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                                        
tGCACGGCGGCCAGCGGGTTTTCCAGGTTACGCTCGGCCGACAGGTCGCGATTCTGCTCTTCACCATGCTTGGCAGGCTCGGCCACCAGGTCACCCTGCCCCGGCGGCTGCGCCTTGACCTGATCCTTGCCGTAACGGGTATCACCCCCTAGCCAGACCTTTTCCGACCCCCAGTACACGTCTTCGTCCGGCTCCCACACATCGG‑CACGCCCGCCGGAGAAGCCGAAGGTCTTGAAGCCCATGGATTCGAGGGCAACATTACCCGTGAGTACGATCAGGTCGGCCCAGGAAATCttg‑t‑tcc                                                                                                                                                                                                                                                                                                   >  1:566694/1‑299 (MQ=255)
    cggcggCCACCGGGTTTTCCAGGTTACGCTCGTCCGACAGGTCTCGATTCTGCTCTTCACCATGCTTGGCAGGCTCGGCCACCAGGTCACCCTGCCCCGGCGGCTGCGCCTTGACCTGATCCTTGCCGTAACGGGTATCACCCCCTAGCCAGACCTTTTCCGACCCCCAGTACACGTCTTCGTCCGGCTCCCACACATCGG‑CACGCCCGCCGGAGAAGCCGAAGGTCTTGAAGCCCATGGATTCGAGGGCAACATTACCCGTGAGTACGATCAGGTCGGCCCAGGAAATCTTG‑T‑TGCCGTa                                                                                                                                                                                                                                                                                               <  2:459124/301‑1 (MQ=255)
        ggcCAGCGGGTTTTCCAGGTTACGCTCGGCCGACAGGTCGCGATTCTGCTCTTCACCATGCTTGGCAGGCTCGGCCACCAGGTCACCCTGCCCCGGCGGCTGCGCCTTGACCTGATCCTTGCCGTAACGGGTATCACCCCCTAGCCAGACCTTTTCCGACCCCCAGTACACGTCTTCGTCCGGCTCCCACACATCGG‑CACGCCCGCCGGAGAAGCCGAAGGTCTTGAAGCCCATGGATTCGAGGGCAACATTACCCGTGAGTACGATCAGGTCGGCCCAGGAAATCTCG‑T‑TGGCCTACttt                                                                                                                                                                                                                                                                                           >  1:445897/1‑300 (MQ=255)
          ccatctggTTTTCCAGGTTACGCTCGGCCGACAGGTCGCGATTCTGCGCTTCACCATGCTTGGCAGGCTCGGCCACCAGGTCACCCTGCCCCGGCGGCTGCGCCTTGACCTGATCCTTGCCGTAACGGGTATCACCCCCTAGCCAGACCTTTTCCGACCCCCAGTACACGTCTTCGTCCGGCTCCCACACATCGG‑CACGCCCGCCGGAGAAGCCGAAGGTCTTGAAGCCCATGGATTCGAGGGCAACATTACCCGTGAGTACGATCAGGTCGGCCCAGGAAATCTTG‑T‑TGCCGTACTTCTg                                                                                                                                                                                                                                                                                         <  1:158676/295‑1 (MQ=255)
                           ttACGCTCGGCCGACAGGTCGCGATTCTGCTCTTCACCATGCTTGGCAGGCTCGGCCACCAGGTCACCCTGCCCCGGCGGCTGCGCCTTGACCTGATCCTTGCCGTAACGGGTATCACCCCCTAGCCAGACCTTTTCCGACCCCCAGTACACGTCTTCGTCCGGCTCCCACACATCGG‑CACGCCCGCCGGAGAAGCCGAAGGTCTTGAAGCCCATGGATTCGAGGGCAACATTACCCGTGAGTACGATCAGGTCGGCCCAGGAAATCTTG‑T‑TTGCGTACTTCTGCTTGATCGGCCAcagca                                                                                                                                                                                                                                                                        >  1:91211/1‑301 (MQ=255)
                                cTCGGCCGACAGGTCGCGATTCTGCTCTTCACCATGCTTGGCAGGCTCGGCCACCAGGTCACCCTGCCCCGGCGGCTGCGCCTTGACCTGATCCTTGCCGTAACGGGTATCACCCCCTAGCCAGACCTTTTCCGACCCCCAGTACACGTCTTCGTCCGGCTCCCACACATCGG‑CACGCCCGCCGGAGAAGCCGAAGGTCTTGAAGCCCATGGATTCGAGGGCAACATTACCCGTGAGTACGATCAGGCCGGCCCAGGAAATCTTG‑T‑TGCCGTACTTCTGCTTGATAGGCCACAGCCGgcgc                                                                                                                                                                                                                                                                   >  1:296144/1‑301 (MQ=255)
                                  cGGCCGACAGGTCGCGATTCTGCTCTTCACCATGCTTGGCAGGCTCGGCCACCAGGTCACCCTGCCCCGGCGGCTGCGCCTTGACCTGATCCTTGCCGTAACGGGTATCACCCCCTAGCCAGACCTTTTCCGACCCCCAGTACACGTCTTCGTCCGGCTCCCACACATCGG‑CACGCCCGCCGGAGAAGCCGAAGGTCTTGAAGCCCATGGATTCGAGGGCAACATTACCCGTGAGTACGATCAGGTCGGCCCAGGAAATCTCG‑TATGCCGT‑CTTCTGCTTGaaaaaaaacagaaggcgcaa                                                                                                                                                                                                                                                                 >  1:350041/1‑282 (MQ=255)
                                  cGGCCGACAGGTCGCGATTCTGCTCTTCACCATGCTTGGCAGGCTCGGCCACCAGGTCACCCTGCCCCGGCGGCTGCGCCTTGACCTGATCCTTGCCGTAACGGGTATCACCCCCTAGCCAGACCTTTTCCGACCCCCAGTACACGTCTTCGTCCGGCTCCCACACATCGG‑CACGCCCGCCGGAGAAGCCGAAGGTCTTGAAGCCAATGGATTCGAGGGCAACATTACCCGTGAGTACCACCAGGTCGGCACAGGAAATCTTG‑T‑TGCCGTACTTCTGCTTGACCGGCCACAGAAGGCGCCg                                                                                                                                                                                                                                                                 >  2:180593/1‑301 (MQ=255)
                                    gCCGACAGGTCGCGATTCTGCTCTTCACCATGCTTGGCAGGCTCGGCCACCAGGTCACCCTGCCCCGGCGGCTGCGCCTTGACCTGATCCTTGCCGTAACGGGTATCACCCCCTAGCCAGACCTTTTCCGACCCCCAGTACACGTCTTCGTCCGGCTCCCACACATCGG‑CACGCCCGCCGGAGAAGCCGAAGGTCTTGAAGCCCATGGATTCGAGGGCAACATTACCCGTGAGTACGATCAGGTCGGCCCAGGAAATCTTG‑T‑TGCCGTACTTCTGCTTGATCGGCCACAGCAGGCGCCggg                                                                                                                                                                                                                                                               >  1:274629/1‑301 (MQ=255)
                                        aCAGGTCGCGATTCTGCTCTTCACCATGCTTGGCAGGCTCGGCCACCAGGTCACCCTGCCCCGGCGGCTGCGCCTTGACCTGATCCTTGCCGTAACGGGTATCACCCCCTAGCCAGACCTTTTCCGACCCCCAGTACACGTCTTCGTCCGGCTCCCACACATCGG‑CACGCCCGCCGGAGAAGCCGAAGGTCTTGAAGCCCATGGATTCGAGGGCAACATTACCCGTGAGTACGATCAGGTCGGCCCAGGAAATCTTG‑T‑TGCCGTACTTCTGCTTGATCGGCCACAGCAGGCGCCGGGCCtt                                                                                                                                                                                                                                                           <  1:522412/301‑1 (MQ=255)
                                                 gATTCTGCTCTTCACCATGCTTGGCAGTCTCGGCCACCAGGTCACCCTGCCCCGGCGGCTGCGCCTTGACCTGATCCTTGCCGTAACGGGTATCACCCCCTAGCCAGACCTTTTCCGACCCCCAGTACACGTCTTCGTCCGACTCCCACACATCGG‑CACGCCCGCCGGAGAAGCCGAAGGTCTTGAAGCCCATGGATTCGACGGCAACATTACCCGTGAGTACGATCAGGTAGGCCCAGGAAACATTG‑T‑TGCCGTACTTCTGCTTGATCGGCCACAGCAAGCGCCGGCCCTTGTACAGGCt                                                                                                                                                                                                                                                  >  2:24078/1‑301 (MQ=255)
                                                   ttCTGCTCTTCACCATGCTTGGCAGGCTCGGCCACCAGGTCACCCTGCCCCGGCGGCTGCGCCTTGACCTGATCCTTGCCGTAACGGGTATCACCCCCTAGCCAGACCTTTTCCGACCCCCAGTACACGTCTTCGTCCGGCTCCCACACATCGG‑CACGCCCGCCGGAGAAGCCGAAGGTCTTGAAGCCCATGGATTCGAGGGCAACATTACCCGTGAGTACGATCAGGTCGGCCCAGGAAATCTCG‑T‑TGGCGTACTTCTGCTTGATAGGCCACAGCAGGCCACGGGCCTTGTCCAgcctct                                                                                                                                                                                                                                                >  1:446933/1‑296 (MQ=255)
                                                              aCCATGCTTGGCAGGCTCGGCCACCAGGTCACCCTGCCCCGGCGGCTGCGCCTTGACCTGATCCTTGCCGTAACGGGTATCACCCCCTAGCCAGACCTTTTCCGACCCCCAGTACACGTCTTCGTCCGGCTCCCACACATCGG‑CACGCCCGCCGGAGAAGCCGAAGGTCTTGAAGCCCATGGATTCGAGGGCAACATTACCCTTGAGTACGATCAGGTCGGCCCAGGAAATCTTG‑T‑TGCCGTACTTCTGCTTGATCGGCCACAGCAGGCGCCGGGCCTTGTCCAGGCTGAAGTTGTCCGGc                                                                                                                                                                                                                                     >  2:235069/1‑301 (MQ=255)
                                                                    gttGGCAGGCTCGGCCACCCGGTCACCCTGCCCCGGCGGCTGCGCCTTGACCTGATCCTTGCCGTAACGGGTATCACCCCCTAGCCAGAACTTTTCCGACCCCCAGTACACGTCTTCGTCCGGCTCCCACACATCGG‑CACGCCCGCCGGAGAAGCCGAAGGTCTTGAAGCCCATGGATTCGAGGGCAACATTACCCGTGAGTACGATCAGGTCGGCCCAGGAAATCTTG‑T‑TGCCGTACTTCTGCTTGATCGGCCACAGCAGGCGCCGGGCCTTGTCCAGGCTGACGTTGTCCGGCCAGCTg                                                                                                                                                                                                                               <  2:699012/300‑1 (MQ=255)
                                                                    cggtgCAGGCTCGGCCACCAGGTGACCCTGCCCCGGCGGCTGCGCCTTGACCTGATCCCTGCCGTAACGGGTATCACCCCCTAGCCAGACCTTTTCCGACCCCCAGTACACGTCTTCGTCCGGCTCCCACACATCGG‑CACGCCCGCCGGAGAAGCCGAAGGTCTTGAAGCCCATGGATTCGAGGGAAACATTACCCGTGAGTACGATCAGGTCGGCCCAGGAAATCTTG‑T‑TGCCGTACTTCTGCTTGATCGGCCACAGCAGGCGCCGGGCCTTGTCCAGGCTGACGTTGTCCGGCCAGCTg                                                                                                                                                                                                                               <  2:189026/297‑1 (MQ=255)
                                                                                ggCCACCAGGTCACCCTGCCCCGGCGGCTGCGCCTTGACCTGATCCTTGCCGTAACGGGTATCACCCCCTAGCCAGACCTTTTCCGACCCCCAGTACACGTCTTCGTCCGGCTCCCACACATCGG‑CACGCCCGCCGGAGAAGCCGAAGGTCTTGAAGCCCATGGATTCGAGGGCAACATTACCCGTGAGTACGATCAGGTCGGCCCAGGAAATCTTG‑T‑TTCCGTACTTCTGCTTGAACGACAAAAGCAGGCGCCGGGCCTTGTCCAGGCTGACGTTGTCCGGCCAGCTGTTGAGCgggtcg                                                                                                                                                                                                                   >  1:153350/1‑298 (MQ=255)
                                                                                     ccaGGTCACCCTGCCCCGGCGGCTGCGCCTTGACCTGATCCTTGCCGTAACGGGTATCACCCCCTAGCCAGACCTTTTCCGACCCCCAGTACACGTCTTCGTCCGGCTCCCACACATCGG‑CACGCCCGCCGGAGAAGCCGAAGGTCTTGAAGCCCATGGATTCGAGGGCAACATTACCCGTGAGTACGATCAGGTCGGCCCAGGAAATCTTG‑T‑TTCCGTACTTCTGCTTGATCGGCCACAGCAGGCGCCGGGCCTTGTCCAGGCTGACGTTGTCCGGCCAGCTGTTGAGCGGGGCGAAGCg                                                                                                                                                                                                              >  1:267800/1‑301 (MQ=255)
                                                                                          ccACCATGCCTCGGCGGCTGCGCCTTGACCTGATCCTTGCCTTAACGGGTATCACCCCCTAGCCAGCCCTTTTCCGACCCCCAGTACACGTCTTCGTCCGGCTCCCACACATCGG‑CACGCCCGCCGGAGAAGCCGAAGGTCTTGAAGCCCATGGATTCGAGGGCAACATTACCCGTGAGTACGATCAGGTCGGCCCAGGAAATCTTG‑T‑TGCCGTACTTCTGCTTGATCGGCCACAGCAGGCGCCGGGCCTTGTCCAGGCTGACGTTGTCCGGCCAGCTGTTGAGCGGGGCGAAGCGTTGCt                                                                                                                                                                                                         <  1:398573/300‑1 (MQ=255)
                                                                                               cTGCCCAGGCGGCGGCGCCTTGACCTGATCCTTGACGTAACGGGGTTCACCCCCTAGCCAGACCTTTTCCGACCCCCAGTACACGTCTTCGTCCGGCTCCCACACATCGG‑CACGCCCGCCGGAGAAGCCGAAGGTCTTGAAGACCATGGATTCGAGGGCAACATTACCCGTGAGTACGATCAGGTCGGCCCAGGAAATCTTG‑T‑TGCCGTACTTCTGCTTGATCGGCCACAGCAGGCGCCGGGCCTTGTCCAGGCTGACGTTGTCCGGCCAGCTGTTGAGCGGGGCGAAGCGTTGCTGGCCa                                                                                                                                                                                                    <  2:91173/301‑1 (MQ=255)
                                                                                                      ggcggcTGCGCCTTGACCTGATCCTTGCCGTAACGGGTATCACCCCCTAGCCAGACCTTTTCCGACCCCCAGTACACGTCTTCGTCCGGCTCCCACACATCGG‑CACGCCCGCCGGAGAAGCCGAAGGTCTTGAAGCCCATGGATTCGAGGGCAACATTACCCGTGAGTACGATCAGGTCGGCCCAGGAAATCTTG‑T‑AGCCGTACTTCTGCTTGATCAACCACAGCAGGCGCCGGGCCTTGTCCAGGCTGACGTTGTCCGGCCAGCTGTTGACCGGGGCGAAGCGTTGCTGGCCAGCGcccc                                                                                                                                                                                             >  1:191669/1‑300 (MQ=255)
                                                                                                                     aCCTGATCCTTGCCGTAACGGGTATCACCCCCTAGCCAGACCTTTTCCGACCCCCAGTACACGTCTTCGTCCGGCTCCCACACATCGG‑CACGCCCGCCGGAGAAGCCGAAGGTCTTGAAGCCCATGGATTCGAGGGCAACATTACCCGTGAGTACGATCAGGTCGGCCCAGGAAATCTTG‑T‑TGCCGTACTTCTGCTTGATCGGCCACAGCAGGCGCCGTGCCTTGTCCAGGCTGACGTTGTCCGGCCAGCTGTTGAGCGGGCCGAAGCGTCGCTGGCCACAGCCCGCACCACCGCGGCCac                                                                                                                                                                              >  2:450055/1‑300 (MQ=255)
                                                                                                                       cTGATCCTTGCCGTAACGGGTATCACCCCCTAGCCAGACCTTTTCCGACCCCCAGTACACGTCTTCGTCCGGCTCCCACACATCGG‑CACGCCCGCCGGAGAAGCCGAAGGTCTTGAAGCCCATGGATTCGAGGGCAACATTACCCGTGAGTACGATCAGGTCGGCCCAGGAAATCTCG‑TATCCCG‑ACTTCTGCTTGATCGGCCACAGCAGGCGCCGGGCCTTCTCCAGGCTGACGTTGTCCGGCCAGCTGTTGAGCCGGGCGAAGCGTTGCTGGCCAGAGCCCGCACCACCGCGGCCATcc                                                                                                                                                                            >  1:91095/1‑300 (MQ=255)
                                                                                                                                                      aGCCAGACCTTTTCCGACCCCCAGTACACGTCTTCGTCCGGCTCCCACACATCGG‑CACGCCCGCCGGAGAAGCCGAAGGTCTTGAAGCCCATGGATTCGAGGGCAACATTACCCGTGAGTACGATCAGGTCGGCCCAGGAAATCTCG‑TATGCCGT‑CTTCTGCTTGAAAAAAAACAGCAGGCGCCATGCCTTGTACAGGCTGACGTTGTCCGGCCAGCTGTTGAGCTGGGCGAAGCTTTGCTTGCCATAGCCCGCACCACTGCGGCCATCACCAATACGGTAGGATCCGgcgccatccccgg                                                                                                                                             >  1:439086/1‑292 (MQ=255)
                                                                                                                                                                  caagacccacagtcaccGTCTTCGGCCGGCTCCCACACATCGGCCAC‑CCCGCCTGACAAGCCGAAGGTCTTGAATCCCATGGCTTGGAGGGCAACATTACCCGTGGGTGCGATCAGGTCGGCCCAGGAGATCTTG‑T‑TGCCGTACTTCTGCTTGATCGGCCACAGCAGGCGCCGGGCCTTGTCCAGGCTGACGTTGTCCGGCCAGCTGTTGAGCGGGGCGAAGCGTTGCTGGCCAGAGCCCGCACCACCGCGGCCATCACCAATGCGGTAGGTGCCGGCGCTATGCCAGGCCATGCGGATGa                                                                                                                                 <  2:566449/285‑1 (MQ=255)
                                                                                                                                                                     gACCCCCAGTACACGTCTTCGTCCGGCTCCCACACATCGG‑CACGCCCGCCGGAGAAGCCGAAGGTCTTGAAGCCCATGGATTCGAGGGCAACATTACCCGTGAGTACGATCAGGTCGGCCCAGGAAATCTCGTT‑TGCCG‑ACTTCTGCTTGATAAGCCAAAGCAGGCGCCGGGCCTTGTCCAGGCTGACGTTGTCCGGCCAGCTGTTGAGCGGGGCGAAGAGTTGCTGGCCAGAGCCCGCACCACCGCGGCCATCACCAATGCGGTATGTTCCGGCTCTATGCCAGGCCATGCGGATGcaca                                                                                                                              >  1:61657/1‑301 (MQ=255)
                                                                                                                                                                            ggTACAGGTCTTCGTCCGGCTCCCACACATCGG‑CACGCCCGCCGGAGAAGCCGAAGGTCTTGAAGCCCATGGATTCGAGGGCAACATTACCCGTGAGTACGATCAGGTCGGCCCAGGAAATCTTG‑T‑TGCCGTACTTCTGCTTGATCGGCCACAGCAGGCGCCGGGCCTTGTCCAGGCTGACGTTGTCCGGCCAGCTGTTGAGCGGGGCGAAGCGTTGCTGGCCAGAGCCCGCACCACCGCGGCCATCACCAATGCGGTAGGTGCCGGCGCTATGCCAGGCCATGCGGATGAACAGCGGGcc                                                                                                                       <  1:450245/300‑1 (MQ=255)
                                                                                                                                                                               acatgtgttCGTCCGGATCACACACATCGT‑GACGCCCGCCGGAGAAGCCGAAGTTCTTGAAGCCCATGGAGTCGAGGGCAACATTACCCGTGAGTACGAGCAGGTCGGCCCAGGAAATCTTG‑T‑TGCCGTACTTCTGCTTGAGCGGCCACAGCAGGCGCCGGGCCTTGTCCAGGCTGACGTTGTCCGGCCAGCTGTTGAGCGGGGCGAAGCGTTGCTGGCCAGAGCCCGCACCACCGCGGCCATCACCAATGCGGTAGGTGCCGGCGCTATGCCAGGCCATGCGGATGAACAGCGGGCCGTa                                                                                                                    <  1:235159/294‑1 (MQ=255)
                                                                                                                                                                                              gCTCCCACACATCGG‑CACGCCCGCCGGAGAAGCCGAAGGTCTTGAAGCCCATGGATTCGAGGGCAACATTACCCGTGAGTACGATCAGGTCGGCCCAGGAAATCTCG‑T‑TTCCGGTCTTCTGCTTGATTGAACACAGCAGGCGCCGGGCCTCGTCCAGGCTGACGTTGCCCGGCCAGCTGTCGAGCGGGGCGAAGCGTTGCTGGCCAGAGCCCGCAGCACCGCGGCAATCCCCAATGCGGTAAGTGCCGGCGCTATTCCAGGCCATGTGGCTGAACAGCCGTCCGTAGGGGCCGAAGTctgg                                                                                                     >  1:536611/1‑298 (MQ=255)
                                                                                                                                                                                                                            aaGCCGAAGGTCTTGAAGCCCATGGATTCGAGGGCAACATTACCCGTGAGTACGATCAGGTCGGCCCAGGAAATCTTG‑T‑TGCCGTACTTCTGCTTGATCGGCCACAGCAGGCGCCGGGCCTTGTCCAGGCTGACGTTGTCCGGCCAGCTGTTGAGCGGGGCGAAGCGTTGCTGGCCAGAGCCCGCACCACCGCGGCCATCACCAATGCGGTAGGTGCCGGCGCTATGCCAGGCCATGCGGATGAACAGCGGGCCGTAGTGGCCGAAGTCGGCCGGCCACCAATCCTGCGAGTGGGtccgca                                                                        >  2:498200/1‑297 (MQ=255)
                                                                                                                                                                                                                               ccGAAGGTCTTGAAGCCCATGGATTCGAGGGCAACATTACCCGTGAGTACGATCAGGTCGGCCCAGGAAATCTCG‑TATG‑CGTACTTCTGCTTGAAAGACAAAAGCAGGCGCCGGGCCTTGTCCAGGCTGACGTTGTCCGGCCAGCTGTTGAGCGGGGCGAAGCGTTGCTGGCCAGAGCCCGCACCACCGCGGCCATCACCAATGCGTTAGGTGCCGGCGCTATCCCAGGCCATGCGGATGCACAGCGGGCCGTAGTGGTCGCATTCGTCCTGCCCCCCATCCTGCGAGCCGGtcatcatgc                                                                     >  1:61880/1‑298 (MQ=255)
                                                                                                                                                                                                                                   aGGTCTTGAAGCCCATGGATTCGAGGGCAACATTACCCGTGAGTACGATCAGGTCGGCCCAGGAAATCTTG‑T‑TGCCGTACTTCTGCTTGATCGGCCACAGCAGGCGCCGGGCCTTGTCCAGGCTGACGTTGTCCGGCCAGCTGTTGAGCGGGGCGAAGCGTTGCTGGCCAGAGCCAGCACCACCGCGGCCATCACCAATGCGGTAGGTGCCGGCGCTATGCCAGGCCATGCGGATGAACAGCGGGCCGTAGTGGCCGAAGTCGGCCGGCCACCAATCCTGCGAGTCGGTCATCAGGGCggg                                                                 >  1:306502/1‑300 (MQ=255)
                                                                                                                                                                                                                                      tCTTGAAGCCCATGGATTCGAGGGCAACATTACCCGTGAGTACGATCAGGTCGGCCCAGGAAATCTGG‑T‑TGCCGTACTTCTGCTTGATCGGCCACAGCAGGCGCCGGGCCTTGTCCAGGCTGACGTTGTCCGGCCAGCTGTTGAGCGGGGCGAAGCGTTGCTGGCCAGAGCCAGCACCACCGCGGCCATCACCAATGCGGTAGGTGCCGGCGCTATGCCAGGCCATGCGGATGAACAGCGGGCCGTAGTGGCCGAAGTCGGCCGGCCACCAATCCTGCGAGTCGGTCATCAGGGCGGTCAg                                                              <  2:306383/301‑1 (MQ=255)
                                                                                                                                                                                                                                       cTTGAAGCCCATGGATTCGAGGGCAACATTACCCGTGAGTACGATCAGGTCGGCCCAGGAAATCTCG‑TATG‑CGTACTTCTGCTTGAACGGACAAAGCAGGCGCCGGGCCTTGTCCAGGCTGACGTTGTCCGGCCAGCTGTTGAGCGGGGCGAAGCGTTGCTGGCCAGAGCCCGCACCACCGCGGCCATCACCAATGCGGTAGGTGCCGCCGCTATGCCAGGCCATGCGGATGAACAGCGGGCCGCAGTGGCCGCAGTCGGCCGGCTACCAATCCTGCGAGTCGGCCATCAGGCCGGTCAgg                                                             >  1:249278/1‑301 (MQ=255)
                                                                                                                                                                                                                                            aGCCCATGGATTCGAGGGCAACATTACCCGTGAGTACGATCAGGTCGGCCCAGGAAATCTCG‑TATGCCGT‑CTTCTGCTTGAAAAAAAAAAGCAGGCGCCGGGCCTTGTCCAGTCTGACGTTGTCCGGCAAGCTGTTGAGCGGGGCGAAGCGTTGCTGGCCAGAGCCCGCACCACCGCGGCCATCACCAATGCGGTAGGTGCCGGCGCTATGCCAGGCCATGCGGATGAACAGCGGGCCGTAGTTGACGAAGTCGGCCGGCCAACAATACTGCGAGTCGGTCATTCGGGCGGTCAGGtcttt                                                        >  1:352658/1‑301 (MQ=255)
                                                                                                                                                                                                                                                   ggagtCGAGGGAAACATTACCCAGAATTACGATCAGGTCGGCCCAGGAAATCTTG‑T‑TGCCGTACTCCTGCTGGATCGGCCACGGCAGGCGCCGGGCCTTGTCCAGGCTGACGTTGTCCGGCCAGCTGTTAGGCGGGGCGAAGCGTTGCTGGCCAGAGCCCGCACCACCGCGGCCATCACCAATGCGGTAGGTGCCGGCGCTATGCCAGGCCATGCGGATGAACAGCGGGCCGTAGTGGCCGAAGTCGGCCGGCCACCAATCCTGCGAGTCGGTCATCAGGGCGGTCAGGTCTTTCTTCAgg                                                 <  2:267695/297‑1 (MQ=255)
                                                                                                                                                                                                                                                   gaaTTCGAGGGCAACATTACCCGTGAGTACGATCAGGTCGGCCCAGGAAATCTTG‑T‑TGCCGTACTTCTGCTTGATCGGCCACAGCAGGCGCCGGGCCTTGTCCAGGCTGACGTTGTCCGGCCAGCTGTTGAGCGGGGCGAAGCGTTGCTGGCCAGAGCCCGCACCACCGCGGCCATCACCAATGCGGTAGGTGCCGGCGCTATGCCAGGCCATGCGGATGAACAGCGGGCCGTAGTGGCCGAAGTCGGCCGGCCACCAATCCTGCGAGTCGGTCATCAGGGCGGTCAGGTCTTTCTTCAgg                                                 <  1:498405/299‑1 (MQ=255)
                                                                                                                                                                                                                                                    gATTCGAGGGCAACATTACCCGTGAGTACGATCAGGTCGGCCCAGGAAATCTTG‑T‑TGCCGTACTTCTGCTTGATAGGCCACAGCAGGCGCCGGGCCTTGTCCAGGCTGACGTTGTCCGGCCAGCTGTTGAGCGGGGCGAAGCGTTGCTGGCCAGAGCCCGCACCACCGCGGCTATCACCAATGCGGTAGGTGCCGGCGCTATGCCAGGCCATGCGGATGAACAGCGGGCCGTAGTGTCCGAAGTCGGCCGGCAACCAATCCTGCGAGGCGGTCATCAGGGCGGACAGGTCTTTCTTCAggg                                                >  1:389802/1‑301 (MQ=255)
                                                                                                                                                                                                                                                     aTTCGAGGGCAACATTACCCGTGAGTACGATCAGGTCGGCCCAGGAAATCTCGTT‑TGCCGT‑CTTCTGCTTGATCGACCACAGCAGGCGCCGGGCCTTGTCCAGGCTGACGTTGTCCGGCCAGCTGTTGAGCGGGGCGAAGCGTTGCTGGCCAGAGCCCGCACCACCGCGGCCATCACCAATGCGGTAGGTGCCGGCGCTATGCCAGTCCATGCGGATGAACAGCGGGCCGTAGTGGCCGAAGTCGGCGGGCCACCAATCCTGCGAGTCCACCATCACGGCGGTCGGTTCTTTTTTCCGGGc                                               >  1:264625/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                            aGTACGATCAGGTCGGCCCAGGAAATCTCG‑T‑TTGCGTACTTCTGCTTGATCGGCCACAGCAGGCGCCGGGCCTTGTCCAGGCTGACGTTGTCCGGCCAGCTGTTGAGCGGGGCGAAGCGTTGCTGGCCAGAGCCCGCACCACCGCGGCCATCACCAATGCGGTAGGTGCCGGCGCTATGCCAGGCCATGCGGATGAACAGCGGGCCGTAGTGGCCGAAGTCGGCCGGCCACCAATCCTGCGAGTCGGTCATCAGGGCGGTCAGGTCTTTCTTCAGGGCCTGAAAGTCGAGGATCTTGAACg                        >  1:564326/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                             gTACGATCAGGTCGGCCCAGGAAATCTTG‑T‑TGCCGTACTTCTGCTTGATCGGCCACAGCAGGCGCCGGGCCTTGTCCAGGCTGACGTTGTCCGGTCAGCTGTTGAGCGGGGCGAAGCGTTGCTGGCCAGAGCCCGCACCACCGCGGCCATCACCAATGCGGTAGGTGCCGGCGCTATGCCAGGCCATGCGGATGAACAGCGGGCCGTAGTGGCCGAAGTCGGCCGGCCACCAATCCTGCGAGTCGGTCATCAGGGCGGTCAGGTCTTTCTTCAGGGCCTGAAAGTCGAGGCTCTTGAACGc                       <  1:522953/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                  aaaTCTTG‑T‑TGCCGTACTTCTGCTTGATCGGCCACAGCAGGCGCCGGGCCTTGTCCAGGCTGACGTTGTCCGGCCAGCTGTTGAGCGGGGCGAAGCGTTGCTGGCCAGAGCCCGCACCACCGCGGCCATCACCAATGCGGTAGGTGCCGGCGCTATGCCAGGCCATGCGGATGAACAGCGGGCCGTAGTGGCCGAAGTCGGCCGGCCACCAATCCTGCGAGTCGGTCATCAGGGCGGTCAGGTCTTTCTTCAGGGCCTGAAAGTCGAGGCTCTTGAACGCCTTGTCTTAGTCGACGTCCgg  >  1:9338/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                                        
TGCACGGCGGCCAGCGGGTTTTCCAGGTTACGCTCGGCCGACAGGTCGCGATTCTGCTCTTCACCATGCTTGGCAGGCTCGGCCACCAGGTCACCCTGCCCCGGCGGCTGCGCCTTGACCTGATCCTTGCCGTAACGGGTATCACCCCCTAGCCAGACCTTTTCCGACCCCCAGTACACGTCTTCGTCCGGCTCCCACACATCGG‑CACGCCCGCCGGAGAAGCCGAAGGTCTTGAAGCCCATGGATTCGAGGGCAACATTACCCGTGAGTACGATCAGGTCGGCCCAGGAAATCTTG‑T‑TGCCGTACTTCTGCTTGATCGGCCACAGCAGGCGCCGGGCCTTGTCCAGGCTGACGTTGTCCGGCCAGCTGTTGAGCGGGGCGAAGCGTTGCTGGCCAGAGCCCGCACCACCGCGGCCATCACCAATGCGGTAGGTGCCGGCGCTATGCCAGGCCATGCGGATGAACAGCGGGCCGTAGTGGCCGAAGTCGGCCGGCCACCAATCCTGCGAGTCGGTCATCAGGGCGGTCAGGTCTTTCTTCAGGGCCTGAAAGTCGAGGCTCTTGAACGCCTTGGCATAGTCGAAGTCCGG  >  AE015451/4169177‑4169766

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