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% 162.2 / NA 51*531* (TAG→TAAlsrKautoinducer‑2 (AI‑2) kinase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (26/25);  total (26/25)

CAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAACTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTCCGGGTCTGGTGTG  >  NC_000913/1598516‑1598746
                                                                                                     |                                                                                                                                 
cAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCGGGCGCTTTCCATAACGACGTCGTCAGTc                                                                                              <  2:1134377/139‑1 (MQ=255)
     tGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGa                                                                                         <  2:599350/139‑1 (MQ=255)
                 aGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGcccc                                                                             >  1:1582882/1‑139 (MQ=255)
                                 agGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGAcgtcgt                                                                                                   <  2:463284/101‑1 (MQ=255)
                                 agGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGAcgtcgt                                                                                                   >  1:463284/1‑101 (MQ=255)
                                  gggAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGAcgtcgt                                                                                                   >  2:569441/1‑100 (MQ=255)
                                  gggAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGAcgtcgt                                                                                                   <  1:569441/100‑1 (MQ=255)
                                      aaTAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGcc                                                                               <  2:35344/116‑1 (MQ=255)
                                      aaTAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGcc                                                                               >  1:35344/1‑116 (MQ=255)
                                          aGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGgcg                                                                                                                      <  1:1559542/73‑1 (MQ=255)
                                          aGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGgcg                                                                                                                      >  2:1559542/1‑73 (MQ=255)
                                            aaTAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGt                                                                                               >  2:975341/1‑94 (MQ=255)
                                            aaTAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGt                                                                                               <  1:975341/94‑1 (MQ=255)
                                             aTAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGgcg                                                                                                                      >  2:1123733/1‑70 (MQ=255)
                                             aTAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGgcg                                                                                                                      >  2:944805/1‑70 (MQ=255)
                                             aTAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGgcg                                                                                                                      <  1:944805/70‑1 (MQ=255)
                                             aTAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGgcg                                                                                                                      <  1:1123733/70‑1 (MQ=255)
                                                    taCTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCaa                                                                                     <  1:15031/96‑1 (MQ=255)
                                                    taCTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCaa                                                                                     >  2:15031/1‑96 (MQ=255)
                                                    taCTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATcgcg                                          <  2:1582882/139‑1 (MQ=255)
                                                      cTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGAcgtcgt                                                                                                   <  1:1267236/80‑1 (MQ=255)
                                                      cTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGAcgtcgt                                                                                                   >  2:1267236/1‑80 (MQ=255)
                                                      cTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATa                                                               <  1:431926/116‑1 (MQ=255)
                                                      cTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATa                                                               >  2:431926/1‑116 (MQ=255)
                                                      cTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTg                                        >  1:227141/1‑139 (MQ=255)
                                                          aaaaCTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGaa                                      >  1:383991/1‑137 (MQ=255)
                                                          aaaaCTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGaa                                      <  2:383991/137‑1 (MQ=255)
                                                              cTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGgcg                                                                                                                      >  2:1486361/1‑53 (MQ=255)
                                                              cTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGgcg                                                                                                                      <  1:1486361/53‑1 (MQ=255)
                                                              cTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTg                                        >  1:727674/1‑131 (MQ=255)
                                                              cTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTg                                        <  2:727674/131‑1 (MQ=255)
                                                                 aCAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGAcgtcgt                                                                                                   >  1:1084621/1‑69 (MQ=255)
                                                                 aCAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGAcgtcgt                                                                                                   <  2:1084621/69‑1 (MQ=255)
                                                                         ggTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGgcg                                                                                                                      <  1:899276/42‑1 (MQ=255)
                                                                         ggTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGgcg                                                                                                                      >  2:899276/1‑42 (MQ=255)
                                                                         ggTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGAcgtcgt                                                                                                   <  1:1772483/61‑1 (MQ=255)
                                                                         ggTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGAcgtcgt                                                                                                   >  2:1772483/1‑61 (MQ=255)
                                                                               catGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTc               >  2:147424/1‑139 (MQ=255)
                                                                                    tgtgGCTGTCAATGAAATTATAACCCAGGCGCttt                                                                                                                  <  2:1266014/35‑1 (MQ=255)
                                                                                    tgtgGCTGTCAATGAAATTATAACCCAGGCGCttt                                                                                                                  >  1:1266014/1‑35 (MQ=255)
                                                                                     gtgGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAacgac                                                                                                        >  2:301515/1‑44 (MQ=255)
                                                                                     gtgGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAacgac                                                                                                        <  1:301515/44‑1 (MQ=255)
                                                                                     gtgGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTCCGGGTc         >  1:1830443/1‑139 (MQ=255)
                                                                                       ggCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGcc                                                                               >  2:1196382/1‑67 (MQ=255)
                                                                                       ggCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGcc                                                                               <  1:1196382/67‑1 (MQ=255)
                                                                                        gCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGa                                                                       <  2:543215/74‑1 (MQ=255)
                                                                                        gCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGa                                                                       >  1:543215/1‑74 (MQ=255)
                                                                                            tCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTCCGGGTCTGgtgtg  >  2:401291/1‑139 (MQ=255)
                                                                                            tCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTCCGGGTCTGgtgtg  <  2:769991/139‑1 (MQ=255)
                                                                                               aTGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTc                                                                                              >  2:631467/1‑44 (MQ=255)
                                                                                               aTGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTc                                                                                              <  1:631467/44‑1 (MQ=255)
                                                                                                     |                                                                                                                                 
CAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAACTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTCCGGGTCTGGTGTG  >  NC_000913/1598516‑1598746

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