Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,065,588 A→G N263D (AAT→GAT)  scpC → propionyl‑CoA:succinate CoA transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,065,5880AG100.0% 202.0 / NA 57N263D (AAT→GAT) scpCpropionyl‑CoA:succinate CoA transferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (30/27);  total (30/27)

AACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCAATATCAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGTTTATGATGTATTCGGAAGTGCTACAGGAATCGGTGGTGCATTTACTGGAAACCGGCAAAATCAGCGGGGCCAGCGCCTCCAGCCTGACAATCT  >  NC_000913/3065440‑3065735
                                                                                                                                                    |                                                                                                                                                   
aaCCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCGat                                                                                                                                                   >  2:249412/1‑151 (MQ=255)
   cAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCGATATc                                                                                                                                                <  1:1003753/151‑1 (MQ=255)
        tGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCGATATCaata                                                                                                                                            >  1:918065/1‑150 (MQ=255)
         gCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCGATATCaataat                                                                                                                                          <  1:438466/151‑1 (MQ=255)
           ccGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCGATATCAATAATGc                                                                                                                                        <  1:617660/151‑1 (MQ=255)
            cGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCGATATCAATAATGCg                                                                                                                                       >  1:503855/1‑151 (MQ=255)
            agACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCGATATCAATAATGCg                                                                                                                                       >  2:677268/2‑151 (MQ=255)
             gACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCGATATCAATAATGCgg                                                                                                                                      <  1:1047899/151‑1 (MQ=255)
               cGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGATGCAAAGTGGCGTGGGCGATATCAATAATGCGGTa                                                                                                                                    <  1:558950/151‑1 (MQ=255)
                         atGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCGATATCAATAATGCGGTAATg                                                                                                                                 >  2:880639/1‑144 (MQ=255)
                         atGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCGATATCAATAATGCGGTAATGGCGCGTc                                                                                                                          <  2:310245/151‑1 (MQ=255)
                         atGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCGATATCAATAATGCGGTAATGGCGCGTc                                                                                                                          >  2:712631/1‑151 (MQ=255)
                           gCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCGATATCAATAATGCGGTAATGGCGCGTCt                                                                                                                         <  1:507176/150‑1 (MQ=255)
                            cTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCGATATCAATAATGCGGTAATGGCGCGTCTgg                                                                                                                       >  2:950718/1‑151 (MQ=255)
                                  aaGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCGATATCAATAATGCGGTAATGGCGCGTCTGGGGGaaa                                                                                                                 <  1:825768/151‑1 (MQ=255)
                                  aaGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCGATATCAATAATGCGGTAATGGCGCGGCTgggggca                                                                                                                  >  2:524221/1‑148 (MQ=255)
                                    gCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCGATATCAATAATGCGGTAATGGCGCGTCTGGGGGAAAAc                                                                                                               >  1:126019/1‑151 (MQ=255)
                                      aaaaTCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCGATATCAATAATGCGGTAATGGCGCGTCTGGGGGAAAAccc                                                                                                             <  2:496655/151‑1 (MQ=255)
                                        aaTCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCGATATCAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCgg                                                                                                           <  2:771514/151‑1 (MQ=255)
                                        aaTCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCGATATCAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCgg                                                                                                           <  1:637940/151‑1 (MQ=255)
                                           cccATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCGATATCAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTaa                                                                                                        <  2:120729/151‑1 (MQ=255)
                                             cATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCGATATCAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAAtt                                                                                                      >  1:823011/1‑151 (MQ=255)
                                              aTGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCGATATCAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATtc                                                                                                     >  1:143977/1‑151 (MQ=255)
                                                  gCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCGATATCAATAATGCGGTAATGGCGCGGCTGGGGGAAAACCCGGTAATtcctcc                                                                                                 >  2:595518/1‑151 (MQ=255)
                                                    cagcagATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCGATATCAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGt                                                                                               >  1:1023320/1‑151 (MQ=255)
                                                     agcagATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCGATATCAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGtt                                                                                              >  2:432600/1‑151 (MQ=255)
                                                     agcagATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCGATATCAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGtt                                                                                              >  2:540876/1‑151 (MQ=255)
                                                       cagATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCGATATCAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGTTTa                                                                                            <  2:250244/151‑1 (MQ=255)
                                                        agATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCATGGGCGATATCAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGTTTat                                                                                           >  2:30629/1‑151 (MQ=255)
                                                          aTTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCGATATCAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGTTTatg                                                                                          >  1:125965/1‑150 (MQ=255)
                                                          aTTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGCGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCGATATCAATAATGCGGTAATGGCGCGTCTGGGGGAAAAACCGGTAATTCCTCCGtt                                                                                              >  2:13818/1‑146 (MQ=255)
                                                              ccGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCGATATCAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGTTTATGATGt                                                                                      >  1:1016480/1‑150 (MQ=255)
                                                                 attaCGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCGATATCCATAATGCGGTAATGCCGCCGCTGGGGGAAAACCCGGTAATTCCTCCGCTTATGATGTATTc                                                                                  >  2:1052122/4‑151 (MQ=255)
                                                                        ggTCGCGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCGATATCAATAATGCGATAATGTCGCGTCTGGGAGAAAACCCGGTAATTCCTCCGTTTATGATGTATTCGGAAGt                                                                            <  2:425081/150‑1 (MQ=255)
                                                                        ggTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCGATATCAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGTTTATGATGTATTCGGAAGTg                                                                           >  1:111070/1‑151 (MQ=255)
                                                                          tCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCTGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCGATATCAATAATGCGGTAATGGCGCGTTCGGGGGAAAACCCGGCAATTCCTCCGTTTTTGATGTATTCGGCGGTGCt                                                                         >  1:131112/1‑151 (MQ=255)
                                                                              gTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCGATATCAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGTTTATGATGTATTCGGAAGTGCTACAg                                                                     <  1:76581/151‑1 (MQ=255)
                                                                               ttCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGTTGCAAAGTGGCGTGGGCGCTCTCAATAATGCGGTAAAGGCGCGTCTGGGGGAAAAACCGGTAATTCCTCCGTTTATCTTGTATTCGGAAGGGCGAc                                                                       >  2:177736/1‑148 (MQ=255)
                                                                                  ttATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCGATATCAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGGAATTCCTCCGGTTATGATGTATTCGGAAGTGCTACAGGaa                                                                  >  2:836147/1‑150 (MQ=255)
                                                                                   tATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCGATATCAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGTTTATGATGTATTCGGAAGTGCTACAGGAAt                                                                 <  1:1063408/150‑1 (MQ=255)
                                                                                   tATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCGATATCAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGTTTATGATGTATTCGGAAGTGCTACAGGAAt                                                                 <  1:907642/150‑1 (MQ=255)
                                                                                   tATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCGATATCAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGTTTATGATGTATTCGGAAGTGCTACAGGAATc                                                                >  2:829466/1‑151 (MQ=255)
                                                                                         aGGACATGGCGCATGGGCGTATTCCGCCGGAATTTCTGACGCTGCAAAGTGGCGTGGGCGGTATCAATAATGCGGTAATGGTGCGTCCGGGGGAAAACCCGGTGATTCCTCCGGTTATGATGTATTCGGAAGTGGTACAGGAAGCggtgg                                                           >  2:42084/1‑150 (MQ=255)
                                                                                          ggAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCGATATCAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGTTTATGATGTATTCGGAAGTGCTACAGGAATCGgtggtg                                                         <  1:742536/151‑1 (MQ=255)
                                                                                          ggAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCGATATCAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGTTTATGATGTATTCGGAAGTGCTACAGGAATCGgtggtg                                                         <  1:184008/151‑1 (MQ=255)
                                                                                           gAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCGATATCAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGTTTATGATGTATTCGGAAGTGCTACAGGAATCGGTGGTGc                                                        <  2:187070/151‑1 (MQ=255)
                                                                                                 gcgcATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCGATATCAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGTTTATGATGTATTCGGAAGTGCTACAGGAATCGGTGGTGCATTTAc                                                  <  1:238443/151‑1 (MQ=255)
                                                                                                   gcATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCGATATCAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGTTTATGATGTATTCGGAAGTGCTACAGGAATCGGTGGTGCATTTACTg                                                <  2:1014662/151‑1 (MQ=255)
                                                                                                               tCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCGATATCAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGTTTATGATGTATTCGGAAGTGCTACAGGAATCGGTGGTGCATTTACTGGAAACCGGCaa                                     >  1:939482/1‑150 (MQ=255)
                                                                                                                   ccgGAATTTCTGCCGCTGCAAAGTGGCGTGGGCGATATCAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGTTTATGATGTATTCGGAAGTGCTACAGGAATCGGTGGTGCATTTACTGGAAACCGGCAAAATCa                                <  2:71717/151‑1 (MQ=255)
                                                                                                                                    gCAAAGTGGCGTGGGCGATATCAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGTTTATGATGTATTCGGAAGTGCTACAGGAATCGGTGGTGCATTTACTGGAAACCGGCAAAATCAGCGGGGCCAGCGCCt                 <  2:337150/149‑1 (MQ=255)
                                                                                                                                     cAAAGTGGCGTGGGCGATATCAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGTTTATGATGTATTCGGAAGTGCTACAGGAATCGGTGGTGCATTTACTGGAAACCGGCAAAATCAGCGGGGCCAGCGCCTCCa              <  2:706978/151‑1 (MQ=255)
                                                                                                                                       aaGTGGCGTGGGCGATATCAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGTTTATGATGTATTCGGAAGTGCTACAGGAATCGGTGGTGCATTTACTGGAAACCGGCAAAATCAGCGGGGCCAGCGCCTCCAGc            >  1:469871/1‑151 (MQ=255)
                                                                                                                                         gTGGCGTGGGCGATATCAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGTTTATGATGTATTCGGAAGTGCTACAGGAATCGGTGGTGCATTTACTGGAAACCGGCAAAATCAGCGGGGCCAGCGCCTCCAGCCt          >  1:561409/1‑151 (MQ=255)
                                                                                                                                                gggCGATATCAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGTTTATGATGTATTCGGAAGTGCTACAGGAATCGGTGGTGCATTTACTGGAAACCGGCAAAATCAGCGGGGCCAGCGCCTCCAGCCTGACAAtc   >  1:332037/1‑151 (MQ=255)
                                                                                                                                                 ggCGATATCAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGTTTATGATGTATTCGGAAGTGCTACAGGAATCGGTGGTGCATTTACTGGAAACCGGCAAAATCAGCGGGGCCAGCGCCTCCAGCCTGACAAtct  <  1:939916/151‑1 (MQ=255)
                                                                                                                                                 ggCGATATCAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGTTTATGATGTATTCGGAAGTGCTACAGGAATCGGTGGTGCATTTACTGGAAACCGGCAAAATCAGCGGGGCCAGCGCCTCCAGCCTGACAAtca  <  2:33024/151‑2 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AACCAACTTGCCCGACGCCGGTAATATGCTGGATAAGCAAAATCCCATGTGCCAGCAGATTGCCGATAACGTGGTCACGTTCTTATTGCAGGAAATGGCGCATGGGCGTATTCCGCCGGAATTTCTGCCGCTGCAAAGTGGCGTGGGCAATATCAATAATGCGGTAATGGCGCGTCTGGGGGAAAACCCGGTAATTCCTCCGTTTATGATGTATTCGGAAGTGCTACAGGAATCGGTGGTGCATTTACTGGAAACCGGCAAAATCAGCGGGGCCAGCGCCTCCAGCCTGACAATCT  >  NC_000913/3065440‑3065735

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