Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,624,622 C→T V56V (GTG→GTA ytjB ← putative inner membrane protein Smp

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,624,6220CT91.4% 90.9 / 0.9 35V56V (GTG→GTAytjBputative inner membrane protein Smp
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base T (16/16);  total (18/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.04e-01

GTCCGTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAAATTTGAGACTTAATGTTGCCAGAAGCAATGGATACAAGGTAGCCTCATGCGTT  >  NC_000913/4624378‑4624858
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                             
gTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTgg                                                                                                                                                                                                                                         >  1:259552/1‑250 (MQ=255)
gTCCGTTTTTCCCGGCAATTGGTTCGACAACCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATCCACGCCCGCGTCGAGGATACGGCTTTCATCCGTTATCTGATCGAGGATCGCCTGAATGCGTTGTTCATCCGTCGAGTCGGTACGCATCAGCGGTGCAACTTTCAGCGTTACCTGa                                                                                                                                                                                                                                         <  1:331406/251‑1 (MQ=255)
          cccggCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAgg                                                                                                                                                                                                                                <  1:130340/246‑1 (MQ=255)
                 aTTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGc                                                                                                                                                                                                                        >  1:96470/1‑251 (MQ=255)
                       tCGACAATCTGCGGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGtt                                                                                                                                                                                                                   <  1:58309/250‑1 (MQ=255)
                         gACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCAAGCCGGTCGCGCACTTCGACGCTTTCGCACGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGAAACGGCTTTCATACGTTAACTGCTCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGCTCAGCGTTACCTGACGCGCCAGGGTGAGGGCCAGttctt                                                                                                                                                                                                                >  1:182916/1‑251 (MQ=255)
                                       tAAAATAGCCGCCCGCTTTTTTACCGTCTAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGCGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGGGCAAAGTTCAGCGTTACCTGACGCGCCAGGGGGCGGGCCCGTTCTTCCGGCTGTgggta                                                                                                                                                                                                   >  2:284025/1‑247 (MQ=255)
                                        aaaaTAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTGTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTa                                                                                                                                                                                                  <  1:284025/250‑1 (MQ=255)
                                            tAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGCTAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACgttg                                                                                                                                                                                             >  1:279220/1‑248 (MQ=255)
                                              gcctcccgCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGACGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTc                                                                                                                                                                                           <  1:80294/244‑1 (MQ=255)
                                                        tttttACCGTCGGGCCCCAGCCGGTCGCGCACTTCGACGCTTTCGCCCGAACGTGCGATAAGACCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTggtggt                                                                                                                                                                                  <  2:96470/250‑1 (MQ=255)
                                                                  cGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAAc                                                                                                                                                                       <  1:295830/251‑1 (MQ=255)
                                                                      cgcCAGCCGGTCGCGAACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTCCCTGACGCGCCAGGGTGCGGGCCAGTTCTTCAAGCTGTGGATTACGTTGTCGCTGGTGGTTTTAACTAAACAATg                                                                                                                                                                  >  1:187333/1‑251 (MQ=255)
                                                                                  gcgcACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGGTCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGc                                                                                                                                                       >  1:43477/1‑250 (MQ=255)
                                                                                                 ttCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAAc                                                                                                                                        <  1:64123/250‑1 (MQ=255)
                                                                                                    gCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGCGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGAGCCAGGGAGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCAGCAACGCCCCTAACaggg                                                                                                                                    >  1:55798/1‑248 (MQ=255)
                                                                                                          acgggCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGAGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACACGGCAAGAc                                                                                                                              >  1:58416/5‑251 (MQ=255)
                                                                                                           cgggCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGCTCGAGGATCGCCTGAATGCGTTTTTCATACGGCGAGTCGGTACGCATCAGCGGGGCAACGTTCAGCGTTACCTGACGCCCCAGGGGGCGGGCCCGTTCATCCAGCTGTGGATGACGTTTTCGCTGGTGGTTTTGCATAACCCATTGACCTCCCAGCATCCGCGCCCCAACCAAGGCCAGAc                                                                                                                              >  2:195673/4‑250 (MQ=255)
                                                                                                           cgggCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACa                                                                                                                             <  1:332829/248‑1 (MQ=255)
                                                                                                           cgggCGATAAGATAGACTTGTTCGTCCTACACGCCCGCGACGAGGATACGGCTTTCATACGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCCGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCAAGCTGTGGATTACGTTGTAGCTGGTGGTTTTGCCTAAACCATGACGCTCCCTGCATAAGCGACACTACTAAGGCAAgcc                                                                                                                              >  2:321267/4‑248 (MQ=255)
                                                                                                                 agaaGATCGCTTTTTTCGTCATACACGCCCGCGGCTATGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCACCCGGCGAGTCGGTACGCATCAGCGTTGCAACTTTCAGCGTTACCTGACGCGTCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAg                                                                                                                       <  2:324571/249‑1 (MQ=255)
                                                                                                                                tCGTCATACTCTCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACTTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCGGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACACGGCAATACAGAACAGGACAATCTCTGccc                                                                                                         <  2:307536/250‑1 (MQ=255)
                                                                                                                                 cGTCATACACACCCGCGTCGAGGATACGTCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCTTTTTTCATCCGGTGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACCAACAAGGCAAGACAGAACAGGACGATCACTGCCCgg                                                                                                       >  1:41237/1‑250 (MQ=255)
                                                                                                                                                                ttCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGcc                                                                         >  1:258307/1‑250 (MQ=255)
                                                                                                                                                                     ccGTTAACTGAGCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGTTTCAACGTTCAGCGTTACCTGACGCGGCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAGCAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCtt                                                                    <  1:31330/250‑1 (MQ=255)
                                                                                                                                                                            gcggtgaaGGATCGCTTGAATGCGTTTTACATCCGGCGATTCGTGAGTTGTCAGCGGTTCCACGTTGAGCTTTACCTGACGGGCCTGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTGTTGACTAACCCATGACGCTCCCTGCAGCAGCGCCACTAAAAAGGCAAGTCAGAACAGGACAATCACTGCCCGATGCATCCGGAATGTCAGTTTTGTGCGAGCTATCTTCCAACCtt                                                            <  2:245741/244‑1 (MQ=255)
                                                                                                                                                                              gATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTg                                                          <  1:211432/251‑1 (MQ=255)
                                                                                                                                                                                    ggATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAAAtt                                                    <  1:142685/251‑1 (MQ=255)
                                                                                                                                                                                            tGAATGCGTTTTTCATCCGGCGAGTCGGTACGCCTCAGCGGTGCAACGTTCAGCGTTACCTGAAGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGAGCGAGCCATCTTCCACCCTTTGAAAATTTGAGACtt                                            >  1:324834/1‑251 (MQ=255)
                                                                                                                                                                                             gAATGCGTTTTTAATCCAACAAGTCGGTACGCATCAACGGTGCAACGTGCAGCGTTACCTGGCGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTACACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACCGGACAATAACTGCCCGCTGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAAATTTgagaattt                                           >  2:257117/1‑247 (MQ=255)
                                                                                                                                                                                                       ttCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAAATTTGAGACTTAATGTTGCCa                                  <  1:195673/250‑1 (MQ=255)
                                                                                                                                                                                                                           gCATCAGCGGTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGCTGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAAATTTGAGACTTAATGTTGCCAGAAGCAATGGATACAAGGt               >  1:363087/1‑249 (MQ=255)
                                                                                                                                                                                                                                    gTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCgagatcggaagagc                                                                                  >  2:167567/1‑160 (MQ=255)
                                                                                                                                                                                                                                    gTGCAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCgagatcggaagagc                                                                                  <  1:167567/173‑14 (MQ=255)
                                                                                                                                                                                                                                       cAACGTTCAGCGTTACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGTTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAAATTTGAGACTTAATGTTGCCAGAAGCAATGGATACAAGGTAGCCTCATGCGtt  >  1:202867/1‑250 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                             
GTCCGTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAAATTTGAGACTTAATGTTGCCAGAAGCAATGGATACAAGGTAGCCTCATGCGTT  >  NC_000913/4624378‑4624858

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