Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,815,917 Δ1 bp pseudogene (633/669 nt) rph ← ribonuclease PH (defective),enzyme, Degradation of RNA, RNase PH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,815,9170A.100.0% 286.8 / NA 70pseudogene (633/669 nt)rphribonuclease PH (defective),enzyme, Degradation of RNA, RNase PH
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (38/32);  total (38/32)

TTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTC  >  NC_000913/3815789‑3816042
                                                                                                                                  |                                                                                                                             
ttttCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATT‑‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCTGGGTg                                                                                                                         <  1:1349485/134‑1 (MQ=255)
ttttCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATT‑‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCTGGGTg                                                                                                                         >  2:1349485/1‑134 (MQ=255)
   tCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATTCCCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCTGGGTGAACGGc                                                                                                                   >  2:1294925/1‑139 (MQ=255)
     tACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATT‑‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCTGGGTGAACGGCTCGc                                                                                                               >  2:977323/1‑139 (MQ=255)
            aaaaaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATT‑‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCTGGGTGAACGGCTCGCCTTCTGc                                                                                                        >  2:448736/1‑139 (MQ=255)
                     gCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATT‑‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCTGGGTGAACGGCTCGCCTTCTGCCGTCCCCTg                                                                                               <  1:448736/139‑1 (MQ=255)
                        aCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATT‑‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCTGGGTGAACGGCTCGCCTTCTGCCGTCCCCTg                                                                                               >  1:1507430/1‑136 (MQ=255)
                        aCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATT‑‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCTGGGTGAACGGCTCGCCTTCTGCCGTCCCCTg                                                                                               <  2:1507430/136‑1 (MQ=255)
                        aCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATT‑‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCTGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCAc                                                                                            >  1:43908/1‑139 (MQ=255)
                                   gCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATT‑‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCTGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATgc                                                                                 >  1:1331955/1‑139 (MQ=255)
                                    ccTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATT‑‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCTGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATgcg                                                                                >  2:923962/1‑139 (MQ=255)
                                      ttAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATT‑‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCTGGGTGAACGGCTCGCCTTCTg                                                                                                         <  1:1161407/112‑1 (MQ=255)
                                      ttAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATT‑‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCTGGGTGAACGGCTCGCCTTCTg                                                                                                         >  2:1161407/1‑112 (MQ=255)
                                                  ttGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATT‑‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCTGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGtca                                                                  <  2:43908/139‑1 (MQ=255)
                                                  ttGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATT‑‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCTGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGtca                                                                  <  2:1331955/139‑1 (MQ=255)
                                                   tGCCAGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATT‑‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCTGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGtcat                                                                 >  1:602551/1‑139 (MQ=255)
                                                       aGCGCCGCCTTCTGCGTCGTTACAATGGATTCGATT‑‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCTGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACt                                                             <  2:452632/139‑1 (MQ=255)
                                                         cgccgcCTTCTGCGTCGTTACAATGGATTCGATT‑‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCTGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGcc                                                                              >  2:314162/1‑120 (MQ=255)
                                                         cgccgcCTTCTGCGTCGTTACAATGGATTCGATT‑‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCTGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGcc                                                                              <  1:314162/120‑1 (MQ=255)
                                                              ccTTCTGCGTCGTTACAATGGATTCGATT‑‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCTGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCa                                                      >  2:846107/1‑139 (MQ=255)
                                                                   tGCGTCGTTACAATGGATTCGATT‑‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCTGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCg                                                 <  1:977323/139‑1 (MQ=255)
                                                                    gCGTCGTTACAATGGATTCGATT‑‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCTGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCgg                                                >  1:674986/1‑139 (MQ=255)
                                                                                  gattcgattccCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCTGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCag                                    <  1:1294925/139‑1 (MQ=255)
                                                                                    ttcgatt‑‑cCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCTGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGt                                                                    >  2:1661625/1‑103 (MQ=255)
                                                                                    ttcgatt‑‑cCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCTGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGt                                                                    <  1:1661625/103‑1 (MQ=255)
                                                                                    ttcgatt‑‑cCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCTGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTc                                >  1:1249130/1‑139 (MQ=255)
                                                                                    ttcgatt‑‑cCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCTGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTc                                >  1:86946/1‑139 (MQ=255)
                                                                                    ttcgatt‑‑cCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCTGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTc                                >  1:761142/1‑139 (MQ=255)
                                                                                     tcgatt‑‑cCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCTGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAatgat                                                                                   >  1:1480216/1‑87 (MQ=255)
                                                                                     tcgatt‑‑cCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCTGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAatgat                                                                                   <  2:1480216/87‑1 (MQ=255)
                                                                                     tcgatt‑‑cCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCTGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGTTGCGCCCGTCTTCGGtcat                                                                 <  2:495481/105‑1 (MQ=255)
                                                                                     tcgatt‑‑cCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCTGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGTTGCGCCCGTCTTCGGtcat                                                                 >  1:495481/1‑105 (MQ=255)
                                                                                     tcgatt‑‑cCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCTGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGt                                                                    <  2:524501/102‑1 (MQ=255)
                                                                                     tcgatt‑‑cCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCTGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGt                                                                    >  1:524501/1‑102 (MQ=255)
                                                                                     tcgatt‑‑cCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCTGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCAc                                                              >  2:1516892/1‑108 (MQ=255)
                                                                                     tcgatt‑‑cCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCTGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCAc                                                              <  1:1516892/108‑1 (MQ=255)
                                                                                     tcgatt‑‑cCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCTGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGt                                                         >  2:162392/1‑113 (MQ=255)
                                                                                     tcgatt‑‑cCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCTGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGt                                                         <  1:162392/113‑1 (MQ=255)
                                                                                      cgatt‑‑cCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCTGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCtt                              <  2:1249130/139‑1 (MQ=255)
                                                                                        att‑‑cCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCTGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCa                            <  2:86946/139‑1 (MQ=255)
                                                                                          t‑‑cCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCTGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAAc                          >  1:1260469/1‑139 (MQ=255)
                                                                                              cccTCGGGCCAGAGCCAACAAGATGGGTAGCTCTTCTGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCAt                                                     <  2:1597000/110‑1 (MQ=255)
                                                                                              cccTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCTGGGTGAACGGCTCGCCTTCTGCCGTccc                                                                                                  <  2:1578252/65‑1 (MQ=255)
                                                                                              cccTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCTGGGTGAACGGCTCGCCTTCTGCCGTccc                                                                                                  >  1:1578252/1‑65 (MQ=255)
                                                                                              cccTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCTGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCAt                                                     >  1:1597000/1‑110 (MQ=255)
                                                                                              cccTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCTGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGtct                                             <  2:1486409/118‑1 (MQ=255)
                                                                                              cccTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCTGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGtct                                             >  1:1486409/1‑118 (MQ=255)
                                                                                               ccTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCTGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTa                       >  1:443498/1‑139 (MQ=255)
                                                                                                cTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCTGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGtt                                                        <  2:272936/105‑1 (MQ=255)
                                                                                                cTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCTGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGtt                                                        >  1:272936/1‑105 (MQ=255)
                                                                                                   gggCCAGAGCCAACAAGATGAGTAGCTCTTCTGGGTGAACGGCTCGcc                                                                                                              >  1:1599670/1‑48 (MQ=255)
                                                                                                   gggCCAGAGCCAACAAGATGAGTAGCTCTTCTGGGTGAACGGCTCGcc                                                                                                              <  2:1599670/48‑1 (MQ=255)
                                                                                                   gggCCAGAGCCAACAAGATGAGTAGCTCTTCTGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTc                                                                                         <  1:69368/69‑1 (MQ=255)
                                                                                                   gggCCAGAGCCAACAAGATGAGTAGCTCTTCTGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTc                                                                                         >  2:69368/1‑69 (MQ=255)
                                                                                                    ggCCAGAGCCAACAAGATGAGTAGCTCTTCTGGGTGAACGGCTCGCCTTCTGCCGTccc                                                                                                  >  2:1403446/1‑59 (MQ=255)
                                                                                                    ggCCAGAGCCAACAAGATGAGTAGCTCTTCTGGGTGAACGGCTCGCCTTCTGCCGTccc                                                                                                  <  1:1403446/59‑1 (MQ=255)
                                                                                                        agagCCAACAAGATGAGTAGCTCTTCTGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATc              >  2:131870/1‑139 (MQ=255)
                                                                                                          agCCAACAAGATGAGTAGCTCTTCTGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCAc                                                              <  2:92775/89‑1 (MQ=255)
                                                                                                          agCCAACAAGATGAGTAGCTCTTCTGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCAc                                                              >  1:92775/1‑89 (MQ=255)
                                                                                                               acaaGATGAGTAGCTCTTCTGGGTGAACGGCTCGCCTTCTGCCGTccc                                                                                                  <  1:669710/48‑1 (MQ=255)
                                                                                                               acaaGATGAGTAGCTCTTCTGGGTGAACGGCTCGCCTTCTGCCGTccc                                                                                                  >  2:669710/1‑48 (MQ=255)
                                                                                                               acaaGATGAGTAGCTCTTCTGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAAcc       >  2:525472/1‑139 (MQ=255)
                                                                                                                caaGATGAGTAGCTCTTCTGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCg      <  2:674986/139‑1 (MQ=255)
                                                                                                                caaGATGAGTAGCTCTTCTGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCg      <  2:1260469/139‑1 (MQ=255)
                                                                                                                  aGATGAGTAGCTCTTCTGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCt    <  2:761142/139‑1 (MQ=255)
                                                                                                                    aTGAGTAGCTCTTCTGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTc  >  2:717209/1‑139 (MQ=255)
                                                                                                                         tAGCTCTTCTGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCAc                                                              >  2:866139/1‑74 (MQ=255)
                                                                                                                         tAGCTCTTCTGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCAc                                                              <  1:866139/74‑1 (MQ=255)
                                                                                                                            ctctTCTGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAatgatg                                                                                  <  2:645076/51‑1 (MQ=255)
                                                                                                                            ctctTCTGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAatgatg                                                                                  >  1:645076/1‑51 (MQ=255)
                                                                                                                                  |                                                                                                                             
TTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTC  >  NC_000913/3815789‑3816042

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