Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,741,011 G→T 7.2% L308M (CTG→ATG) ‡ PP_3305 ← TerC family membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,741,0110GT 7.2% 162.0 / 3.8 56L308M (CTG→ATG) ‡PP_3305TerC family membrane protein
Reads supporting (aligned to +/- strand):  ref base G (26/26);  new base T (2/2);  total (28/28)
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.95e-01

CAGCGGAGGCGACCGCTTTCACGCCTTGAATCCTAGGCGATCCGGTCATCGCGCTCACGCAGCAGCCAGCGGCCGTAATTATTCAAGAAGCATTAGCGGACGCTCGGGTTTTCCAAAGCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTACCCACAATGTCATGAAGGAAAATTTTCGAGCCGATAAAGACAAGCACCAGAGCCAGCGCATATTTAAGGTAAACAAAGCGACTGATCATGGCAGCCAATGCAA  >  NC_002947/3740867‑3741154
                                                                                                                                                |                                                                                                                                               
cagcGGAGGCGACCGCTTTCACGCCTTGAATCCTAGGCGATCCGGTCATCGCGCTCACGCAGCAGCCAGCGGCCGTAATTATTCAAGAAGCATTAGCGGACGCTCGGGTTTTCCAAAGCGAAATCAACACGCCGCCAATCAGAAgg                                                                                                                                                >  1:220786/1‑146 (MQ=255)
       ggCGACCGCTTTCACGCCTTGAATCCTAGGCGATCCGGTCATCGCGCTCACGCAGCAGCCAGCGGCCGTAATTATTCAAGAAGCATTAGCGGACGCTCGGGTTTTCCAAAGCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTaa                                                                                                                                      >  7:129809/1‑149 (MQ=255)
       ggCGACCGCTTTCACGCCTTGAATCCTAGGCGATCCGGTCATCGCGCTCACGCAGCAGCCAGCGGCCGTAATTATTCAAGAAGCATTAGCGGACGCTCGGGTTTTCCAAAGCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTaa                                                                                                                                      <  7:225299/149‑1 (MQ=255)
          gACCGCTTTCACGCCTTGAATCCTAGGCGATCCGGTCATCGCGCTCACGCAGCAGCCAGCGGCCGTAATTATTCAAGAAGCATTAGCGGACGCTCGGGTTTTCCAAAGCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTAATGc                                                                                                                                   >  8:256538/1‑149 (MQ=255)
          gACCGCTTTCACGCCTTGAATCCTAGGCGATCCGGTCATCGCGCTCACGCAGCAGCCAGCGGCCGTAATTATTCAAGAAGCATTAGCGGACGCTCGGGTTTTCCAAAGCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTAATGc                                                                                                                                   >  2:241992/1‑149 (MQ=255)
             cGCTTTCACGCCTTGAATCCTAGGCGATCCGGTCATCGCGCTCACGCAGCAGCCAGCGGCCGTAATTATTCAAGAAGCATTAGCGGACGCTCGTGTTTTCCAAAGCGAAAGCAACACGCCGCCAATCGGCCTGCTGATTggcg                                                                                                                                      <  8:270505/143‑4 (MQ=255)
             cGCTTTCACGCCTTGAATCCTAGGCGATCCGGTCATCGCGCTCACGCAGCAGCCAGCGGCCGTAATTATTCAAGAAGCATTAGCGGACGCTCGGGTTTTCCAAAGCGAAAGCAACACGCCGCCAATCGGCCTGCTGATTggcggaga                                                                                                                                  >  7:270505/1‑140 (MQ=255)
                tttCACGCCTTGAATCCTAGGCGATCCGGTCATCGCGCTCACGCAGCAGCCAGCGGCCGTAATTATTCAAGAAGCATTAGCGGACGCTCGGGTTTTCCAAAGCGAAAGCAACACGCCGCCAATCGGCCTGCTGATTggcg                                                                                                                                      >  1:130186/1‑137 (MQ=255)
                tttCACGCCTTGAATCCTAGGCGATCCGGTCATCGCGCTCACGCAGCAGCCAGCGGCCGTAATTATTCAAGAAGCATTAGCGGACGCTCGGGTTTTCCAAAGCGAAAGCAACACGCCGCCAATCGGCCTGCTGATTggcg                                                                                                                                      <  2:130186/140‑4 (MQ=255)
                  tCACGCCTTGAATCCTAGGCGATCCGGTCATCGCGCTCACGCAGCAGCCAGCGGCCGTAATTATTCAAGAAGCATTAGCGGACGCTCGGGTTTTCCAAAGCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGa                                                                                                                            >  8:133070/1‑148 (MQ=255)
                      gCCTTGAATCCTAGGCGATCCGGTCATCGCGCTCACGCAGCAGCCAGCGGCCGTAATTATTCAAGAAGCATTAGCGGACGCTCGGGTTTTCCAAAGCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTAAt                                                                                                                                     >  7:36333/1‑135 (MQ=255)
                      gCCTTGAATCCTAGGCGATCCGGTCATCGCGCTCACGCAGCAGCCAGCGGCCGTAATTATTCAAGAAGCATTAGCGGACGCTCGGGTTTTCCAAAGCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTAAt                                                                                                                                     <  8:36333/135‑1 (MQ=255)
                           gAATCCTAGGCGATCCGGTCATCGCGCTCACGCAGCAGCCAGCGGCCGTAATTATTCAAGAAGCATTAGCGGCCCCTCGGGTTTTCCAAAGCGAAAGCAACACGCCGCCATACAGCAGGCCGATTGTAATCATCAACGATATTTCgcat                                                                                                                  >  2:10146/1‑146 (MQ=255)
                           gAATCCTAGGCGATCCGGTCATCGCGCTCACGCAGCAGCCAGCGGCCGTAATTATTCAAGAAGCATTAGCGGACGCTCGGGTTTTCCAAAGCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCgg                                                                                                                    <  6:15221/147‑1 (MQ=255)
                             aTCCTAGGCGATCCGGTCATCGCGCTCACGCAGCAGCCAGCGGCCGTAATTATTCAAGAAGCATTAGCGGACGCTCGGGTTTTCCAAAGCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTg                                                                                                                 <  7:133070/148‑1 (MQ=255)
                                   ggCGATCCGGTCATCGCGCTCACGCAGCAGCCAGCGGCCGTAATTATTCAAGAAGCATTAGCGGACGCTCGGGTTTTCCAAAGCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACtt                                                                                                           <  3:55893/148‑1 (MQ=255)
                                      gATCCGGTCATCGCGCTCACGCAGCAGCCAGCGGCCGTAATTATTCAAGAAGCATTAGCGGACGCTCGGGTTTTCCAAAGCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTAcc                                                                                                       <  8:129809/149‑1 (MQ=255)
                                      gATCCGGTCATCGCGCTCACGCAGCAGCCAGCGGCCGTAATTATTCAAGAAGCATTAGCGGACGCTCGGGTTTTCCAAAGCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTAcc                                                                                                       <  3:144402/149‑1 (MQ=255)
                                            gTCATCGCGCTCACGCAGCAGCCAGCGGCCGTAATTATTCAAGAAGCATTAGCGGACGCTCGGGTTTTCCAAAGCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTACCCACAAt                                                                                                 >  6:283592/1‑149 (MQ=255)
                                                 cgcgCTCACGCAGCAGCCAGCGGCCGTAATTATTCAAGAAGCATTAGCGGACGCTCGGGTTTTCCAAAGCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTACCca                                                                                                     <  4:158371/140‑1 (MQ=255)
                                                 cgcgCTCACGCAGCAGCCAGCGGCCGTAATTATTCAAGAAGCATTAGCGGACGCTCGGGTTTTCCAAAGCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTACCca                                                                                                     >  3:158371/1‑140 (MQ=255)
                                                  gcgcTCACGCAGCAGCCAGCGGCCGTAATTATTCAAGAAGCATTAGCGGACGCTCGGGTTTTCCAAAGCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTACCCACAATGTCAt                                                                                            >  4:251052/1‑148 (MQ=255)
                                                    gcTCACGCAGCAGCCAGCGGCCGTAATTATTCAAGAAGCATTAGCGGACGCTCGGGTTTTCCAAAGCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTACCCACAATGTCATGaa                                                                                         <  3:251052/149‑1 (MQ=255)
                                                          gcagcaGCCAGCGGCCGTAATTATTCAAGAAGCATTAGCGGACGCTCGGGTTTTCCAAAGCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTACACACAATGTAATGAAGGaaaa                                                                                   >  8:268308/1‑149 (MQ=255)
                                                               agccagaGTCCGTAATTATTCAAGAACCATTAGCGGAAGCTCGGGTTTTCCAAAGCGAAAGCAACACGCCGCCAATCATCAGGCCGATTGTAATCCTCAACTATATTTCGGCTGGCACTTTACCCACAATGTCATGAAGGAAAATTTTc                                                                              <  2:220786/149‑1 (MQ=255)
                                                                 ccagcGGCCGTAATTATTCAAGAAGCATTAGCGGACGCTCGGGTTTTCCAAAGCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTGTACCCACAATGTCATGAAGGAAAATTTTc                                                                              >  6:35627/1‑147 (MQ=255)
                                                                    gcGGCCGTAATTATTCAAGAAGCATTAGCGGACGCTCGGGTTTTCCAAAGCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTACCCACAATGTCATGAAGGAAAATTTTCGAGcc                                                                         >  1:189589/1‑149 (MQ=255)
                                                                       gCCGTAATTATTCAAGAAGCATTAGCGGACGCTCGGGTTTTCCAAAGCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTACCCACAATGTCATGAAGGAAAATTTTCGAGCCGAt                                                                      <  5:31028/149‑1 (MQ=255)
                                                                       gCCGTAATTATTCAAGAAGCATTAGCGGACGCTCGGGTTTTCCAAAGCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTACCCACAATGTCATGAAGGAAAATTTTCGAGCCGAt                                                                      <  5:35755/149‑1 (MQ=255)
                                                                              ttattCAAGAAGCATTAGCGGACGCTCGGGTTTTCCAAAGCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTACCCACAATGTCATGAAGGAAAATTTTCGAGCCGATAAAGACa                                                               <  7:256538/149‑1 (MQ=255)
                                                                                    aagaagCATTAGCGGACGCTCGGGTTTTCCAAAGCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTACCCACAATGTCATGAAGGAAAATTTTCGAGCCGATAAAGACAAGCAcc                                                         <  2:61361/149‑1 (MQ=255)
                                                                                    aagaagCATTAGCGGACGCTCGGGTTTTCCAAAGCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTACCCACAATGTCATGAAGGAAAATTTTCGAGCCGATAAAGACAAGCAc                                                          >  3:275309/1‑148 (MQ=255)
                                                                                    aagaagCATTAGCGGACGCTAGGGTTTTCCAAAGCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTACCCACAATGTCATGAAGGAAAATTTTCGAGCCGATAAAGACAAGCAcc                                                         <  7:39615/149‑1 (MQ=255)
                                                                                            ttAGCGGACGCTCGGGTTTTCCAAAGCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTACCCACAATGTCATGAAGGAAAATTTTCGAGCCGATAAAGACAAGCACCATCGCCAg                                                 >  7:124818/1‑149 (MQ=255)
                                                                                                  gACGCTCGGGTTTTCCAAAGCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTACCCACAATGTCATGAAGGAAAATTTTCGAGCCGATAAAGACAAGCACCAGAGCCAGCGCata                                           <  1:10146/149‑1 (MQ=255)
                                                                                                   aCGCTCGGGTTTTCCAAAGCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTACCCACAATGTCATGAAGGAAAATTTTCGAGCCGATAAAGACAAGCACCAGAGCCAGCGCatat                                          <  4:275309/149‑1 (MQ=255)
                                                                                                   aCGCTCGGGTTTTCCAAAGCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTACCCACAATGTCATGAAGGAAAATTTTCGAGCCGATAAAGACAAGCACCAGAGCCAGCGCatat                                          >  6:250089/1‑149 (MQ=255)
                                                                                                       tCGGGTTTTCCAAAGCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTACCCACAATGTCATGAAGGAAAATTTTCGAGCCGATAAAGACAAGCACCAGAGCCAGCGCATAttt                                        >  5:293595/1‑147 (MQ=255)
                                                                                                       tCGGGTTTTCCAAAGCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTACCCACAATGTCATGAAGGAAAATTTTCGAGCCGATAAAGACAAGCACCAGAGCCAGCGAATAttt                                        >  8:129565/1‑147 (MQ=255)
                                                                                                         gggTTTTCCAAAGCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTACCCACAATGTCATGAAGGAAAATTTTCGAGCCGATAAAGACAAGCACCAGAGCCAGCGCATATTTAAgg                                    >  1:225319/1‑149 (MQ=255)
                                                                                                          ggTTTTCCAAAGCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTACCCACAATGTCATGAAGGAAAATTTTCGAGCCGATAAAGACAAGCACCAGAGCCAGCGCATATTTAAGGt                                   <  5:283592/149‑1 (MQ=255)
                                                                                                           gatttCCAAACCGAAAGCAACACCCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTACCCACAATGTCATGAAGGAAAATTTTCGAGCCGATAAAGACAAGCACCAGAGCCAGCGCATATTTAAGGTa                                  <  8:124818/147‑1 (MQ=255)
                                                                                                           gTTTTCCAAAGCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTACCCACAATGTCATGAAGGAAAATTTTCGAGCCGATAAAGACAAGCACCAGAGCCAGCGCATATTTAAGGTa                                  <  6:293595/149‑1 (MQ=255)
                                                                                                            ttttCCAAAGCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTACCCACAATGTCATGAAGGAAAATTTTCGAGCCGATAAAGACAAGCACCAGAGCCAGCGCATATTTAAGGTa                                  >  2:10479/1‑148 (MQ=255)
                                                                                                                  aaaGCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTACCCACAATGTCATGAAGGAAAATTTTCGAGCCGATAAAGACAAGCACCAGAGCCAGCGCATATTTAAGGTAAACAAAg                           >  2:232451/1‑149 (MQ=255)
                                                                                                                     gCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTACCCACAATGTCATGAAGGAAAATTTTCGAGCCGATAAAGACAAGCa                                                           <  8:197256/114‑1 (MQ=255)
                                                                                                                     gCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTACCCACAATGTCATGAAGGAAAATTTTCGAGCCGATAAAGACAAGCa                                                           >  7:197256/1‑114 (MQ=255)
                                                                                                                       gaaagaaaCAAGCCGCCAATCAGCAGGCCTATTGTAATGCTCAACGATATTTCGGCTGGCACTTTACCCACAATGTCATGAAGTAAAATTTTCGAGCCGATAAAGACAAGCACCAGAGCCAGCGCATATTTAAGGTAAACAAAGCGAc                       <  7:268308/148‑1 (MQ=255)
                                                                                                                       gAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTACCCACAATGTCATGAAGGAAAATTTTCGAGCCGATAAAGACAAGCACCAGAGCCAGCGCATATTTAAGGTAAACAAAGCGACt                      >  3:258913/1‑149 (MQ=255)
                                                                                                                          aGCAACACGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTACCCACAATGTCATGAAGGAAAATTTTCGAGCCGATAAAGACAAGCACCAGAGCCAGCGCATATTTAAGGTAAACAAAGCGACTGa                    >  3:261978/1‑148 (MQ=255)
                                                                                                                             aaCACGCCTCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGTCACTTTACCCACAATGTCATGAAGGAAAATTTTCGAGCCGATAAAGACAAGCACCAGAGCCAGCGCATATTTAAgg                                    <  5:75539/129‑1 (MQ=255)
                                                                                                                             aaCACGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTACCCACAATGTCATGAAGGAAAATTTTCGAGCCGATAAAGACAAGCACCAGAGCCAGCGCATATTTAAGGTAAACAAAGCGACTGATCAt                <  6:24108/149‑1 (MQ=255)
                                                                                                                             aaCACGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTACCCACAATGTCATGAAGGAAAATTTTCGAGCCGATAAAGACAAGCACCAGAGCCAGCGCATATTTAAGGTAAACAAAGCGACTGATCAt                >  8:82898/1‑149 (MQ=255)
                                                                                                                              acacGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTACCCACAATGTCATGAAGGAAAATTTTCGAGCCGATAAAGACAAGCACCAGAGCCAGCGCATATTTAAGGTAAACAAAGCGACTGATCATg               >  5:108888/1‑149 (MQ=255)
                                                                                                                               cacGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTACCCACAATGTCATGAAGGAAAATTTTCGAGCCGATAAAGACAAGCACCAGAGCCAGCGCATATTTAAGGTAAAc                               >  6:51486/1‑132 (MQ=255)
                                                                                                                               cacGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTACCCACAATGTCATGAAGGAAAATTTTCGAGCCGATAAAGACAAGCACCAGAGCCAGCGCATATTTAAGGTAAAc                               <  5:51486/132‑1 (MQ=255)
                                                                                                                                       cAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTACCCACAATGTCATGAAGGAAAATTTTCGAGCCGATAAAGACAAGCACCAGAGCCAGCGCATATTTAAGGTAAACAAAGCGACTGATCATGGCAGCCAAt      <  2:225319/149‑1 (MQ=255)
                                                                                                                                            agcagGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTACCCACAATGTCATGAAGGAAAATTTTCGAGCCGATAAAGACAAGCACCAGAGCCAGCGCATATTTAAGGTAAACAAAg                           <  6:15762/123‑1 (MQ=255)
                                                                                                                                            agcagGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTACCCACAATGTCATGAAGGAAAATTTTCGAGCCGATAAAGACAAGCACCAGAGCCAGCGCATATTTAAGGTAAACAAAGCGACTGATCATGGCAGCCAATGCaa  <  6:108888/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                               
CAGCGGAGGCGACCGCTTTCACGCCTTGAATCCTAGGCGATCCGGTCATCGCGCTCACGCAGCAGCCAGCGGCCGTAATTATTCAAGAAGCATTAGCGGACGCTCGGGTTTTCCAAAGCGAAAGCAACACGCCGCCAATCAGCAGGCCGATTGTAATGCTCAACGATATTTCGGCTGGCACTTTACCCACAATGTCATGAAGGAAAATTTTCGAGCCGATAAAGACAAGCACCAGAGCCAGCGCATATTTAAGGTAAACAAAGCGACTGATCATGGCAGCCAATGCAA  >  NC_002947/3740867‑3741154

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