Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,579,530 C→T G104E (GGA→GAA)  symE ← nucleoid‑associated protein SymE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,579,5300CT94.1% 92.6 / ‑2.9 34G104E (GGA→GAA) symEnucleoid‑associated protein SymE
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (22/10);  total (23/11)
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

TTGTTGCTGTGCTTGCTTAATAAATTTTTCAATCCAGGGTTGAATAGATTCCATAAGATATGCCTTCCTCATTGCTAAGCCTCTATTATCGCTTTCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCT‑GCACCTGCTTTTGTTT‑ACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGTG  >  NC_000913/4579284‑4579772
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    
ttgttgCTGTGCTTGCTTAATAAATTTTTCAATCCAGGGTTGAATAGATTCCATAAGATATGCCTTCCTCATTGCTAAGCCTCTATTATCGCTTTCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGTGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGATGTTACGCGACTTTCTGTTTACCGGCAATCACTTCAAt                                                                                                                                                                                                                                                  >  1:297846/1‑251 (MQ=255)
       tgtgCTTGCTTAATAAATTTTTCAATCCAGGGTTGAATAGATTCCATAAGATATGCCTTCCTCATTGCTAAGCCTCTATTATCGCTTTCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTTAAATAAACGcc                                                                                                                                                                                                                                           >  2:1317/1‑251 (MQ=255)
             tgcttAATAAATTTTTCAATCCAGGGTTGAATAGATTCCATAAGATATGCCTTCCTCATTGCTAAGCCTCTATTATCGCTTTCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTTCAATAAACGCCT‑GCAc                                                                                                                                                                                                                                     <  1:488663/250‑1 (MQ=255)
                        tttttCATTCCAGGGTTGAATAGATTCCCTAAGATATACCTACCTCATTGATAAGCATCTATTATCGCATTCGAAACGAACTCAAACAATAGACTTTTCCTGCAAAATCAGACTGGTAAATATTTAATGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAAACGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTTCAATAAACGCCT‑GCACCTGATTTTGttt                                                                                                                                                                                                                         >  2:285522/1‑251 (MQ=255)
                                       ttGAATAGATTCAATACGATACCCCTTCCTCATTGCTAAGTCTCTCTTATCGATTTCCCAACGTACTGAAACAATACACTTTTACTGCAACATCAGACTGCTAAATATTTACTGAGGGGGCAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGCGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTTCAATAAACGCCT‑GCACCTGCTTTTGTTT‑ACGCGCCGACAGttc                                                                                                                                                                                                         >  2:616301/1‑250 (MQ=255)
                                                ttCCATAAGATATGCCTTCCTCATTGCTAAGCCTCTATTATCGCTTTCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTTCAATAAACGCCT‑GCACCTGCTTTTGTTT‑ACGCGCCGACAGTTTGCACACCTg                                                                                                                                                                                                >  2:500652/1‑251 (MQ=255)
                                                    aTAAGATATGCCTTCCTCATTGCTAAGCCTCTATTATCGCTTTCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGTAAATATTTACGGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTTCAATAAACGCCG‑GCACCTGCTTTTGTTT‑ACGCGCCGACAGTTTGCACACCGGGCGt                                                                                                                                                                                            <  1:660861/251‑1 (MQ=255)
                                                         atatGCCTTCCTCATTGCTAAGCCTCTATTATCGCTTTCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTTCAATAAACGCCT‑GCACCTGCTTTTGTTT‑ACGCGCCGACAGTTTGCACACCTGGCGTAGCg                                                                                                                                                                                        >  1:136668/1‑250 (MQ=255)
                                                                                            tttCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTTCAATAAACGCCT‑GCACCTGCTTTTGTTT‑ACGCGCCGAAAGTTTGCAAACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGGGGtt                                                                                                                                                    >  2:470041/1‑251 (MQ=255)
                                                                                               cGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTTCAATAAACGCCT‑GCACCTGCTTTTGTTT‑ACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTggg                                                                                                                                                 >  1:233951/1‑251 (MQ=255)
                                                                                                 cAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTTCAATAAACGCCT‑GCACCTGCTTTTGTTT‑ACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCg                                                                                                                                               >  2:319283/1‑251 (MQ=255)
                                                                                                  aaCGTACTGACACAATAGCTTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTTCAATAAACGCCT‑GCACCTGCTTTTGTTT‑ACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGAGGGTGGTTGGGCgg                                                                                                                                              >  2:316060/1‑251 (MQ=255)
                                                                                                     gTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGCGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGCCTTTCTGTTTACCGGCAATCACTTCAATAACCGCCT‑GAACCTGCTTTTGTTT‑ACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTAGGCGGCGGGTGGTTGGGCGGTg                                                                                                                                            >  1:311806/1‑250 (MQ=255)
                                                                                                          gAAACAATAGATTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTTCAATAAACGCCT‑GCACCTGCTTTTGTTT‑ACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGCGAGGa                                                                                                                                        <  2:561963/249‑1 (MQ=255)
                                                                                                            aaccaTAGATTTTTAATGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGTTGGTTAACCGAGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTTCAATAAACGCCT‑GCACCTGCTTTTGTTT‑ACGCGCCGACACTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGAAAAt                                                                                                                                    <  2:516083/247‑1 (MQ=255)
                                                                                                                  tagattttACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAACGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCT‑GCACCTGCTTTTGTTT‑ACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGcc                                                                                                                              >  1:87428/5‑251 (MQ=255)
                                                                                                                    aTTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTTCAATAAACGCCT‑GCACCTGCTTTTGTTT‑ACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCtt                                                                                                                            <  2:233951/251‑1 (MQ=255)
                                                                                                                        ttACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTTCAATAAACGCCT‑GCACCTGCTTTTGTTT‑ACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCa                                                                                                                         <  1:71078/250‑1 (MQ=255)
                                                                                                                                                tatTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGATCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGAAATCACTTCAATAAACGCCT‑GCACCTGCTTTTGTTT‑ACGCGCCGAAAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGGGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGcc                                                                                                >  2:480495/1‑251 (MQ=255)
                                                                                                                                                    tACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTTCAATAAACGCCT‑GCACCTGCTTTTGTTT‑ACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTg                                                                                            >  1:33528/1‑251 (MQ=255)
                                                                                                                                                           ggggAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTTCAATAAACGCCT‑GCACCTGCTTTTGTTT‑ACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGAGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAAc                                                                                     >  2:498048/1‑251 (MQ=255)
                                                                                                                                                                  gTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTTCAATAAACGCCT‑GCACCTGCTTTTGTTT‑ACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGCCGGCGGGTGGTTGGGCGGTGAGGACAATCCAGCCTTCCATCACTTTGACATCTACCGCCGTGCCCATGGCAAAACCggcggc                                                                              >  2:441074/1‑251 (MQ=255)
                                                                                                                                                                          tGAGTCAGAGGAAGGCTCCCGGTGGTTAACCGGGAGTACACCCTGTTACGCGACTTTCTGTTTACCGGCAATCACTTCAATAAACGCCT‑GCACCTGCTTTTGTTT‑ACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACCATACAGCCTTCCATCACTTTGACATCTACTGCCGTGCCAGTGGCAAAACCGTCGGCTTCCAGcc                                                                      >  2:623051/1‑251 (MQ=255)
                                                                                                                                                                                aGTGGACGGCGCCCGGGGGTAAACGGTGAGTAAGTGCGGTTACGCGACTTTCTGTTGCCCGCCAGTCACTCCAATAGACGCCTCGC‑CCTGCTTTTGTTT‑CCGCGCCGACATTTTGCACACCTGGCGTAGCGCCTGCAGCATTCCGCTCGCCTCGGCGGCGGGTGGTAGGGCGGGCAGGACAATACAGCCTCCCTTCACTTTGACATGTTCCGCCGTGCCATTGGCAAAACCGGCGGCTTCCAGCCACTGAc                                                                <  1:285522/251‑1 (MQ=255)
                                                                                                                                                                                aGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTTCAATAAACGCCT‑GCACCTGCTTTTGTTT‑ACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTAccgcc                                                                                                     >  2:5902/1‑214 (MQ=255)
                                                                                                                                                                                aGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTTCAATAAACGCCT‑GCACCTGCTTTTGTTT‑ACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTAccgcc                                                                                                     <  1:5902/214‑1 (MQ=255)
                                                                                                                                                                                           cccGGTGGTTAACCGAGAGAAAACGCTGTTACGCGACTTTCTGTTTCCCGGCAATCACTTCAATAAACGCCT‑GCACCTGCTTTTGTTT‑ACGCGCCGACAGTTTGCACACATGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGCCATCTACCGCCGTGCCAGTGGCAAAACCGCCGCCTTCCAGCCCCTGCCCTTTCAGGGTg                                                     >  2:82061/1‑251 (MQ=255)
                                                                                                                                                                                                            aGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTTCAATAAACGCCT‑GCACCTGCTTTTGTTT‑ACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGcc                                    >  2:59018/1‑250 (MQ=255)
                                                                                                                                                                                                                aaCGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTTCAATCAACGCCT‑GCACCTGCTTTTGTTT‑ACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCAATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTa                                 >  1:210396/1‑250 (MQ=255)
                                                                                                                                                                                                                  cGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTTCAATAAACGCCT‑GCACCTGCTTTTGTTT‑ACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGCATACGGCTGTAATc                              >  2:237985/1‑251 (MQ=255)
                                                                                                                                                                                                                           gcgccTTTCTGTTTCCCGGCAATCACTTCAATAAACCCCT‑GCCCCTGCTTTTTTTT‑ACGCGCCGCAAGTTTGCCCACCGGGCGTGGCGACTGCATCAGTCCGCTCTCCTCGGCGGGGGGTGTTTGGGCGGGGAGGACAATACAGCTTTCCATCATTTGGACATCTCCCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCg                     <  2:33528/247‑1 (MQ=255)
                                                                                                                                                                                                                                    tGTTTACCGGCAATCACTTCAATAAACGCCT‑GCACCTGCTTTTGTTT‑ACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATTCGGGTAGCCACTCGCAt             >  2:178483/1‑250 (MQ=255)
                                                                                                                                                                                                                                          cgggCAATAACTTCAATAACCGCCT‑CCACAT‑CCTTTGTTTGACGCGCCGACAGTGTGCCCGCCTGGCGTAGCGATTGCACCGGTTCGCTTTCCTCGGCGGGGGGTGGTTTGGCGGTGGGGACAAAACAGCTTTCATTCACTTTGCCAGCTTCCGCCGTGCGAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGGAATCCGGGTAGCGACTCGCATAACCGAc      <  2:51373/249‑1 (MQ=255)
                                                                                                                                                                                                                                              cAATAACTTCAATAAACGCCT‑GCCCCTGCTTTTTTTT‑ACGCGCGGACAGTTTGCACCCCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATTCAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGTg  <  1:470041/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    
TTGTTGCTGTGCTTGCTTAATAAATTTTTCAATCCAGGGTTGAATAGATTCCATAAGATATGCCTTCCTCATTGCTAAGCCTCTATTATCGCTTTCGCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCT‑GCACCTGCTTTTGTTT‑ACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGTG  >  NC_000913/4579284‑4579772

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