Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 803,662 C→A L54I (CTT→ATT)  ybhJ → putative hydratase YbhJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913803,6620CA97.1% 108.1 / ‑6.3 35L54I (CTT→ATT) ybhJputative hydratase YbhJ
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base A (20/14);  minor base T (1/0);  total (21/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCTGCAACTCGAATTATTTAGAGTATATCCATTTATTATCTTTCTGCGCACTTCACGGTGCGCAGATATCTGGAGCATTTGATGATCAAGTTATCTGAAAAAGGCGTGTTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAACTTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATAT  >  NC_000913/803421‑803903
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                 
tCCTGCAACTCGAATTATTTAGAGTATATCCATTTATTATCTTTCTGCTCACTTCACGGTGCGCAGATATCTGGAGCAATTGATGATCAAGTTATCTGAAAAAGGCGTGTTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATaaaattaaaatt                                                                                                                                                                                                                                           <  2:55842/250‑1 (MQ=255)
tCCTGCAACTCGAATTATTTAGAGTATATCCATTTATTATCTTTCTGCGCACTTCACGGTGCGCAGATATCTGGAGCATTTGATGATCAAGTTATCTGTCAACGGCGTGTTTCTCGCCAGTGATAACGAAATAATTGCCGAAGAACATTTCACCGGCGATATTAAACAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTATATTCTCTCTTCGCATAATACGTCTGGCAATATGGATAACATTAAAAtt                                                                                                                                                                                                                                           >  2:497834/1‑250 (MQ=255)
                     gagTATATCCATTTATTATCTTTCTGCGCACTTCACGGTGCGCAGATATCTGGAGCATTTGATGATCAAGTTATCTGAAAAAGGCGTGTTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCtctc                                                                                                                                                                                                                     <  2:160657/251‑1 (MQ=255)
                       gTATATCCATTTATTATCTTTCTGCGCACTTCACGGTGCGCAGATATCTGGAGCATTTGATGATCAAGTTATCTGAAAAAGGCGTGTTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTGTGATTCATTAGCCTCTCAc                                                                                                                                                                                                                   <  2:66007/251‑1 (MQ=255)
                         atatCCATTTATTATCTTTCTGCGCACTTCACGGTGCGCAGATATCTGGAGCATTTGATGATCAAGTTATCTGAAAAAGGCGTGTTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACg                                                                                                                                                                                                                  <  1:28987/250‑1 (MQ=255)
                                                       aCGGTACGCAGATATCTGGAGCATTTGACCATCACGTTATCCGAAAAAGCCGAGTTTCACGCCAGAAATCACGAAACACTTGCCGCAGCATATTTCACTGCCGAAATTAAAAAAGACGCAGCCAAAAAAAGCACTATTGCATGGTCTATTCTCCCTTCGCCTCATAAGTCCGGAAATATGGATAAATTTAAAATTAAGTTTGATTCATTAGACTCTCACGATATTACCTTTGTTGGTATTGTCCAGACCGc                                                                                                                                                                                   >  2:531108/1‑251 (MQ=255)
                                                            gcgcAGATATCTGGAGCATTTGATGATCAAGTTATCTGAAAAAGGCGTGTTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGCAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTaa                                                                                                                                                                                >  1:305034/1‑249 (MQ=255)
                                                                      cTGGAGCATTTGATGATCAAGTTATCTGAAAAAGGCGTGTTTCTCGCCAGTAATAACGAACTAATTGCCGCAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTCTTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTATAGACCGCTAAAGCGTCCGGTAt                                                                                                                                                                    >  2:366287/1‑251 (MQ=255)
                                                                                tgatgaTCAAGTTATCTGAAAAAGGCGTGTTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGttt                                                                                                                                                           <  2:454766/250‑1 (MQ=255)
                                                                                 gatgatCAAGTTATCTGAAAAAGGCGTGTTTCTCTCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCGCTATTGCCTGGTCTATTCTCTCTTCTCATAATACGGCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCGCACGATATTCCCTTTGTTGGTAGTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTc                                                                                                                                                          <  1:341521/250‑1 (MQ=255)
                                                                                     atCAAGTTATCTGAAAAAGGCGTGTTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCTCGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCt                                                                                                                                                     <  1:331878/251‑1 (MQ=255)
                                                                                           ttATCTGAAAAAGGCGTGTTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCAATATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAACTTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGt                                                                                                                                                >  1:526106/1‑250 (MQ=255)
                                                                                                  aaaaaGGAGTGTTTCTAGCCAGTAATAACGAAATAATTGCAGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTATTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTg                                                                                                                                        >  2:483426/1‑251 (MQ=255)
                                                                                                           tgtTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGc                                                                                                                               >  2:12067/1‑251 (MQ=255)
                                                                                                           tgtTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGAACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGc                                                                                                                               >  2:221051/1‑251 (MQ=255)
                                                                                                            gtTTCGCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTTTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGc                                                                                                                               <  2:282433/250‑1 (MQ=255)
                                                                                                                            aaCGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTgcg                                                                                                               <  1:12067/250‑1 (MQ=255)
                                                                                                                              cGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAACATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCg                                                                                                            >  1:236387/1‑251 (MQ=255)
                                                                                                                                    aatTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTcggc                                                                                                       <  1:12743/250‑1 (MQ=255)
                                                                                                                                            aagaaCATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTAtt                                                                                              >  2:197099/1‑251 (MQ=255)
                                                                                                                                               aaCATTTAACCGGCGAATTTAAAAAAGAAGAAACCAAAAAAGGCACTATTGCCTGGTCTATTCTCGCTTCGCATAATACGTCCGGAAATATGAAGAAAAGTAAAATTAAGGTTGATTCATTAGCCGCTCACGATATTACCGTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTACCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAAc                                                                                           <  2:458172/251‑1 (MQ=255)
                                                                                                                                                      cACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGtgat                                                                                     <  1:172900/250‑1 (MQ=255)
                                                                                                                                                         cGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGAcc                                                                                 >  1:560333/1‑251 (MQ=255)
                                                                                                                                                                         gaagaagCCCAAACACGCAATATTGCCTGGTCTCTTATCTCTTCGAATAATACGTCCGGAAATATGGATCAAATTAAAATTAAGTTTGATTCATTAGCCTCTCCCGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAATGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGTCGGCACTATTAACGGTGATGACAATGTTTTTGGTTTATc                                                                 >  2:291422/1‑251 (MQ=255)
                                                                                                                                                                          aagaagCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTATTAACATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATcg                                                                >  1:511885/1‑251 (MQ=255)
                                                                                                                                                                                  aaaaaaGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGcccc                                                        >  1:44895/1‑250 (MQ=255)
                                                                                                                                                                                  aaaaaaGGCAATATTGCATGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCATCTCACGATCTTAACTTTGTTGGTATTGTACAGACAGCTAAAGCGCCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCCCTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGcccc                                                        >  2:382024/1‑250 (MQ=255)
                                                                                                                                                                                   aaaaaGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGCAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCAt                                                                                                                            >  1:97085/1‑182 (MQ=255)
                                                                                                                                                                                   aaaaaGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCGGACCAACTGCCAt                                                                                                                            <  2:97085/182‑1 (MQ=255)
                                                                                                                                                                                    aaaaGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGc                                                      >  2:272507/1‑251 (MQ=255)
                                                                                                                                                                                    aaaaGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGc                                                      >  1:219759/1‑251 (MQ=255)
                                                                                                                                                                                              aTTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCAGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCgg                                            >  1:366187/1‑251 (MQ=255)
                                                                                                                                                                                              aTTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTGGATTCATTAGCCTCTCACGATATTGCCTTTGTTGGGATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACACTGCGCCGTCGGCGGCACTATTAACGTTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGc                                              <  1:382024/249‑1 (MQ=255)
                                                                                                                                                                                                             ctctctTCGCATAATACGTCCGGAAATATGGATAAAATTAACTTTAAGTTTGATTCATTAGCCTCCCACGATATTAACTTTGATGGTATTGTAAAGACCGCTAAAGCGTCCGGCATGGAACGTTTCCAGCTGCCGTATGTGCTGACCAACTGCCCTAACTCACTCTGCGCCGTCGGCGGCACTATTACCGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGGTATgggcgggtctttttttcttc                             >  2:228539/1‑233 (MQ=255)
                                                                                                                                                                                                                  ttCGCATAATACGCCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCCGCGATTTGGCGGGAttttttttgcttctcctt                        >  2:125784/1‑238 (MQ=255)
                                                                                                                                                                                                                                        atGGATAAAATGAAAATTAAGTTTGATTCATTAGCCTCTCACTATATGACCTTTGTGGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCGGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGGCGGCGGCACTATTAACGGTGATGACCATGGTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAatatat  <  2:305034/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                 
TCCTGCAACTCGAATTATTTAGAGTATATCCATTTATTATCTTTCTGCGCACTTCACGGTGCGCAGATATCTGGAGCATTTGATGATCAAGTTATCTGAAAAAGGCGTGTTTCTCGCCAGTAATAACGAAATAATTGCCGAAGAACATTTCACCGGCGAAATTAAAAAAGAAGAAGCCAAAAAAGGCACTATTGCCTGGTCTATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAACTTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATAT  >  NC_000913/803421‑803903

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