Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,680,766 A→G C660C (TGT→TGC yphG ← DUF5107 domain‑containing protein YphG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,680,7660AG100.0% 213.6 / NA 57C660C (TGT→TGCyphGDUF5107 domain‑containing protein YphG
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (32/25);  total (32/25)

GGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCATCTCTGCCGCCAGGGTGGGCATTCCCCAGTTAATCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGACACAGACCGCGCCACTGCGCGAGAGCTGATTCGCTACGACCATCGAACCAGTTCAGCCACCACAGAGTGGCGTTCAGCGGATAATCGCGCAGCAGTTTCTCGCGCTGCACACGCGCGTTGTCCTGACGACCAC  >  NC_000913/2680629‑2680899
                                                                                                                                         |                                                                                                                                     
ggCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCATCTCTGCCGCCAGGGTGGGCATTCCCCAGTTAATCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGgca                                                                                                                                     <  1:1254991/140‑1 (MQ=255)
  cagcaAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCATCTCTGCCGCCAGGGTGGGCATTCCCCAGTTAATCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGGcaca                                                                                                                                   <  2:1383273/140‑1 (MQ=255)
          tGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCATCTCTGCCGCCAGGGTGGGCATTCCCCAGTTAATCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGGCACAGACCgcgc                                                                                                                           >  2:836533/1‑140 (MQ=255)
          tGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCATCTCTGCCGCCAGGGTGGGCATTCCCCAGTTAATCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGGCACAGACCgcgc                                                                                                                           <  1:1443221/140‑1 (MQ=255)
             cTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCATCTCTGCCGCCAGGGTGGGCATTCCCCAGTTAATCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGGCACAGACCGCGCCAc                                                                                                                        <  2:903144/140‑1 (MQ=255)
                 cAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCATCTCTGCCGCGAGGGTGGGCATTCCCCAGTTAATCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGGCACAGACCGCGCCACTgcg                                                                                                                    >  2:1118387/1‑140 (MQ=255)
                 cAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCATCTCTGCCGCCAGGGTGGGCATTCCCCAGTTAATCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGGCACAGACCGCGCCACTgcg                                                                                                                    >  1:601838/1‑140 (MQ=255)
                  aGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCATCTCTGCCGCCAGGGTGGGCATTCCCCAGTTAATCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGGCACAGACCGCGCCACTgcgc                                                                                                                   >  2:484644/1‑140 (MQ=255)
                   ggTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCATCTCTGCCGCCAGGGTGGGCATTCCCCAGTTAATCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGGCACAGACCGCGCCACTGcgcg                                                                                                                  <  2:53433/140‑1 (MQ=255)
                      aaaGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCATCTCTGCCGCCAGGGTGGGCATTCCCCAGTTAATCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGGCACAGACCGCGCCACTGCGCgaga                                                                                                               >  2:981531/1‑140 (MQ=255)
                           ggCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCATCTCTGCCGCCAGGGTGGGCATTCCCCAGTTAATCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGGCACAGACCGCGCCACTGCGCGAGAGCTGa                                                                                                          <  1:1461121/140‑1 (MQ=255)
                            gCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCATCTCTGCCACCAGGGTGGGCATTCCCCAGTTAATCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGGCACAGACCGCGCCACTGCGCGAGAGCTGAt                                                                                                         >  1:180944/1‑140 (MQ=255)
                               gCGTGCGCTGGCAGTCCAGTGCGTTCAGCATCTCTGCCGCCAGGGTGGGCATTCCCCAGTTAATCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGGCACAGACCGCGCCACTGCGCGAGAGCTGATTCg                                                                                                      >  2:1009739/1‑140 (MQ=255)
                                 gTGCGCTGGCAGTCCAGTGCGTTCAGCATCTCTGCCGCCAGGGTGGGCATTCCCCAGTTAATCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGGCACAGACCGCGCCACTGCGCGAGAGCTGATTCGCt                                                                                                    >  1:1482712/1‑140 (MQ=255)
                                  tGCGCTGGCAGTCCAGTGCGTTCAGCATCTCTGCCGCCAGGGTGGGCATTCCCCAGTTAATCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGGCACAGACCGCGCCACTGCGCGAGAGCTGATTCGCTa                                                                                                   >  1:1460890/1‑140 (MQ=255)
                                  tGCGCTGGCAGTCCAGTGCGTTCAGCATCTCTGCCGCCAGGGTGGGCATTCCCCAGTTAATCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGGCACAGACCGCGCCACTGCGCGAGAGCTGATTCGCTa                                                                                                   <  2:211420/140‑1 (MQ=255)
                                  tGCGCTGGCAGTCCAGTGCGTTCAGCATCTCTGCCGCCAGGGTGGGCATTCCCCAGTTAATCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGGCACAGACCGCCCCACTGCGCGAGAGCTGATTCGCTa                                                                                                   >  2:893502/1‑140 (MQ=255)
                                   gcgcTGGCAGTCCAGTGCGTTCAGCATCTCTGCCGCCAGGGTGGGCATTCCCCAGTTAATCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGGCACAGACCGCGCCACTGCGCGAGAGCTGATTCGCTAc                                                                                                  <  1:913304/140‑1 (MQ=255)
                                         gCAGTCCAGAGCGTTCAGCATCTCTGCCGCCAGGGAGGGCATTCCCCAGTTAATCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGGCACAGACCGCGCCACTGCGCGAGAGCTGATTCGCTACGACCAt                                                                                            >  2:232224/1‑140 (MQ=255)
                                             tCCAGTGCGTTCAGCATCTCTGCCGCCAGGGTGGGCATTCCCCAGTTAATCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGGCACAGACCGCGCCACTGCGCGAGAGCTGATTCGCTACGACCATCGaa                                                                                        <  1:316279/140‑1 (MQ=255)
                                                 gTGCGTTCAGCATCTCTGCCGCCAGGGTGGGCATTCCCCAGTTAATCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGGCACAGACCGCGCCACTGCGCGAGAGCTGATTCGCTACGACCATCGAAccgg                                                                                    >  2:1125164/1‑138 (MQ=255)
                                                  tGCGTTCAGCATCTCTGCCGCCAGGGTGGGCATTCCCCAGTTAATCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGGCACAGACCGCGCCACTGCGCGAGAGCTGATTCGCTACGACCATCGAACCAGt                                                                                   <  2:1043037/140‑1 (MQ=255)
                                                   gCGTTCAGCATCTCTGCCGCCAGGGTGGGCATTCCCCAGTTAATCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGGCACAGACCGCGCCACTGCGCGAGAGCTGATTCGCTACGACCATCGAACCAGtt                                                                                  >  2:1504734/1‑140 (MQ=255)
                                                      ttCAGCATCTCTGCCGCCAGGGTGGGCATTCCCCAGTTAATCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGGCACAGACCGCGCCACTGCGCGAGAGCTGATTCGCTACGACCATCGAACCAGTTCAg                                                                               <  1:637812/140‑1 (MQ=255)
                                                       tCAGCATCTCTGCCGCCAGGGTGGGCATTCCCCAGTTAATCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGGCACAGACCGCGCCACTGCGCGAGAGCTGATTCGCTACGACCATCGAACCAGTTCAGc                                                                              <  1:484644/140‑1 (MQ=255)
                                                            aTCTCTGCCGCCAGGGTGGGCATTCCCCAGTTAATCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGGCACAGACCGCGCCACTGCGCGAGAGCTGATTCGCTACGACCATCGAACCAGTTCAGccacca                                                                         >  2:753375/1‑140 (MQ=255)
                                                                ctGCCGCCAGGGTGGGCATTCCCCAGTTAATCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGGCACAGACCGCGCCACTGCGCGAGAGCTGATTCGCTACGACCATCGAACCAGTTCAGCCACCACaga                                                                     >  1:442063/1‑140 (MQ=255)
                                                                ctGCCGCCAGGGTGGGCATTCCCCAGTTAATCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGGCACAGACCGCGCCACTGCGCGAGAGCTGATTCGCTACGACCATCGAACCAGTTCAGCCACCACaga                                                                     >  2:1215101/1‑140 (MQ=255)
                                                                 tGCCGCCAGGGTGGGCATTCCCCAGTTAATCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGGCACAGACCGCGCCACTGCGCGAGAGCTGATTCGCTACGACCATCGAACCAGTTCAGCCACCACagag                                                                    >  2:285345/1‑140 (MQ=255)
                                                                     gccAGGGTGGGCATTCCCCAGTTAATCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGGCACAGACCGCGCCACTGCGCGAGAGCTGATTCGCTACGACCATCGAACCAGTTCAGCCACCACAGAGTGGc                                                                >  1:91856/1‑140 (MQ=255)
                                                                        aGGGTGGGCATTCCCCAGTTAATCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGGCACAGACCGCGCCACTGCGCGAGAGCTGATTCGCTACGACCATCGAACCAGTTCAGCCACCACAGAGTGGCGtt                                                             >  1:215232/1‑140 (MQ=255)
                                                                        aGGGTGGGCATTCCCCAGTTAATCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGGCACAGACCGCGCCACTGCGCGAGAGCTGATTCGCTACGACCATCGAACCAGTTCAGCCACCACAGAGTGGCGtt                                                             >  2:1526610/1‑140 (MQ=255)
                                                                         gggTGGGCATTCCCCAGTTAATCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGGCACAGACCGCGCCACTGCGCGAGAGCTGATTCGCTACGACCATCGAACCAGTTCAGCCACCACAGAGTGGCGTTc                                                            <  2:957402/140‑1 (MQ=255)
                                                                          ggcgggCATTCCCCAGTTAATCAGCTGCCCGGCGGGCATGGGGGCGTTCACGTCCCTACCGCCCCGCAGACTGCGCCACTGCGCGAGAGCTGAGTCGCGACTCCCATCGACCCAGTACGGCGCCCACAGCGCGGCGTTCa                                                           >  2:94035/4‑140 (MQ=255)
                                                                             gggCATTCCCCAGTTAATCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGGCACAGACCGCGCCACTGCGCGAGAGCTGATTCGCTACGACCATCGAACCAGTTCAGCCACCACAGAGTGGCGTTCAGCg                                                        >  1:951981/1‑140 (MQ=255)
                                                                               gCATTCCCCAGTTAATCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGGCACAGACCGCGCCACTGCGCGAGAGCTGATTCGCTACGACCATCGAACCAGTTCAGCCACCACAGAGTGGCGTTCAGCGGa                                                      >  1:1104664/1‑140 (MQ=255)
                                                                                     cccAGTTAATCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGGCACAGACCGCGCCACTGCGCGAGAGCTGATTCGCTACGACCATCGAACCAGTTCAGCCACCACAGAGTGGCGTTCAGCGGATAATcg                                                >  1:779880/1‑140 (MQ=255)
                                                                                         gTTAATCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGGCACAGACCGCGCCACTGCGCGAGAGCTGATTCGCTACGACCATCGAACCAGTTCAGCCACCACAGAGTGGCGTTCAGCGGATAATCGCgca                                            <  1:1118387/140‑1 (MQ=255)
                                                                                             aTCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGGCACAGACCGCGCCACTGCGCGAGAGCTGATTCGCTACGACCATCGAACCAGTTCAGCCACCACAGAGTGGCGTTCAGCGGATAATCGCGcagcag                                        >  2:1425700/1‑140 (MQ=255)
                                                                                               cAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGGCACAGACCGCGCCACTGCGCGAGAGCTGATTCGCTACGACCATCGAACCAGTTCAGCCACCACAGAGTGGCGTTCAGCGGATAATCGCGCAGCAGtt                                      >  2:484481/1‑140 (MQ=255)
                                                                                                 gTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGGCACAGACCGCGCCACTGCGCGAGAGCTGATTCGCTACGACCATCGAACCAGTTCAGCCACCACAGAGTGGCGTTCAGCGGATAATCGCGCAGCAGTTtc                                    <  1:772326/140‑1 (MQ=255)
                                                                                                  ttGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGGCACAGACCGCGCCACTGCGCGAGAGCTGATTCGCTACGACCATCGAACCGGTTCAGCCACCACAGAGTGGCGTTCAGCGGATAATCGCGCAGCAGTTtct                                   <  1:1125164/140‑1 (MQ=255)
                                                                                                  ttGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGGCACAGACCGCGCCACTGCGCGAGAGCTGATTCGCTACGACCATCGAACCAGTTCAGCCACCACAGAGTGGCGTTCAGCGGATAATCGCGCAGCAGTTtct                                   >  2:964008/1‑140 (MQ=255)
                                                                                                    gCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGGCACAGACCGCGCCACTGCGCGAGAGCTGATTCGCTACGACCATCGAACCAGTTCAGCCACCACAGAGTGGCGTTCAGCGGATAATCGCGCAGCAGTTTCTcg                                 >  1:974972/1‑140 (MQ=255)
                                                                                                    gCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGGCACAGACCGCGCCACTGCGCGAGAGCTGATTCGCTACGACCATCGAACCAGTTCAGCCACCACAGAGTGGCGTTCAGCGGATAATCGCGCAGCAGTTTCTcg                                 <  1:893502/140‑1 (MQ=255)
                                                                                                    gCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGGCACAGACCGCGCCACTGCGCGAGAGCTGATTCGCTACGACCATCGAACCAGTTCAGCCACCACAGAGTGGCGTTCAGCGGATAATCGCGCAGCAGTTTCTcg                                 >  1:1196464/1‑140 (MQ=255)
                                                                                                    gCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGGCACAGACCGCGCCACTGCGCGAGAGCTGATTCGCTACGACCATCGAACCAGTTCAACCACCACAGAGTGGCGTTCAGCGGATAATCGCGCAGCAGTTTCTcg                                 <  2:323188/140‑1 (MQ=255)
                                                                                                      ccGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGGCACAGACCGCGCCACTGCGCGAGAGCTGATTCGCTACGACCATCGAACCAGTTCAGCCACCACAGAGTGGCGTTCAGCGGATAATCGCGCAGCAGTTTCTcgcg                               <  2:91856/140‑1 (MQ=255)
                                                                                                        ggcggTCATCAGGGCGTTAACGTCGCGTCCCTGGCACAGACCGCGCCACTGCGCGAGAGCTGATTCGCTACGACCATCGAACCAGTTCAGCCACCACAGAGTGGCGTTCAGCGGATAATCGCGCAGCAGTTTCTCGCGCt                             <  1:131486/140‑1 (MQ=255)
                                                                                                         gcggTCATCAGGGCGTTAACGTCGCGTCCCTGGCACAGACCGCGCCACTGCGCGAGAGCTGATTCGCTACGACCATCGAACCAGTTCAGCCACCACAGAGTGGCGTTCAGCGGATAATCGCGCAGCAGTTTCTCGCGCTg                            <  2:121231/140‑1 (MQ=255)
                                                                                                            gTCATCAGGGCGTTAACGTCGCGTCCCTGGCACAGACCGCGCCACTGCGCGAGAGCTGATTCGCTACGACCATCGAACCAGTTCAGCCACCACAGAGTGGCGTTCAGCGGATAATCGCGCAGCAGTTTCTCGCGCTGcac                         >  1:413511/1‑140 (MQ=255)
                                                                                                                 caGGGCGTTAACGTCGCGTCCCTGGCACAGACCGCGCCACTGCGCGAGAGCTGATTCGCTACGACCATCGAACCAGTTCAGCCACCACAGAGTGGCGTTCAGCGGATAATCGCGCAGCAGTTTCTCGCGCTGCACAcgcg                    <  1:285345/140‑1 (MQ=255)
                                                                                                                   gggCGTTAACGTCGCGTCCCTGGCACAGACCGCGCCACTGCGCGAGAGCTGATTCGCTACGACCATCGAACCAGTTCAGCCACCACAGAGTGGCGTTCAGCGGATAATCGCGCAGCAGTTTCTCGCGCTGCACAcgcgcg                  <  1:556167/140‑1 (MQ=255)
                                                                                                                    ggCGTTAACGTCGCGTCCCTGGCACAGACCGCGCCACTGCGCGAGAGCTGATTCGCTACGACCATCGAACCAGTTCAGCCACCACAGAGTGGCGTTCAGCGGATAATCGCGCAGCAGTTTCTCGCGCTGCACACGCGCGt                 >  2:663030/1‑140 (MQ=255)
                                                                                                                        ttAACGTCGCGTCCCTGGCACAGACCGCGCCACTGCGCGAGAGCTGATTCGCTACGACCATCGAACCAGTTCAGCCACCACAGAGTGGCGTTCAGCGGATAATCGCGCAGCAGTTTCTCGCGCTGCACACGCGCGTTGTc             >  1:512602/1‑140 (MQ=255)
                                                                                                                         tAACGTCGCGTCCCTGGCACAGACCGCGCCACTGCGCGAGAGCTGATTCGCTACGACCATCGAACCAGTTCAGCCACCACAGAGTGGCGTTCAGCGGATAATCGCGCAGCAGTTTCTCGCGCTGCACACGCGCGTTGTcc            >  2:267885/1‑140 (MQ=255)
                                                                                                                               cgcgTCCCTGGCACAGACCGCGCCACTGCGCGAGAGCTGATTCGCTACGACCATCGAACCAGTTCAGCCACCACAGAGTGGCGTTCAGCGGATAATCGCGCAGCAGTTTCTCGCGCTGCACACGCGCGTTGTCCTgacga      >  2:861593/1‑140 (MQ=255)
                                                                                                                               cgcgTCCCTGGCACAGACCGCGCCACTGCGCGAGAGCTGATTCGCTACGACCATCGAACCAGTTCAGCCACCACAGAGTGGCGTTCAGCGGATAATCGCGCAGCAGTTTCTCGCGCTGCACACGCGCGTTGTCCTgacga      <  1:978421/140‑1 (MQ=255)
                                                                                                                               cgcgTCCCTGGCACAGACCGCGCCACTGCGCGAGAGCTGATTCGCTACGACCATCGAACCAGTTCAGCCACCACAGAGTGGCGTTCAGCGGATAATCGCGCAGCAGTTTCTCGCGCTGCACACGCGCGTTGTCCTgacga      <  1:1526610/140‑1 (MQ=255)
                                                                                                                                   tCCCTGGCACAGACCGCGCCACTGCGCGAGAGCTGATTCGCTACGACCATCGAACCAGTTCAGCCACCACAGAGTGGCGTTCAGCGGATAATCGCGCAGCAGTTTCTCGCGCTGCACACGCGCGTTGTCCTGACGACCAc  <  2:696916/140‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                     
GGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCATCTCTGCCGCCAGGGTGGGCATTCCCCAGTTAATCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGACACAGACCGCGCCACTGCGCGAGAGCTGATTCGCTACGACCATCGAACCAGTTCAGCCACCACAGAGTGGCGTTCAGCGGATAATCGCGCAGCAGTTTCTCGCGCTGCACACGCGCGTTGTCCTGACGACCAC  >  NC_000913/2680629‑2680899

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