Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,340,588 T→C L105S (TTA→TCA)  yrbG → putative calcium/sodium:proton antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,340,5880TC100.0% 258.7 / NA 72L105S (TTA→TCA) yrbGputative calcium/sodium:proton antiporter
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (34/38);  total (34/38)

GCGTCTCTGCACGAACAACGCGATTTAGCCGTTGGTACAGCCCTCGGCTCAAACATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATTACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTGGCTGTGCTATGGCTGCTGTTCATTGTTAAACTTGCACGTC  >  NC_000913/3340452‑3340719
                                                                                                                                        |                                                                                                                                   
gcgTCTCTGCACGAACAACGCGATTTAGCCGTTGGTACAGCCCTCGGCTCAAACATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCAc                                                                                                                                   <  1:648584/139‑1 (MQ=255)
 cgTCTCTGCACGAACAACGCGATTTAGCCGTTGGTACAGCCCTCGGCTCAAACATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCAcc                                                                                                                                  <  1:1078316/139‑1 (MQ=255)
 cgTCTCTGCACGAACAACGCGATTTAGCCGTTGGTACAGCCCTCGGCTCAAACATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCAcc                                                                                                                                  <  2:597744/139‑1 (MQ=255)
              acaacGCGATTTAGCCGTTGGTACAGCCCTCGGCTCAAACATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGttgttg                                                                                                                     >  1:724219/1‑139 (MQ=255)
                        ttAGCCGTTGGTACAGCTCTCGGCTCAAACATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCgtgg                                                                                                           >  2:296401/1‑139 (MQ=255)
                              gTTGGTACAGCCCTCGGCTCAAACATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCg                                                                                                     >  1:851383/1‑139 (MQ=255)
                                   tACAGCCCTCGGCTCAAACATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAACGTGGTGGCCGGTTcc                                                                                                <  1:70067/139‑1 (MQ=255)
                                         ccTCGGCTCAAACATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTActc                                                                                          >  1:433704/1‑139 (MQ=255)
                                             ggCTCAAACATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCAcc                                                                                                                                  <  1:386718/95‑1 (MQ=255)
                                             ggCTCAAACATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCAcc                                                                                                                                  >  2:386718/1‑95 (MQ=255)
                                             ggCTCAAACATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGc                                                                                                       >  1:847990/1‑122 (MQ=255)
                                             ggCTCAAACATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTAct                                                                                           <  1:674182/134‑1 (MQ=255)
                                             ggCTCAAACATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTAct                                                                                           >  2:674182/1‑134 (MQ=255)
                                             ggCTCAAACATTATCAATATATGGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGc                                                                                                       <  2:847990/122‑1 (MQ=255)
                                                  aaaCATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCg                                                                                               >  1:233159/1‑125 (MQ=255)
                                                  aaaCATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCg                                                                                               <  2:233159/125‑1 (MQ=255)
                                                  aaaCATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTAt                                                                                       >  1:564505/1‑133 (MQ=255)
                                                  aaaCATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTAt                                                                                       <  2:564505/133‑1 (MQ=255)
                                                          tCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACgg                                                                                  <  1:274299/130‑1 (MQ=255)
                                                          tCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACgg                                                                                  >  2:274299/1‑130 (MQ=255)
                                                                  ttGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCgtg                                                                                                            <  2:596205/96‑1 (MQ=255)
                                                                  ttGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCgtg                                                                                                            >  1:596205/1‑96 (MQ=255)
                                                                  ttGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCg                                                                 <  2:293874/139‑1 (MQ=255)
                                                                   tGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGa                                                                <  2:851383/139‑1 (MQ=255)
                                                                    gCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCg                                                                      <  2:1101831/132‑1 (MQ=255)
                                                                    gCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCg                                                                      >  1:1101831/1‑132 (MQ=255)
                                                                        aTCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTc                                                                                                                  >  2:12079/1‑84 (MQ=255)
                                                                        aTCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTc                                                                                                                  <  1:12079/84‑1 (MQ=255)
                                                                        aTCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGTGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTa                                                           >  1:125154/1‑139 (MQ=255)
                                                                        aTCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCg                                                                                               >  2:521272/1‑103 (MQ=255)
                                                                        aTCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCg                                                                                               <  1:521272/103‑1 (MQ=255)
                                                                        aTCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTa                                                           >  2:865491/1‑139 (MQ=255)
                                                                             cGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCttt                                                      >  2:394833/1‑139 (MQ=255)
                                                                               gTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTtctc                                                   <  1:190751/139‑1 (MQ=255)
                                                                                tCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTtctc                                                   <  1:394833/139‑1 (MQ=255)
                                                                                    gCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGc                                                                                                       <  2:636551/83‑1 (MQ=255)
                                                                                    gCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGc                                                                                                       >  1:636551/1‑83 (MQ=255)
                                                                                    gCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCctct                                               >  1:782217/1‑139 (MQ=255)
                                                                                      tGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCttt                                             <  1:296401/139‑1 (MQ=255)
                                                                                       gcgcTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTc                                            <  2:928775/139‑1 (MQ=255)
                                                                                        cgcTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCa                                                                    >  2:948318/1‑114 (MQ=255)
                                                                                        cgcTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCa                                                                    <  1:948318/114‑1 (MQ=255)
                                                                                          cTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTgg                                         >  2:1057996/1‑139 (MQ=255)
                                                                                          cTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTgg                                         <  1:865491/139‑1 (MQ=255)
                                                                                           tGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTgtc                                        <  1:1057996/139‑3 (MQ=255)
                                                                                            ggTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGc                                                                                                       <  2:305500/75‑1 (MQ=255)
                                                                                            ggTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGc                                                                                                       >  1:305500/1‑75 (MQ=255)
                                                                                                  tCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGa                                                                                 <  2:27706/91‑1 (MQ=255)
                                                                                                  tCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGa                                                                                 >  1:27706/1‑91 (MQ=255)
                                                                                                   ccTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTGGCTGTGCTAt                                <  2:433704/139‑1 (MQ=255)
                                                                                                       ttACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGc                                                                                                       >  1:830431/1‑64 (MQ=255)
                                                                                                       ttACCGTCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGc                                                                                                       <  2:830431/64‑1 (MQ=255)
                                                                                                             tCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGc                                                                                                       >  1:905216/1‑58 (MQ=255)
                                                                                                             tCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGc                                                                                                       <  2:905216/58‑1 (MQ=255)
                                                                                                             tCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCg                                                                                                     <  2:301927/60‑1 (MQ=255)
                                                                                                             tCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCg                                                                                                     >  1:301927/1‑60 (MQ=255)
                                                                                                             tCCATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTGGCTGTGCTATGGCTGCTGtt                      >  2:515811/1‑139 (MQ=255)
                                                                                                               cATTCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTGGCTGTGCTATGGCTGCTGTTCa                    >  2:589108/1‑139 (MQ=255)
                                                                                                                  tCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGa                                                                >  1:686856/1‑92 (MQ=255)
                                                                                                                  tCTGATGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGa                                                                <  2:686856/92‑1 (MQ=255)
                                                                                                                      aTGTTCTACGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTGGCTGTGCTATGGCTGCTGTTCATTGTTaa             <  2:724219/139‑1 (MQ=255)
                                                                                                                              cGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGc                                                                                                       <  2:980318/41‑1 (MQ=255)
                                                                                                                              cGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGc                                                                                                       >  1:980318/1‑41 (MQ=255)
                                                                                                                              cGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTAct                                                                                           >  2:534373/1‑53 (MQ=255)
                                                                                                                              cGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTAct                                                                                           <  1:534373/53‑1 (MQ=255)
                                                                                                                              cGCCGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTGGCTGTGCTATGGCTGCTGTTCATTGTTAAACTTGCAc     >  1:252418/1‑139 (MQ=255)
                                                                                                                                ccGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTGGCTGTGCTATGGCTGCTGTTCATTGTTAAACTTGCACGt   <  1:589108/139‑1 (MQ=255)
                                                                                                                                 cGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGc                                                                  <  2:796866/75‑1 (MQ=255)
                                                                                                                                 cGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGc                                                                  >  1:796866/1‑75 (MQ=255)
                                                                                                                                 cGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTa                                                           >  2:1022948/1‑82 (MQ=255)
                                                                                                                                 cGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTa                                                           <  1:1022948/82‑1 (MQ=255)
                                                                                                                                 cGTGAATCACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTGGCTGTGCTATGGCTGCTGTTCATTGTTAAACTTGCACGTc  <  2:95023/139‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                   
GCGTCTCTGCACGAACAACGCGATTTAGCCGTTGGTACAGCCCTCGGCTCAAACATTATCAATATATTGCTGATCCTCGGTCTGGCTGCGCTGGTTCGTCCTTTTACCGTCCATTCTGATGTTCTACGCCGTGAATTACCCTTAATGTTGTTGGTCAGCGTGGTGGCCGGTTCCGTACTCTATGACGGACAACTTAGTCGCAGCGATGGTATCTTTCTCCTCTTTCTGGCTGTGCTATGGCTGCTGTTCATTGTTAAACTTGCACGTC  >  NC_000913/3340452‑3340719

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