Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,632,369 T→C intergenic (‑84/‑703) ydfJ ← / → ydfK pseudogene, MFS transporter family; interrupted by Qin prophage;Phage or Prophage Related; putative transport protein/cold shock protein, function unknown, Qin prophage

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,632,3690TC98.2% 203.1 / NA 56intergenic (‑84/‑703)ydfJ/ydfKpseudogene, MFS transporter family; interrupted by Qin prophage;Phage or Prophage Related; putative transport protein/cold shock protein, function unknown, Qin prophage
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (33/22);  minor base A (1/0);  total (34/22)
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

GAAAAAATATTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCATCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCATATATGCCCCAAAAAAACTCCCCAACAGATAAGTAGTTTTTTCATGGATTTATGCGTAAAATCAAGAACGGCTGGAAATCATTCAATACTCACACT  >  NC_000913/1632228‑1632510
                                                                                                                                             |                                                                                                                                             
gAAAAAATATTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCACCCTAAGTTc                                                                                                                                      >  2:400121/1‑151 (MQ=255)
        atTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCACCCTAAGTTCGAtctc                                                                                                                                >  1:20746/1‑149 (MQ=255)
        atTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCACCCTAAGTTCGATCTCAg                                                                                                                              <  1:514070/151‑1 (MQ=255)
            cATGAAACACTAACGCTGCGCCGAGCGAATATAACGGGAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCACCCTAAGTTCGATCTCAGTGTc                                                                                                                          >  2:246661/1‑151 (MQ=255)
                 aaCACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCACCCTAAGTTCGATCTCAGTGTCTATCt                                                                                                                     <  1:154669/151‑1 (MQ=255)
                     cTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCACCCTAAGTTCGATCTCAGTGTCTATCTgggg                                                                                                                 >  2:346563/1‑151 (MQ=255)
                      tAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCACCCTAAGTTCGATCTCAGTGTCTATCTGGGGc                                                                                                                >  2:793345/1‑151 (MQ=255)
                                  gAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCGAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCACCCTAAGTTCGCTCTCAGTGTCTATCTGGCGCCTATTTCTGTc                                                                                                     <  2:14512/150‑1 (MQ=255)
                                   aGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCACCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTccc                                                                                                   <  2:323137/151‑1 (MQ=255)
                                           taACTGGAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCACCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCATATATGc                                                                                           >  1:409925/1‑151 (MQ=255)
                                           taACTGGAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCACCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCATATATGc                                                                                           >  2:18697/1‑151 (MQ=255)
                                            aaCTGGAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCACCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCATATATGcc                                                                                          >  2:750826/1‑151 (MQ=255)
                                                ggAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCACCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCATATATGcccc                                                                                        <  1:303150/149‑1 (MQ=255)
                                                ggAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCACCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCATATATGCCCCa                                                                                       >  1:384162/1‑150 (MQ=255)
                                                      ccATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCACCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCATATATGCCCCAAAAAAAc                                                                                >  2:278160/1‑151 (MQ=255)
                                                      ccATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCACCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCATATATGCCCCAAAAAAAc                                                                                >  1:829641/1‑151 (MQ=255)
                                                                     gtTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCACCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCATATATGCCCCAAAAAAACTCCCCAACAGATa                                                                   >  1:799543/1‑149 (MQ=255)
                                                                     gtTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCACCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCATATATGCCCCAAAAAAACTCCCCAACAGATa                                                                   >  1:690592/1‑149 (MQ=255)
                                                                         tttCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCACCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCATATATGCCCCAAAAAAACTCCCCAACAGATAagtagt                                                             >  2:122588/1‑151 (MQ=255)
                                                                         tttCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCACCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCATATATGCCCCAAAAAAACTCCCCAACAGATAagtag                                                              <  1:166492/150‑1 (MQ=255)
                                                                           tCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCACCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCATATATGCCCCAAAAAAACTCCCCAACAGATAAGTAGttt                                                           >  1:507139/1‑151 (MQ=255)
                                                                             ccATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCACCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCATATATGCCCCAAAAAAACTCCCCAACAGATAAGTAGtttt                                                          <  2:200053/150‑1 (MQ=255)
                                                                                 ccTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCACCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCATATATGCCCCAAAAAAACTCCCCAACAGATAAGTAGTTTTTTCa                                                      <  1:644445/150‑1 (MQ=255)
                                                                                   tgtgGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCACCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCATATATGCCCCAAAAAAACTCCCCAACAGATAAGTAGTTTTTTCATgg                                                   >  2:900870/1‑151 (MQ=255)
                                                                                    gtgGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCACCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCATATATGCCCCAAAAAAACTCCCCAACAGATAAGTAGTTTTTTCATg                                                    <  1:394518/149‑1 (MQ=255)
                                                                                      ggTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCGCCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCATATATGCCCCAAAAAAACTCCCCAACAGATAAGTAGTTTTTTCATGGAtt                                                <  2:650866/151‑1 (MQ=255)
                                                                                      ggTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCACCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCATATATGCCCCAAAAAAACTCCCCAACAGATAAGTAGTTTTTTCATGGAtt                                                >  2:635587/1‑151 (MQ=255)
                                                                                      ggTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCACCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCATATATGCCCCAAAAAAACTCCCCAACAGATAAGTAGTTTTTTCATGGAtt                                                <  2:339969/151‑1 (MQ=255)
                                                                                      ggTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCACCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCATATATGCCCCAAAAAAACTCCCCAACAGATAAGTAGTTTTTTCATGGAt                                                 <  2:835293/150‑1 (MQ=255)
                                                                                          tCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCACCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCATATATGCCCCAAAAAAACTCCCCAACAGATAAGTAGTTTTTTCATGGATTTAt                                             <  1:293637/150‑1 (MQ=255)
                                                                                              tGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCACCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCATATATGCCCCAAAAAAACTCCCCAACAGATAAGTAGTTTTTTCATGGATTTATGCGt                                         <  1:523658/150‑1 (MQ=255)
                                                                                                gCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCACCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCATATATGCCCCAAAAAAACTCCCCAACAGATAAGTAGTTTTTTCATGGATTTATGCGTaaa                                      <  1:531901/151‑1 (MQ=255)
                                                                                                 cGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCACCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCATATATGCCCCAAAAAAACTCCCCAACAGATAAGTAGTTTTTTCATGGATTTATGCGTaaaa                                     <  2:748273/151‑1 (MQ=255)
                                                                                                  gTTTTCTAGCTTTTTTCAGCTAGTTGCATTTTTTTAAAAAGCAACCTAAGTACGATCTCAGTGACTATCTGCGGCCTATTTCTGTCCCATATATGCCCCTCAAAAACTCCCCAACCGATAAGTAGTTTTTTCATGGATTTATGCGCaaaa                                     >  2:246449/1‑150 (MQ=255)
                                                                                                   ttttCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCACCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCATATATGCCCCAAAAAAACTCCCCAACAGATAAGTAGTTTTTTCATGGATTTATGCGTAAAATc                                   >  1:664536/1‑151 (MQ=255)
                                                                                                    tttCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCACCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCATATATGCCCCAAAAAAACTCCCCAACAGATAAGTAGTTTTTTCATGGATTTATGCGTAAAATc                                   <  1:599901/150‑1 (MQ=255)
                                                                                                    tttCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCACCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCATATATGCCCCAAAAAAACTCCCCAACAGATAAGTAGTTTTTTCATGGATTTATGCGTAAAATc                                   >  2:858092/1‑150 (MQ=255)
                                                                                                    tttCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCACCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCATATATGCCCCAAAAAAACTCCCCAACAGATAAGTAGTTTTTTCATGGATTTATGCGTAAAATc                                   >  2:386453/1‑150 (MQ=255)
                                                                                                     ttCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCACCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCATATATGCCCCAAAAAAACTCCCCAACAGATAAGTAGTTTTTTCATGGATTTATGCGTAAAATCa                                  >  1:454045/1‑150 (MQ=255)
                                                                                                       cTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCACCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCATATATGCCCCAAAAAAACTCCCCAACAGATAAGTAGTTTTTTCATGGATTTATGCGTAAAATCa                                  >  1:603404/1‑148 (MQ=255)
                                                                                                            ttttttCAGATAGTTGCATTTTTTTAAAAAGCACCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCATATATGCCCCAAAAAAACTCCCCAACAGATAAGTAGTTTTTTCATGGATTTATGCGTAAAATCAAGAAc                             >  1:491251/1‑148 (MQ=255)
                                                                                                              ttttCAGATAGTTGCATTTTTTTAAAAAGCACCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCATATATGCCCCAAAAAAACTCCCCAACAGATAAGTAGTTTTTTCATGGATTTATGCGTAAAATCAAGAACGGCt                         <  2:799478/150‑1 (MQ=255)
                                                                                                              ttttCAGATAGTTGCATTTTTTTAAAAAGCACCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCATATATGCCCCAAAAAAACTCCCCAACAGATAAGTAGTTTTTTCATGGATTTATGCGTAAAATCAAGAACGGCTg                        >  2:241977/1‑151 (MQ=255)
                                                                                                                ttCAGATAGTTGCATTTTTTTAAAAAGCACCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCATATATGCCCCAAAAAAACTCCCCAACAGATAAGTAGTTTTTTCATGGATTTATGCGTAAAATCAAGAACGGCTgg                       >  1:202888/1‑150 (MQ=255)
                                                                                                                 tCAGATAGTTGCATTTTTTTAAAAAGCACCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCATATATGCCCCAAAAAAACTCCCCAACAGATAAGTAGTTTTTTCATGGATTTATGCGTAAAATCAAGAACGGCTGGa                      >  1:140194/1‑150 (MQ=255)
                                                                                                                 tCAGATAGTTGCATTTTTTTAAAAAGCACCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCATATATGCCCCAAAAAAACTCCCCAACAGATAAGTAGTTTTTTCATGGATTTATGCGTAAAATCAAGAACGGCTGGa                      <  1:212284/150‑1 (MQ=255)
                                                                                                                  cAGATAGTTGCATTTTTTTAAAAAGCACCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCATATATGCCCCAAAAAAACTCCCCAACAGATAAGTAGTTTTTTCATGGCTTTATGCGTAAAATCAAGAACGGCTGGaa                     >  1:173609/1‑150 (MQ=255)
                                                                                                                  cAGATAGTTGCATTTTTTTAAAAAGCACCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCATATATGCCCCAAAAAAACTCCCCAACAGATAAGTAGTTTTTTCATGGATTTATGCGTAAAATCAAGAACGGCTGGaa                     <  2:899368/150‑1 (MQ=255)
                                                                                                                       aGTTGCATTTTTTTAAAAAGCACCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCATATATGCCCCAAAAAAACTCCCCAACAGATAAGTAGTTTTTTCATGGATTTATGCGTAAAATCAAGAACGGCTGGAAATCAtt               >  1:54589/1‑151 (MQ=255)
                                                                                                                       aGTTGCATTTTTTTAAAAAGCACCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCATATATGCCCCAAAAAAACTCCCCAACAGATAAGTAGTTTTTTCATGGATTTATGCGTAAAATCAAGAACGGCTGGAAATCAt                <  1:123982/150‑1 (MQ=255)
                                                                                                                         ttGCATTTTTTTAAAAAGCACCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCATATATGCCCCAAAAAAACTCCCCAACAGATAAGTAGTTTTTTCATGGATTTATGCGTAAAATCAAGAACGGCTGGAAATCATTc              >  1:291771/1‑150 (MQ=255)
                                                                                                                          tGCATTTTTTTAAAAAGCACCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCATATATGCCCCAAAAAAACTCCCCAACAGATAAGTAGTTTTTTCATGGATTTATGCGTAAAATCAAGAACGGCTGGAAATCATTCa             >  2:97455/1‑150 (MQ=255)
                                                                                                                                tttttAAAAAGCACCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCATATATGCCCCAAAAAAACTCCCCAACAGATAAGTAGTTTTTTCATGGATTTATGCGTAAAATCAAGAACGGCTGGAAATCATTCAATACTc       >  1:794742/1‑150 (MQ=255)
                                                                                                                                  tttAAAAAGCACCCTAATTTCGATCTCAGTGTCTATCTGAGGCCTATTTCTGTCCCATATATGCCCCAAAAAAACTCCCCAACAGATAAGTAGTTTTTTCATGGATTTATGCGTAAAATCAAGAACGGCTGGAAATCATTCAATAcgcgc     <  2:539453/150‑5 (MQ=255)
                                                                                                                                     aaaaaGCACCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCATATATGCCCCAAAAAAACTCCCCAACAGATAAGTAGTTTTTTCATGGATTTATGCGTAAAATCAAGAACGGCTGGAAATCATTCAATACTCACACt  >  2:899652/1‑150 (MQ=255)
                                                                                                                                      aaaaCCACCCTAAGTTCGATCTGAGTGTCTATCTGGGGCCTATTTCTGTCCCATATATGCCCCAAAAAAACTCCCCAACAGATAAGTAGTTTTTTCATGGATTTATGCGTAAAATCAAGAACGGCTGGAAATCATTCAATACTCACACt  >  2:699501/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                             
GAAAAAATATTTCATGAAACACTAACGCTGCGCCGAGCGAATATAACTGGAAATCCATAATTCACAGGTGTTTTTTCCCATCCTGTGGTTTCCTTGGCGTTTTCTAGGTTTTTTCAGATAGTTGCATTTTTTTAAAAAGCATCCTAAGTTCGATCTCAGTGTCTATCTGGGGCCTATTTCTGTCCCATATATGCCCCAAAAAAACTCCCCAACAGATAAGTAGTTTTTTCATGGATTTATGCGTAAAATCAAGAACGGCTGGAAATCATTCAATACTCACACT  >  NC_000913/1632228‑1632510

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