Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*1487,7240CA42.7% 140.7 / 30.9 75intergenic (‑13558/‑62969)unknown/unknowntRNA‑SeC(p)‑TCA/SSU rRNA ## 16S rRNA, small subunit ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base C (39/4);  new base A (0/32);  total (39/36)
Fisher's exact test for biased strand distribution p-value = 3.77e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GGTTTACAACAACCCTGAATATGCCGGCGAACGCATGATGTTGCTGTTTACCAATCCGACGGACGTATTGCGTCTTTCCGATGCGGGCGTTGATTTT‑AAAT‑CCATCAATGTG‑GGTGGCATGGCGTATCGCGATGGTAAAAAAATGATTACCAA‑TGCGGTATCCGTGGACGACAAAGATATTG‑CCGCCTTCAA‑A‑GCCCTTGGTGCC‑AAA‑GGCATTGAGCTGGACGTGC‑GTAAAGTTTCAAATGACAACCGCCAGTACATGATGGATTTATTGAAGAAAAATAACTTACTTTAAAGGAAATGATTATGGAAATTACCACTTTACAAATTATTCTGGTGTTTCTGGTCGCCTGTATTTGCGGCGCAGGATCGATTCTTGACGAATGGCAAACCCAACGCCCGTTAATTGCCTGTACCTTAATCGGCTTAGTGTTGGGCGATATTCAAACCGGTATCATTATCGGCGGGACCTTGGAA  >  1/487488‑487961
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                             
ggTTTACAACAACCCTGAATATGCCGGCGAACGCATGATGTTGCTGTTTACCAATCCGACGGACGTATTGCGTCTTTCCGATGCGGGCGTTGATTTT‑AAAT‑CCATCAATGTG‑GGTGGCATGGCGTATCGCGATGGTAAAAAAATGATTACCAA‑TGCGGTATCCGTGGACGACAAAGATATTG‑CCGCCTTCAA‑A‑GCCCTTGGTGCC‑AAA‑GGCATTGAGCTGGACGTGC‑GTAAAGTTTCAAATGACAACCGc                                                                                                                                                                                                                                  >  1:165869/1‑250 (MQ=255)
ggTTTACAACAACCCTGAATATGCCGGCGAACGCATGATGTTGCTGTTTACCAATCCGACGGACGTATTGCGTCTTTCCGATGCGGGCGTTGATTTT‑AAAT‑CCATCAATGTG‑GGTGGCATGGCGTATCGCGATGGTAAAAAAATGATTACCAA‑TGCGGTATCCGTGGACGACAAAGATATTG‑CCGCCTTCAA‑A‑GCCCTTGGTGCC‑AAA‑GGCATTGAGCTGGACGTGC‑GTAAAGTTTCAAATGACAACCGc                                                                                                                                                                                                                                  >  1:583523/1‑250 (MQ=255)
           acCCTGAATATGCCGGCGAACGCATGATGTTGCTGTTTACCAATCCGACGGACGTATTGCGTCTTTCCGATGCGGGCGTTGATTTT‑AAAT‑CCATCAATGTG‑GGTGGCATGGCGTATCGCGATGGTAAAAAAATGATTACCAA‑TGCGGTATCCGTGGACGACAAAGATATTG‑CCGCCTTCAA‑A‑GCCCTTGGTGCC‑AAA‑GGCATTGAGCTGGACGTGC‑GTAAAGTTTCAAATGACAACCGCCAGTACatgatg                                                                                                                                                                                                                      >  1:191954/1‑251 (MQ=255)
           acCCTGAATATGCCGGCGAACGCATGATGTTGCTGTTTACCAATCCGACGGACGTATTGCGTCTTTCCGATGCGGGCGTTGATTTT‑AAAT‑CCATCAATGTG‑GGTGGCATGGCGTATCGCGATGGTAAAAAAATGATTACCAA‑TGCGGTATCCGTGGACGACAAAGATATTG‑CCGCCTTCAA‑A‑GCCCTTGGTGCC‑AAA‑GGCATTGAGCTGGACGTGC‑GTAAAGTTTCAAATGACAACCGCCAGTACatgatg                                                                                                                                                                                                                      >  1:757812/1‑251 (MQ=255)
           acCCTGAATATGCCGGCGAACGCATGATGTTGCTGTTTACCAATCCGACGGACGTATTGCGTCTTTCCGATGCGGGCGTTGATTTT‑AAAT‑CCATCAATGTG‑GGTGGCATGGCGTATCGCGATGGTAAAAAAATGATTACCAA‑TGCGGTATCCGTGGACGACAAAGATATTG‑CCGCCTTCAA‑A‑GCCCTTGGTGCC‑AAA‑GGCATTGAGCTGGACGTGC‑GTAAAGTTTCAAATGACAACCGCCAGTACatgatg                                                                                                                                                                                                                      >  1:467221/1‑251 (MQ=255)
           acCCTGAATATGCCGGCGAACGCATGATGTTGCTGTTTACCAATCCGACGGACGTATTGCGTCTTTCCGATGCGGGCGTTGATTTT‑AAAT‑CCATCAATGTG‑GGTGGCATGGCGTATCGCGATGGTAAAAAAATGATTACCAA‑TGCGGTATCCGTGGACGACAAAGATATTG‑CCGCCTTCAA‑A‑GCCCTTGGTGCC‑AAA‑GGCATTGAGCTGGACGTGC‑GTAAAGTTTCAAATGACAACCGCCAGTACatgat                                                                                                                                                                                                                       >  1:662247/1‑250 (MQ=255)
                gAATATGCCGGCGAACGCATGATGTTGCTGTTTACCAATCCGACGGACGTATTGCGTCTTTCCGATGCGGGCGTTGATTTT‑AAAT‑CCATCAATGTG‑GGTGGCATGGCGTATCGCGATGGTAAAAAAATGATTACCAA‑TGCGGTATCCGTGGACGACAAAGATATTG‑CCGCCTTCAA‑A‑GCCCTTGGTGCC‑AAA‑GGCATTGAGCTGGACGTGC‑GTAAAGTTTCAAATGACAACCGCCAGTACATGATGGAttt                                                                                                                                                                                                                 >  1:80463/1‑251 (MQ=255)
                gAATATGCCGGCGAACGCATGATGTTGCTGTTTACCAATCCGACGGACGTATTGCGTCTTTCCGATGCGGGCGTTGATTTT‑AAAT‑CCATCAATGTG‑GGTGGCATGGCGTATCGCGATGGTAAAAAAATGATTACCAA‑TGCGGTATCCGTGGACGACAAAGATATTG‑CCGCCTTCAA‑A‑GCCCTTGGTGCC‑AAA‑GGCATTGAGCTGGACGTGC‑GTAAAGTTTCAAATGACAACCGCCAGTACATGATGGAt                                                                                                                                                                                                                   >  1:487619/1‑249 (MQ=255)
                      gCCGGCGAACGCATGTTGTTGCTTTTTACCAACCCGCCGGACGTATTGCGTCTTTCCGATGGGGGCGTTATTTTT‑AAAT‑CCATCAATGTG‑GGTGGCAGGGCGTATCGCGATGGAAAAAAAATGATTCCCAA‑TGCGGTATCCGTGGACGAAAAAGATATTG‑CCGCCTTAAA‑A‑GCCCTTGGTGCA‑AAA‑GGCATTGAGCTGGCCGTGC‑GTAAAGTTTCAAATGACAACCGCCAGTACATGATGGATTTATTg                                                                                                                                                                                                             <  1:147261/249‑1 (MQ=255)
                          gCGAACGCATGATGTTGCTGTTTACCAATCCGACGGACGTATTGCGTCTTTCCGATGCGGGCGTTGATTTT‑GAAT‑CCATCAATGTG‑GGTGGCATGGCGTATCGCGATGGTAAAAAAATGATTACCAA‑TGCGGTATCCGTGGACGACAAAGATATTG‑CCGCCTTCAA‑A‑GCCCTTGGTGCC‑AAA‑GGCATTGAGCTGGACGTGC‑GTAAAGTTTCAAATGACAACCGCCAGTACATGATGGATTTATTGAAGaaa                                                                                                                                                                                                       >  1:168356/1‑251 (MQ=255)
                             aaCGCTTGATGTTGCTTTTTCCCAATCCGACGGGCTTTTTGCGTCTTGCGGATGGGGGGTTTGTTTTT‑AAAT‑CCATAAATGGG‑GGTGGCATGGCGTATCGCGATGAAAAAAAAATGATTACAAA‑GGCGGTATCCGGGGACGAAAAAGATATTC‑CCGCCTTAAA‑A‑GCCCTGGGTGCA‑AAG‑GGCTTTGGGCGGGACGTGC‑GTAAATTTTAAAATGACAACCGCCAGTACATGATGGATTTATTGAAGAAAAAt                                                                                                                                                                                                    <  1:826796/251‑1 (MQ=255)
                             aaCGAATGATGTTGTTGTTTACCAATCCGACGGACGTATTGCGTCTTTCGGAGGGGGGAGTTAATTTT‑AAAT‑CCATCAATGTG‑GGTGGCATGGCGTAGCGCGATGGAAAAAAAATGATTCCAAA‑TGCGGTACCCGGGGCCAAAAAAGATATTG‑CCGCTTACAA‑A‑CCCCTTGGTGCA‑AAA‑GGCATTGAGCTGGCCGTGC‑GAAAATTTTAAAATGACAACCGCCAGTACATGATGGATTTATTGAAGaaaaa                                                                                                                                                                                                     <  1:106893/250‑1 (MQ=255)
                               cGCATGATGTTGCTGTTTACCAATCCGACGGACGTATTGCGTCTTTCCGATGCGGGCGTTGATTTT‑AAAT‑CCATCAATGTG‑GGTGGCATGGCGTATCGCGATGGTAAAAAAATGATTACCAA‑TGCGGTATCCGTGGACGACAAAGATATTG‑CCGCCTTCAA‑A‑GCCCTTGGTGCC‑AAA‑GGCATTGAGCTGGACGTGC‑GTAAAGTTTCAAATGACAACCGCCAGTACATGATGGATTTATTGAAGAAAAATaa                                                                                                                                                                                                  >  1:329788/1‑251 (MQ=255)
                                   tgatgTTGCTGTTTACCAATCCGACGGACGTATTGCGTCTTTCCGATGCGGGCGTTGATTTT‑AAAT‑CCATCAATGTG‑GGTGGCATGGCGTATCGCGATGGTAAAAAAATGATTACCAA‑TGCGGTATCCGTGGACGACAAAGATATTG‑CCGCCTTCAA‑A‑GCCCTTGGTGCC‑AAA‑GGCATTGAGCTGGACGTGC‑GTAAAGTTTCAAATGACAACCGCCAGTAc                                                                                                                                                                                                                            >  1:282805/1‑221 (MQ=255)
                                   tgatgTTGCTGTTTACCAATCCGACGGACGTATTGCGTCTTTCCGATGCGGGCGTTGATTTT‑AAAT‑CCATCAATGTG‑GGTGGCATGGCGTATCGCGATGGTAAAAAAATGATTACCAA‑TGCGGTATCCGTGGACGACAAAGATATTG‑CCGCCTTCAA‑A‑GCCCTTGGTGCC‑AAA‑GGCATTGAGCTGGACGTGC‑GTAAAGTTTCAAATGACAACCGCCAGTACATGATGGAt                                                                                                                                                                                                                   >  1:667154/1‑230 (MQ=255)
                                     atgTTGCTGTTTACCAATCCGACGGACGTATTGCGTCTTTCCGATGCGGGCGTTGATTTT‑AAAT‑CCATCAATGTG‑GGTGGCATGGCGTATCGCGATGGTAAAAAAATGATTACCAA‑TGCGGTATCCGTGGACGACAAAGATATTG‑CCGCCTTCAA‑A‑GCCCTTGGTGCC‑AAA‑GGCATTGAGCTGGACGTGC‑GTAAAGTTTCAAATGACAACCGCCAGTACATGATGGATTTATTGAAGAAAAATAactta                                                                                                                                                                                              >  1:712032/1‑249 (MQ=255)
                                           ttGTTTCCCAATCCGACGGACGTTTTGCTTCTTTCCGTGGGGGGCGTTGTTTTT‑AAAC‑CCATCAATGGG‑GGTGGCAGGGCTTATCGCGATGGAAAAAAAATGATTACCAA‑TGCGGTATCCGTGGCCAAAAAAGATATTG‑CCGCCTACAA‑A‑CCCCTTGGTGCA‑AAA‑GGCATTGAGCTGGACGGGC‑GAAAAGTTTAAAATGACAACCGCCAGTACATGATGGATTTATTGAAGAAAAATAACTTACTTTAAAg                                                                                                                                                                                      <  1:694173/250‑1 (MQ=255)
                                           cTGTTTACCAATCCGACGGACGTATTGCGTCTTTCCGATGCGGGCGTTGATTTT‑AAAT‑CCATCAATGTG‑GGTGGCATGGCGTATCGCGATGGTAAAAAAATGATTACCAA‑TGCGGTATCCGTGGACGACAAAGATATTG‑CCGCCTTCAA‑A‑GCCCTTGGTGCC‑AAA‑GGCATTGAGCTGGACGTGC‑GTAAAGTTTCAAATGACAACCGCCAGTACATGATGGATTTATTGAAGAAAAATAACTTACTTTAAAg                                                                                                                                                                                      >  1:803843/1‑251 (MQ=255)
                                            tgtttccaaacccGACGGACGTATGGCTTCTTTCGGAGGGGGGGGTTTATTTA‑AAAC‑CCATAATTGTG‑GGTGGCATGGCTTATCGCGATGGAAAAAAAATGATTACCAA‑TGCGGTATCCTTGGACAAAAAAGATATTG‑CCGCTTTAAA‑A‑GCCCTTGGTCCC‑AAA‑GGAATTGAGCTGTACGTGC‑GTAAAGTTTAAAATGACAACCGCCAGTACATGATGGATTTATTGAAGAAAAATAACTTACTTTAAAgg                                                                                                                                                                                     <  1:20690/240‑1 (MQ=255)
                                            tGTTTACCAACCCGACGGGCGTTTGGCGTC‑TTCCGGAGGGGGGTTTGATTTTAAAAC‑CCATAAAGGTG‑GGTGGCAGGGGGTATCGCGAGGGAAAAAAAATGATTCCAAA‑GGCGGTACCGGGGGACAAAAAAAAT‑TTGCCCGCTTTAAA‑A‑CCCCTTGG‑GCCAAAA‑GGTTTTGGGCGGGAGGGGG‑GAAAAGTTAAAAAGGAAAACCGCCAGTACATGATGGATTTATTGAAGAAAAATAACTTACTTTAAAgg                                                                                                                                                                                     <  1:476585/251‑1 (MQ=255)
                                             gTTTACCAATCCGACGGACGTATTGCGTCTTTCCGATGCGGGCGTTGATTTT‑AAAT‑CCATCAATGTG‑GGTGGCATGGCGTATCGCGATGGTAAAAAAATGATTACCAA‑TGCGGTATCCGTGGACGACAAAGATATTG‑CCGCCTTCAA‑A‑GCCCTTGGTGCC‑AAA‑GGCATTGAGCTGGACGTGC‑GTAAAGTTTCAAATGACAACCGCCAGTACATGATGGATTTATTGAAGAAAAATAACTTACTTTAAAGGa                                                                                                                                                                                    >  1:77622/1‑251 (MQ=255)
                                               ttACCAATCCGACGGACGTATTGCGTCTTTCCGATGCGGGCGTTGATTTT‑AAAT‑CCATCAATGTG‑GGTGGCATGGCGTATCGCGATGGTAAAAAAATGATTACCAA‑TGCGGTATCCGTGGACGACAAAGATATTG‑CCGCCTTCAA‑A‑GCCCTTGGTGCC‑AAA‑GGCATTGAGCTGGACGTGC‑GTAAAGTTTCAAATGACAACCGCCAGTACATGATGGATTTATTGAAGAAAAATAACTTACTTTAAAGGaaa                                                                                                                                                                                  >  1:216543/1‑251 (MQ=255)
                                                 aCCAATCCGACGGACGTATTGCGTCTTTCCGATGCGGGCGTTGA‑TTT‑AAAT‑CCATCAATGTG‑GGTGGCATGGCGTATCGCGATGGTAAAAAAATGATTACCAA‑TGCGGTATCCGTGGACGACAAAGATATTG‑CCGCCTTCAA‑A‑GCCCTTGGTGCC‑AAA‑GGCATTGAGCTGGACGTGC‑GTAAAGTTTCAAATGACAACCGCCAGTACATGATGGATTTATTGAAGAAAAATAACTTACTTTAAAGGAAATGa                                                                                                                                                                               >  1:214806/1‑251 (MQ=255)
                                                        cgccggacgTATTGCTTCTTCCGGAGGGGGGTGTTGTTTTA‑AAAT‑CCATAATTGTG‑GGTGGCGTGGCGTATCGCGATGAAAAAAAAATGATTACCAA‑TGCGTTACCCGTGGACGAAAAAGATATTC‑CCGCCTTAAA‑A‑GCCTTTGGTCCC‑AAA‑GGCATTGAGCTGGACGTGG‑GAAAAGTTTAAAATGAAAACCGCCAGTACATGATGGATTTATTGAAGAAAAATAACTTACTTTAAAGGAAATGATTATgg                                                                                                                                                                         <  1:442049/248‑1 (MQ=255)
                                                        cgacggcggttttgcttttttCCGATGCGGGCGTTAATTTT‑AAAC‑CCATAAATGGG‑GGTGGCGGGGCGTATCGCGATGGAAAAAAAATGATTACCAA‑TGGGGTACCCGTGGACGAAAAAGATATTG‑CCGCCTTAAA‑A‑GCCCTTGGTGCA‑AAA‑GGCATTGAGCGGGACGTGC‑GTAAAGTTTAAAATGACAACCGCCAGTACATGATGGATTTATTGAAGAAAAATAACTTACTTTAAAGGAAATGATTATgg                                                                                                                                                                         <  1:265588/233‑1 (MQ=255)
                                                             gacgTATTGCGTCTTTCCGATGCGGGCGTTGATTTT‑AAAT‑CCATCAATGTG‑GGTGGCATGGCGTATCGCGATGGTAAAAAAATGATTACCAA‑TGCGGTATCCGTGGACGACAAAGATATTG‑CCGCCTTCAA‑A‑GCCCTTGGTGCC‑AAA‑GGCATTGAGCTGGACGTGC‑GTAAAGTTTCAAATGACAACCGCCAGTACATGATGGATTTATTGAAGAAAAATAACTTACTTTAAAGGAAATGATTATGGAAAtt                                                                                                                                                                    >  1:308380/1‑251 (MQ=255)
                                                               cgTATTGCGTCTTTCCGATGCGGGCGTTGATTTT‑AAAT‑CCATCAATGTG‑GGTGGCATGGCGTATCGCGATGGTAAAAAAATGATTACCAA‑TGCGGTATCCGTGGACGACAAAGATATTG‑CCGCCTTCAA‑A‑GCCCTTGGTGCC‑AAA‑GGCATTGAGCTGGACGTGC‑GTAAAGTTTCAAATGACAACCGCCAGTAc                                                                                                                                                                                                                            >  1:720946/1‑193 (MQ=255)
                                                                 tATTGCGTCTTTCCGATGCGGGCGTTGATTTT‑AAAT‑CCATCAATGTG‑GGTGGCATGGCGTATCGCGATGGTAAAAAAATGATTACCAA‑TGCGGTATCCGTGGACGACAAAGATATTG‑CCGCCTTCAA‑A‑GCCCTTGGTGCC‑AAA‑GGCATTGAGCTGGACGTGC‑GTAAAGTTTCAAATGACAACCGCCAGTACATGATGGATTTATTGAAGAAAAATAACTTACTTTAAAGGAAATGATTATGGAAATTACCa                                                                                                                                                                >  1:861264/1‑251 (MQ=255)
                                                                     gCGTCTTTCCGATGCGGGCGTTGATTTT‑AAAT‑CCATCAATGTG‑GGTGGCATGGCGTATCGCGATGGTAAAAAAATGATTACCAA‑TGCGGTATCCGTGGACGACAAAGATATTG‑CCGCCTTCAA‑A‑GCCCTTGGTGCC‑AAA‑GGCATTGAGCTGGACGTGC‑GTAAAGTTTCAAATGACa                                                                                                                                                                                                                                       >  1:35544/1‑176 (MQ=255)
                                                                                 tgcgggtttttttttt‑AAAT‑CCAT‑AATGGGTGGTGGAATGGCGTATCGGGATGGTAAAAAAAGTATTCCCAT‑TGCGGTTTCCGTAGACGCAAAAGATATTG‑CCGCTTTAAA‑G‑GCCTTTGTTGCC‑AAA‑GGCTTTGAGTTGGATG‑GCTATAATGTTTAAAATGAAAACCGCCAGTACATGATGGATTTATTGAAGAAAAATAACTTACTTTAAAGGAAATGATTATGGAAATTACCACTTTACAAATTATTc                                                                                                                                                 <  1:861954/238‑1 (MQ=255)
                                                                                      gctttgttttt‑AAAT‑CAATAAATGTG‑GGTGGCATGGCGTATCGCGATGGAAAAAAAATGTTTACCAA‑TGCGTTACCCGTGGACAAAAAAGATATTG‑CCCCTTTAAA‑A‑GCCCTTGGTGCC‑AAA‑GGCATTGAGCTGGGGGTGC‑GTAAATTTTAAAATGACAACCGCCAGTACATGATGGATTTATTGAAGAAAAATAACTTACTTTAAAGGAAATGATTATGGAAATTACCACTTTACAAATTATTCTGgtg                                                                                                                                            <  1:506709/243‑1 (MQ=255)
                                                                                          ttaattt‑AAATCCCATCATTGTG‑GGTGGCTGGGCGTATGGGGATGGAAAAAAAATGATTCCAAA‑TGCGGTATCCGGGAACGAAAAAAATATTG‑CCGCCTTAAA‑A‑GCCCTTGGTCCA‑AAA‑GGCTTTGAGTTGGAGGTGT‑AAAAAGTTTAAAATGAAAGCCGCCAGTACATGATGGATTTATTGAAGAAAAATAACTTACTTTAAAGGAAATGATTATGGAAATTACCACTTTACAAATTATTCTGGTGTTTc                                                                                                                                        <  1:779071/247‑1 (MQ=255)
                                                                                            ttttttaaaacccATCATTGTG‑GGTGGCAGGGCGTATCGGGATGGAAAAAAAATGATTACCAA‑GGCGGTATCCGTGGACAAAAAAGATATTG‑CCGCCTTCAA‑A‑CCCTTTGGTCCA‑AAA‑GGCATTGAGCGGGACGTGC‑GAAAATTTTAAAATGACAACCGCCAGTACATGATGGATTTATTGAAGAAAAATAACTTACTTTAAAGGAAATGATTATGGAAATTACCACTTTACAAATTATTCTGGTGTTTCt                                                                                                                                       <  1:592857/240‑1 (MQ=255)
                                                                                                        cATCAATGTG‑GGTGGCATGGCGTATCGCGATGGTAAAAAAATGATTACCAA‑TGCGGTATCCGTGGACGACAAAGATATTG‑CCGCCTTCAA‑A‑GCCCTTGGTGCC‑AAA‑GGCATTGAGCTGGACGTGC‑GTAAAGTTTCAAATGACAACCGCCAGTACATGATGGATTTATTGAAGAAAAATAACTTACTTTAAAgg                                                                                                                                                                                     >  1:607268/1‑193 (MQ=255)
                                                                                                          tCAATGTG‑GGTGGCATGGCGTATCGCGATGGTAAAAAAATGATTACCAA‑TGCGGTATCCGTGGACGACAAAGATATTG‑CCGCCTTCAA‑A‑GCCCTTGGTGCC‑AAA‑GGCATTGAGCTGGACGTGC‑GTAAAGTTTCAAATGACAACCGCCAGTACATGATGGATTTATTGAAGAAAAATAACTTACTTTAAAGGAAATGATTATGGAAATTAcc                                                                                                                                                                 >  1:47196/1‑211 (MQ=255)
                                                                                                             aTGTG‑GGTGGCATGGCGTATCGCGATGGTAAAAAAATGATTACCAA‑TGCGGTATCCGTGGACGACAAAGATATTG‑CCGCCTTCAA‑A‑GCCCTTGGTGCC‑AAA‑GGCATTGAGCTGGACGTGC‑GTAAAGTTTCAAATGACAACCGCCAGTACATGATGGATTTATTGAAGAAAAATAACTTACTTTAAAGGAAATGATTATGGAAATTACCACTTTACAAATTATTCTGGTGTTTCTGGTCGCCTGTATTTgcg                                                                                                                      >  1:662126/1‑251 (MQ=255)
                                                                                                                       gCATGGCGTATCGCGATGGTAAAAAAATGATTACCAA‑TGCGGTATCCGTGGACGACAAAGATATTG‑CCGCCTTCAA‑A‑GCCCTTGGTGCC‑AAA‑GGCATTGAGCTGGACGTGC‑GTAAAGTTTCAAATGACAACCGCCAGTACATGATGGATTTATTGAAGAAAAATAACTTACTTTAAAGGAAATGATTATGGAAATTACCACTTTACAAATTATTCTGGTGTTTCTGGTCGCCTGTATTTGCGGCGCAGGAt                                                                                                             >  1:3587/1‑251 (MQ=255)
                                                                                                                        cagggCGTATCGCGATGGAAAAAAAATGATTCCCAA‑GGCGGTATCCGGGGACGAAAAAAATATTG‑CCGCCTTAAA‑A‑GCCCTGGGTGCA‑AAA‑GGAATTGAGCTGGCCGTGC‑GAAAATTTTAAAATAACAACCGCCAGTACATGATGGATTTATTGAAGAAAAATAACTTACTTTAAAGGAAATGATTATGGAAATTACCACTTTACAAATTATTCTGGTGTTTCTGGTCGCCTGTATTTGCGGCGCAGGAt                                                                                                             <  1:571178/247‑1 (MQ=255)
                                                                                                                        cATGGCTTATCGCTATGGAAAAAAAATGATTACCAA‑TGCGGTATCCGTGGACAAAAAAGATATTC‑CCCCCTAAAA‑C‑CCCCTTGGTCCA‑AAA‑GGCTTGGAGCTGGACGTGC‑GTAAATTTTCAAATGACAACCGCCAGTACATGATGGATTTATTGAAGAAAAATAACTTACTTTAAAGGAAATGATTATGGAAATTACCACTTTACAAATTATTCTGGTGTTTCTGGTCGCCTGTATTTGCGGCGCAGGAt                                                                                                             <  1:328125/250‑1 (MQ=255)
                                                                                                                        aagggcgttcgagtaTGGAAAAAAAATGATT‑CCAATTGCGGTCTCCGTGGACACCAAAGTTATTC‑CCG‑CTTCAACA‑CCCTTTGGTGCC‑AAA‑GGCTTTGTGCGGGACGTGG‑GTAAAGTTCAAAAAGACACCCGCCAGTACATGATGGATTTATTGAAGAAAAATAACTTACTTTAAAGGAAATGATTATGGAAATTACCACTTTACAAATTATTCTGGTGTTTCTGGTCGCCTGTATTTGCGGCGCAGGAt                                                                                                             <  1:228552/236‑1 (MQ=255)
                                                                                                                            gCGTATCGCGATGGTAAAAAAATGATTACCAA‑TGCGGTATCCGTGGACGACAAAGATATTG‑CCGCCTTCAA‑A‑GCCCTTGGTGCC‑AAA‑GGCATTGAGCTGGACGTGC‑GTAAAGTTTCAAATGACAACCGCCAGTACATGATGGATTTATTGAAGAAAAATAACTTACTTTAAAGGAAATGATTATGGAAATTACCACTTTACAAATTATTCTGGTGTTTCTGGTCGCCTGTATTTGCGGCGCAGGATCGAtt                                                                                                        >  1:339529/1‑251 (MQ=255)
                                                                                                                                        ggaaaaaaaaTGATTCCCAA‑TGCGGTACCCGGGGACGAAAAAGATATTC‑CCGCCTTCAA‑A‑CCCCTTGGTGCC‑AAA‑GGCATTGAGCTGGACGTGC‑ATAAATTTAAAAATGACAACCGCCAGTACATGATGGATTTATTGAAGAAAAATAACTTACTTTAAAGGAAATGATTATGGAAATTACCACTTTACAAATTATTCTGGTGTTTCTGGTCGCCTGTAttt                                                                                                                         <  1:719571/219‑1 (MQ=255)
                                                                                                                                        ggTAAAAAAATGATTACCAA‑TGCGGTATCCGTGGACGACAAAGATATTG‑CCGCCTTCAA‑A‑GCCCTTGGTGCC‑AAA‑GGCATTGAGCTGGACGTGC‑GTAAAGTTTCAAATGACAACCGCCAGTACATGATGGATTTATTGAAGAAAAATAACTTACTTTAAAGGAAATGATTATGGAAATTACCACTTTACAAATTATTCTGGTGTTTCTGGTCGCCTGTATTTGCGGCGCAGGATCGATTCTTGACGAATgg                                                                                            >  1:222706/1‑251 (MQ=255)
                                                                                                                                         aaaaaaaaaTGATTACCAA‑GGCGGTATCCGTGGACGAAAAAGATATTG‑CCGCCTAAAA‑A‑GCCCTGGGTGCA‑AAG‑GGCATGGAGCGGGACGTGA‑GAAAAGTTTAAAATGACAACCGCCAGTACATGATGGATTTATTGAAGAAAAATAACTTACTTTAAAGGAAATGATTATGGAAATTACCACTTTACAAATTATTCTGGTGTTTCTGGTCGCCTGTATTTGCGGCGCAGGATCGATTCTTGACGAATGGc                                                                                           <  1:192174/249‑1 (MQ=255)
                                                                                                                                                 aTGATTACCAA‑TGCGGTATCCGTGGACGACAAAGATATTG‑CCGCCTTCAA‑A‑GCCCTTGGTGCC‑AAA‑GGCATTGAGCTGGACGTGC‑GTAAAGTTTCAAATGACAACCGCCAGTACATGATGGATTTATTGAAGAAAAATAACTTACTTTAAAGGAAATGATTATgg                                                                                                                                                                         >  1:313333/1‑165 (MQ=255)
                                                                                                                                                        ccAA‑TGCGGTATCCGTGGACAACAAAGATATTG‑CCGCCTTAAA‑A‑CCCCTTGGTGCA‑AAA‑GGCTTTGAGCTGGACGTGC‑GTAAAGTTTAAAATGACAACCGCCAGTACATGATGGATTTATTGAAGAAAAATAACTTACTTTAAAgg                                                                                                                                                                                     <  1:429530/146‑1 (MQ=255)
                                                                                                                                                              ggggTACCCGTGGACGAAAAAGATATTC‑CCGCCTTCAA‑A‑CCCCTTGGTGCA‑AAG‑GGCATTGAGCTGGACGTGC‑GAAAATTTAAAAATGACAACCGCCAGTACATGATGGATTTATTGAAGAAAAATAACTTACTTTAAAGGAAATGATTATGGAAATTACCACTTTACAAATTATTCTGGTGTTTCTGGTCGCCTGTATTTGCGGCGCAGGATCGATTCTTGACGAATGGCAAACCCAACGCCCGTTAAt                                                                        <  1:666495/248‑1 (MQ=255)
                                                                                                                                                                ggTATCCGTGGACGACAAAGATATTC‑CCGCTTTAAA‑A‑CCCCTTGGTGCA‑AAA‑GGCATTGGGCTGGACGTGC‑GTAAATTTTCAAATGACAACCGCCAGTACATGATGGATTTATTGAAGAAAAATAACTTACTTTAAAGGAAATGATTATGGAAATTACCACTTTACAAATTATTCTGGTGTTTCTGGTCGCCTGTATTTGCGGCGCAGGATCGATTCTTGACGAATGGCAAACCCAACGCCCGTTAATTg                                                                      <  1:240580/250‑1 (MQ=255)
                                                                                                                                                                 gTATCCGTGGACGACAAAGATATTG‑CCGCCTTCAA‑A‑GCCCTTGGTGCC‑AAA‑GGCATTGAGCTGGACGTGC‑GTAAAGTTTCAAATGACAACCGCCAGTACATGATGGATTTATTGAAGAAAAATAACTTACTTTAAAGGAAATGATTATGGAAATTACCACTTTACAAATTATTCTGGTGTTTCTGGTCGCCTGTATTTGCGGCGCAGGATCGATTCTTGACGAATGGCAAACCCAACGCCCGTTAATTGc                                                                     >  1:427022/1‑250 (MQ=255)
                                                                                                                                                                   aTCCGTGGACGACAAAGATATTG‑CCGCCTTCAA‑A‑GCCCTTGGTGCC‑AAA‑GGCATTGAGCTGGACGTGC‑GTAAAGTTTCAAATGACAACCGCCAGTACATGATGGATTTATTGAAGAAAAATAACTTACTTTAAAGGAAATGATTATGGAAATTACCACTTTACAAATTATTCTGGTGTTTCTGGTCGCCTGTATTTGCGGCGCAGGATCGATTCTTGACGAATGGCAAACCCAACGCCCGTTAATTGc                                                                     >  1:253436/1‑248 (MQ=255)
                                                                                                                                                                      cGTGGACGACAAAGATATTG‑CC‑CATTTAA‑ACGCCCTTGGTGCA‑AAA‑GGTGTTGAGCGGGCCGTGC‑GTAAAGTTTAAAATAACAACCGCCAGTACATGATGGATTTATTGAAGAAAAATAACTTACTTTAAAGGAAATGATTATGGAAATTACCACTTTACAAATTATTCTGGTGTTTCTGGTCGCCTGTATTTGCGGCGCAGGATCGATTCTTGACGAATGGCAAACCCAACGCCCGTTAATTGCCTGTAc                                                               <  1:483973/251‑1 (MQ=255)
                                                                                                                                                                            cgaaaaaGATTTTG‑CCGCCTTAAA‑A‑GCCTTTGGTGCA‑AAG‑GGCATTGAGCTGGAGGTGC‑GAAAAGTTAAAAATGACAACCGCCAGTACATGATGGATTTATTGAAGAAAAATAACTTACTTTAAAGGAAATGATTATGGAAATTACCACTTTACAAATTATTCTGGTGTTTCTGGTCGCCTGTATTTGCGGCGCAGGATCGATTCTTGACGAATGGCAAACCCAACGCCCGTTAATTGCCTGTACCTTAAt                                                         <  1:351601/247‑1 (MQ=255)
                                                                                                                                                                                   gATATTG‑CCGCCTTCAA‑A‑GCCCTTGGTGCC‑AAA‑GGCATTGAGCTGGACGTGC‑GTAAAGTTTCAAATGACAACCGCCAGTACATGATGGATTTATTGAAGAAAAATAACTTACTTTAAAGGAAATGATTATGGAAATTACCACTTTACAAATTATTCTGGTGTTTCTGGTCGCCTGTATTTGCGGCGCAGGATCGATTCTTGACGAATGGCAAACCCAACGCCCGTTAATTGCCTGTACCTTAATCGGCTTa                                                  >  1:642285/1‑251 (MQ=255)
                                                                                                                                                                                    atatTG‑CCGCCTTCAA‑A‑GCCCTTGGTGCC‑AAA‑GGCATTGAGCTGGACGTGC‑GTAAAGTTTCAAATGACAACCGCCAGTACATGATGGATTTATTGAAGAAAAATAACTTACTTTAAAGGAAATGATTATGGAAATTACCACTTTACAAATTATTCTGGTGTTTCTGGTCGCCTGTATTTGCGGCGCAGGATCGATTCTTGACGAATGGCAAACCCAACGCCCGTTAATTGCCTGTACCTTAATCGGCTTAg                                                 >  1:563754/1‑251 (MQ=255)
                                                                                                                                                                                       ttG‑CCGCCTACAA‑A‑GCCTTGGGTGCA‑AAA‑GGCTTTGAGCGGGACGTGC‑GAAAAGTTTAAAATGACAACCGCCAGTACATGATGGATTTATTGAAGAAAAATAACTTACTTTAAAGGAAATGATTATGGAAATTACCACTTTACAAATTATTCTGGTGTTTCTGGTCGCCTGTATTTGCGGCGCAGGATCGATTCTTGACGAATGGCAAACCCAACGCCCGTTAATTGCCTGTACCTTAATCGGCTTAgtgt                                              <  1:510497/251‑1 (MQ=255)
                                                                                                                                                                                        ttccccctttCAA‑G‑GCCTTTGGTGCA‑AAA‑GGAATTGGGCTGGACGTGC‑ATAAATTTTAAAATGAAAACCGCCAGTACATGATGGATTTATTGAAGAAAAATAACTTACTTTAAAGGAAATGATTATGGAAATTACCACTTTACAAATTATTCTGGTGTTTCTGGTCGCCTGTATTTGCGGCGCAGGATCGATTCTTGACGAATGGCAAACCCAACGCCCGTTAATTGCCTGTACCTTAATCGGCTTAgtgt                                              <  1:28804/243‑1 (MQ=255)
                                                                                                                                                                                        atcccgcaTTCAA‑A‑CCCCTTGGTGCA‑AAA‑GGCTTGGAGCTGGACGTGG‑GTAAATTTTAAAATGACAACCGCCAGTACATGATGGATTTATTGAAGAAAAATAACTTACTTTAAAGGAAATGATTATGGAAATTACCACTTTACAAATTATTCTGGTGTTTCTGGTCGCCTGTATTTGCGGCGCAGGATCGATTCTTGACGAATGGCAAACCCAACGCCCGTTAATTGCCTGTACCTTAATCGGCTTAgtgt                                              <  1:483406/248‑1 (MQ=255)
                                                                                                                                                                                            agacatttaaaggcctttGGTCCC‑AAA‑GGCATTGAGCTGGACGGGC‑AAAAATTTTAAAATGACAACCGCCAGTACATGATGGATTTATTGAAGAAAAATAACTTACTTTAAAGGAAATGATTATGGAAATTACCACTTTACAAATTATTCTGGTGTTTCTGGTCGCCTGTATTTGCGGCGCAGGATCGATTCTTGACGAATGGCAAACCCAACGCCCGTTAATTGCCTGTACCTTAATCGGCTTAgtgt                                              <  1:642360/233‑1 (MQ=255)
                                                                                                                                                                                              ctttaaa‑a‑GCCCTGGGTGCA‑AAA‑GGCATTGAGCTGGACGTGT‑ATAAAGTTTAAAATGACAACCGCCAGTACATGATGGATTTATTGAAGAAAAATAACTTACTTTAAAGGAAATGATTATGGAAATTACCACTTTACAAATTATTCTGGTGTTTCTGGTCGCCTGTATTTGCGGCGCAGGATCGATTCTTGACGAATGGCAAACCCAACGCCCGTTAATTGCCTGTACCTTAATCGGCTTAGTGTTGGGCg                                        <  1:751973/246‑1 (MQ=255)
                                                                                                                                                                                              ccTTCAA‑A‑GCCCTTGGTGCC‑AAA‑GGCATTGAGCTGGACGTGC‑GTAAAGTTTCAAATGACAACCGCCAGTACATGATGGATTTATTgaag                                                                                                                                                                                                          >  1:132284/1‑89 (MQ=255)
                                                                                                                                                                                                ttCAA‑A‑GCCCTTGGTGCC‑AAA‑GGCATTGAGCTGGACGTGC‑GTAAAGTTTCAAATGACAACCGCCAGTACATGATGGATTTATTGAAGAAAAATAACTTACTTTAAAGGAAATGATTATGGAAATTACCACTTTACAAATTATTCTGGTGTTTCTGGTCGCCTGTATTTGCGGCGCAGGATCGATTCTTGACGAATGGCAAACCCAACGCCCGTTAATTGCCTGTACCt                                                             >  1:832973/1‑228 (MQ=255)
                                                                                                                                                                                                    a‑a‑GCCCTTGGTGCA‑AAG‑GGCTTTGAGCGGGACGTGC‑AAAAATTTTAAAATGACAACCGCCAGTACATGATGGATTTATTGAAGAAAAATAACTTACTTTAAAGGAAATGATTATGGAAATTACCACTTTACAAATTATTCTGGTGTTTCTGGTCGCCTGTATTTGCGGCGCAGGATCGATTCTTGACGAATGGCAAACCCAACGCCCGTTAATTGCCTGTACCTTAATCGGCTTAGTGTTGGGCGATATTc                                  <  1:75751/251‑1 (MQ=255)
                                                                                                                                                                                                      aaccccTTGGT‑CC‑AAAGGGCTTGGAGCTGGCCGTGC‑AAAAATTTTAAAATGACAACCGCCAGTACATGATGGATTTATTGAAGAAAAATAACTTACTTTAAAGGAAATGATTATGGAAATTACCACTTTACAAATTATTCTGGTGTTTCTGGTCGCCTGTATTTGCGGCGCAGGATCGATTCTTGACGAATGGCAAACCCAACGCCCGTTAATTGCCTGTACCTTAATCGGCTTAGTGTTGGGCGATATTc                                  <  1:519172/248‑1 (MQ=255)
                                                                                                                                                                                                          ccTTGGTGCC‑AAA‑GGCTTTGAGCTGGACGTGC‑GTAAATTTTAAAATGACAACCGCCAGTACATGATGGATTTATTGAAGAAAAATAACTTACTTTAAAGGAAATGATTATGGAAATTACCACTTTACAAATTATTCTGGTGTTTCTGGTCGCCTGTATTTGCGGCGCAGGATCGATTCTTGAc                                                                                                  <  1:657456/183‑1 (MQ=255)
                                                                                                                                                                                                           cTTGGTGCA‑AAA‑GGCATTGAGCTGGACGTGC‑GAAAAGTTAAAAATGACAACCGCCAGTACATGATGGATTTATTGAAGAAAAATAACTTACTTTAAAGGAAATGATTATGGAAATTACCACTTTACAAATTATTCTGGTGTTTCTGGTCGCCTGTATTTGCGGCGCAGGATCGATTCTTGACGAATGGCAAACCCAACGCCCGTTAATTGCCTGTACCTTAATCGGCTTAGTGTTGGGCGATATTCAAAcc                             <  1:696229/251‑1 (MQ=255)
                                                                                                                                                                                                                  gcaaaa‑GGCTTTGAGCTGGACGTGC‑ATAAATTTTAAAATGACAACCGCCAGTACATGATGGATTTATTGAAGAAAAATAACTTACTTTAAAGGAAATGATTATGGAAATTACCACTTTACAAATTATTCTGGTGTTTCTGGTCGCCTGTATTTGCGGCGCAGGATCGATTCTTGACGAATGGCAAACCCAACGCCCGTTAATTGCCTGTACCTTAATCGGCTTAGTGTTGGGCGATATTCAAACCGGTAt                        <  1:844237/247‑1 (MQ=255)
                                                                                                                                                                                                                      aa‑GGCATGGAGCTGGCCGTGC‑GTAAATTTAAAAATGAAAACCGCCAGTACATGATGGATTTATTGAAGAAAAATAACTTACTTTAAAGGAAATGATTATGGAAATTACCACTTTACAAATTATTCTGGTGTTTCTGGTCGCCTGTATTTGCGGCGCAGGATCGATTCTTGACGAATGGCAAACCCAACGCCCGTTAATTGCCTGTACCTTAATCGGCTTAGTGTTGGGCGATATTCAAACCGGTATCATTa                   <  1:364023/251‑1 (MQ=255)
                                                                                                                                                                                                                               gAGCTGGACGTGC‑GTAAAGTTTCAAATGACAACCGCCAGTACATGATGGATTTATTGAAGAAAAATAACTTACTTTAAAGGAAATGATTATGGAAATTACCACTTTACAAATTATTCTGGTGTTTCTGGTCGCCTGTATTTGCGGCGCAGGATCGATTCTTGACGAATGGCAAACCCAACGCCCGTTAATTg                                                                      >  1:4830/1‑192 (MQ=255)
                                                                                                                                                                                                                               gAGCTGGACGTGC‑GTAAAGTTTCAAATGACAACCGCCAGTACATGATGGATTTATTGAAGAAAAATAACTTACTTTAAAGGAAATGATTATGGAAATTACCACTTTACAAATTATTCTGGTGTTTCTGGTCGCCTGTATTTGCGGCGCAGGATCGATTCTTGACGAATGGCAAACCCAACGCCCGTTAATTGCCTGTACCTTAATCGGCTTAGTGTTGGGCGATATTCAAACCGGTATCATTATCGGCggg           >  1:408363/1‑251 (MQ=255)
                                                                                                                                                                                                                                aGCTGGACGTGC‑GTAAAGTTTCAAATGACAACCGCCAGTACATGATGGATTTATTGAAGAAAAATAACTTACTTTAAAGGAAATGATTATGGAAATTACCACTTTACAAATTATTCTGGTGTTTCTGGTCGCCTGTATTTGCGGCGCAGGATCGATTCTTGACGAATGGCAAACCCAACGCCCGTTAATTGCCTGTACCTTAATc                                                        >  1:74648/1‑205 (MQ=255)
                                                                                                                                                                                                                                aGCTGGACGTGC‑GTAAAGTTTCAAATGACAACCGCCAGTACATGATGGATTTATTGAAGAAAAATAACTTACTTTAAAGGAAATGATTATGGAAATTACCACTTTACAAATTATTCTGGTGTTTCTGGTCGCCTGTATTTGCGGCGCAGGATCGATTCTTGACGAATGGCAAACCCAACGCCCGTTAATTGCCTGTACCTTAATCGGCTTAGTGTTGGGCGATATTCAAACCGGTATCATTATCGGCgggg          >  1:311837/1‑250 (MQ=255)
                                                                                                                                                                                                                                  cgggACGTGC‑GTAAAGTTTAAAATAACAACCGCCAGTACATGATGGATTTATTGAAGAAAAATAACTTACTTTAAAGGAAATGATTATGGAAATTACCACTTTACAAATTATTCTGGTGTTTCTGGTCGCCTGTATTTGCGGCGCAGGATCGATTCTTGACGAATGGCAAACCCAACGCCCGTTAATTGCCTGTACCTTAATCGGCTTAGTGTTGGGCGATATTCAAACCGGTATCATTATCGGCGGGAcc        <  1:39874/249‑1 (MQ=255)
                                                                                                                                                                                                                                    ggACGTGC‑GTAAAGTTTCAAATGACAACCGCCAGTACATGATGGATTTATTGAAGAAAAATAACTTACTTTAAAGGAAATGATTATGGAAATTACCACTTTACAAATTATTCTGGTGTTTCTGGTCGCCTGTATTTGCGGCGCAGGATCGATTCTTGACGAATGGCAAACCCAACGCCCGTTAATTGCCTGTACCTTAATCGGCTTAGTGTTGGGCGATATTCAAACCGGTATCATTATCGGCGGGACCtt      >  1:772292/1‑251 (MQ=255)
                                                                                                                                                                                                                                         tGC‑GTAAATTTTAAAATGACAACCGCCAGTACATGATGGATTTATTGAAGAAAAATAACTTACTTTAAAGGAAATGATTATGGAAATTACCACTTTACAAATTATTCTGGTGTTTCTGGTCGCCTGTATTTGCGGCGCAGGATCGATTCTTGACGAATGGCAAACCCAACGCCCGTTAATTGCCTGTACCTTAATc                                                        <  1:473761/196‑1 (MQ=255)
                                                                                                                                                                                                                                          tactaaaaGTTTCAAAGCCCAACCGCCAGTACATGATGGATTTATTGAAGAAAAATAACTTACTTTAAAGGAAATGATTATGGAAATTACCACTTTACAAATTATTCTGGTGTTTCTGGTCGCCTGTATTTGCGGCGCAGGATCGATTCTTGACGAATGGCAAACCCAACGCCCGTTAATTGCCTGTACCTTAATCGGCTTAGTGTTGGGCGATATTCAAACCGGTATCATTATCGGCGGGACCTTGGaa  <  1:332513/245‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                             
GGTTTACAACAACCCTGAATATGCCGGCGAACGCATGATGTTGCTGTTTACCAATCCGACGGACGTATTGCGTCTTTCCGATGCGGGCGTTGATTTT‑AAAT‑CCATCAATGTG‑GGTGGCATGGCGTATCGCGATGGTAAAAAAATGATTACCAA‑TGCGGTATCCGTGGACGACAAAGATATTG‑CCGCCTTCAA‑A‑GCCCTTGGTGCC‑AAA‑GGCATTGAGCTGGACGTGC‑GTAAAGTTTCAAATGACAACCGCCAGTACATGATGGATTTATTGAAGAAAAATAACTTACTTTAAAGGAAATGATTATGGAAATTACCACTTTACAAATTATTCTGGTGTTTCTGGTCGCCTGTATTTGCGGCGCAGGATCGATTCTTGACGAATGGCAAACCCAACGCCCGTTAATTGCCTGTACCTTAATCGGCTTAGTGTTGGGCGATATTCAAACCGGTATCATTATCGGCGGGACCTTGGAA  >  1/487488‑487961

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑