Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,624,433 T→C 7.7% K119K (AAA→AAG)  ytjB ← protein Smp

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,624,4330TC 7.7% 125.9 / 2.1 53K119K (AAA→AAG) ytjBprotein Smp
Reads supporting (aligned to +/- strand):  ref base T (26/22);  new base C (2/2);  total (29/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.54e-01

GCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTT  >  NC_000913/4624292‑4624575
                                                                                                                                             |                                                                                                                                              
gCGTAAAATGTTAGTGGTGTTATCCGCCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGGTTTTACCCGCAATTGGCTCGACAATCCGCTGGTTGCAATAGCCCGCTGCtt                                                                                                                                               >  2:932033/1‑143 (MQ=255)
gCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCtt                                                                                                                                               >  2:948908/1‑143 (MQ=255)
gCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCCAAGATAGCCGAGCGGTCCGTTGTTACCCGCAATTGGCGCGACAACCTGCTGGCTAAAATAGCCGCCTGCtc                                                                                                                                               >  2:398599/1‑142 (MQ=255)
   tAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCttttt                                                                                                                                            >  1:815960/1‑143 (MQ=255)
     aaaTGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATCTGCTGGCTAAAATAGCCGCCTGCTTTTTTa                                                                                                                                          >  2:515473/1‑143 (MQ=255)
       aggTTTGTGGTGTTGTCCGCCTGTTGGGCTTCGGTGGCGAGCGTGTGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTAcc                                                                                                                                        <  2:236461/141‑1 (MQ=255)
         gTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGt                                                                                                                                      <  1:670910/143‑1 (MQ=255)
          ttAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTACAATAGCCGCCTGCTTTTTTACCGTc                                                                                                                                     >  1:1084146/1‑143 (MQ=255)
          ttAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGACGTGCGCAAACAGGGGTCCCTCGGTGAAGACTCCGCCGTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTc                                                                                                                                     >  1:1130271/1‑143 (MQ=255)
             ggggtggtATCCACCGGTGGGGCTTCGGTGGCGAGCGTATGTGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAg                                                                                                                                  <  2:458567/141‑1 (MQ=255)
                         aCCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGGGCGGTCCGTTTTTACGCGCACTTGGCTCGACAATCTGCTGGTTAACATAGCCGCCTGGTCTTTTACCGTCGAGCGCCAGGCgggc                                                                                                                      >  1:678480/1‑141 (MQ=255)
                         aCCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATCTGATGGTTAAAACAGCCCCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTc                                                                                                                      >  2:496261/1‑143 (MQ=255)
                            tGTTGGGCTTCGGTGGCGAGCGTATGTGTGTCGATTGTCAGGCGCATTTTTTTGTGCGGTTCTTTTTTAACCGGAATTGGCGCGACAATCGGCGGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTcgcg                                                                                                                   <  1:40812/143‑1 (MQ=255)
                             gTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCGCCATTGGCTCGACAATCTGCTCGTTAAAATAGCCGCCTGCTTTCTTACCGTCTAGCGCCAACCCGTcgcca                                                                                                                  >  1:224595/1‑141 (MQ=255)
                             gTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCATTTTTACCGTCGAGCGCCAGCAGGTCgcgc                                                                                                                  >  2:374686/1‑143 (MQ=255)
                              ttGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCa                                                                                                                 <  1:685483/143‑1 (MQ=255)
                                 ggCTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGACGCCGGCCGTTTTACCGTCGAGCGCCAGCCGGTCGCGCACtt                                                                                                              >  2:124385/1‑143 (MQ=255)
                                    gtCGTTGGCGAGCGTATGGGGGCCGGGTGTCAGGCGCAGACAGCCGAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTCACCGTCGAGCGCCAGCCGGTCGCGCACTTCGa                                                                                                           <  2:591553/142‑1 (MQ=255)
                                        gTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCt                                                                                                       <  1:788882/143‑1 (MQ=255)
                                          ggCGAGCGTAGGGGTGGCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCttt                                                                                                     <  2:96989/143‑1 (MQ=255)
                                             gAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCAAACCGGACGCGCACTTCGACGCTATCGc                                                                                                  >  1:624094/1‑143 (MQ=255)
                                                     gggTGGCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGt                                                                                          <  1:655525/143‑1 (MQ=255)
                                                      ggTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCTAGCGCCAGCCGGTCGCGCACTTCGCCGCTTTCGCCATAACGTg                                                                                         >  1:352303/1‑143 (MQ=255)
                                                      ggTGTCGAGTGTCAGGCGCAGATAGCCCAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGGCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTg                                                                                         <  1:416759/143‑1 (MQ=255)
                                                            gAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTCCCGTCGAGCGCCAGCCGGTCCCGCACTTCGACGCTGTCGCCCGAACATGCGATac                                                                                   >  2:619181/1‑142 (MQ=255)
                                                            gAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTAGCGTCGAGCGCCAGCCGGCCGCGCCCTTCGACGCTTTCGCCAGAACATACGATaa                                                                                   >  1:575556/1‑143 (MQ=255)
                                                            gAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCCCTTCGACGCTTTCGCCAGAAAGTGCGATaa                                                                                   >  1:752053/1‑143 (MQ=255)
                                                               tgtCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCCCTTCGACGCTTTCGCCAGAACGTGCGATAAgcc                                                                                >  2:148300/1‑141 (MQ=255)
                                                                     gcgcAGATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGCGCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCtt                                                                          <  1:635195/143‑1 (MQ=255)
                                                                       gcAGATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGt                                                                        <  1:204001/143‑1 (MQ=255)
                                                                         aGATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTc                                                                      <  1:639466/143‑1 (MQ=255)
                                                                         aGATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGCACCTGCGATAAGACCGCCCTGCTc                                                                      >  2:415728/1‑143 (MQ=255)
                                                                          gATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATCTCCTGGTTAAAATAGCCGCCTGCTTTTTCAACGTCGATTGCTTCTCGGTCGCGCACTCCGCAGCTTTCGCCAGGCCGCGCGATAATATCGCCCTGCTCg                                                                     >  2:279213/1‑143 (MQ=255)
                                                                               ccttgtggTCCGTTTTTACCCGCAATTGGCTCGACAAGCTGCAGGTTAAAATAGCCGCCTGCCTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATa                                                                <  1:124385/137‑1 (MQ=255)
                                                                                 gAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCCTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATc                                                                               >  2:1017225/1‑126 (MQ=255)
                                                                                 gAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCCTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATc                                                                               <  1:1017225/126‑1 (MQ=255)
                                                                                    cGGTCCGTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGCGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATccatac                                                           >  2:967618/1‑143 (MQ=255)
                                                                                         cGTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGCCGCCCACTCCGGCGCTTTCGCCAGACCCTCCGATACGATCGCCTCGTTCGTCCTCCCCACCcgcg                                                      >  2:129841/1‑143 (MQ=255)
                                                                                             tttACCCTCAAGTGGCTCGACACTCTGCTCGTTAAAATAGCCGCCGGCTTTTTTACCGTAGAGCGCCAGCTGCTCGCGAACTTCGACGCTTTTGCCAGCTCATGCGTTAATATCGCCTTGTTCGTCATACACACCCGCGTCGa                                                  <  1:687379/143‑1 (MQ=255)
                                                                                                 cccGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGCCGACCCCCAGCCGGTCGCGCACTTCGCCGCTCTCGCCACACCGTCCGATAAGACCGCCCTGCTCCTCATACACCACCGCGTCGAGGac                                              >  2:124939/1‑142 (MQ=255)
                                                                                                    gCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACg                                           <  1:948908/143‑1 (MQ=255)
                                                                                                     cAATTGCTTCGACAATCTGCTAAGGAAACTATCCGCCTGCTTTTTTTCCGTCGAGCGCCCGCCGGTCGCTGACTGCTACGCTTTCGACAGAACGGGCGAGAAGATCGCCTTGTTCTTCATACACACCCGCGTCGAGGATACgg                                          <  2:599742/143‑1 (MQ=255)
                                                                                                         tGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCTAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCttt                                      <  2:524126/143‑1 (MQ=255)
                                                                                                                aCAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGt                               <  2:815960/143‑1 (MQ=255)
                                                                                                                   aTCGGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGCCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATTAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTaa                            <  1:712664/143‑1 (MQ=255)
                                                                                                                      tgctgGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGGACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTGATCCGTTAACTg                         >  1:586244/1‑143 (MQ=255)
                                                                                                                         tgGTTAATATGGCCGCCTGCTTTTTTAGCGTCTGGCGCCAGCCGGTCGCGGACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATc                      <  2:104515/143‑1 (MQ=255)
                                                                                                                           gTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACTCACCCGCTTCGGGGCTCCGGCTTTCATCCGTTACCTGATcca                    >  1:99532/1‑141 (MQ=255)
                                                                                                                           gTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTAATAGa                    >  2:306651/1‑143 (MQ=255)
                                                                                                                               aaaTAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCCTCCGTTAACTGATCGAGGAt                >  1:628002/1‑143 (MQ=255)
                                                                                                                                 aTAGCCGCCTGCTTTTTTACCGTCGTGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTGTCATGCATTACCTGGTTGAGGGTCg              >  2:990843/1‑143 (MQ=255)
                                                                                                                                   aGCCGCCTGCTTTTGTACCTTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGcc            <  2:471280/143‑1 (MQ=255)
                                                                                                                                    gccgccGGCTTTTTTACCGTCGAGTCCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGAGCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCt           <  2:1084146/143‑1 (MQ=255)
                                                                                                                                     ccgccTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACCCCCGCGTCGAGGATACGGCTTTCATCCATTAACTCATCGAGGATCGccgg          >  1:742599/1‑141 (MQ=255)
                                                                                                                                        ccTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTACCTGATCGAGGATCGCCTGAAt       >  1:119488/1‑143 (MQ=255)
                                                                                                                                          tGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATGCGTTAACTGATAGAGGAACGCCAGAATGc     >  1:43827/1‑143 (MQ=255)
                                                                                                                                          tGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGc     >  1:834566/1‑143 (MQ=255)
                                                                                                                                             ttttttACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTGCATCCGTTAACTGATCGAGGATCGCCTGAATGCGtg  >  1:741855/1‑142 (MQ=255)
                                                                                                                                             |                                                                                                                                              
GCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTT  >  NC_000913/4624292‑4624575

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