Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,077,514 Δ1 bp intergenic (‑165/+149) cmtB ← / ← tktA PTS system mannitol‑specific enzyme IIA component/transketolase 1 isozyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,077,5130C.100.0% 775.5 / NA 169intergenic (‑164/+150)cmtB/tktAPTS system mannitol‑specific enzyme IIA component/transketolase 1 isozyme
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (112/57);  total (112/57)

Only 100 of 175 total aligned reads displayed.

GCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACCGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCG  >  NC_000913/3077421‑3077600
                                                                                            |                                                                                       
gCAATTGAGCAGGGAAGCCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATAc                                                                                         <  2:1523791/93‑1 (MQ=255)
  aaTTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATAcgc                                                                                       >  2:3248173/1‑91 (MQ=255)
   aTTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATAcgca                                                                                      <  2:781841/93‑4 (MQ=255)
   aTTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATAcgca                                                                                      >  1:5000689/1‑90 (MQ=255)
    ttGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATAcgcat                                                                                     >  1:7045540/1‑89 (MQ=255)
     tGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCata                                                                                   <  2:5194387/93‑1 (MQ=255)
       aGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATc                                                                                 >  2:3452071/1‑93 (MQ=255)
               aaCCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCa                                                                         <  1:7784593/93‑1 (MQ=255)
                 ccGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCata                                                                       >  1:6700505/1‑93 (MQ=255)
                 ccGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCata                                                                       <  2:3098843/93‑1 (MQ=255)
                 ccGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCata                                                                       <  1:7637983/93‑1 (MQ=255)
                 ccGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATATCAATACGCATATCGTGATGCata                                                                       <  1:1112420/93‑1 (MQ=255)
                  cGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCatat                                                                      <  1:6330898/93‑1 (MQ=255)
                      gtTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACa                                                                  <  2:7547432/93‑1 (MQ=255)
                       ttttGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACAt                                                                 <  2:5534179/93‑1 (MQ=255)
                        tttGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATa                                                                >  1:5633258/1‑93 (MQ=255)
                        tttGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATa                                                                <  2:5365144/93‑1 (MQ=255)
                          tGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAAt                                                              >  2:2319513/1‑93 (MQ=255)
                           gAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAAtt                                                             <  1:2601039/93‑1 (MQ=255)
                           gAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAAtt                                                             <  2:5887954/93‑1 (MQ=255)
                              aaCAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGat                                                          <  2:4295054/93‑1 (MQ=255)
                               aCAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGata                                                         >  1:205682/1‑93 (MQ=255)
                               aCAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGata                                                         >  2:5205404/1‑93 (MQ=255)
                                 aGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGAtatt                                                       <  2:2826292/93‑1 (MQ=255)
                                 aGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGAtatt                                                       <  2:238562/93‑1 (MQ=255)
                                 aGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGAtatt                                                       <  1:1037303/93‑1 (MQ=255)
                                 aGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGAtatt                                                       >  1:34602/1‑93 (MQ=255)
                                  gTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGAtattt                                                      >  2:3936326/1‑93 (MQ=255)
                                  gTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGAtattt                                                      >  2:6352177/1‑93 (MQ=255)
                                   tGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGAtatttt                                                     >  2:3025991/1‑92 (MQ=255)
                                   tGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGAtattta                                                     <  2:3554784/93‑1 (MQ=255)
                                   tGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGAtattta                                                     >  2:7460091/1‑93 (MQ=255)
                                   tGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGAtattta                                                     <  1:654482/93‑1 (MQ=255)
                                     aaTTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGAtatttatt                                                   >  1:6736244/1‑93 (MQ=255)
                                       ttattaATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGc                                                 >  1:2885240/1‑93 (MQ=255)
                                        tattaATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCa                                                >  2:7642124/1‑93 (MQ=255)
                                        tattaATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCa                                                >  1:7185532/1‑93 (MQ=255)
                                          ttaATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCata                                              <  2:6878867/93‑1 (MQ=255)
                                          ttaATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCata                                              >  1:2226218/1‑93 (MQ=255)
                                           taATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCatat                                             >  1:1338250/1‑93 (MQ=255)
                                            aaTTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCAtata                                            <  1:241076/93‑1 (MQ=255)
                                            aaTTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCAtata                                            <  1:4630196/93‑1 (MQ=255)
                                             aTTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATaa                                           >  1:1177161/1‑93 (MQ=255)
                                              tttCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATaaa                                          >  2:2213628/1‑93 (MQ=255)
                                              tttCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATaaa                                          <  2:3724475/93‑1 (MQ=255)
                                              tttCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATaaa                                          <  2:2626280/93‑1 (MQ=255)
                                              tttCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATaaa                                          <  2:6722741/93‑1 (MQ=255)
                                              tttCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATaaa                                          <  1:1543362/93‑1 (MQ=255)
                                              tttCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATaaa                                          >  1:2197370/1‑93 (MQ=255)
                                              tttCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGCGATGCATATTACATAATTGATATTTATTGCATATaaa                                          >  2:7070948/1‑93 (MQ=255)
                                               ttCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAat                                         <  2:1174622/93‑1 (MQ=255)
                                                tCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAatt                                        >  1:2308317/1‑92 (MQ=255)
                                                tCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAata                                        >  2:2902074/1‑93 (MQ=255)
                                                tCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAata                                        <  1:5409576/93‑1 (MQ=255)
                                                tCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAata                                        >  1:4848366/1‑93 (MQ=255)
                                                tCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAata                                        <  1:4931572/93‑1 (MQ=255)
                                                  ttATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATGTAAATAtt                                      >  1:3046355/1‑93 (MQ=255)
                                                  ttATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATAtt                                      >  2:209034/1‑93 (MQ=255)
                                                  ttATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATAtt                                      <  2:672601/93‑1 (MQ=255)
                                                   tatataACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATAttt                                     <  2:48320/93‑1 (MQ=255)
                                                    atataACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTtg                                    <  1:879481/93‑1 (MQ=255)
                                                     tataACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTtgt                                   >  1:7397844/1‑93 (MQ=255)
                                                        aaCATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGAt                                <  2:4690219/93‑1 (MQ=255)
                                                         aCATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATc                               >  1:4996270/1‑93 (MQ=255)
                                                         aCATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATc                               >  1:5145418/1‑93 (MQ=255)
                                                         aCATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATc                               >  2:574419/1‑93 (MQ=255)
                                                          cATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCt                              <  2:7333131/93‑1 (MQ=255)
                                                           aTTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTa                             >  1:1786540/1‑93 (MQ=255)
                                                           aTTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTa                             >  1:5434999/1‑93 (MQ=255)
                                                            ttGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTac                            >  2:3534250/1‑93 (MQ=255)
                                                            ttGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTac                            >  2:1944000/1‑93 (MQ=255)
                                                              gAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTacaa                          >  2:1701629/1‑93 (MQ=255)
                                                               aGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTacaac                         >  1:2538653/1‑93 (MQ=255)
                                                               aGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTacaac                         >  1:458697/1‑93 (MQ=255)
                                                                 ttATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAacac                       >  1:1417658/1‑93 (MQ=255)
                                                                 ttATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAacac                       >  2:4727275/1‑93 (MQ=255)
                                                                 ttATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAacac                       >  2:3813372/1‑93 (MQ=255)
                                                                   ataGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGc                     >  2:4104950/1‑93 (MQ=255)
                                                                   ataGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGc                     >  2:4965083/1‑93 (MQ=255)
                                                                    taGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGcc                    <  2:7578068/93‑1 (MQ=255)
                                                                    taGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGcc                    >  1:4924376/1‑93 (MQ=255)
                                                                    taGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGcc                    >  2:5895519/1‑93 (MQ=255)
                                                                      gATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCtt                  >  1:5804590/1‑93 (MQ=255)
                                                                      gATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCtt                  <  1:6427060/93‑1 (MQ=255)
                                                                        tataACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTAt                >  2:7091684/1‑93 (MQ=255)
                                                                        tataACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTAt                >  1:4157528/1‑93 (MQ=255)
                                                                        tataACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTAt                >  1:6207215/1‑93 (MQ=255)
                                                                        tataACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTAt                >  1:2901081/1‑93 (MQ=255)
                                                                          taACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCt              >  2:3451251/1‑93 (MQ=255)
                                                                           aacaacAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTa             >  2:6639207/1‑93 (MQ=255)
                                                                           aacaacAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTa             >  2:4898187/1‑93 (MQ=255)
                                                                              aacaaATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTg          >  1:5946699/1‑93 (MQ=255)
                                                                              aacaaATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTg          >  2:6940858/1‑93 (MQ=255)
                                                                                 aaaTGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCtt       <  1:4068455/93‑1 (MQ=255)
                                                                                 aaaTGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCtt       >  1:4215053/1‑93 (MQ=255)
                                                                                  aaTGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTg      >  1:385598/1‑93 (MQ=255)
                                                                                  aaTGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTg      >  1:5342228/1‑93 (MQ=255)
                                                                                    tGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTc    >  2:3709315/1‑93 (MQ=255)
                                                                                     gTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTcc   >  2:3364677/1‑93 (MQ=255)
                                                                                      tCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCg  >  1:4345671/1‑93 (MQ=255)
                                                                                            |                                                                                       
GCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACCGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCG  >  NC_000913/3077421‑3077600

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