Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,667,460 G→A Q369* (CAG→TAG)  mlc ← DNA‑binding transcriptional repressor Mlc

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,667,4600GA83.3% 81.6 / 7.9 36Q369* (CAG→TAG) mlcDNA‑binding transcriptional repressor Mlc
Reads supporting (aligned to +/- strand):  ref base G (4/2);  new base A (18/12);  total (22/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 = 8.37e-01

TCACAAATATAAACCAGGAAAATAATTAACCTTGAAAGTCTAAGTTATGCTTTCCTGGCCCAAATTGAGATAGCGCAAATTTTGGTAGAACAGTTAAAAAATGTTAACCCTGCAACAGACGAATCAACAAAGAACCGTTATACATCGCGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAACTGAGTACTCTCAACGCTGATGTGCTGACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGATATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGAATGCGCCCGACATGCGCGCCCACCCCGGTAATGATGTCTTTTGCCAGTAGATCGCCGCGCAATGCCGCCTGACACAATGAGTCCACGGTTAACGGTTGTCCATGTAACATCGAGCTC  >  NC_000913/1667241‑1667707
                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                       
tcacaaaaatacacagggtcaataatTAACCTTGAAAGGCTACGATCTGCTTTCCTGGCCCAAATTGAGATAGCGCAAATTTTGGTAGAATAGTTAAAAAATGATAACCCTGCAACAGAAGAATCAATAAAGACCGGTTATACACTGTGTCTTTTCCCCGTGCAACGCCAGCCATCGTGACCTGGTTCGAAAACTGAGTACTCTCAACGCTGATGTGCTAACTATACGCCGTAAGGTCCTGCTGACgaatc                                                                                                                                                                                                                          >  2:308883/21‑247 (MQ=255)
        ataAACCAGGAAAATAATTAACCTTGAAAGTCTAAGTTATGCTTTCCTGGCCCAAATTGAGATAGCGCAAATTTTGGTAGAACAGTTAAAAAATGTTAACCCTGCAACAGACGAATCAACAAAGAACCGTTATACATCGCGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAACTGAGTACTCTCAACGCTGATGTGCTAACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTga                                                                                                                                                                                                                  <  2:218287/251‑1 (MQ=255)
              cAGGAAACTCATTAACCTTGAAAGTCTAAGTTATGCTTTCCTGGCCCAAATTGAGATAGCGCAAATTTTGGTAGAACAGTTAAAAAATGTTAACCCTACAACAGACGAATCAACAAAGAACCGTTATACATCGCGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAACATGAGTACTCTAAACGCTGATGTGCTGACTATACGCAGGAAGGGCCTGATGACGGATGCTGTCTGCGANGAc                                                                                                                                                                                                            >  2:441863/1‑251 (MQ=255)
                  aaaaTAATTAACCTTGAAAGTCTAAGTTATGCTTTCCTGGCCCAAATTGAGATAGCGCAAATTTTGGTAGAACAGTTAAAAAATGTTAACCCTGCAACAGACGAATCAACAAAGAACCGTTATACATCGCGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAACTGAGTACTCTCAACGCTGATGTGCTAACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCggg                                                                                                                                                                                                        >  2:424145/1‑251 (MQ=255)
                                     gTCTAAGTTATGCTTTCCTGGCCCAAATTGAGATAGCGCAAATTTTGGTAGAACAGTTAAAAAATGTTAACCCTGCAACAGACGAATCAACAAAGAACCGTTATACATCGCGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAACTGAGTACTCTCAACGCTGATGTGCTAACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGATATCTGCCGCtt                                                                                                                                                                                     <  2:189989/251‑1 (MQ=255)
                                       cTAAGTTATGCTTTCCTGGCCCAAATTGAGATAGCGCAAATTTTGGTGGAACAGTTAAAAAATGTTAACCCTGCAACAGACGAATCAACAAAGAACCGTTATACATCGCGTCTTTTACAAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAACTGAGTACTCTCAACGCTGATGTGCTAACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGATATCTGACGCttt                                                                                                                                                                                    >  1:43686/1‑250 (MQ=255)
                                            ttATGCTTTCCTGGCCCAAATTGAGATAGCGCAAATTTTGGTAGAACAGTTAAAAAATGTTAACCCTGCAACAGACGAATCAACAAAGAACCGTTATACATCGCGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAACTGAGTACTCTCAACGCTGATGTGCTAACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGATATCTGCCGCTTTACTTa                                                                                                                                                                               >  1:454706/1‑250 (MQ=255)
                                                gCTTTCCTGGCCCAAATTGAGATAGCGCAAATTTTGGTAGAACAGTTAAAAAATGTTAACCCTGCAACAGACGAATCAACAAAGAACCGTTATACATCGCGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAACTGAGTACTCTCAACGCTGATGTGCTAACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGATATCTGCCGCTTTACTTAACGGt                                                                                                                                                                          <  1:383745/251‑1 (MQ=255)
                                                      cTGGCCCAAATTGAGATAGCGCAAATTTTGGTAGAACAGTTAAAAAATGTTAACCCTGCAACAGACGAATCAACAAAGAACCGTTATACATCGCGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAACTGAGTACTCTCAACGCTGATGTGCTAACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGATATCTGCCGCTTTACTTAACGGTGAGcc                                                                                                                                                                     >  2:222706/1‑250 (MQ=255)
                                                           ccAAATTGAGATAGCGCAAATTTTGGTAGAACAGTTAAAAAATGTTAACCCTGCAACAGACGAATCAACAAAGAACCGTTATACATCGCGGCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAACTGAGTACTCTCAACGCTGATGTGCTAACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGATATCTGCCGCTTTACTTACCGGTGAGCCAATCAg                                                                                                                                                               <  2:180009/251‑1 (MQ=255)
                                                                     aTAGCGCAAATTTTGGTAGAACAGTTAAAAAATGTTAACCCTGCAACAGACGAATCAACAAAGAACCGTTATACATCGCGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAACTGAGTACTCTCAACGCTGATGTGCTAACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGATATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTtgt                                                                                                                                                     >  1:270960/1‑251 (MQ=255)
                                                                       aGCGCAAATGTTGGTAGAACAGTTAAAAAATGTTAACCCTGCAACAGACGAATCAACAAAGAACCGTTATACATCGCGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAACTGAGTACTCTCAACGCTGATGTGCTAACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGATATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTgg                                                                                                                                                   <  2:43686/251‑1 (MQ=255)
                                                                               ttttGGTAGAACAGTTAAAAAATGTTAACCCTGCAACAGACGAATCAACAAAGAACCGTTATACATCGCGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAACTGAGTACTCTCAACGCTGATGTGCTAACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGATATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTtaaat                                                                                                                                            >  1:93584/1‑250 (MQ=255)
                                                                                          cAGTTAAAAAATGTTAACCCTGCAACAGACGAATCAACAAAGAACCGTTATACATCGCGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAACTGAGTACTCTCAACGCTGATGTGCTAACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGATATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCatg                                                                                                                                <  1:250892/251‑1 (MQ=255)
                                                                                                  aaaTGTTAACCCTGCAACAGACGAATCAACAAAGAACCGTTATACATCGCGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAACTGAGTACTCTCAACGCTGATGTGCTGACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGATATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCa                                                                                                                          >  2:346773/1‑249 (MQ=255)
                                                                                                   aaTGTTAACCCTGCAACAGACGAATCAACAAAGAACCGTTATACATCGCGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAACTGAGTACTCTCAACGCTGATGTGCTAACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGATATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCa                                                                                                                          >  1:310892/1‑248 (MQ=255)
                                                                                                              tGCAACAGACGAATCAACAAAGAACCGTTATACATCGCGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAACTGAGTACTCTCAACGCTGATGTGCTAACTATACGCGGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGCAGAGGATATCTGCCGCTTTACTTAACGGTGAGCCAAt                                                                                                                                                                  >  1:241610/1‑197 (MQ=255)
                                                                                                              tGCAACAGACGAATCAACAAAGAACCGTTATACATCGCGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAACTGAGTACTCTCAACGCTGATGTGCTAACTATACGCGGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGATATCTGCCGCGTTACTTAACGGTGAGCCAAt                                                                                                                                                                  <  2:241610/197‑1 (MQ=255)
                                                                                                                cAACAGACGAATCAACAAAGAACCGTTATACATCGCGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAACTGAGTACTCTCAACGCTGATGTGCTAACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGATATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGAATGCGCCCGACAt                                                                                                          >  1:97703/1‑251 (MQ=255)
                                                                                                                            caacaaAGAACCGTTATACATCGCGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAACTGAGTACTCTCAACGCTGATGTGCTAACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGATATCTGCCGCTTTACTTAACGGTGTGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGAATGCGCCCGACATGCGCGCCCAccc                                                                                              >  1:369702/1‑251 (MQ=255)
                                                                                                                             aacaaAGAACCGTTATACATCGCGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAACTGAGTACTCTCAACGCTGATGTGCTGACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGATAGCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGAATGCGCCCGACATGCTCGCCCAcccc                                                                                             <  2:326772/251‑1 (MQ=255)
                                                                                                                                  aGAACCGTTATACATCGCGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAACTGAGTACTCTCAACGCTGATGTGCTGACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGATATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGAATGCGCCCGACATGCGCGCCCACCCCGGTa                                                                                         >  1:311578/1‑250 (MQ=255)
                                                                                                                                            taCATCGCGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAACTGAGTACTCTCAACGCTGATGTGCTAACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGATATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGAATGCGCCCGACATGCGCGCCCACCCCGGTAATGATGTCttt                                                                              >  1:409465/1‑251 (MQ=255)
                                                                                                                                                   cgTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAACTGAGTACTCTCAACGCTGATGTGCTAACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGATATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGAATGCGCCCGACATGCGCGCCCACCCCGGTAATGATGTCTTTTGCCAGt                                                                       <  1:320771/251‑1 (MQ=255)
                                                                                                                                                        tttACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAACTGAGTACTCTCAACGCTGATGTGCTAACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGATATCTGCCGCTTTACTTAACGGTTAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGAATGCGCCCGACATGCGCGCCCACCCCGGTAATGATGTCTTTTGCCAGTAGAt                                                                   <  2:369702/250‑1 (MQ=255)
                                                                                                                                                          tACCAGTGCAGCGCCTGCCATCGCGCCCTGGTTAGAAAACTGCGTACTCTCAACGCTGATGTGCTAACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGATATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGAATGCGCCCGACATGCGCGCCCACCCCGGTAATGATGTCTTTTGCCAGTAGATcgc                                                                >  1:148984/1‑251 (MQ=255)
                                                                                                                                                                     cgcCTGCCATCGTGCCCTGGTTAGAAAACTGAGTACTCTCAACGCTGATGTGCTGACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGATATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGAATGCGCCCGACATGCGCGCCCACCCCGGTAATGATGTCTTTTGCCAGTAGATCGCCGCGCAATgcc                                                     <  1:335665/251‑1 (MQ=255)
                                                                                                                                                                        cTGCCATCGTGCCCTGGTTAGAAAACTGAGTACTCTCAACGCTGATGTGCTAACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGATATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGAATGCGCCCGACATGCGCGCCCACCCCGGTAATGATGTCTTTTGCCAGTAGATCGCCGCGCAATgccgcc                                                  >  2:324485/1‑251 (MQ=255)
                                                                                                                                                                          gCCATCGTGCCCTGGTTAGAAAACTGAGTACTCTCAACGCTGATGTGCTAACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGATATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTTGGTTAAATAAATTCACCATGAGGGCAAGAATGCGCCCGACATGCGCGCCCACCCCGGTAATGATGTCTTTTGCCAGTAGATCTCCGCGCAATGCCTCCTg                                                <  2:454706/251‑1 (MQ=255)
                                                                                                                                                                            cATCGTGCCCTGGTTAGAAAACTGAGTACTCTCAACGCTGGTGTGCTAACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGATATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGAATGCGCCCGACATGCGCGCCCACCCCGGTAATGATGTCTTTTGCCAGTAGATCGCCGCGCAATGCCGCCTGac                                              >  2:339719/1‑251 (MQ=255)
                                                                                                                                                                                    ccTGGTTAGAAAACTGAGTACTCTCAACGCTGATGTGCTAACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGATATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGAATGCGCCCGACATGCGCGCCCACCCCGGTAATGATGTCTTTTGCCAGTAGATCGCCGCGCAATGCCGCCTGACACAATg                                        >  1:348949/1‑249 (MQ=255)
                                                                                                                                                                                                aCTGAGTACTCTCAACGCTGATGTGCTAACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGATATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGAATGCGCCCGACATGCGCGCCCACCCCGGTAATGATGTCTTTTGCCAGTAGATCGCCGCGCAATGCCGCCTGACACAATGAGTCCACGGTTAAc                          <  2:409465/251‑1 (MQ=255)
                                                                                                                                                                                                        ctctCAACGCTGATGTGCTAACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGATATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGAATGCGCCCGACATGCGCGCCCACCCCGGTAATGATGTCTTTTGCCAGTAGATCGCCGCGCAATGCCGCCTGACACAATGAGTCCACGGTTAACGGTTGTcc                  >  1:4744/1‑251 (MQ=255)
                                                                                                                                                                                                            cAACGCTGATGTTCTAACTAGACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGATATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGAATGCGCCCGACATGCGCGCCCACCCCGGTAATGATGTCTTTTGCGAGTAGATc                                                                  <  1:115131/199‑1 (MQ=255)
                                                                                                                                                                                                            cAACGCTGATGTGCTAACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGATATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGAATGCGCCCGACATGCGCGCCCACCCCGGTAATGATGTCTTTTGCCAGTAGATc                                                                  >  2:115131/1‑199 (MQ=255)
                                                                                                                                                                                                                         cTGACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGATCTCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGAATGCGCCCGACATGCGCGCCCACCCAGGTAATGATGTCTTTTGCCAGTAGATCGCCGCGCAATGCCGCCTGACACTATGAGTCCCCGGTTAACGGTTGTCCATGTAACATCGAGCTc  >  2:266695/1‑250 (MQ=255)
                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                       
TCACAAATATAAACCAGGAAAATAATTAACCTTGAAAGTCTAAGTTATGCTTTCCTGGCCCAAATTGAGATAGCGCAAATTTTGGTAGAACAGTTAAAAAATGTTAACCCTGCAACAGACGAATCAACAAAGAACCGTTATACATCGCGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAACTGAGTACTCTCAACGCTGATGTGCTGACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGATATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGAATGCGCCCGACATGCGCGCCCACCCCGGTAATGATGTCTTTTGCCAGTAGATCGCCGCGCAATGCCGCCTGACACAATGAGTCCACGGTTAACGGTTGTCCATGTAACATCGAGCTC  >  NC_000913/1667241‑1667707

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: