Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,878,3310TG73.7% 56.3 / 20.8 39I206R (ATA→AGA) yeaPdiguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base T (4/6);  new base G (13/15);  total (18/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.63e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CCAGCAGTGAGAAGCGCCAGTGGAGTGAACGCGCGGAACAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCGAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTAAACAACGAGAATGCGGATATTTCGTCGCTGCGAGAACGCATTCAGCAGCAAATACGTGGAGAA  >  NC_000913/1878097‑1878563
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                           
ccAGCAGTGAGAAGCGCCAGTGGAGTGAACGCGCGGAACAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCGAATCGGCGGGCGCTTTTTGAAAATCTGACGCCACTGTTTTCCCTCTCCCGGCATCTTACCCATAAGATaaat                                                                                                                                                                                                                                         >  2:1115516/1‑237 (MQ=255)
                     ggAGTGAACGCGCGGAACAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCGAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGAGAATGATCGCGTTTAt                                                                                                                                                                                                                              <  1:468645/229‑1 (MQ=255)
                         tGAACGCGCGGAACAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCGAATCGGCGGGCGATTTTTTAAAATCTTACGACACTGTTTTCCCTCGCCCGGCCTCTTAACCCTAAGAGAATGATCGCGTTTTTCGCTCTTGATAACt                                                                                                                                                                                                                >  2:560251/1‑239 (MQ=255)
                                         gTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCGAATCGGCGGGGGGTTTTTTAAAATCTTACGACACTGTTTTTCCTCCCCCCGGCGCCTTACCACAAAAAAATGATCCCGGTTTTTGTTCCTGTTAACTTACAATTAATCAAtaa                                                                                                                                                                                                >  2:20972/1‑237 (MQ=255)
                                              aCAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCGAATCGGCGGGCGGTTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCCTCTTTAACATAAGAGAATGATCCCGGTTATCCATCTGGATAAATTCCAATTTAACAATTATCGtt                                                                                                                                                                                           >  2:143428/1‑239 (MQ=255)
                                                   tATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCGAATCGGCGGGCGATTTTTGGAAATCTGACGACACTGTTTTCCCTCGCCCGGTCTGTTGACCAAGAGGTATTTAACGGGTTTGTCGATATGGAAAACTTCAAATTAATCAATCATCGtt                                                                                                                                                                                           <  2:415333/234‑1 (MQ=255)
                                                   tATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCGAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGtt                                                                                                                                                                                           <  1:166932/234‑1 (MQ=255)
                                                   tATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCGAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGAGAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGtttt                                                                                                                                                                                         <  1:143428/236‑1 (MQ=255)
                                                   tATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCGAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGAGAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGt                                                                                                                                                                                      <  1:560251/239‑1 (MQ=255)
                                                   tATTCGCCGGACTGATTGCAAAATTTATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATGGCGCAATCGTATACCGACTCGTTAACCGGGCTTCCGAAGCGGCAGGCTATTTTTGAAAATCTAACGACACTGTTGTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGAGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGt                                                                                                                                                                                      <  2:1181748/239‑1 (MQ=255)
                                                          cGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCGAATCGGCGGGCGATTTTTGAAAATCTGACGACCCTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATTGATCTGGATAACTTCAAATTTATCAATGATCGTTTTGGTCataata                                                                                                                                                                               >  1:1175822/1‑239 (MQ=255)
                                                                        aaTATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCGAATCGGCTGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGAGAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGttt                                                                                                                                                                 <  1:120380/239‑1 (MQ=255)
                                                                             atTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCGAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCCCCCGGCCTCTTAACCATAAGAGAATGATCGCGTTTTTCGATCTGGATAACTTTAAATTAATCAATGATCGTTTTGGTCATAATATGGGCGATATGttttct                                                                                                                                                              >  2:73330/1‑234 (MQ=255)
                                                                               tCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCGAATCGGCGGGCGATTTTTTAAAAACTGACGAACCTGTTTTTCCTCCCCCGGCATCTTAACCAAAAAAGAATGATCGCGTTTATCGATCTGGGTAACTTCCAAATAAACAATGATCGTTTTTGGTATAATAGAGGCGa                                                                                                                                                                        >  2:1160699/1‑225 (MQ=255)
                                                                                       agGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCGAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCCGCATCTTAACCATAAGAGAATGATCGCGGTTATCGATCTGGGTAAATTCAAATTAATCAATGAACGTTTTGGTCATAATAATGGCGATCTGTTTTTCCTTTCAGGtggc                                                                                                                                                  >  2:1117459/1‑236 (MQ=255)
                                                                                           aCTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCGAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGAGAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAg                                                                                                                                              <  1:139929/238‑1 (MQ=255)
                                                                                             gggTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCGAATCGGCGGGCGATTGTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGAGAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGtt                                                                                                                                                                                           <  1:446639/191‑1 (MQ=255)
                                                                                                ttGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCGAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGAGAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCt                                                                                                                                         <  2:977595/238‑1 (MQ=255)
                                                                                                       gCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCGAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCGTCTTAACCATAAGAGAATGATCGCGTTTTTAGATCTGGGTTAATTCAAATTAATAAATGATTGTTTTGGGTataa                                                                                                                                                                                 >  2:446639/1‑192 (MQ=255)
                                                                                                         tGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCGAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTTCCTCGCCCGGCATCTTAACCATAAGAGAATGATCGCGTTTTTCGATCTGGATAACTTCCAATTAATCAATGATTGTTTTGGTCATAATAGTGGCGATCTGTTTTTCATTTAGGGTGGCGAGGGCCTTAATACGcc                                                                                                                               >  2:1145780/1‑237 (MQ=255)
                                                                                                         tGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCGAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGAGAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTc                                                                                                                              <  2:535878/239‑1 (MQ=255)
                                                                                                                      tgctgcGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCGAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGAGAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCg                                                                                                                 <  1:1117459/239‑1 (MQ=255)
                                                                                                                      tgctgcGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCGAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGAGAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCg                                                                                                                 <  1:1145780/239‑1 (MQ=255)
                                                                                                                                     gcAATCGTATACCGACTCGTTAACCGGGCTACCGAATCGGCGGGCGATTTTTTAAAATCTAACGACACTGTTTTCCCTCGCCCGGCATCTTTACCATAAGAGAATGATCGCGTTTTTCGAACTGGATAACTTCAACATAATCCATTATTGTTTTGGTCATAATAGTGGCGAACTGTTTTTCATTCCGGTTTGCGAGGGCCTTTATACGCTCCAGCAAAAAGGCGGAGTTTTTGgtccg                                                                                                  >  2:131783/1‑236 (MQ=255)
                                                                                                                                         tCGTATACCGACTCGTTAACCGGGCTACCGAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGAGAATGATCGTGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGgt                                                                                                     <  1:571082/231‑1 (MQ=255)
                                                                                                                                              taCCGACTCGTTAACCGGGCTACCGAATCGGCGGGGGATTTTTTAAAAACTGACGACACTGTTTTCCCTCGCCCGGAATCTTAACCAAAAGAGAATGATCGCGTTTATAGATCTGGGTAACCTTAAATTAATTAATGATTGGTTTGGTCATAATAATGGTGAGTTGTTTCTCATTAAGGTTGGGGAGGGGCTTCATACGCTCCAGGAAAATGGATAAGATATTGgtcggttcgcggggg                                                                                        >  2:153503/1‑228 (MQ=255)
                                                                                                                                                   aCTCGTTAACCGGGCTACCGAATCGGCGGGCGGTTTTTTAAAAACTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGAGAATGATCGGGTTTATCGATCTGGATAAATTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTTAGGTTGGCGAGCGCCTTTATACGCCCCAGCAAAATGGCGAAATTTTTGGTCGGTTcggtg                                                                                           >  2:571082/1‑229 (MQ=255)
                                                                                                                                                   aCTCGTTAACCGGGCTACCGAATCGGCGGGCGATTTTTTAAAAACTGACGACACCGTTTTTCCCCCCCCCGCCTCTTTACCATAAAAGAATGATCGCGTTTATCGCTCTTGGTAACTTCAAATTAATCAATGGTCCGTTTTGTCATAATAATGGCGGTTTGTTTTTCATTCCTGGTTGCGGACGCCTTTATAACGCCCCACCAAATGGTGGAGTTTTTTGGTGGTCCCGGGGGtggggg                                                                                    >  2:1152102/1‑235 (MQ=255)
                                                                                                                                                   aCTCGTTAACCGGGCTACCGAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGAGAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCt                                                                                                                                         >  2:562583/1‑187 (MQ=255)
                                                                                                                                                   aCTCGTTAACCGGGCTACCGAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGAGAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCt                                                                                                                                         <  1:562583/187‑1 (MQ=255)
                                                                                                                                                               ggCTACCGAACCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAAt                                                                                                                                                                                                          >  1:942868/1‑111 (MQ=255)
                                                                                                                                                               ggCTACCGAACCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAAt                                                                                                                                                                                                          <  2:942868/111‑1 (MQ=255)
                                                                                                                                                                     cGAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTa                                                                 <  2:557502/239‑1 (MQ=255)
                                                                                                                                                                           ggcgggcgATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTAaacaac                                                           <  2:230605/239‑1 (MQ=255)
                                                                                                                                                                                                                  tCGCCCGGCATCTTAACCATAAGAGAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTAAACAACGAGAATGCGGATATTTCGTCGCTGCGAGAACGCATTcag                    <  1:979650/239‑1 (MQ=255)
                                                                                                                                                                                                                       ccGGCATCTTAACCATAAGAGAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTAAACAACGAGAATGCGGATATTTCGTCGCTg                                   <  1:991417/220‑1 (MQ=255)
                                                                                                                                                                                                                                tAACCATAAGAGAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTAAACAACGAGAATGCGGATATTTCGTCGCTGCGAGAACGCATTCAGCAGCAAATACGTg       <  2:94660/239‑1 (MQ=255)
                                                                                                                                                                                                                                 aaCCATAAGAGAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGGTTTTAGTTGTTTTACTAAAAAACGAGAATGCGGATATTTTGTTGCTGCGAGAAAGCa                         >  2:991417/1‑220 (MQ=255)
                                                                                                                                                                                                                                     aTAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGGTGGGGTTTTTGTTTTTTTACTAAAAAAACAGAAATGGGAGATTTTGGTGCTGGGGGGACGGATTTCACAGCAAAAACGGGGGGaa  >  2:1019300/1‑239 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                           
CCAGCAGTGAGAAGCGCCAGTGGAGTGAACGCGCGGAACAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCGAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCTCATTCAGGTTGGCGAGCGCCTTAATACGCTCCAGCAAAATGGCGAAGTTATTGGTCGTCTCGGCGGTGATGAGTTTTTAGTTGTTTCACTAAACAACGAGAATGCGGATATTTCGTCGCTGCGAGAACGCATTCAGCAGCAAATACGTGGAGAA  >  NC_000913/1878097‑1878563

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