Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,173,363 Δ2 bp intergenic (‑490/+918) gatD ← / ← gatB galactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,173,3630C.85.3% 106.2 / 15.4 34intergenic (‑490/+919)gatD/gatBgalactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base C (0/5);  new base . (17/12);  total (17/17)
Fisher's exact test for biased strand distribution p-value = 4.45e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01
*NC_0009132,173,3640C.85.3% 107.7 / 13.9 34intergenic (‑491/+918)gatD/gatBgalactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base C (0/5);  new base . (17/12);  total (17/17)
Fisher's exact test for biased strand distribution p-value = 4.45e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

CGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCA‑TCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCAT  >  NC_000913/2173131‑2173588
                                                                                                                                                                                                                                          ||                                                                                                                                                                                                                                
cgccAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGCTGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCA‑TCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCagcagc                                                                                                                                                                                                                >  2:190635/1‑250 (MQ=255)
         ggTGTACAAACCAATAGGTTGCGTCGCGACCCACAGGGTGATGCTCATAAGGATGACCCCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCA‑TCGCCACGAAGAAGCCGATGGTGGCAAGATCGCAAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGAt                                                                                                                                                                                                      <  2:12737/251‑1 (MQ=255)
                        gtgttGCGTCTCGATCCACAGGGTGATTCGCATAACGTGGACCCCCGTATTTAAGTTGCGAAACAGATTTCCACGAGGCCCGCCGACCTCCA‑TCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGTAAGCACCTGATTCCCCGGCACACAAACAGCCATTAAAATCGGT‑AGGGTGATAAAAATCAGGCTTGCCGATAC‑‑CGCCGCCGTATGTCCCAGCAGCATCTCCGGATCAAGGCCAATCAGGa                                                                                                                                                                                       <  2:21130/248‑1 (MQ=255)
                            tGCGTCGCGATCCACAGGGGGATGCTCATAATGATGCCAACCGAGCTTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCA‑TCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGAAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGATTGCCGATACCACC‑‑GCCGTATGTcc                                                                                                                                                                                                                       >  2:459364/1‑215 (MQ=255)
                            ggCGTCGCGAACCGCAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTCCCACGAGGCACGTCGCCCGCCA‑TCGCCACGAAGAAGCCGATGGTTGCAATATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTcc                                                                                                                                                                                                                       <  1:459364/214‑1 (MQ=255)
                                  gcgATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCA‑TCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAcc                                                                                                                                                                             >  1:276771/1‑251 (MQ=255)
                                     aTCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCA‑TCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccgcc                                                                                                                                                                          >  1:165926/1‑251 (MQ=255)
                                           ggggTGTTGCTCAAAAGGGTGACACCCGATATTAAGGTGCGTAACATATTTCCGCGATGC‑CGGGGGCCGCCAATCGCCACGAAGAAGTGAATGGTGGCAAGATCGCCAATTGGCAGCACCGGATTCCCCGGCCCAAAAACAGCAATTAAAAT‑GGTGTGTGGAATAAAAACCAGGCTTGCCGAT‑CC‑CTCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAAGCAGGAGCTCCGGACCGCCTAACtt                                                                                                                                                                    <  2:429325/250‑1 (MQ=255)
                                                   gCTCATAAGTATGACACCCGAGATTAAGGGGCGGAACAGATTTCCGCGATGCCCGGCGACCGCCA‑TCACCGCGAAGAAGCCGATGGTGGAAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTa                                                                                                                                                            <  1:332190/251‑1 (MQ=255)
                                                      cATAATGATGCCCCCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCA‑TCGCCACGAAGAAGCCGATGGTGGCAAGAGCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCAGTAAAc                                                                                                                                                         <  1:131077/251‑1 (MQ=255)
                                                                acacCCGAGATTAAGGTGCGGAACAGATTTACACGATGCACGGCGCCCGCCA‑TCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCCATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGc                                                                                                                                               >  1:208115/1‑251 (MQ=255)
                                                                     cGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGGGACCGCCA‑TCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATGACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGt                                                                                                                                                             <  1:26075/232‑1 (MQ=255)
                                                                     cGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCA‑TCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGt                                                                                                                                                             >  2:26075/1‑232 (MQ=255)
                                                                               gTGCGGAACAGATTTCCACGATGCACGGCGAACGCCA‑TCGCCACGAAGAAGCCGATGGTGGCACGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGACCCAGCAGCAACGCCGGATCACGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCAGTAAACAACTACGAGCTTGCTTAACGAtagg                                                                                                                                >  2:265506/1‑248 (MQ=255)
                                                                                    gAACAGATTTCCACGAGGCACGGCGCCCGCCA‑TCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTa                                                                                                                           <  1:405204/251‑1 (MQ=255)
                                                                                              tttccgatcccGGCGCCCGCCA‑TCGCCCCGAAGAAGCCGAGGGTGGCAAGATCGCAAAACGTCAGCACCTGATTCCCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGTTTGCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAAccatcca                                                                                                                  <  1:313598/240‑1 (MQ=255)
                                                                                                        cGGCGACCGCCA‑TCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGAAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAAccatcca                                                                                                                  >  2:313598/1‑240 (MQ=255)
                                                                                                            gACCGCCA‑TCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGCTTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGACCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGACCTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATc                                                                                                    >  2:235095/1‑250 (MQ=255)
                                                                                                            gACCGCCA‑TCGCCACGAAGAAGCAGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGCGCGTTAAACCATCCATGATGGCTTTAATc                                                                                                    >  2:381359/1‑250 (MQ=255)
                                                                                                              ccGCCA‑TCGCCAAGAAACACCCGATGCTGCCCAGACCGCCAAACGTCCGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGACTGGGATAACAATCAGACTTGACGCAACCACA‑‑TCCGTATGTACCAACAGCAACGCCGGATCAAGGCCAATCCCAAACTCCGGACCGCAGAACTTCGCCTGCACACGCCTACGCACCTCCTTAGCTCTGGGTGTTACAccagccagacacgtgttacccacc                                                                                                 >  2:182586/1‑234 (MQ=255)
                                                                                                                  cA‑TCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCa                                                                                                                                                                                              >  1:1528/1‑154 (MQ=255)
                                                                                                                  cA‑TCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCa                                                                                                                                                                                              <  2:1528/154‑1 (MQ=255)
                                                                                                                             aagaagCCGATGGGGGCAAAATCGCCAAACGGCAGCCCCTGATTACCCGGCACACAAACAGCAATAAAAAT‑GGTGAGTGGGATAAAAATCGGGCTGGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTagcag                                                                                   <  2:378106/251‑1 (MQ=255)
                                                                                                                                   gcgaggggggcaaaaccGCCAAACGGCAGCACCGGATTCCCCGGCAAACAAACAGCAATTAAAAG‑GGGGAGTGGGATAAAAATCAGTCGTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGGCCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCAc                                                                             <  2:146498/235‑1 (MQ=255)
                                                                                                                                   ccGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCAc                                                                             >  2:133347/1‑251 (MQ=255)
                                                                                                                                                gATCGCCAACCGGCACCACCTGATTCCCCGGCACACAAACAGCAATAAAAAT‑TGTGAGTGGGATAAAAATCAGGCTTGCCGAT‑CC‑CCCCGCCGTATTTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGttttt                                                                <  1:67052/251‑1 (MQ=255)
                                                                                                                                                         aCGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGACCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGc                                                       >  1:102656/1‑251 (MQ=255)
                                                                                                                                                                 cccTGATTCCCCGGCACACAAACAGCAATTAAAAT‑GGGGAGTGGAAAAAAAATCAGGCTTGCCGATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGAAACTCCGGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATac                                               <  2:383149/250‑1 (MQ=255)
                                                                                                                                                                              gCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCg                                  >  1:108535/1‑251 (MQ=255)
                                                                                                                                                                                        aGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCcgcg                         >  1:52934/1‑250 (MQ=255)
                                                                                                                                                                                          cAATTAATAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACA‑‑GTCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTaa                               >  2:39141/1‑242 (MQ=255)
                                                                                                                                                                                          cAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GTCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTaa                               <  1:39141/242‑1 (MQ=255)
                                                                                                                                                                                            aTTAAAAT‑GGGGAGTGGGTTAAAAATCAGGTTTGCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACCCCGGATCAAGGCCAATCAGGAACGCCTGACCGCCGAACTTCGCCTGTAAACGTCTACGTGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGAt                    <  1:260547/251‑1 (MQ=255)
                                                                                                                                                                                                          gTGGGATAAAAATCAGGCTTGCCGATCCCACC‑‑GCCTTATTTCCCAGCAGCATCGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGTTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTTCCGCCAGCTGCAATACACTTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGc       <  1:454895/251‑1 (MQ=255)
                                                                                                                                                                                                                 ataaatataagacttgccgatccaaccgccgTTTGACGCAGCAGAAACGCCGAATCAAGGCCAATCAGGAAATACTGACCGACGAACTTCGCCTGTAAACGACTACGAGCTTGATTAGCGAGCAGCGTTAAACAATCCATGATGGGTTTACTCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTCCCGCCAGCTGCCATACACCTTTTCCATCGTAACCCGCGCGGATGCCGATAATCAGCCCCAt  >  2:260496/26‑251 (MQ=255)
                                                                                                                                                                                                                                          ||                                                                                                                                                                                                                                
CGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCA‑TCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCAT  >  NC_000913/2173131‑2173588

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