Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 76,576 G→A L242L (CTG→TTG)  sgrR ← DNA‑binding transcriptional dual regulator SgrR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091376,5760GA100.0% 119.4 / NA 35L242L (CTG→TTG) sgrRDNA‑binding transcriptional dual regulator SgrR
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (16/19);  total (16/19)

CCCTGTGGACCTTTTAGCATCAGCCCTCCGGCTGGCTCGTCGGCAATTTCCGGCAGAACCCAGACGTTAACTTCGTCGATTAATGCCCGGTAACCGAAGAAGTCATCGAATGCCTGAATTTTCAGTTGATTGGTGCTGTTGCGAATCACCGCATACGGACCGGTGCCGATGGGATGGCTGGCAAAGTTACTGAGGGTTTCCCATTCGCGCGGCAGGATCATCGCCGGAACTTGCCCCAGCAGTAACGGTAACCAGCGGTCCGGTTGCGTGAGATGGATATCCAGCGTCCAGGGCGTCGGCGACACAATGTCAGCAATATGCGAATAGAGCGGCAGCGTATTGATTCGTTTTAAAGAGGCGATCACATCGTCCATTTCCAGTTCACGACCATGGTGAAAATGGACTCCTGGACGCAAAAAGAAACGCCAGTGAAGCGGTGAAATTTGCTGCCAGTGGTGGGCGATGTCTGCTTCCAGTTCCCCAT  >  NC_000913/76335‑76818
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                  
cccTGTGGACCTTTTAGCATCAGACCTCCAGCTGGCTCGTCGGCAATTTCCGGCAGAACCCAGACGTTAACTTCGTCGATTAATGCCCGGTAACCGAAGAAGTCATCGAATGCCTGAATTTTCAGTTGATTGGTGCTGTTGCGGATCACCGCATACGGACCGGTGCCGATGGGATGGCTGGCAAAGTTACTGAGGGTTTCCCATTCGCGCGGCAGGATCATCGCCGGAACTTGCCCCAGCAATAACGGTa                                                                                                                                                                                                                                            >  1:145114/1‑250 (MQ=255)
                    cagaccTCCAGCTGGCTCGTCGGCAATTGCCGGCAGAACCCAGACGTTAACTTCGTCGATTAATGCCCGGTAACCGAAGAAGTCATCGAATGCCTGAATTTTCAGTTGATTGGTGCTGTTGCGGATCACCGCATACGGACCGGTGCCGATGGGATGGCTGGCAAAGTTACTTAGGGTTTCCCATTCGCGCGGCAGGATCATCGCCGGAACTTGCCCCAGCAATAACGGTAACCAGCGGTCCGGTTGCGTg                                                                                                                                                                                                                        <  1:292674/246‑1 (MQ=255)
                           ccagctggctCGTCGGCAATTTCCGGCAGAACCCAGACGTTAACTTCGTCGATTAATGCCCGGTAACCGAAGAAGTCATCGAATGCCTGAATTTTCAGTTGATTGGTGCTGTTGCGGATCACCGCATACGGACCGGTGCCGATGGGATGGCTGGCAAAGTTACTGAGGGTTTCCCATTCGCGCGGCAGGATCATCGCCGGAACTTGCCCCAGCAATAACGGTAACCAGCGGTCCGGTTGCGTGAGATGGa                                                                                                                                                                                                                 <  2:254326/247‑1 (MQ=255)
                            cagctggctCGTCGGCAATTTCCGGCAGAACCCAGACGTTAACTTCGTCGATTAATGCCCGGTAACCGAAGAAGTCATCGAATGCCTGAATTTTCAGTTGATTGGTGCTGTTGCGGATCACCGCATACGGACCGGTGCCGATGGGATGGCTGGCAAAGTTACTGAGGGTTTCCCATTCGCGCGGCAGGATCATCGCCGGAACTTGCCCCAGCAATAACGGTAACCAGCGGTCCGGTTGCGTGAGATGGat                                                                                                                                                                                                                >  1:179296/3‑250 (MQ=255)
                                 ggctCGTCGGCAATTTCCGGCAGAACCCAGACGTTAACTTCGTCGATTAATGCCCGGTAACCGAAGAAGTCATCGAATGCCTGAATTTTCAGTTGATTGGTGCTGTTGCGGATCACCGCATACGGACCGGTGCCGATGGGATGGCTGGCAAAGTTACTGAGGGTTTCCCATTCGCGCGGCAGGATCATCGCCGGAACTTGCCCCAGCAATAACGGTAACCAGCGGTCCGGTTGCGTGAGATGGATATCCa                                                                                                                                                                                                           <  2:287855/250‑1 (MQ=255)
                                     cgtcgGCAATTTCCGGCAGAACCCAGACGTTAACTTCGTCGATTAATGCCCGGTAACCGAAGAAGTCATCGAATGCCTGAATTTTCAGTTGATTGGTGCTGTTGCGGATCACCGCATACGGACCGGTGCCGATGGGATGGCTGGCAAAGTTACTGAGGGTTTCCCATTCGCGCGGCAGGATCATCGCCGGAACTTGCCCCAGCAATAACGGTAACCAGCGGTCCGGTTGCGTGAGATGGATATCCAGCg                                                                                                                                                                                                        <  2:20096/249‑1 (MQ=255)
                                         ggCAATTTCCGGCAGAACCCAGACGTTAACTTCGTCGATTAATGCCCGGTAACCGAAGAAGTCATCGAATGCCTGAATTTTCAGTTGATTGGTGCTGTTGCGGATCACCACATACGGACCGGTGCCGATGGGATGGCTGGCAAAGTTACTGAGGGTTTCCCATTCGCGCGGCAGGATCATCGCCGGAACTTGCCCCAGCAATAACGGTAACCAGCGGTCCGGTTGCGTGAGATGGATATCCAGCGTCCAgg                                                                                                                                                                                                  >  2:74038/1‑251 (MQ=255)
                                         agCAATTTTCGTCAGAACCCAGACGTTAAGTTCGTCGATTAATGCCCGGTAAACGAAGACGTCATACAATGCCTGAATTTTCAGTTGATTGGTGCTGTTGCGGATCACCGCATACGGACCGTTGCCGATGGGATGGCTGGCAAAGTTACTGAGGGTTTCCCATTCGCGCGGCAGGATCATCGCCGGATCTTGCCCCAGCAATCACGGTCACCAGCGGTCCGGTTGCGTGGGATGGATATCCAGCGTCCAgg                                                                                                                                                                                                  >  2:179784/2‑251 (MQ=255)
                                                         aCCCAGACGTTAACTGCGTCGATTAATGCCCGGTAACCGAAGAAGTCATCGAATGCCTGAATTTTCATTTGATTGGTGCTGTTGCGGATCACCGCATACGGACCGGTGCCGATGGGATGGCTGGCAAAGTTACTGAGTGTTTCCCATTCGCGCGGCAGGATCATCGCCGGAACTTGCCCCAGCAATAACGGTAACCAGCGGTCCGGTTGCGTGAGATGGATATCCAGCGTCCAGGGCGTCGGCGACACAAt                                                                                                                                                                                  <  1:5874/251‑1 (MQ=255)
                                                                cGTTAACTTCGTCGATTAATGCCCGGTAACCGAAGAAGTCATCGAATGCCTGAATTTTCAGTTGATTGGTGCTGTTGCGGATCACCGCATACGGACCGGTGCCGATGGGATGGCTGGCAAAGTTACTGAGGGTTTCCCATTCGCGCGGCAGGATCATCGCCGGAACTTGCCCCAGCAATAACGGTAACCAGCGGTCCGGTTGCGTGAGATGGATATCCAGCGTCCAGGGCGTCGGCGACACAATGTCAGc                                                                                                                                                                            >  2:314896/1‑250 (MQ=255)
                                                                           tcgATTAATGCCCGGTAACCGAAGAAGTCATCGAATGCCTGAATTTTCAGTTGATTGGTGCTGTTGCGGATCACCGCATACGGACCGGTGCCGATGGGATGGCTGGCAAAGTTACTGAGGGTTTCCCATTCGCGCGGCAGGATCATCGCCGGAACTTGCCCCAGCAATAACGGTAACCAGCGGTCCGGTTGCGTGAGATGGATATCCAGCGTCCAGGGCGTCGGCGACACAATGTCAGCAATATGCGAAt                                                                                                                                                                 >  2:91833/1‑250 (MQ=255)
                                                                                    gCCCGGTAACCGAAGAAGTCATCGAATGCCTGAATTTTCAGTTGATTGGTGCTGTTGCGGATCACCGCATACGGACCGGTGCCGATGGGATGGCTGGCAAAGTTACTGAGGGTTTCCCATTCGCGCGGCAGGATCATCGCCGGAACTTGCCCCAGCAATAACGGTAACCAGCGGTCCGGTTTCGTGAGATGGATATCCAGCGTCCAGGGCGTCGGCGACACAATGTCTGCAATATGCGAATAGAGCGGCa                                                                                                                                                        <  2:145114/250‑1 (MQ=255)
                                                                                     cccGGTAACCGAAGACGACATCGAATGCCTGAATTTTCAGTTGCTTGGTGCTGTTGCGGATCACCGCATACGGACCAGTGCCGATGGGATGGCTGGCAAAGTTACTGAGGGTTTCCCATTCGCGCGGCAGGATCATCGCCGGAACTTGCCCCCGCAATAACGGTAACCAGCGGGCCGGTTGCGTGAGATGGATATCCAGCGTCCAGGGCGTCGGCGACACAATGTCAGCAATATGCGAATAGCGCGGctgc                                                                                                                                                      >  2:159550/1‑248 (MQ=255)
                                                                                       cGGTAACCGAAGAAGTCATCGAATGCCTGAATTTTCAGTTGATTGGTGCTGTTGCGGATCACCGCATACGGACCGGTGCCGATGGGATGGCTGGCAAAGTTACTGAGGGTTTCCCATTCGCGCGGCAGGATCATCGCCGGAACTTGCCCCAGCAATAACGGTAACCAGCGGTCCGGTTGCGTGAGATGGATATCCAGCGTCCAGGGCGTCGGCGACACAATGTCAGCAATATGCGAATAGAGCGGCAGCGt                                                                                                                                                    >  1:239733/1‑251 (MQ=255)
                                                                                            aCCGAAGAAGTCATCGAATGCCTGAATTTTCAGTTGATTGGTGCTGTTGCGGCTCACCGCATAAGGAACGGTGCCGATGGGATAGCTGGCACAGTTACTGCGGGTTTCCCATTCGCGCGGAAGGATCATCGCCGGAACTTGCCCCAGCAATAACGGTAACCAGCGGTCCGGTTGCGTGAGATGGATATCCAGCGTCCAGGGCGTCGGCGACACAATGTCAGCAATATGCGACTAGAGCGGCAGCGTATTg                                                                                                                                                >  1:27800/1‑250 (MQ=255)
                                                                                            aCCGAAGAAGTCATCGAAGGCCTGAATTTTCCGTTGATTGGTGATGTTGCGGATCACCGCATACGGACAGGTGCCGATGGGATGGCTGGCAAAGTTACTGAGGGTTTCCCATTCGCGCGGCAGGATCATCGCCGGAACTTGCCCCAGCAATAACGGTAACCATCGGTCCGGTTGCGTGAGATGGAGATCCAGCGTCCAGGGCGTCGGCGACACAATGTCAGCAATATGCGAATAGAGCGGCAGCGTATTGa                                                                                                                                               <  2:29475/251‑1 (MQ=255)
                                                                                                     gTCATCGAATGCCTGAATTTTCAGTTGATTGGTGCTGTTGCGGATCACCGCATACGGACCGGTGCCGATGGGATGGCTGGCAAAGTTACTGAGGGTTTCCCATTCGCGCGGCAGGATCATCGCCGGAACTTGCCCCAGCAATAACGGTAACCAGCGGTCCGGTTGCGTGAGATGGATATCCAGCGTCCAGGGCGTCGGCGACACAATGTCAGCAATATGCGAATAGAGCGGCAGCGTATTGATTCGtttt                                                                                                                                       >  1:100212/1‑250 (MQ=255)
                                                                                                     gTCATCGAATGCCTGAATTTTCAGTTGATTGGTGCTGTTGCGGATCACCGCATACGGACCGGTGCCGATGGGATGGCTGGCAAAGTTACTGAGGGTTTCCCATTCGCGCGGCAGGATCATCGCCGGAACTTGCCCCAGCAATAACGGTAACCAGCGGTCCGGTTGCGTGAGATGGATATCCAGCGTCCAGGGCGTCGGCGACACAATGTCAGCAATATGCGAATAGAGCGGCAGCGTATTGATTCGTTTTa                                                                                                                                      <  2:268216/251‑1 (MQ=255)
                                                                                                          cGAATGCCTGAATTTTCAGTTGATTGGTGCTGTTGCGGATCACCGCATACGGACCGGTGCCGATGGGATGGCTGGCAAAGTTACTGAGGGTTTCCCATTCGCGCGGCAGGATCATCGCCGGAACTTGCCCCAGCAATAACGGTAACCAGCGGTCCGGTTGCGTGAGATGGATATCCAGCGTCCAGGGCGTCGGCGACACAATGTCAGCAATATGCGAATAGAGCGGCAGCGTATTGATTCGTTTTAAaga                                                                                                                                  <  2:19087/250‑1 (MQ=255)
                                                                                                                ccTGAATTTTCAGTTGATTGGTGCTGTTGCGGATCACCGCATACGGACCGGTGCCGATGGGATGGCTGGCAAAGTTACTGAGGGTTTCCCATTCGCGCGGCAGGATCATCGCCGGAACTTGCCCCAGCAATAACGGTAACCAGCGGTCCGGTTGCGTGAGATGGATATCCAGCGTCCAGGGCGTCGGCGACACAATGTCAGCAATATGCGAATAGAGCGGCAGCGTATTGATTCGTTTTAAAGAGGCGATc                                                                                                                           <  1:91833/251‑1 (MQ=255)
                                                                                                                        ttCAGTTGATTGGTGCTGTTGCGGATCACCGCATACGGACCGGTGCCGATGGGATGGCTGGCAAAGTTACTGAGGGTTTCCCATTCGCGCGGCAGGATCATCGCCGGAACTTGCCCCAGCAATAACGGTAACCAGCGGTCCGGTTGCGTGAGATGGATATCCAGCGTCCAGGGCGTCGGCGACACAATGTCAGCAATATGCGAATAGAGCGGCAGCGTATTGATTCGTTTTAAAGAGGCGATCACgtcgtc                                                                                                                   >  1:207964/1‑251 (MQ=255)
                                                                                                                               gATTGGTGCTGTTGCGGATCACCGCATACGGACCGGTGCCGATGGTATGGCTGGCAAAGTTACTGAGGGTTTCCCATGCGCGCGGCAGGATCATCGCCGGAACTTGCCCTAGCAATAACGGTATCCAGCGGTCCGGTTGCGTGAGATGGATAGCCAGCGTCCATTGCGGCGGCTACACAATGTCAGCAATATGCGAGTAGAGCAGCAGCGTAGGGTTTCGTTTTAAAGAGGCgctcgcgttttccatttcc                                                                                                            <  2:293769/251‑19 (MQ=255)
                                                                                                                                aTTGGTGCTGTTGCGGATCACCGCATACGGACCGGTGCCGATGGGATGGCTGGCAAAGTTACTGCGGGTTTCCCATTCGCGCGGCAGGATCATCGCCGGAACTTGCCCCAGCAATAACGGTAACTAGCGGTCCGGTTGCGTGAGATGGATATCCAGCGTCCAGGGCGTCGGCGACACAATGTCAGCAATATGCGAATAGAGCGGCAGCGTATTGATTCGTTTTAAAGAGGCGATCACGTCGTCCATTTcc                                                                                                            >  2:214942/1‑250 (MQ=255)
                                                                                                                                aTTGGTGCTGTTGCGGATCACCGCATACGGACCGGTGCCGATGGGATGGCTGGCAAAGTTACTGAGGGTTTCCCATTCGCGCGGCAGGATCATCGCCGGAACTTGCCCCAGCAATAACGGTAACCAGCGGTCCGGTTGCGTGAGATGGATATCCAGCGTCCAGGGCGTCGGCGACACAATGTCAGCAATATGCGAATAGAGCGGCAGCGTATTGATTCGTTTTAAAGAGGCGATCACGTCGTCCATTTCCa                                                                                                           <  1:213757/251‑1 (MQ=255)
                                                                                                                                                    ccGCATACGGACCGGTGCCGATGGGATGGCTGGCAAAGTTACTGAGGGTTTCCCATTCGCGCGGCAGGATCATCGCCGGAACTTGCCCCAGCAATAACGGTAACCAGCGGTCCGGTTGCGTGAGATGGATATCCAGCGTCCAGGGCGTCGGCGACACAATGTCAGCAATATGCGAATAGAGCGGCAGCGTATTGATTCGTTTTAAAGAGGCGATCACGTCGTCCATTTCCAGTTCACGACCATGATGaaaa                                                                                       <  2:80356/251‑1 (MQ=255)
                                                                                                                                                     cGCATACGGACCGGTGCCGATGGGATGGCTGGCAAAGTTACTGAGGGTTTCCCATTCGCGCGGCAGGATCATCGCCGGAACTTGCCCCAGCAATAACGGTAACCAGCGGTCCGGTTGCGTGAGATGGATATCCAGCGTCCAGGGCGTCGGCGACACAATGTCAGCAATATGCGAATAGAGCGGCAGCGTATTGATTCGTTTTAAAGAGGCGATCACGTCGTCCATTTCCAGTTCACGACCATGATGAAAAt                                                                                      >  1:174224/1‑251 (MQ=255)
                                                                                                                                                          aCGGACCGGTGCCGATGGGATGGCTGGCAAAGTTACTGAGGGTTTCCCATTCGCGCGGCAGGATCATCGCCGGAACTTGCCCCAGCAATAACGGTAACCAGCGGTCCGGTTGCGTGAGATGGATATCCAGCGTCCAGGGCGTCGGCGACACAATGTCAGCAATATGCGAATAGAGCGGCAGCGTATTGATTCGTTTTAAAGAGGCGATCACGTCGTCCATTTCCAGTTCACGACCATGATGAAAATGGACt                                                                                 <  2:127266/251‑1 (MQ=255)
                                                                                                                                                          aCGGACCGGTGCCGATGGGATGGCTGGCAAAGTTACTGAGGGTTTCCCATTCGCGCGGCAGGATCATCGCCGGAACTTGCCCCAGCAATAACGGTAACCAGCGGTCCGGTTGCGTGAGATGGATATCCAGCGTCCAGGGCGTCGGCGACACAATGTCAGCAATATGCGAATAGAGCGGCAGCGTATTGATTCGTTTTAAAGAGGCGATCACGTCGTCCATTTCCAGTTCACGACCATGATGAAAATGGACt                                                                                 >  1:87305/1‑251 (MQ=255)
                                                                                                                                                                               ggctggcAAAGTTACTGAGGGTTTCCCATTCGCGCGGCAGGATCATCGCCGGAACTTGCCCCAGCAATAACGGTAACCAGCGGTCCGGTTGCGTGAGATGGATATCCAGCGTCCAGGGCGTCGGCGACACAATGTCAGCAATATGCGAATAGAGCGGCAGCGTATTGATTCGTTTTAAAGAGGCGATCACGTCGTCCATTTCCAGTTCACTACCATGATGAAAATGGACTCCTGGACTCAAAAAGAAACg                                                             <  2:192818/250‑1 (MQ=255)
                                                                                                                                                                                  gggcAAAGTTACTGAGGGTTTTCCATTCGCGCGGCAGGATCATCGCCGGAACTTGCCCCAGCAATAACGGTAACCAGCGGTCGGGTTGCGTGAGATGGATATCCAGCGTCGAGGGCGTCGGCGACACAATGTCAGCAATATTCTAATAGAGCGGCAGCGTATTGATTCGTTTTAAAGAGGCGATCACGTCGTCCATTTCCAGTTCACGACGGTGATGAAAATGGAGTTCTGGACGCAAAAAGAAACGCCAg                                                         <  2:113146/250‑1 (MQ=255)
                                                                                                                                                                                      aaaattACGGAGGGTTTCCCATTCGCGCGGCACGAGCATCGCCGGAACTTGCCCCTGCAATAACGGTAACCAGCGGTCCGGTTGCGTGAGAAGGATATCCAGCGTCCAGTGCGTCGGCGACACAATGTCAGCAATATGCGAATAGAGCGGCAGCGTATTGATTCGTTTTAAAGAGGCGATCACTTCTTCCATATCCAGTTCACGATCATGATGAAAATGGACTCCTGGACGCAAAAAGAAACGCCAGTGaa                                                     <  2:231027/247‑1 (MQ=255)
                                                                                                                                                                                                 gggTTTCCCATTCGCGCGGCAGGATCATCGCCGGAACTTGCCCCAGCAATAACGGTAACCAGCGGTCCGGTTGCGTGAGATGGATATCCAGCGTCCAGGGCGTCGGCGACACAATGTCAGCAATATGCGAATAGAGCGGCAGCGTATTGATTCGTTTTAAAGAGGCGATCACGTCGTCCATTTCCAGTTCACGACCATGATGAAAATGGACTCCTGGACGCAAAAAGAAACGCCAGTGAAGCGGTGATAtt                                          >  1:48144/1‑251 (MQ=255)
                                                                                                                                                                                                      tCCCATTCGCGCGGCAGGATCATCGCCGGAACTTGCCCCAGCAATAACGGTAACCAGCGGTCCGGTTGCGTGAGATGGATATCCAGCGTCCAGGGCGTCGGCGACACAATGTCAGCAATATGCGAATAGAGCGGCAGCGTATTGATTCGTTTTAAAGAGGCGATCACGTCGTCCATTTCCAGTTCACGACCATGATGAAAATGGACTCCTGGACGCAAAAAGAAACGCCAGTGAAGCGGTGATATTtgct                                      <  1:12568/250‑1 (MQ=255)
                                                                                                                                                                                                      tCCCATTCGCGCGGCAGGATCATCGCCGGAACTTGCCCCAGCAATAACGGTAACCAGCGGTCCGGTTGCGTGAGATGGATATCCAGCGTCCAGGGCGTCGGCGACACAATGTCAGCAATATGCGAATAGAGCGGCAGCGTATTGATTCGTTTTAAAGAGGCGATCACGTCGTCCATTTCCAGTTCACGACCATGATGAAAATGGACTCCTGGACGCAAAAAGAAACGCCAGTGAAGCGGTGATATTtgct                                      <  1:11031/250‑1 (MQ=255)
                                                                                                                                                                                                                                         ccccAGCAATAACGGTAACCAGCGGTCCGGTTGCGTGAGATGGATATCCAGCGTCCAGGGCGTCGGCGACACAATGTCAGCAATATGCGAATAGAGCGGCAGCGTATTGATTCGTTTTAAAGAGGCGATCACGTCGTCCATTTCCAGTTCACGACCATGATGAAAATGGACTCCTGGACGCAAAAAGAAACGCCAGTGAAGCGGTGATATTTGCTGCCAGTGATGGGCGATGTCTGCTTCCAGTTCCCCAt  >  2:107178/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                  
CCCTGTGGACCTTTTAGCATCAGCCCTCCGGCTGGCTCGTCGGCAATTTCCGGCAGAACCCAGACGTTAACTTCGTCGATTAATGCCCGGTAACCGAAGAAGTCATCGAATGCCTGAATTTTCAGTTGATTGGTGCTGTTGCGAATCACCGCATACGGACCGGTGCCGATGGGATGGCTGGCAAAGTTACTGAGGGTTTCCCATTCGCGCGGCAGGATCATCGCCGGAACTTGCCCCAGCAGTAACGGTAACCAGCGGTCCGGTTGCGTGAGATGGATATCCAGCGTCCAGGGCGTCGGCGACACAATGTCAGCAATATGCGAATAGAGCGGCAGCGTATTGATTCGTTTTAAAGAGGCGATCACATCGTCCATTTCCAGTTCACGACCATGGTGAAAATGGACTCCTGGACGCAAAAAGAAACGCCAGTGAAGCGGTGAAATTTGCTGCCAGTGGTGGGCGATGTCTGCTTCCAGTTCCCCAT  >  NC_000913/76335‑76818

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