Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,208,534 Δ1 bp coding (403/414 nt);coding (586/603 nt) tfaP →
tfaE ←
e14 prophage, uncharacterized protein
e14 prophage, putative tail fiber assembly protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,208,5340C.100.0% 262.1 / NA 73coding (403/414 nt)
coding (586/603 nt)
tfaP
tfaE
e14 prophage, uncharacterized protein
e14 prophage, putative tail fiber assembly protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (41/32);  total (41/32)

AGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCCTACGCCTCCGGCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCCTGCAATTTCAGTTGAATC  >  NC_000913/1208403‑1208664
                                                                                                                                   |                                                                                                                                  
aGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGa                                                                                                                            >  1:664877/1‑139 (MQ=255)
            gAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGgcg                                                                                                                <  1:641010/139‑1 (MQ=255)
              aGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGcgcg                                                                                                              <  2:155709/139‑1 (MQ=255)
                        aagaTGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTAt                                                                                                    <  2:664877/139‑1 (MQ=255)
                          gaTGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGgcg                                                                                                                >  1:469754/1‑125 (MQ=255)
                          gaTGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGgcg                                                                                                                <  2:469754/125‑1 (MQ=255)
                          gaTGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCt                                                                                                  >  2:1212/1‑139 (MQ=255)
                          gaTGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCACTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCt                                                                                                  >  1:891207/1‑139 (MQ=255)
                               gCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGc                                                                                             >  1:370516/1‑139 (MQ=255)
                                    aaCAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCa                                                                                        >  1:537871/1‑139 (MQ=255)
                                            aaTTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCa                                                                                <  1:1212/139‑1 (MQ=255)
                                             aTTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCaa                                                                               >  2:941430/1‑139 (MQ=255)
                                                 gtgGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGt                                                                                  <  2:541372/132‑1 (MQ=255)
                                                 gtgGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGt                                                                                  >  1:541372/1‑132 (MQ=255)
                                                 gtgGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTa                                                                           >  2:1303620/1‑139 (MQ=255)
                                                   ggCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAAt                                                                         >  2:43418/1‑139 (MQ=255)
                                                   ggCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAAt                                                                         >  1:792514/1‑139 (MQ=255)
                                                                   aGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAg                                                         >  2:424853/1‑139 (MQ=255)
                                                                    gCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGt                                                        <  1:424853/139‑1 (MQ=255)
                                                                    gCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGt                                                        <  1:1092160/139‑1 (MQ=255)
                                                                       cTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGcc                                                                                                                               <  2:1196990/65‑1 (MQ=17)
                                                                       cTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGcc                                                                                                                               >  1:1196990/1‑65 (MQ=17)
                                                                       cTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGgcgct                                                                                                              <  2:893951/82‑2 (MQ=17)
                                                                       cTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGcgcg                                                                                                              >  1:893951/1‑82 (MQ=17)
                                                                       cTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCa                                                                                        <  2:529161/104‑1 (MQ=255)
                                                                       cTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCa                                                                                        >  1:529161/1‑104 (MQ=255)
                                                                        tGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTgg                                                    >  2:1263162/1‑139 (MQ=255)
                                                                        tGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTgg                                                    >  2:518066/1‑139 (MQ=255)
                                                                                gTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGcgcg                                                                                                              >  1:1279114/1‑73 (MQ=17)
                                                                                gTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGcgcg                                                                                                              <  2:1279114/73‑1 (MQ=17)
                                                                                   gATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTctgcc                                         <  1:941430/139‑1 (MQ=255)
                                                                                        aTCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGt                                    <  1:518066/139‑1 (MQ=255)
                                                                                          ccAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCgg                                                                                                             <  2:296177/64‑1 (MQ=11)
                                                                                          ccAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCgg                                                                                                             >  1:296177/1‑64 (MQ=11)
                                                                                          ccAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCa                                  >  1:885854/1‑139 (MQ=255)
                                                                                          ccAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCa                                  >  1:1238020/1‑139 (MQ=255)
                                                                                           cAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTgg                                                    <  1:1253252/120‑1 (MQ=255)
                                                                                           cAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTgg                                                    >  2:1253252/1‑120 (MQ=255)
                                                                                             gTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGtt                               <  2:370516/139‑1 (MQ=255)
                                                                                              tGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACcgcg                                                                                   >  1:710785/1‑86 (MQ=255)
                                                                                              tGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACcgcg                                                                                   <  2:710785/86‑1 (MQ=255)
                                                                                               gcgcCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGcgcg                                                                                                              >  1:806803/1‑58 (MQ=11)
                                                                                               gcgcCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGcgcg                                                                                                              <  2:806803/58‑1 (MQ=11)
                                                                                                  ccAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGtt                                                            >  1:899485/1‑105 (MQ=255)
                                                                                                  ccAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGtt                                                            <  2:899485/105‑1 (MQ=255)
                                                                                                   cAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTc                         >  1:1065176/1‑139 (MQ=255)
                                                                                                          tGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAGGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCCTg                  >  2:900179/1‑139 (MQ=255)
                                                                                                            aaTGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGt                                                                                  >  1:1030465/1‑73 (MQ=11)
                                                                                                            aaTGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGt                                                                                  <  2:1030465/73‑1 (MQ=11)
                                                                                                            aaTGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCCTGCa                >  1:1211789/1‑139 (MQ=255)
                                                                                                               ggcttACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGt                                    >  1:1081001/5‑116 (MQ=255)
                                                                                                               ggcttACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGt                                    <  2:1081001/112‑1 (MQ=255)
                                                                                                                gcttACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGc                                                                                                          <  2:984659/42‑1 (MQ=14)
                                                                                                                gcttACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGc                                                                                                          >  1:984659/4‑45 (MQ=11)
                                                                                                                gcttACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGt                                                                                  >  1:1116771/4‑69 (MQ=11)
                                                                                                                gcttACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGt                                                                                  <  2:1116771/66‑1 (MQ=11)
                                                                                                                gcttACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATg                                                                             >  2:726541/4‑74 (MQ=11)
                                                                                                                gcttACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATg                                                                             <  1:726541/71‑1 (MQ=11)
                                                                                                                 cttACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGt                                    >  1:106490/3‑114 (MQ=255)
                                                                                                                 cttACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGt                                    <  2:106490/112‑1 (MQ=255)
                                                                                                                 cttACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCCTGCAATTTc           >  2:1003940/3‑139 (MQ=255)
                                                                                                                  ttACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGt                                    <  1:438125/112‑1 (MQ=255)
                                                                                                                  ttACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGt                                    >  2:438125/2‑113 (MQ=255)
                                                                                                                    aCGCCTCCGGCAGTTAGGCCAGATGACATCCGgcg                                                                                                                >  1:222265/1‑35 (MQ=14)
                                                                                                                    aCGCCTCCGGCAGTTAGGCCAGATGACATCCGgcg                                                                                                                <  2:222265/35‑1 (MQ=14)
                                                                                                                     cGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCAc                                                                                       >  2:22548/1‑59 (MQ=255)
                                                                                                                     cGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCAc                                                                                       <  1:22548/59‑1 (MQ=255)
                                                                                                                         tCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCCTGCAATTTCAGTTGAAt   <  1:1263162/139‑1 (MQ=255)
                                                                                                                          ccGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCCTGCAATTTCAGTTGAATc  >  2:958260/1‑139 (MQ=255)
                                                                                                                           cGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGt                                    >  1:93270/1‑104 (MQ=255)
                                                                                                                           cGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGt                                    <  2:93270/104‑1 (MQ=255)
                                                                                                                           cGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCCTGCAATTTCAGt        <  1:1417490/132‑1 (MQ=255)
                                                                                                                           cGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCCTGCAATTTCAGt        >  2:1417490/1‑132 (MQ=255)
                                                                                                                                   |                                                                                                                                  
AGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCCTACGCCTCCGGCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCCTGCAATTTCAGTTGAATC  >  NC_000913/1208403‑1208664

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