Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,637,154 A→C 6.9% V86G (GTC→GGC)  yhiN ← putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,637,1540AC 6.9% 189.3 / 3.6 56V86G (GTC→GGC) yhiNputative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (24/28);  new base C (2/2);  total (26/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

GCCACACTCAGCACTTCGCTACGCAATCTGAAGGTCACATTGCCCTTCTCGCACTCATCCACCAGCATGTCGACAATCTGCTGCGCGGAGTCATCGCAGAAGAGTTGCCCTAACGTTTTCTCGTGCCAGGCGATGCCGTGTTTATTGACCAGATCAATGAAATCCCACTGGGTAAAACGTGCGAGTGCAGACTTACAAAAATGCGGATTCTGGCTCAGATAAGCGCCTGGTTCGACATAAAGGTTGGTAAAGTTGCAGCGCCCACCGCCAGACATAAGGATTTTGCGCCCTGGTTT  >  NC_000913/3637007‑3637302
                                                                                                                                                   |                                                                                                                                                    
gCCACACTCAGCACTTCGCTACGCAATCTGCAGGTCACATTGCCCTTCTCGCACTCATCCACCAGCATGTCGACAATCTGCTGCGTGGAGTCATCGCAGAAGAGTTGCCCTAACGTTTTCTCGTGCCAGGCGATGCCGTGTTTATTGAc                                                                                                                                                     <  4:12173/149‑1 (MQ=255)
gCCACACTCAGCACTTCGCTACGCAATCTGAAGGTCACATTGCCCTTCTCGCACTCATCCACCAGCATGTCGACAATCTGCTGCGCGGAGTCATCGCAGAAGAGTTGCCCTAACGTTTTCTCGTGCCAGGCGATGCCGTGTTTATTGAc                                                                                                                                                     <  4:273418/149‑1 (MQ=255)
  cacaCTCAGCACTTCGCTACGCAATCTGAAGGTCACATTGCCCTTCTCGCACTCATCCACCAGCATGTCGACAATCTGCTGCGCGGAGTCATCGCAGAAGAGTTGCCCTAACGTTTTCTCGTGCCAGGCGATGCCGTGTTTATTGACCa                                                                                                                                                   <  1:214006/149‑1 (MQ=255)
  cacaCTCAGCACTTCGCTACGCAATCTGAAGGTCACATTGCCCTTCTCGCACTCATCCACCAGCATGTCGACAATCTGCTGCGCGGAGTCATCGCAGAAGAGTTGCCCTAACGTTTTCTCGTGCCAGGCGATGCCGTGTTTATTGACCa                                                                                                                                                   >  2:379898/1‑149 (MQ=255)
   acacTCAGCACTTCGCTACGCAATCTGAAGGTCACATTGCCCTTCTCGCACTCATCCACCAGCATGTCGACAATCTGCTGCGCGGAGTCATCGCAGAAGAGTTGCCCTAACGTTTTCTCGTGCCAGGCGATGCCGTGTTTATTGACCAg                                                                                                                                                  >  1:395504/1‑149 (MQ=255)
       tCAGCACTTCGCTACGCAATCTGAAGGTCACATTGCCCTTCTCGCACTCATCCACCAGCATGTCGACAATCTGCTGCGCGGAGTCATCGCAGAAGAGTTGCCCTAACGTTTTCTCGTGCCAGGCGATGCCGTGGTTATTGACCAGATCa                                                                                                                                              >  4:249067/1‑149 (MQ=255)
             cTTCGCTACGCAATCTGAAGGTCACATTGCCCTTCTCGCACTCATCCACCAGCATGTCGACAATCTGCTGCGCGGAGTCATCGCAGAAGAGTTGCCCTAACGTTTTCTCGTGCCAGGCGATGCCGTGTTTATTGACCAGATCAATGaaa                                                                                                                                        <  1:4944/149‑1 (MQ=255)
              ttCGCTACGCAATCTGAAGGTCACATTGCCCTTCTCGCACTCATCCACCAGCATGTCGACAATCTGCTGCGCGGAGTCATCGCAGAAGAGTTGCCCTAACGTTTTCTCGTGCCAGGCGATGCCGTGTTTATTGa                                                                                                                                                      >  3:356255/1‑134 (MQ=255)
              ttCGCTACGCAATCTGAAGGTCACATTGCCCTTCTCGCACTCATCCACCAGCATGTCGACAATCTGCTGCGCGGAGTCATCGCAGAAGAGTTGCCCTAACGTTTTCTCGTGCCAGGCGATGCCGTGTTTATTGa                                                                                                                                                      <  4:356255/134‑1 (MQ=255)
                    aCGCAATCTGAAGGTCACATTGCCCTTCTCGCACTCATCCACCAGCATGTCGACAATCTGCTGCGCGGAGTCATCGCAGAAGAGTTGCCCTAACGTTTTCTCGTGCCAGGCGATGCCGTGTTTATTGACCAGATCAATGAAATCCCACt                                                                                                                                 <  2:1036/149‑1 (MQ=255)
                             gAAGGTCACATTGCCCTTCTCGCACTCATCCACCAGCATGTCGACAATCTGCTGCGCGGAGTCATCGCAGAAGAGTTGCCCTAACGTTTTCTCGTGCCAGGCGATGCCGTGTTTATTGACCAGATCAATGAAATCCCACTGGGTAAAAc                                                                                                                        <  1:379898/149‑1 (MQ=255)
                              aaGGTCACATTGCCCTTCTCGCACTCATCCACCAGCATGTCGACAATCTGCTGCGCGGAGTCATCGCAGAAGAGTTGCCCTAACGTTTTCTCGTGCCAGGCGATGCCGTGTTTATTGACCAGATCAATGAAATCCCACTGGGTAAAACg                                                                                                                       >  4:274784/1‑149 (MQ=255)
                                ggTCACATTGCCCTTCTCGCACTCATCCACCAGCATGTCGACAATCTGCTGCGCGGAGTCATCGCAGAAGAGTTGCCCTAACGTTTTCTCGTGCCAGGCGATGCCGTGTTTATTGACCAGATCAATGAAATCCCACTGGGTAAAACGTg                                                                                                                     <  1:219537/149‑1 (MQ=255)
                                 gTCACATTGCCCTTCTCGCACTCATCCACCAGCATGTCGACAATCTGCTGCGCGGAGTCATCGCAGAAGAGTTGCCCTAACGTTTTCTCGTGCCAGGCGATGCCGTGTTTATTGCCCAGATCAATGAAATCCCACTGGGTAAAACGTGc                                                                                                                    >  1:222405/1‑149 (MQ=255)
                                 gTCACATTGCCCTTCTCGCACTCATCCACCAGCATGTCGACAATCTGCTGCGCGGAGTCATCGCAGAAGAGTTGCCCTAACGTTTTCTCGTGCCAGGCGATGCCGTGTTTATTGACCAGATCAAt                                                                                                                                            >  1:15946/1‑125 (MQ=255)
                                 gTCACATTGCCCTTCTCGCACTCATCCACCAGCATGTCGACAATCTGCTGCGCGGAGTCATCGCAGAAGAGTTGCCCTAACGTTTTCTCGTGCCAGGCGATGCCGTGTTTATTGACCAGATCAAt                                                                                                                                            <  2:15946/125‑1 (MQ=255)
                                 gTCACATTGCCCTTCTCGCACTCATCCACCAGCATGTCGACAATCTGCTGCGCGGAGTCATCGCAGAAGAGTTGCCCTAACGTTTTATCGTGCCAGGCGATGCCGTGTTTATTGCCCAGATCAATGAAATCCCACTGGGTAAAACGTGc                                                                                                                    >  3:167208/1‑149 (MQ=255)
                                  tCACATTGCCCTTCTCGCACTCATCCACCAGCATGTCGACAATCTGCTGCGCGGAGTCATCGCAGAAGAGTTGCCCTAACGTTTTCTCGTGCCAGGCGATGCCGTGTTTATTGACCAGATCAATGAAATCCCACTGGGTAAAACGTGCg                                                                                                                   >  3:365388/1‑149 (MQ=255)
                                      aTTGCCCTTCTCGCACTCATCCACCAGCATGTCGACAATCTGCTGCGCGGAGTCATCGCAGAAGAGTTGCCCTAACGTTTTCTCGTGCCAGGCGATGCCGTGTTTATTGACCAGATCAATGAAATCCCACTGGGTAAAACGTGCGAGTg                                                                                                               >  2:7013/1‑149 (MQ=255)
                                          cccTTCTCGCACTCATCCACCAGCATGTCGACAATCTGCTGCGCGGAGTCATCGCAGAAGAGTTGCCCTAACGTTTTCTCGTGCCAGGCGATGCCGTGTTTATTGACCAGATCAATGAAATCCCACTGGGTAAAACGTGCGAGTGCAGa                                                                                                           <  1:58226/149‑1 (MQ=255)
                                           ccatctcGCACTCATCCACCAGCATGTCGACAATCTGCTGCGCGGAGTCATCGCAGAAGAGTTGCCCTAACGTTTTCTCGTGCCAGGCGATGCCGTGTTTATTGACCAGATCAATGAAATCCCACTGGGTAAAACGTGCGAGTGCAGAc                                                                                                          <  4:34333/146‑1 (MQ=255)
                                           ccTTCTCGCACTCATCCACCAGCATGTCGACAATCTGCTGCGCGGAGTCATCGCAGAAGAGTTGCCCTAACGTTTTCTCGTGCCAGGCGATGCCGTGTTTATTGACCAGATCAATGAAATCCCACTGGGTAAAACGTGc                                                                                                                    >  4:337526/1‑139 (MQ=255)
                                           ccTTCTCGCACTCATCCACCAGCATGTCGACAATCTGCTGCGCGGAGTCATCGCAGAAGAGTTGCCCTAACGTTTTCTCGTGCCAGGCGATGCCGTGTTTATTGACCAGATCAATGAAATCCCACTGGGTAAAACGTGc                                                                                                                    <  3:337526/139‑1 (MQ=255)
                                            cTTCTCGCACTCATCCACCAGCATGTCGACAATCTGCTGCGCGGAGTCATCGCAGAAGAGTTGCCCTAACGTTTTCTCGTGCCAGGCGATGCCGTGTTTATTGACCAGATCAATGAAATCCCACTGGGTAAAACGTGCGAGTGCAGACt                                                                                                         >  4:248795/1‑149 (MQ=255)
                                             ttCTCGCACTCATCCACCAGCATGTCGACAATCTGCTGCGCGGAGTCATCGCAGAAGAGTTGCCCTAACGTTTTCTCGTGCCAGGCGATGCCGTGTTTATTGACCAGATCAATGAAATCCCACTGGGTAAAACGTGCGAGTGCAGACtt                                                                                                        >  2:188662/1‑149 (MQ=255)
                                                tcGCACTCATCCACCAGCATGTCGACAATCTGCTGCGCGGAGTCATCGCAGAAGAGTTGCCCTAACGTTTTCTCGTGCCAGGCGATGCCGTGTTTATTGACCAGATCAATGAAATCCCACTGGGTAAAACGTGCGAGTGCAGACTTACa                                                                                                     <  1:202273/149‑1 (MQ=255)
                                                    aCTCATCCACCAGCATTTCGACAATCTGCTGCGCGGAGTCATCGCAGAAGAGTTGCCCTAACGTTTTCTCGTGCCAGGCGATGCCGTGTTTATTGACCAGATCAATGAAATCCCACTGGGTAAAACGTGCGAGTGCAGACTTACaaaaa                                                                                                 <  2:331897/149‑1 (MQ=255)
                                                       cATCCACCAGCATGTCGACAATCTGCTGCGCGGAGTCATCGCAGAAGAGTTGCCCTAACGTTTTCTCGTGCCAGGCGATGCCGTGTTTATTGACCAGATCAATGAAATCCCACTGGGTAAAACGTGCGAGTGCAGACTTACAAAAATGc                                                                                              <  2:99773/149‑1 (MQ=255)
                                                        aTCCACCAGCATGTCGACAATCTGCTGCGCGGAGTCATCGCAGAAGAGTTGCCCTAACGTTTTCTCGTGCCAGGCGATGCCGTGTTTATTGACCAGATCAATGAAATCCCACTGGGTAAAACGTGCGAGTGCAGACTTACAAAAATGCg                                                                                             >  2:123400/1‑149 (MQ=255)
                                                          ccaccaGCATGTCGACAATCTGCTGCGCGGAGTCATCGCAGAAGAGTTGCCCTAACGTTTTCTCGTGCCAGGCGATGCCGTGTTTATTGCCCAGATCAATGAAATCCCACTGGGTAAAACGTGCGAGTGCAGACTTACAAAAATGCGGa                                                                                           <  4:167208/149‑1 (MQ=255)
                                                          ccaccaGCATGTCGACAATCTGCTGCGCGGAGTCATCGCAGAAGAGTTGCCCTAACGTTTTCTCGTGCCAGGCGATGCCGTGTTTATTGCCCAGATCAATGAAATCCCACTGGGTAAAACGTGCGAGTGCAGACTTACAAAAATGCGGa                                                                                           <  2:222405/149‑1 (MQ=255)
                                                              caGCATGTCGACAATCTGCTGCGCGGAGTCATCGCAGAAGAGTTGCCCTAACGTTTTCTCGTGCCAGGCGATGCCGTGTTTATTGACCAGATCAATGAAATCCCACTGGGTAAAACGTGCGAGTGCAGACTTACAAAAATGCGGATTCt                                                                                       <  4:245739/149‑1 (MQ=255)
                                                                gCATGTCGACAACCTGCTGCGCGGAGTCATCGCAGAAGAGTTGCCCTAACGTTTTCTCGTGCCAGGCGATGCCGTGTTTGTTGACCAGATCAATGAAATCCCACTGGGTAAAACGTGCGAGTCCAGACTTACAAAAATGCGGATTCTgg                                                                                     <  1:256820/149‑1 (MQ=255)
                                                                     tCGACAATCTGCTGCGCGGAGTCATCGCAGAAGAGTTGCCCTAACGTTTTCTCGTGCCAGGCGATGCCGTGTTTATTGACCAGATCAATGAAATCCCACTgg                                                                                                                               >  3:373349/1‑102 (MQ=255)
                                                                     tCGACAATCTGCTGCGCGGAGTCATCGCAGAAGAGTTGCCCTAACGTTTTCTCGTGCCAGGCGATGCCGTGTTTATTGACCAGATCAATGAAATCCCACTgg                                                                                                                               <  4:373349/102‑1 (MQ=255)
                                                                       gACAATCTGCTGCGCGGAGTCATCGCAGAAGAGTTGCCCTAACGTTTTCTCGTGCCAGGCGATGCCGTGTTTATTGACCAGATCAATGAAATCCCACTGGGTAAAACGTGCGAGTGCAGACTTACAAAAATGCGGATTCTGGCTCAGAt                                                                              >  2:435654/1‑149 (MQ=255)
                                                                           aTCTGCTGCGCGGAGTCATCGCAGAAGAGTTGCCCTAACGTTTTCTCGTGCCAGGCGATGCCGTGTTTATTGACCAGATCAATGAAATCCCACTGGGTAAAACGTGCGAGTGCAGACTTACAAAAATGCGGATTCTGGCTCAGATAAgc                                                                          <  1:188662/149‑1 (MQ=255)
                                                                              tgctgcGCGGAGTCATCGCAGAAGAGTTGCCCTAACGTTTTCTCGTGCCAGGCGATGCCGTGTTTATTGACCAGATCAATGAAATCCCACTGGGTAAAACGTGCGAGTGCAGACTTACAAAAATGCGGATTCTGGCTCAGATAAGCGcc                                                                       <  3:248795/149‑1 (MQ=255)
                                                                                ctgcGCGGAGTCATCGCAGAAGAGTTGCCCTAACGTTTTCTCGTGCCAGGCGATGCCGTGTTTATTGACCAGATCAAt                                                                                                                                            >  3:107419/1‑78 (MQ=255)
                                                                                ctgcGCGGAGTCATCGCAGAAGAGTTGCCCTAACGTTTTCTCGTGCCAGGCGATGCCGTGTTTATTGACCAGATCAAt                                                                                                                                            <  4:107419/78‑1 (MQ=255)
                                                                                ctgcGCGGAGTCATCGCAGAAGAGTTGCCCTAACGTTTTCTCGTGCCAGGCGATGCCGTGTTTATTGACCAGATCAATGAAATCCCACTGGGTAAAACGTGCGAGTGCAGACTTACAAAAATGCGGATTCTGGCTCAGATAAGCGCCTg                                                                     <  1:435654/149‑1 (MQ=255)
                                                                                  gcgcGGAGTCATCGCAGAAGAGTTGCCCTAACGTTTTCTCGTGCCAGGCGATGCCGTGTTTATTGACCAGATCAATGAAATCCCACTGGGTAAAACGTGCGAGTGCAGACTTACAAAAATGCGGATTCTGGCTCAGATAAGCGCCTGGt                                                                   >  3:133960/1‑149 (MQ=255)
                                                                                      ggAGTCATCGCAGAAGAGTTGCCCTAACGTTTTCTCGTGCCAGGCGATGCCGTGTTTATTGACCAGATCAATGAAATCCCACTGGGTAAAACGTGCGAGTGCAGACTTACAAAAATGCGGATTCTGGCTCAGATAAGCGCCTGGTTCGa                                                               <  3:7752/149‑1 (MQ=255)
                                                                                          tCATCGCAGAAGAGTTGCCCTAACGTTTTCTCGTGCCAGGCGATGCCGTGTTTATTGACCAGATCAATGAAATCCCACTGGGTAAAACGTGCGAGTGCAGACTTACAAAAATGCGGATTCTGGCTCAGATAAGCGCCTGGTTCGACATa                                                           >  1:326255/1‑149 (MQ=255)
                                                                                             tCGCAGAAGAGTTGCCCTAACGTTTTCTCGTGCCAGGCGATGCCGTGTTTATTGACCAGATCAATGAAATCCCACTGGGTAAAACGTGCGAGTGCAGACTTACAAAAATGCGGATTCTGGCTCAGATAAGCGCCTGGTTCGACATAAAg                                                        <  3:165905/149‑1 (MQ=255)
                                                                                                         tGCCCTAACGTTTTCTCGTGCCAGGCGATGCCGTGTTTATTGACCAGATCAATGAAATCCCACTGGGTAAAACGTGCGAGTGCAGACTTACAAAAATGCGGATTCTGGCTCAGATAAGCGCCTGGTTCGACATAAAGGTTGGTAAAGtt                                            <  3:468133/149‑1 (MQ=255)
                                                                                                          gCCCTAACGTTTTCTCGTGCCAGGCGATGCCGTGTTTATTGACCAGATCAATGAAATCCCACTGGGTAAAACGTGCGAGTGCAGACTTACAAAAATGCGGATTCTGGCTCAGATAAGCGCCTGGTTCGACATAAAGGTTGGTAAAGTTg                                           <  2:186883/149‑1 (MQ=255)
                                                                                                          gCCCTAACGTTTTCTCGTGCCAGGCGATGCCGTGTTTATTGACCAGATCAATGAAATCCCACTGGGTAAAACGTGCGAGTGCAGACTTACAAAAATGCGGATTCTGGCTCAGATAAGCGCCTGGTTCGACATAAAGGTTGGTAAAGTTg                                           <  4:146048/149‑1 (MQ=255)
                                                                                                           cccTAACGTTTTCTCGTGCCAGGCGATGCCGTGTTTATTGACCAGATCAATGAAATCCCACTGGGTAAAACGTGCGAGTGCAGACTTACAAAAATGCGGATTCTGGCTCAGATAAGCGCCTGGTTCGACATAAAGGTTGGTAAAGTTGc                                          >  3:123876/1‑149 (MQ=255)
                                                                                                            ccTAACGTTTTCTCGTGCCAGGCGATGCCGTGTTTATTGACCAGATCAATGAAATCCCACTGGGTAAAACGTGCGAGTGCAGACTTACAAAAATGCGGATTCTGGCTCAGATAAGCGCCTGGTTCGACATAAAGGTTGGTAAAGTTGCa                                         >  3:180545/1‑149 (MQ=255)
                                                                                                              tAACGTTTTCTCGTGCCAGGCGATGCCGTGTTTATTGACCAGATCAATGAAATCCCACTGGGTAAAACGTGCGAGTGCAGACTTACAAAAATGCGGATTCTCGCTCAGATAAGCGCCCGGTTCGACATAAAGGTTGGTAAAGTTGCAgc                                       >  2:223349/1‑149 (MQ=255)
                                                                                                                    tttCTCGTGCCAGGCGATGCCGTGTTTATTGACCAGATCAATGAAATCCCACTGGGTAAAACGTGCGAGTGCAGACTTACAAAAATGCGGATTCTGGCTCAGATAAGCGCCTGGTTCGACATAAAGGTTGGTAAAGTTGCAGCGCCCAc                                 >  4:102215/1‑149 (MQ=255)
                                                                                                                       ctcGTGCCAGGCGATGCCGTGTTTATTGACCAGATCAATGAAATCCCACTGGGTAAAACGTGCGAGTGCAGACTTACAAAAATGCGGATTCTGGCTCAGATAAGCGCCTGGTTCGACATAAAg                                                        >  2:465260/1‑123 (MQ=255)
                                                                                                                       ctcGTGCCAGGCGATGCCGTGTTTATTGACCAGATCAATGAAATCCCACTGGGTAAAACGTGCGAGTGCAGACTTACAAAAATGCGGATTCTGGCTCAGATAAGCGCCTGGTTCGACATAAAg                                                        <  1:465260/123‑1 (MQ=255)
                                                                                                                        tcGTGCCAGGCGATGCCGTGTTTATTGACCAGATCAATGAAATCCCACTGGGTAAAACGTGCGAGTGCAGACTTACAAAAATGCGGATTCTGGCTCAGATAAGCGCCTGGTTCGACATAAAGGTTGGTAAAGTTGCAGCGCCCACCGcc                             <  2:326255/149‑1 (MQ=255)
                                                                                                                               aGGCGATGCCGTGTTTATTGACCAGATCAATGAAATCCCACTGGGTAAAACGTGCGAGTGCAGACTTACAAAAATGCGGATTCTGGCTCAGATAAGCGCCTGGTTCGACATAAAGGTTGGTAAAGTTGCAGCGCCCACCGCCAGACATa                      >  1:36271/1‑149 (MQ=255)
                                                                                                                                          tgtTTATTGACCAGATCAATGAAATCCCACTGGGTAAAACGTGCGAGTGCAGACTTACAAAAATGCGGATTCTGGCTCAGATAAGCGCCTGGTTCGACATAAAGGTTGGTAAAGTTGCAGCGCCCACCGCCAGACATAAGGATTTTgag           >  4:256516/1‑147 (MQ=255)
                                                                                                                                           gtTTATTGACCAGATCAATGAAATCCCACTGGGTAAAACGTGCGAGTGCAGACTTACAAAAATGCGGATTCTGGCTCAGATAAGCGCCTGGTTCGACATAAAGGTTGGTAAAGTTGCAGCGCCCACCGCCAGACATAAGGATTTTgcgc          >  3:394431/1‑149 (MQ=255)
                                                                                                                                                   aCCAGGTCAATGAAATCCCACTGGGTAAAACGTGCGAGTGCAGACTTACAAAAATGCGGATTCTGGCTCAGATAAGCGCCTGGTTCGACATAAAGGTTGGTAAAGTTGCAGCGCCCACCGCCAGACATAAGGATTTTGCGCCCTGGttt  >  1:365277/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GCCACACTCAGCACTTCGCTACGCAATCTGAAGGTCACATTGCCCTTCTCGCACTCATCCACCAGCATGTCGACAATCTGCTGCGCGGAGTCATCGCAGAAGAGTTGCCCTAACGTTTTCTCGTGCCAGGCGATGCCGTGTTTATTGACCAGATCAATGAAATCCCACTGGGTAAAACGTGCGAGTGCAGACTTACAAAAATGCGGATTCTGGCTCAGATAAGCGCCTGGTTCGACATAAAGGTTGGTAAAGTTGCAGCGCCCACCGCCAGACATAAGGATTTTGCGCCCTGGTTT  >  NC_000913/3637007‑3637302

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