Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,901,435 A→C 7.1% L51V (TTA→GTA)  yoaE ← putative membrane protein/conserved protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,901,4350AC 7.1% 149.9 / 3.7 56L51V (TTA→GTA) yoaEputative membrane protein/conserved protein
Reads supporting (aligned to +/- strand):  ref base A (35/17);  new base C (2/2);  total (37/19)
Fisher's exact test for biased strand distribution p-value = 5.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.52e-01

AACAGCAAGAATATCCCCCCGAACAACATAATCAGGTCGCGTCCGGAGAAGGAGAAATCCATGACGGTAAATAGCGGTTTGGTCAGCGTGACCATCCATGAAATCAGCGACAGCAGCCCCAGACGCATAATCAGCGCCAGTGATAACCCCAGCAAACGCGCTTTATCGCGTTGTTTTGGCGGCAGTTTGTCAGCAAGAATGGCGATGAAGACCAGGT‑TATCGATACCCAGCACAATTTCGAGAACAACAAGCGTGAGTAGCCCCGCCCAAATTGAGGGGTCCATTAAGAATTC  >  NC_000913/1901290‑1901582
                                                                                                                                                 |                                                                                                                                                    
aaCAGCAAGAATATCCCCCCGAACAACATAATCAGGTCGCGTCCGGAGAAGGAGAAATCCATGACGGTAAATAGCGGTTTGGTCAGCGTGACCATCCATGAAATCAGCGACAGCAGCCCCAGACGCATAATCAGCGCCAGTGATAAccc                                                                                                                                                   >  2:1143945/1‑149 (MQ=255)
    gCAAGAATATCCCCCCGAACAACATAATCAGGTCGCGTCCGGAGAAGGAGAAATCCATGACGGTAAATAGCGGTTTGGTCAGCGTGACCATCCATGAAATCAGCGACAGCAGCCCCAGACGCATAATCAGCGCCAGTGATAACCCCAGc                                                                                                                                               <  1:250938/149‑1 (MQ=255)
     cAAGAATATCCCCCCGAACAACATAATCAGGTCGCGTCCGGAGAAGGAGAAATCCATGACGGTAAATAGCGGTTTGGTCAGCGTGACCATCCATGAAATCAGCGACAGCAGCCCCAGACGCATAATCAGCGCCAGTGATAACCCCAGCa                                                                                                                                              >  2:744328/1‑149 (MQ=255)
        gAATATCCCCCCGAACAACATAATCAGGTCGCGTCCGGAGAAGGAGAAATCCATGACGGTAAATAGCGGTTTGGTCAGCGTGACCATCCATGAAATCAGCGACAGCAGCCCCAGACGCATAATCAGCGCCAGTGATAACCCCAGCAAAc                                                                                                                                           >  2:877595/1‑149 (MQ=255)
                 cccGAACAACATAATCAGGTCGCGTCCGGAGAAGGAGAAATCCATGACGGTAAATAGCGGTTTGGTCAGCGTGACCATCCATGAAATCAGCGACAGGAGCCCCAGACGCATAATAAGCGCCAGTGATAACCCCAGCAAACACGCtttct                                                                                                                                  >  2:7488/1‑147 (MQ=255)
                 cccGAACAACATAATCAGGTCGCGTCCGGAGAAGGAGAAATCCATGACGGTAAATAGCGGTTTGGTCAGCGTGACCATCCATGAAATCAGCGACAGCAGCCCCAGACGCATAATCAGCGCCAGTGATAACCCCAGCAAACGCGCTTTAt                                                                                                                                  >  2:330136/1‑149 (MQ=255)
                        aacaTAATCAGGTCGCGTCCGGAGAAGGAGAAATCCATGACGGTAAATAGCGGTTTGGTCAGCGTGACCATCCATGAAATCAGCGACAGCAGCCCCAGACGCATAATCAGCGCCAGTGATAACCCCAGCAAACGCGCTTTATCGCggtg                                                                                                                           >  1:540424/1‑149 (MQ=255)
                        aacaTAATCAGGTCGCGTCCGGAGAAGGAGAAATCCATGACGGTAAATAGCGGTTTGGTCAGCGTGACCATCCATGAAATCAGCGACAGCAGCCCCAGACGAGTAATCAGCGCCAGTGATAACCCCAGCAAACGCGCTTTATCGCgtcg                                                                                                                           >  2:377503/1‑147 (MQ=255)
                                  ggTCGCGTCCGGAGAAGGAGAAATCCATGACGGTAAATAGCGGTTTGGTCAGCGTGACCATCCATGAAATCAGCGACAGCAGCCCCAGACGCATAATCAGCGCCAGTGATACCCCCAGCAAACGCGCTTTATCGCGTTGTTTTggcggc                                                                                                                 >  1:717431/1‑149 (MQ=255)
                                  ggTCGCGTCCGGAGAAGGAGAAATCCATGACGGTAAATAGCGGTTTGGTCAGCGTGACCATCCATGAAATCAGCGACAGCAGCCCCAGACGCATAATCAGCGCCAGTGATACCCCCAGCAAACGCGCTTTATCGCGTTGTTTTggcggc                                                                                                                 >  1:717386/1‑149 (MQ=255)
                                   ggctcgTCCGGAGAAGGAGAAATCCATGACGGTAAATAGCGGTTTGGTCAGCGTGACCATCCATGAAATCAGCGACAGCAGCCCCAGACGCATAATCAGCGCCAGTGATAACCCCAGCAAACGCGCTTTATCGCGTTGTTTTGGCGGCa                                                                                                                <  2:1149041/145‑1 (MQ=255)
                                   gTCGCGTCCGGAGAAGGAGAAATCCATGACGGTAAATAGCGGTTTGGTCAGCGTGCCCATCCATGAAATCAGCGACAGCAGCCCCACACGCATAATCAGCGCCAGTGCTAACCCCAGTACAGGCGATTTTTTGCGTTGTTTTggcggcc                                                                                                                >  1:209373/1‑148 (MQ=255)
                                    tCGCGTCCGGAGAAGGAGAAATCCATGACGGTAAATAGCGGTTTGGTCAGCGTGACCATCCATGAAATCAGCGACAGCAGCCCCAGACGCATAATCAGCGCCAGTGATAACCCCAGCAAACGCGCTTTATCGCGTTGTTTTGGCGGCAg                                                                                                               >  1:589006/1‑149 (MQ=255)
                                      gcgTCCGGAGAAGGAGAAATCCATGACGGTAAATAGCGGTTTGGTCAGCGTGACCATCCATGAAATCAGCGACAGCAGCCCCAGACGCATAATCAGCGCCAGTGATAACCCCAGCAAACGCGCTTTATCGCGTTGTTTTGGCGGCAGtt                                                                                                             >  2:310107/1‑149 (MQ=255)
                                         tCCGGAGAAGGTGAAATCCATGACGGTAAATAGCGGTTTGGTCAGCGTGACCATCCATGAAATCAGCGACAGCAGCCCCAGACGCATACTCAGCGTCAGTGATAACCCCAGCAAACGCGCTTTATCGCGTTGTTTTGGCGGCAGTTTGt                                                                                                          <  2:5114/149‑1 (MQ=255)
                                            ggagaaggagaaATCCATGACGGTAAATAGCGGTTTGGTCAGCGTGACCATCCATGAAATCAGCGACAGCAGCCCCAGACGCATAATCAGCGCCAGTGATAACCCCAGCAAACGCGCTTTATCGCGTTGTTTTGGCGGCAGTTTGTCAg                                                                                                       >  2:1087332/1‑149 (MQ=255)
                                             gagaaggagaaATCGATGACGGTAAATAGCGGTTTGGTCAGCGTGACCATCCATGAAATCAGCGACAGCAGCCCCAGACGCATAATCAGCGCCAGTGATAACCCCAGCAAACGCGCTTTATCGCGTTGTTTTGGCGGCAGTTTGTCAGc                                                                                                      <  1:681037/149‑1 (MQ=255)
                                             gagaaggagaaATCCATGACGGTAAATAGCGGTTTGGTCAGCGTGACCATCCATGAAATCAGCGACAGCAGCCCCAGACGCATAATCAGCGCCAGTGATAACCCCAGCAAACGCGCTTTATCGCGTTGTTTTGGCGGCACTTTGTCAGc                                                                                                      >  2:155798/1‑149 (MQ=255)
                                             gagaaggagaaATCCATGACGGTAAATAGCGGTTTGGTCAGCGTGACCATCCATGAAATCAGCGACAGCAGCCCCAGACGCATAATCAGCGCCAGTGATAACCCCAGCAAACGCGCTTTATCGCGTTGTTTTGGCAGCAGATTGTAAGc                                                                                                      >  2:1039692/1‑149 (MQ=255)
                                             gagaaggagaaATCCATGACGGGAAATAGCGGTTTGGTCGGCGTGAGCATCAATGAAATCAGCGACCGGAGCGCCAGACGCAGAATAAGCGACAGTGCTAACCCCCGCAAACGCGCTTAATCGAGTTGTTTTGGCGGCAGTTTGTCAGc                                                                                                      >  2:421459/1‑149 (MQ=255)
                                                  ggagaaATCCATGACGGTAAATAGCGGTTTGGTCAGCGTGACCATCCATGAAATCAGCGACAGCAGCCCCAGACGCATAATCAGCGCCAGTGATAACCCCAGCAAACGCGCTTTATCGCGTTGTTTTGGCGGCAGTTTGTCAGCAAGaa                                                                                                 >  1:121975/1‑149 (MQ=255)
                                                     gaaACCCATGACGGTAAATAGCGGTTTGGTCAGCGTGACCATCCATGAAATCAGCGACAGCAGCCCCAGACGCATAATCAGCGCCAGTGATAACCCCAGCAAACGCGCTTTATCGCGTTGTTTTGGCGGCAGTTTGTCAGCAAGAATgg                                                                                              <  2:1071711/149‑1 (MQ=255)
                                                      aaaTCCATGACGGTAAATAGCGGTTTGGTCGGCGTGACCATCCATGAAATCAGCGACAGCAGCCCCAGACGCATAATCAGCGCCAGTGATAACCCCAGCAAACGCGCTTTATCGCGTTGTTTTGGCGGCAGTTTGTCAGCAAGAATGGc                                                                                             >  1:598484/1‑149 (MQ=255)
                                                      aaaTCCATGACGGTAAATAGCGGTTTGGTCAGCGTGACCATCCATGAAATCAGCGACAGCAGCCCCAGACGCATAATCAGCGCCAGTGATAACCCCAGCAAACGCGCTTTATCGCGTTGTTTTGGCGGCAGTTTGTCAGCAAGAATGGc                                                                                             >  1:183537/1‑149 (MQ=255)
                                                         tCCATGACGGTAAATAGCGGTTTGGTCAGCGTGACCATCCATGAAATCAGCGACAGCAGCCCCAGACGCATAATCAGCGCCAGTGATACCCCCAGCAAACGCGCTTTATCGCGTTGTTTTGGCGGCAGTTTGTCAGCAAGAATGGCGAt                                                                                          <  2:717431/149‑1 (MQ=255)
                                                         tCCATGACGGTAAATAGCGGTTTGGTCAGCGTGACCATCCATGAAATCAGCGACAGCAGCCCCAGACGCATAATCAGCGCCAGTGATACCCCCAGCAAACGCGCTTTATCGCGCTGTTTTGGCGGCAGTTTGTCAGCAAGAATGGCGAt                                                                                          <  2:717386/149‑1 (MQ=255)
                                                           cATGACGGTAAATAGCGGTTTGGTCAGCGTGACCATCCATGAAATCAGCGACAGCAGCCCCAGACGCATAATCAGCGCCAGTGAGAACCCCAGCAAACGCGCTTTATCGCGTTGTTTTGGCGGAAGTTTGTCAGCAAGAATGGCGAGGa                                                                                        >  1:971958/1‑149 (MQ=255)
                                                            aTGACGGTAAATAGCGGTTTGGTCAGCGTGACCAGCCATGAAATCAGCGACAGTAGCCCCAGACGCATAATCAGCGCCAGTGATAACCCCAGCAAACGCGCTTTATCGCGTTGTTTTGGCGGCAGTTTGTCAGCGAGAATGGCGATGaa                                                                                       >  1:634742/1‑149 (MQ=255)
                                                                cGGGAAATAGCGGTTTGGTCAGCGTGACCATCCATGAAATCAGCGACAGCAGCCCCAGACGCATAATCAGCGCCAGTGATAACCCCAGCAAACGCGCTTTATCGCGTTGTTTTGGCGGCAGTTTGTCAGCAAGAATGGCGATGAAGAcc                                                                                   <  1:1005383/149‑1 (MQ=255)
                                                                  gTAAATAGCGGGTTGGTCAGCGTGACCATCCATGAAATCAGGGACAGCAGCCCCAGACGCATAATCAGCGCTAGTGATAACCCCAGCAAACGCGGTTTATCGCGTTGTTTTGGCGGCAGTTTGTCAGCAAGAATGGCGATGAAGACCAg                                                                                 <  2:904205/149‑1 (MQ=255)
                                                                   tAAATGGCGGTTTGGTCAGCGTGACCATCCATGAAATCAGCGACAGCAGCCCCAGACGCATAATCAGCGCCAGTGATAACCCCAGCAAACGCGCTTTATCGCGTTGTTTTGGCGGCAGTTTGTCAGCAAGAATGGCGATGAAGACCAgg                                                                                <  2:206373/149‑1 (MQ=255)
                                                                       tAGCGGTTTGGTCAGCGTGACCATCCATGAAATCAGCGACAGCAGCCCCAGACGCATAATCAGCGCCAGTGATAACCCCAGCAAACGCGCTTTATCGCGTTGTTTTGGCGGCAGTCTGTCAGCAAGAATGGCGATGAAGACCAGGT‑TAt                                                                           <  1:739808/149‑1 (MQ=255)
                                                                        aGCGGTTTGGTCAGCGTGACCATCCATGAAATCAGCGGCAGCAGCCCCAGACGCATAATCAGCGCCAGTGATAACCCCAGCAAACGCGCTTTATCGCGTTGTTTTGGCGGCAGTTTGTCAGCAAGAATGGCGATGAAGACCAGGT‑TATc                                                                          <  1:957631/149‑1 (MQ=255)
                                                                                ggTCAGCGTGACCATCCATGAAATCAGCGACAGCAGCCCCAGACGCATAATCAGCGCCAGTGATAACCCCAGCAAACGCGCTTTATCGCGTTGTTTTGGCGGCAGTTTGTCAGCAAGAATGGCGATGAAGACCAGGT‑TATCGATACCCa                                                                  <  1:641297/149‑1 (MQ=255)
                                                                                       gTGACCATCCATGAAATCAGCGACAGCAGCCCCAGACGCATAATCAGCGCCAGTGATAACCCCAGCAAACGCGCTTTATCGCGTTGTTTTGGCGGCAGTTTGTCAGCAAGAATGGCGATGAAGACCACGG‑TATAGATACCCAGCACaac                                                           >  1:397703/1‑148 (MQ=255)
                                                                                        tGACCATCCATGAAATCAGCGACAGCAGCCCCAGACGCATAATCAGCGCCAGTGATAACCCCAGCAAACGCGCTTTATCGCGTTGTTTTGGCGGCAGTTTGTCAGCAAGAATGGCGATGAAGACCAGGT‑TATCGATACCCAGCACAAtt                                                          >  2:1061760/1‑149 (MQ=255)
                                                                                          aCCATCCATGAAATCAGCGACAGCAGCCCCAGACGCATAATCAGCGCCAGTGATAACCCCAGCAAACGCGCTTTATCGCGTTGTTTTGGCGGCAGTTTGTCAGCAAGAATGGCGATGAAGACCAGGT‑TCTCGATACCCAGCACAATTTc                                                        >  2:146359/1‑149 (MQ=255)
                                                                                            catccatGAAATCAGCGACAGCAGCCCCAGACGCATAATCAGCGCCAGTGATAACCCCAGCAAACGCGCTTTATCGCGTTGTTTTGGCGGCAGTTTGTCAGCAAGAATGGCGATGAAGACCAGGT‑TATCGATACCCAGCACAATTTCga                                                      >  2:295690/1‑149 (MQ=255)
                                                                                            catccatGAAATCAGCGACAGCAGCCCCAGACGCATAATCAGCGCCAGTGATAACCCCAGCAAACGCGCTTTATCGCGTTGTTTTGGCGGCAGTTTGTCAGCAAGAATGGCGATGAAGACCAGGT‑TATCGATACCCAGCACAATTTCga                                                      <  1:1002348/149‑1 (MQ=255)
                                                                                            cacccatGAAATCCGCGACAGCCGCCCCAGACGCATAATCAGCGCCAGTGATAACCCCAGCAAACGCGCTTTATCGCGTTGTTTTGGCGGCAGTTTGTCAGCAAGAATGGCGATGAAGACCAGGT‑TATCGATACCCAGCACAATTTCga                                                      <  2:236258/146‑1 (MQ=255)
                                                                                              tccatGAAATCAGCGACAGCAGCCCCAGACGCATAATCAGCGCCAGTGATAACCCCAGCAAACGCGCTTTATCGCGTTGTTTTGGCGGCAGTTTGTCAGCAAGAATGGCGATGAAGACCAGGT‑TATCGATACCCAGCACAATTTCgaga                                                    >  2:3957/1‑149 (MQ=255)
                                                                                              tccatGAAATCAGCGACAGCAGCCCCAGACGCATAATCAGCGCCAGGGATAACCCCAGCAAACGCGCTTTATCGCGTTGTTTTGGCGGCAGTTTGTCAGCAAGAATGGCGATGAAGACCAGGT‑TATCGATACCCAGCACAATTTCgaga                                                    <  1:330136/149‑1 (MQ=255)
                                                                                                catGAAATCAGCGACAGCAGCCCCAGACGCATAATCAGCGCCAGTGATAACCCCAGCAAACGCGCTTTATCGCGTTGTTTTGGCGGCAGTTTGTCAGCAAGAATGGCGATGAAGACCAGGT‑TATCGATACCCAGCACAATTTCGAGaac                                                  >  2:961389/1‑149 (MQ=255)
                                                                                                   gAAATCAACGACAGCAGCCCCAGACGCATAATCAGCGCCAGTGATAACCCCAGCAAACGCGCTTTATCGCGTTGCTTTGGCGGCAGTTTGTCAGCAAGAATGGCGATGAAGACCAGGT‑TATCGATACCCAGCACCATATCAAGaaccac                                               >  2:943455/1‑149 (MQ=255)
                                                                                                         aGCGACAGCCGCCCCAGACCCATAATCAGCGCCAGTGATAACCCCAGCAAACGCGCTTTATCGCGTTGTTTTGGCGGCAGTTTGTCAGCAAGAATGGCGATGAAGACCAGGT‑TATCGATACCCAGCACAATTTCGAGAACAACAAGCGt                                         <  1:545332/149‑1 (MQ=255)
                                                                                                          gCGACAGCAGCCCCAGACGCATAATCAGCGCCAGTGATAACCCCAGCAAACGCGCTTTATCGCGTTGTTTTGGCGGCAGTTTGTCAGCAAGAATGGCGATGAAGACCAGGT‑TATCGATACCCAGCAAAATTTCGAGAACAACAAGCGTg                                        >  2:467425/1‑149 (MQ=255)
                                                                                                                    ccccAGACGCATAATCAGCGCCAGTGATAACCCCAGCAAACGCGCTTTATCGCGTTGTTTTGGCGGCAGTTTGTCAGCAAGAATGGCGATGAAGACCAGGT‑TATCGATACCCAGCACAATTTCGAGAACAACAAGCGGCAGTAGcccag                              >  2:97511/1‑147 (MQ=255)
                                                                                                                     cccAGACGCATAACCAGCGCCAGTGATAACCCCAGCAAACGCGCTTTATCGCGTTGTTTTGGCGGCAGTTTGTCAGCAAGAATGGCGATGAAGACCAGGT‑TATCGATACCCAGCACAATTTCGAGAACAACAAGCGTGAGTAGCCCCGc                             >  1:786159/1‑149 (MQ=255)
                                                                                                                     cccAGACGCATAACCAGCGCCAGTGATAACCCCAGCAAACGCGCTTTATCGCGTTGTTTTGGCGGCAGTTTGTCAGCAAGAATGGCGATGAAGACCAGGT‑TATCGATACCCAGCACAATTTCGAGAACAACAAGCGTGAGTAGCCCCGc                             >  1:786398/1‑149 (MQ=255)
                                                                                                                          aCGCATAATCAGCGCCAGTGATAACCCCAGCAAACGCGCTTTATCGCGTTGTTTTGGCGGCAGTTTGTCAGCAAGAATGGCGATGAAGACCAGGT‑TATCGATACCCAGCACAATTTCGAGAACAACAAGCGTGAGTAGCCCCGCCCaaa                        <  1:3957/149‑1 (MQ=255)
                                                                                                                            gCATAATCAGCGCCAGTGATAACCCCAGCAAACGCGCTTTATCGCGTTGTTTTGACGGCAGTTTGTCAGCAAGAATGGCGATAAAGA‑CAGGTATCTCGGGAGTTAGGAGGAGGTTTTGAACAACAAGCGAGAGTAGCCCCGCCCAACtt                      >  2:529430/1‑149 (MQ=255)
                                                                                                                               tAATCAGCGCCAGTGATAACCCCAGCAAACGCGCTTTATCGCGTTGTTTTGGCGGCAGTTTGTCAGCAAGAATGGCGATGAAGACCAGGT‑TATCGATACCCAGCACAATTTCGAGAACAACAAGCGAGAGTAGCCCCGCCCAAATTGAg                   >  2:525090/1‑149 (MQ=255)
                                                                                                                                 aTCAGCGCCATGGATAACCCCAGCAAACGCGCTTTATCGCGTTGTTTTGGCGGCAGTTTGTCAGCAAGAATGGCGATGAAGACCAGGT‑TATCGATACCCAGCACAATTTCGAGAACAACAAGCGTGAGTAGCCCCGCCCAAATTGAggg                 <  1:1039692/149‑1 (MQ=255)
                                                                                                                                      cgcCAGTGATAACCCCAGCAAACGCGCTTTATCGCGTTGTTTTGGCGGCAGTTTGTCAGCAAGAATGGCGATGAAGACCAGGT‑TATCGATACCCAGCACAATTTCGAGAACAACAAGCGCGAGTAGCCCCGCCCAAATTGAGGGGACCa            >  2:651582/1‑149 (MQ=255)
                                                                                                                                           gTGATAACCCCAGCAAACGCGCTTTATCGCGTTGTTTTGGCGGCAGTTTGTCAGCAAGAATGGCGATGAAGACCAGGT‑TATCGATACCCAGCACAATTTCGAGAACAACAAGCGTGAGTAGCCCCGCCCAAATTGAGGGGTCCATTAAg       >  1:899094/1‑149 (MQ=255)
                                                                                                                                               tAACCCCAGAAAACGCGCTTAATCGCGTGGTTTTGGCGGCAGTTGGTCAGCAAAAATGGCGATGAAGACCAGGT‑TAGCGATACCCAGCACAACGTCGAGAACAAATAGCGACGATAGCACCGCCGAGTTAGAGGGGTCCATGAAgcaaa   >  2:525374/1‑145 (MQ=255)
                                                                                                                                                aaCCCCAGCAAACGCGCTTTATCGCGTTGTTTTGGCGGCAGTTTGTCAGCAAGAATGGCGATGAAGACCAGGT‑TATCGATACCCAGCACCATTTCGAGAACAACAAGCGTGAGTAGCCCCGCCCAAATTGAGGGGTCCATTAAGAATTc  <  2:123030/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
AACAGCAAGAATATCCCCCCGAACAACATAATCAGGTCGCGTCCGGAGAAGGAGAAATCCATGACGGTAAATAGCGGTTTGGTCAGCGTGACCATCCATGAAATCAGCGACAGCAGCCCCAGACGCATAATCAGCGCCAGTGATAACCCCAGCAAACGCGCTTTATCGCGTTGTTTTGGCGGCAGTTTGTCAGCAAGAATGGCGATGAAGACCAGGT‑TATCGATACCCAGCACAATTTCGAGAACAACAAGCGTGAGTAGCCCCGCCCAAATTGAGGGGTCCATTAAGAATTC  >  NC_000913/1901290‑1901582

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