Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,521,852 T→C D209G (GAT→GGT)  xapR ← transcriptional activator of xapAB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,521,8520TC100.0% 260.8 / NA 71D209G (GAT→GGT) xapRtranscriptional activator of xapAB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (32/39);  total (32/39)

TGTGCGTAGCTGTCCGCTATCAATGTGATACCAATGCCCATACTGACCATAGCGAGTACCGTTTGCGGTTCATTAACTTCGCGAATAACAACCGGTGAAAATCCCACCTGCTGGCAAACTCGCTGCAAAAAATCCCAGTCAGTGTAAACGGGCGGCATTGTAACAAAATACTCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGCCGATTCATG  >  NC_000913/2521720‑2521976
                                                                                                                                    |                                                                                                                            
tgtgCGTAGCTGTCCGCTATCAATGTGATACCAATGCCCATACTGACCATAGCGAGTACCGTTTGCGGTTCATTAACTTCGCGAATAACAACCGGTGAAAATCCCACCTGCTGGCAAACTCGCTGCAAAAAACCCcagt                                                                                                                        <  1:296851/139‑1 (MQ=255)
 gtgCGTAGCTGTCCGCTATCAATGTGATACCAATGCCCATACTGACCATAGCGAGTACCGTTTGCGGTTCATTAACTTCGCGAATAACAACCGGTGAAAATCCCACCTGCTGGCAAACTCGCTGCAAAAAACCCcagt                                                                                                                        <  2:15122/138‑1 (MQ=255)
 gtgCGTAGCTGTCCGCTATCAATGTGATACCAATGCCCATACTGACCATAGCGAGTACCGTTTGCGGTTCATTAACTTCGCGAATAACAACCGGTGAAAATCCCACCTGCTGGCAAACTCGCTGCAAAAAACCCcagt                                                                                                                        >  1:15122/1‑138 (MQ=255)
           gTCCGCTATCAATGTGATACCAATGCCCATACTGACCATAGCGAGTACCGTTTGCGGTTCATTAACTTCGCGAATAACAACCGGTGAAAATCCCACCTGCTGGCAAACTCGCTGCAAAAAACCCCAGTCAGTGTAAACg                                                                                                             <  1:769009/139‑1 (MQ=255)
                  aTCAATGTGATACCAATGCCCATACTGACCATAGCGAGTACCGTTTGCGGTTCATTAACTTCGCGAATAACAACCGGTGAAAATCCCACCTGCTGGCAAACTCGCTGCAAAAAACCCCAGTCAGTGTAAACGGGCGGCa                                                                                                      <  2:1250615/139‑1 (MQ=255)
                     aaTGTGATACCAATGCCCATACTGACCATAGCGAGTACCGTTTGCGGTTCATTAACTTCGCGAATAACAACCGGTGAAAATCCCACCTGCTGGCAAACTCGCTGCAAAAAACCCCAGTGAGTGTAAACGGGCGGCATTg                                                                                                   <  2:841619/139‑1 (MQ=255)
                     aaTGTGATACCAATGCCCATACTGACCATAGCGAGTACCGTTTGCGGTTCATTAACTTCGCGAATAACAACCGGTGAAAATCCCACCTGCTGGCAAACTCGCTGCAAAAAACCCCAGTCAGTGTAAACGGGCGGCATTg                                                                                                   <  2:1236763/139‑1 (MQ=255)
                           aTACCAATGCCCATACTGACCATAGCGAGTACCGTTTGCGGTTCATTAACTTCGCGAATAACAACCGGTGAAAATCCCACCTGCTGGCAAACTCGCTGCAAAAAACCCCAGTCAGTGTAAACGGGCGGCATTGTAACaa                                                                                             >  2:1098313/1‑139 (MQ=255)
                             aCCAATGCCCATACTGACCATAGCGAGTACCGTTTGCGGTTCATTAACTTCGCGAATAACAACCGGTGAAAATCCCACCTGCTGGCAAACTCGCTGCAAAAAACCCCAGTCAGTGTAAACGGGCGGCATTGTAACaaaa                                                                                           >  2:147561/1‑139 (MQ=255)
                             aCCAATGCCCATACTGACCATAGCGAGTACCGTTTGCGGTTCATTAACTTCGCGAATAACAACCGGTGAAAATCCCACCTGCTGGCAAACTCGCTGCAAAAAACCCCAGTCAGTGTAAACGGGCGGCATTGTAACaaaa                                                                                           >  1:15583/1‑139 (MQ=255)
                              ccAATGCCCATACTGACCATAGCGAGTACCGTTTGCGGTTCATTAACTTCGCGAATAACAACCGGTGAAAATCCCACCTGCTGGCAAACTCGCTGCAAAAAACCCCAGTCAGTGTAAACGGGCGGCATTGTAACAAAAt                                                                                          <  2:15583/139‑1 (MQ=255)
                                        tACTGACCATAGCGAGTACCGTTTGCGGTTCATTAACTTCGCGAATAACAACCGGTGAAAATCCCACCTGCTGGCAAACTCGCTGCAAAAAACCCCAGTCAGTGTAAACGGGCGGCATTGTAACAAAATACTCGTCACg                                                                                >  1:451446/1‑139 (MQ=255)
                                         aCTGACCATAGCGAGTACCGTTTGCGGTTCATTAACTTCGCGAATAACAACCGGTGAAAATCCCACCTGCTGGCAAACTCGCTGCAAAAAACCCca                                                                                                                          >  1:954311/1‑96 (MQ=255)
                                         aCTGACCACAGCGAGTACCGTTTGCGGTTCATTAACTTCGCGAATAACAACCGGTGAAAATCCCACCTGCTGGCAAACTCTCTGCAAAAAACCCca                                                                                                                          <  2:954311/96‑1 (MQ=255)
                                           tGACCATAGCGAGTACCGTTTGCGGTTCATTAACTTCGCGAATAACAACCGGTGAAAATCCCACCTGCTGGCAAACTCGCTGCAAAAAACCCCAGTCAGTGTAAACGGGCGGCATTGTAACAAAATACTCGTCACGTAg                                                                             <  1:147561/139‑1 (MQ=255)
                                             aCCATAGCGAGTACCGTTTGCGGTTCATTAACTTCGCGAATAACAACCGGTGAAAATCCCACCTGCTGGCAAACTCGCTGCAAAAAACCCCAGTCAGTGTAAACggg                                                                                                           >  2:1239832/1‑107 (MQ=255)
                                             aCCATAGCGAGTACCGTTTGCGGTTCATTAACTTCGCGAATAACAACCGGTGAAAATCCCACCTGCTGGCAAACTCGCTGCAAAAAACCCCAGTCAGTGTAAACggg                                                                                                           <  1:1239832/107‑1 (MQ=255)
                                             aCCATAGCGAGTACCGTTTGCGGTTCATTAACTTCGCGAATAACAACCGGTGAAAATCCCACCTGCTGGCAAACTCGCTGCAAAAAACCCCAGTCAGTGTAAACGggcgg                                                                                                        <  1:94455/110‑1 (MQ=255)
                                             aCCATAGCGAGTACCGTTTGCGGTTCATTAACTTCGCGAATAACAACCGGTGAAAATCCCACCTGCTGGCAAACTCGCTGCAAAAAACCCCAGTCAGTGTAAACGggcgg                                                                                                        <  2:1366023/110‑1 (MQ=255)
                                             aCCATAGCGAGTACCGTTTGCGGTTCATTAACTTCGCGAATAACAACCGGTGAAAATCCCACCTGCTGGCAAACTCGCTGCAAAAAACCCCAGTCAGTGTAAACGggcgg                                                                                                        >  1:1366023/1‑110 (MQ=255)
                                             aCCATAGCGAGTACCGTTTGCGGTTCATTAACTTCGCGAATAACAACCGGTGAAAATCCCACCTGCTGGCAAACTCGCTGCAAAAAACCCCAGTCAGTGTAAACGggcgg                                                                                                        >  2:94455/1‑110 (MQ=255)
                                              ccATAGCGAGTACCGTTTGCGGTTCATTAACTTCGCGAATAACAACCGGTGAAAATCCCACCTGCTGGCAAACTCGCTGCAAAAAACCCCAGTCAGTGTAAACGGGCGGCATTGTAACAAAATACTCGTCACGTAgcgc                                                                          >  2:1345452/1‑139 (MQ=255)
                                               cATAGCGAGTACCGTTTGCGGTTCATTAACTTCGCGAATAACAACCGGTGAAAATCCCACCTGCTGGCAAACTCGCTGCAAAAAACCCCAGTCAGTGTAAACGGGCGGCATTGTAACAAAATACTCGTCACGTAGCGCt                                                                         >  2:1383053/1‑139 (MQ=255)
                                                    cGAGTACCGTTTGCGGTTCATTAACTTCGCGAATAACAACCGGTGAAAATCCCACCTGCTGGCAAACTCGCTGCAAAAAACCCCAGTCAGTGTAAACGGGCGGCATTGTAACAAAATACTCGTCACGTAGCGCTTCCAg                                                                    >  2:722318/1‑139 (MQ=255)
                                                       gTACCGTTTGCGGTTCATTAACTTCGCGAATAACAACCGGTGAAAATCCCACCTGCTGGCAAACTCGCTGCAAAAAACCCCAGTCAGTGTaaa                                                                                                               >  1:347351/1‑93 (MQ=255)
                                                       gTACCGTTTGCGGTTCATTAACTTCGCGAATAACAACCGGTGAAAATCCCACCTGCTGGCAAACTCGCTGCAAAAAACCCCAGTCAGTGTaaa                                                                                                               <  2:347351/93‑1 (MQ=255)
                                                       gTACCGTTTGCGGTTCATTAACTTCGCGAATAACAACCGGTGAAAATCCCACCTGCTGGCAAACTCGCTGCAAAAAACCCCAGTCAGTGTAAACGGGCGGCATTGTAACAAAATACTCGTCACGTAGCGCTTCCAGCgg                                                                 >  1:977123/1‑139 (MQ=255)
                                                       gTACCGTTTGCGGTTCATTAACTTCGCGAATAACAACCGGTGAAAATCCCACCTGCTGGCAAACTCGCTGCAAAAAACCCCAGTCAGTGTAAACGGGCGGCATTGTAACAAAATACTCGTCACGTAGCGCTTCCAGCgg                                                                 >  1:10161/1‑139 (MQ=255)
                                                         aCCGTTTGCGGTTCATTAACTTCGCGAATAACAACCGGTGAAAATCCCACCTGCTGGCAAACTCGCTGCAAAAAACCCCAGTCAGTGTAAACGGGCGGCATTGTAACAAAATACTCGTCACGTAGCGCTTCCAGCGGGa                                                               <  1:1098313/139‑1 (MQ=255)
                                                          ccGTTTGCGGTTCATTAACTTCGCGAATAACAACCGGTGAAAATCCCACCTGCTGGCAAACTCGCTGCAAAAAACCCCAGTCAGTGTAAACGGGCGGCATTGTAACAAAATACTCGTCACGTAGCGCTTCCAGCGGGAc                                                              <  1:1275497/139‑1 (MQ=255)
                                                           cGTTTGCGGTTCATTAACTTCGCGAATAACAACCGGTGAAAATCCCACCTGCTGGCAAACTCGCTGCAAAAAACCCCAGTCAGTGTAAACGGGCGGCATTGTAACAAAATACTCGTCACGTAGCGCTTCCAGCGGGACg                                                             <  1:1383053/139‑1 (MQ=255)
                                                            gTTTGCGGTTCATTAACTTCGCGAATAACAACCGGTGAAAATCCCACCTGCTGGCAAACTCGCTGCAAAAAAcccc                                                                                                                           >  2:1262170/1‑76 (MQ=255)
                                                            gTTTGCGGTTCATTAACTTCGCGAATAACAACCGGTGAAAATCCCACCTGCTGGCAAACTCGCTGCAAAAAAcccc                                                                                                                           <  1:1262170/76‑1 (MQ=255)
                                                            gTTTGCGGTTCATTAACTTCGCGAATAACAACCGGTGAAAATCCCACCTGCTGGCAAACTCGCTGCAAAAAACCCCAGTCAGTGTAAACGGGCGGCATTGTAACa                                                                                              >  1:140503/1‑105 (MQ=255)
                                                            gTTTGCGGTTCATTAACTTCGCGAATAACAACCGGTGAAAATCCCACCTGCTGGCAAACTCGCTGCAAAAAACCCCAGTCAGTGTAAACGGGCGGCATTGTAACa                                                                                              <  2:140503/105‑1 (MQ=255)
                                                             tttGCGGTTCATTAACTTCGCGAATAACAACCGGTGAAAATCCCACCTGCTGGCAAACTCGCTGCAAAAAACCCCAGTCAGTGTAAACGGGCGGCATTGTAACAAAATACTCGTCACGTAGCGCTTCCAGCGGGACGGt                                                           <  2:451446/139‑1 (MQ=255)
                                                               tGCGGTTCATTAACTTCGCGAATAACAACCGGTGAAAATCCCACCTGCTGGCAAACTCGCTGCAAAAAACCCCAGTCAGTGTAAACGGGCGGCATTGTAACAAAATACTCGTCACGTAGCGCTTCCAGCGGGACGGTgg                                                         >  1:1173020/1‑139 (MQ=255)
                                                               tGCGGTTCATTAACTTCGCGAATAACAACCGGTGAAAATCCCACCTGCTGGCAAACTCGCTGCAAAAAACCCCAGTCAGTGTAAACGGGCGGCATTGTAACAAAATACTCGTCACGTAGCGCTTCCAGCGGGACGGTgg                                                         <  1:1345452/139‑1 (MQ=255)
                                                                gCGGTTCATTAACTTCGCGAATAACAACCGGTGAAAATCCCACCTGCTGGCAAACTCGCTGCAAAAAACCCCAGTCAGTGTAAACGGGCGGCATTGTAACAAAATACTCGTCACGTAGCGCTTCCAGCGGGACGGTGGa                                                        >  1:193579/1‑139 (MQ=255)
                                                                 cGGTTCATTAACTTCGCGAATAACAACCGGTGAAAATCCCACCTGCTGGCAAACTCGCTGCAAAAAACCCCAGTCAGTGTAAACGGGCGGCATTGTAACAAAATACTCGTCACGTAGCGCTTCCAGCGGGACGGTGGaa                                                       <  2:1252012/139‑1 (MQ=255)
                                                                  ggTTCATTAACTTCGCGAATAACAACCGGTGAAAATCCCACCTGCTGGCAAACTCGCTGCAAAAAACCCCAGTCAGTGTAAACGGGCGGCATTGTAACAAAATACTCGTCACGTAGCGCTTCCAGCGGGACGGTGGaaa                                                      <  2:977123/139‑1 (MQ=255)
                                                                  ggTTCATTAACTTCGCGAATAACAACCGGTGAAAATCCCACCTGCTGGCAAACTCGCTGCAAAAAACCCCAGTCAGTGTAAACGGGCGGCATTGTAACAAAATACTCGTCACGTAGCGCTTCCAGCGGGACGGTGGaaa                                                      <  2:10161/139‑1 (MQ=255)
                                                                     tCATTAACTTCGCGAATAACAACCGGTGAAAATCCCACCTGCTGGCAAACTCGCTGCAAAAAACCCCAGTCAGTGTAAACGGGCGGCATTGTAACAAAATACTCGTCACGTAGCGCTTCCAGCGGGACGGTgg                                                         >  1:286756/1‑133 (MQ=255)
                                                                     tCATTAACTTCGCGAATAACAACCGGTGAAAATCCCAACTGCTGGCAAACTCGCTGCAAAAAACCCCAGTCAGTGTAAACGGGCGGCATTGTAACAAAATACTCGTCACGTAGCGCTTCCAGCGGGACGGTgg                                                         <  2:286756/133‑1 (MQ=255)
                                                                              tCGCGAATAACAACCGGTGAAAATCCCACCTGCTGGCAAACTCGCTGCAAAAAACCCCAGTCAGTGTAAACGGGCGGCATTGTAACAAAATACTCGTCACGTAGCGCTTCCAGCGGGa                                                               >  2:710700/1‑118 (MQ=255)
                                                                              tCGCGAATAACAACCGGTGAAAATCCCACCTGCTGGCAAACTCGCTGCAAAAAACCCCAGTCAGTGTAAACGGGCGGCATTGTAACAAAATACTCGTCACGTAGCGCTTCCAGCGGGa                                                               <  1:710700/118‑1 (MQ=255)
                                                                                gcgAATAACAACCGGTGAAAATCCCACCTGCTGGCAAACTCGCTGCAAAAAACCCCAGTCAGTGTAAACGGGCGGCATTGTAACAAAATACTCGTCACGTAGCGCTTCCAGCGGGACGGTGGaaaa                                                     >  1:940665/1‑126 (MQ=255)
                                                                                gcgAATAACAACCGGTGAAAATCCCACCTGCTGGCAAACTCGCTGCAAAAAACCCCAGTCAGTGTAAACGGGCGGCATTGTAACAAAATACTCGTCACGTAGCGCTTCCAGCGGGACGGTGGaaaa                                                     <  2:940665/126‑1 (MQ=255)
                                                                                gcgAATAACAACCGGTGAAAATCCCACCTGCTGGCAAACTCGCTGCAAAAAACCCCAGTCAGTGTAAACGGGCGGCATTGTAACAAAATACTCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAgatga                                        >  2:1223960/1‑139 (MQ=255)
                                                                                 cgAATAACAACCGGTGAAAATCCCACCTGCTGGCAAACTCGCTGCAAAAAACCCCAGTCAGTGTAAACGGGCGGCa                                                                                                      >  1:264613/1‑76 (MQ=255)
                                                                                 cgAATAACAACCGGTGAAAATCCCACCTGCTGGCAAACTCGCTGCAAAAAACCCCAGTCAGTGTAAACGGGCGGCa                                                                                                      <  2:264613/76‑1 (MQ=255)
                                                                                 cgAATAACAACCGGTGAAAATCCCACCTGCTGGCAAACTCGCTGCAAAAAACCCCAGTCAGTGTAAACGGGCGGCATTGTAACAAAATACTCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAgatgat                                       <  1:1223960/139‑1 (MQ=255)
                                                                                         aacCGGTGAAAATCCCACCTGCTGGCAAACTCGCTGCAAAAAACCCcagtc                                                                                                                       <  2:632136/51‑1 (MQ=255)
                                                                                         aacCGGTGAAAATCCCACCTGCTGGCAAACTCGCTGCAAAAAACCCcagtc                                                                                                                       >  1:632136/1‑51 (MQ=255)
                                                                                              gTGAAAATCCCACCTGCTGGCAAACTCGCTGCAAAAAACCCCAGTCAGTGTAAACGGGCGGCATTGTAACAAAATACTCGTCACGTAGCGCTTCCAGCGGGa                                                               >  2:79088/1‑102 (MQ=255)
                                                                                              gTGAAAATCCCACCTGCTGGCAAACTCGCTGCAAAAAACCCCAGTCAGTGTAAACGGGCGGCATTGTAACAAAATACTCGTCACGTAGCGCTTCCAGCGGGa                                                               <  1:79088/102‑1 (MQ=255)
                                                                                                gAAAATCCCACCTGCTGGCAAACTCGCTGCAAAAAACCCCAGTCAGTGTAAACGGGCGGCATTGTAACAAAATACTCGTCACGTAGCGCTTCCAGCGGGAc                                                              <  2:1367015/101‑1 (MQ=255)
                                                                                                gAAAATCCCACCTGCTGGCAAACTCGCTGCAAAAAACCCCAGTCAGTGTAAACGGGCGGCATTGTAACAAAATACTCGTCACGTAGCGCTTCCAGCGGGAc                                                              >  1:1367015/1‑101 (MQ=255)
                                                                                                   aaTCCCACCTGCTGGCAAACTCGCTGCAAAAAACCCCAGTCAGTGTAAACGGGCGGCATTGTAACAAAATACTCGTGACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAgatg                                         >  2:1026617/1‑119 (MQ=255)
                                                                                                   aaTCCCACCTGCTGGCAAACTCGCTGCAAAAAACCCCAGTCAGTGTAAACGGGCGGCATTGTAACAAAATACTCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAgatg                                         <  1:1026617/119‑1 (MQ=255)
                                                                                                             gctgGCAAACTCGCTGCAAAAAACCCCAGTCAGTGTAAACGGGCGGCAt                                                                                                     <  2:130392/49‑1 (MQ=255)
                                                                                                             gctgGCAAACTCGCTGCAAAAAACCCCAGTCAGTGTAAACGGGCGGCAt                                                                                                     >  1:130392/1‑49 (MQ=255)
                                                                                                                ggCAAACTCGCTGCAAAAAACCCCAGTCAGTGTAAACGGGCGGCATTGTAACAAAATACTCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGa                >  2:320262/1‑131 (MQ=255)
                                                                                                                ggCAAACTCGCTGCAAAAAACCCCAGTCAGTGTAAACGGGCGGCATTGTAACAAAATACTCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGa                <  1:320262/131‑1 (MQ=255)
                                                                                                                     aCTCGCTGCAAAAAACCCCAGTCAGTGTAAACGggcggcc                                                                                                      <  2:893830/40‑2 (MQ=255)
                                                                                                                     aCTCGCTGCAAAAAACCCCAGTCAGTGTAAACGGGCGGCa                                                                                                      >  1:893830/1‑40 (MQ=255)
                                                                                                                     aCTCGCTGCAAAAAACCCCAGTCAGTGTAAACGGGCGGCATTGTAACAAAATACTCGTCACGTAGCGCTTCCAGCGGGACg                                                             <  2:244530/81‑1 (MQ=255)
                                                                                                                     aCTCGCTGCAAAAAACCCCAGTCAGTGTAAACGGGCGGCATTGTAACAAAATACTCGTCACGTAGCGCTTCCAGCGGGACg                                                             >  1:244530/1‑81 (MQ=255)
                                                                                                                      cTCGCTGCAAAAAACCCCAGTCAGTGTAAACGGGCGGCATTGTAACAAAATACTCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGCCGATTCATg  <  2:193579/139‑1 (MQ=255)
                                                                                                                      cTCGCTGCAAAAAACCCCAGTCAGTGTAAACGGGCGGCATTGTAACAAAATACTCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGCCGATTCATg  <  2:1173020/139‑1 (MQ=255)
                                                                                                                           tGCAAAAAACCCCAGTCAGTGTAAACGGGCGGCATTGTAACAAAATACTCGTCACGTAgcg                                                                           <  2:1031487/61‑1 (MQ=255)
                                                                                                                           tGCAAAAAACCCCAGTCAGTGTAAACGGGCGGCATTGTAACAAAATACTCGTCACGTAgcg                                                                           >  1:1031487/1‑61 (MQ=255)
                                                                                                                                 aaaCCCCAGTCAGTGTAAACGGGCGGCATTGTAACAAAATACTCGTCACGTAGCGCTTCCAGCggg                                                                <  2:866062/66‑1 (MQ=255)
                                                                                                                                    |                                                                                                                            
TGTGCGTAGCTGTCCGCTATCAATGTGATACCAATGCCCATACTGACCATAGCGAGTACCGTTTGCGGTTCATTAACTTCGCGAATAACAACCGGTGAAAATCCCACCTGCTGGCAAACTCGCTGCAAAAAATCCCAGTCAGTGTAAACGGGCGGCATTGTAACAAAATACTCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGCCGATTCATG  >  NC_000913/2521720‑2521976

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: