Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,598,617 C→T *531* (TAG→TAA lsrK ← autoinducer‑2 (AI‑2) kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,598,6170CT100.0% 244.3 / NA 74*531* (TAG→TAAlsrKautoinducer‑2 (AI‑2) kinase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (40/34);  total (40/34)

TTGAATATCGTACTGGTGATGGAACGATGAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAACTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTCCGGGTCTGGTGTGT  >  NC_000913/1598481‑1598747
                                                                                                                                        |                                                                                                                                  
ttGAATATCGTACTGGTGATGGAACGATGAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAGCTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATta                                                                                                                                  <  1:1480341/139‑1 (MQ=255)
        cGTACTGGTGATGGAACGAGGAATCTGCAAGCTGGCCTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCCATGAAATTATAAACCAg                                                                                                                          >  1:1391271/1‑139 (MQ=255)
          tACTGGTGATGGAACGATGAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGgc                                                                                                                        >  1:89935/1‑139 (MQ=255)
               gTGATGGAACGATGAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAg                                                                                                                          <  2:301468/132‑1 (MQ=255)
               gTGATGGAACGATGAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAg                                                                                                                          >  1:301468/1‑132 (MQ=255)
                    ggAACGATGAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATa                                                                                                              >  2:969001/1‑139 (MQ=255)
                         gATGAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAacgac                                                                                                         >  1:1730054/1‑139 (MQ=255)
                              aTCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGAcgtcgt                                                                                                    <  2:1730054/139‑1 (MQ=255)
                              aTCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGAcgtcgt                                                                                                    <  2:610410/139‑1 (MQ=255)
                              aTCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGAcgtcgt                                                                                                    <  2:635851/139‑1 (MQ=255)
                                   cAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTc                                                                                               <  1:969001/139‑1 (MQ=255)
                                     aGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTGCATAACGACGTCGTCAGTCCa                                                                                             >  2:1103542/1‑139 (MQ=255)
                                         ggCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGgcg                                                                                                                       >  2:1779557/1‑109 (MQ=255)
                                         ggCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGgcg                                                                                                                       <  1:1779557/109‑1 (MQ=255)
                                           cTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCa                                                                                       >  1:171029/1‑139 (MQ=255)
                                              ttAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAAcc                                                                                    <  2:158885/139‑1 (MQ=255)
                                                 aCAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGc                                                                                 <  2:910993/139‑1 (MQ=255)
                                                    aGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGcccc                                                                              >  2:257571/1‑139 (MQ=255)
                                                     gCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCa                                                                             >  1:1494033/1‑139 (MQ=255)
                                                                    agGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATaaa                                                              >  2:1151762/1‑139 (MQ=255)
                                                                    agGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATaaa                                                              >  1:1603126/1‑139 (MQ=255)
                                                                    agGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATaaa                                                              >  2:1514681/1‑139 (MQ=255)
                                                                    agGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATa                                                                <  1:1198268/137‑1 (MQ=255)
                                                                    agGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATa                                                                >  2:1198268/1‑137 (MQ=255)
                                                                     gggAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAAc                                                             >  1:1709612/1‑139 (MQ=255)
                                                                            aaGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGAcgtcgt                                                                                                    <  1:1577524/93‑1 (MQ=255)
                                                                            aaGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGAcgtcgt                                                                                                    >  2:1577524/1‑93 (MQ=255)
                                                                             aGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGgcg                                                                                                                       <  2:937377/73‑1 (MQ=255)
                                                                             aGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGgcg                                                                                                                       >  1:937377/1‑73 (MQ=255)
                                                                                aTAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAActgc                                                          <  1:155465/131‑1 (MQ=255)
                                                                                aTAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAActgc                                                          >  2:155465/1‑131 (MQ=255)
                                                                                      ataCTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAActacctgc                                                      <  2:1303611/129‑1 (MQ=255)
                                                                                      ataCTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAActacctgc                                                      >  1:1303611/1‑129 (MQ=255)
                                                                                         cTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCtt                                                                                                                    <  1:1093902/64‑1 (MQ=255)
                                                                                         cTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCtt                                                                                                                    >  2:1093902/1‑64 (MQ=255)
                                                                                         cTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTg                                                                         <  1:967171/107‑1 (MQ=255)
                                                                                         cTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTg                                                                         >  2:967171/1‑107 (MQ=255)
                                                                                              aaaCTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATcgcg                                           <  2:1855095/132‑1 (MQ=255)
                                                                                              aaaCTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATcgcg                                           >  1:1855095/1‑132 (MQ=255)
                                                                                                   aaCAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATg                                                                                           <  1:776572/79‑1 (MQ=255)
                                                                                                   aaCAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATg                                                                                           >  2:776572/1‑79 (MQ=255)
                                                                                                    aCAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATaa                              >  1:266240/1‑139 (MQ=255)
                                                                                                      aGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGcc                                                                                <  2:1675280/87‑1 (MQ=255)
                                                                                                      aGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGcc                                                                                >  1:1675280/1‑87 (MQ=255)
                                                                                                      aGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCt                                                                    >  1:370429/1‑99 (MQ=255)
                                                                                                      aGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCt                                                                    <  2:370429/99‑1 (MQ=255)
                                                                                                       gCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGt                           <  1:1151762/139‑1 (MQ=255)
                                                                                                                  catGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTc                <  2:89935/139‑1 (MQ=255)
                                                                                                                  catGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTc                <  1:1514681/139‑1 (MQ=255)
                                                                                                                  catGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTc                <  2:266240/139‑1 (MQ=255)
                                                                                                                     gatgTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGaa                                       >  1:168295/1‑113 (MQ=255)
                                                                                                                     gatgTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGaa                                       <  2:168295/113‑1 (MQ=255)
                                                                                                                        gtgGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCaa                                                                                      >  1:251498/1‑63 (MQ=255)
                                                                                                                        gtgGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCaa                                                                                      <  2:251498/63‑1 (MQ=255)
                                                                                                                        gtgGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATa                                                                <  1:209658/85‑1 (MQ=255)
                                                                                                                        gtgGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATa                                                                >  2:209658/1‑85 (MQ=255)
                                                                                                                        gtgGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTCCGGGTc          >  2:690567/1‑139 (MQ=255)
                                                                                                                        gtgGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTCCGGGTc          >  1:981234/1‑139 (MQ=255)
                                                                                                                        gtgGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTCCGGGTc          >  1:731005/1‑139 (MQ=255)
                                                                                                                          ggCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACgtcgtc                                                                                                   <  2:1442782/48‑1 (MQ=255)
                                                                                                                          ggCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACgtcgtc                                                                                                   >  1:1442782/1‑48 (MQ=255)
                                                                                                                          ggCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTCCGGGTCTg        >  2:1709277/1‑139 (MQ=255)
                                                                                                                          ggCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTCCGGGTCTg        <  2:731005/139‑1 (MQ=255)
                                                                                                                           gCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTAt                                               >  2:800208/1‑99 (MQ=255)
                                                                                                                           gCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTAt                                               <  1:800208/99‑1 (MQ=255)
                                                                                                                            cTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGc                                                                           <  1:1291103/70‑1 (MQ=255)
                                                                                                                            cTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGc                                                                           >  2:1291103/1‑70 (MQ=255)
                                                                                                                               tCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATc                                                                      <  2:821951/72‑1 (MQ=255)
                                                                                                                               tCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATc                                                                      >  1:821951/1‑72 (MQ=255)
                                                                                                                               tCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAActgcct                                                        >  1:1460371/1‑86 (MQ=255)
                                                                                                                               tCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAActgcct                                                        <  2:1460371/86‑1 (MQ=255)
                                                                                                                               tCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTCCGGGTCTGgtgtg   <  1:257571/139‑1 (MQ=255)
                                                                                                                               tCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTCCGGGTCTGgtgtg   >  1:412161/1‑139 (MQ=255)
                                                                                                                                cAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATCAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTCCATGCTTTTCCGGGTCTGgtgtgt  >  1:1345034/1‑139 (MQ=255)
                                                                                                                                        |                                                                                                                                  
TTGAATATCGTACTGGTGATGGAACGATGAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAACTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTCCGGGTCTGGTGTGT  >  NC_000913/1598481‑1598747

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