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

TGGTGATGGAACGATGAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAACTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTCCGGGTCTGGTGTG  >  NC_000913/1598494‑1598746
                                                                                                                           |                                                                                                                                 
tGGTGATGGAACGATGAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCt                                                                                                                    >  2:567477/1‑139 (MQ=255)
     aTGGAACGATGAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATGAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCa                                                                                                               >  1:375196/1‑139 (MQ=255)
            gATGAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAacgac                                                                                                        >  1:1151285/1‑139 (MQ=255)
            gATGAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAacgac                                                                                                        <  1:567477/139‑1 (MQ=255)
                 aTCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATGAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGAcgtcgt                                                                                                   <  2:375196/139‑1 (MQ=255)
                 aTCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGAcgtcgt                                                                                                   <  2:1151285/139‑1 (MQ=255)
                      cAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTc                                                                                              <  1:472878/139‑1 (MQ=255)
                                                       agGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATaaa                                                             >  1:221728/1‑139 (MQ=255)
                                                                     aaCTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTa                                               <  2:221728/139‑1 (MQ=255)
                                                                       cTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGcc                                                                               <  1:162485/105‑1 (MQ=255)
                                                                       cTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGcc                                                                               >  2:162485/1‑105 (MQ=255)
                                                                               aaaaaCTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCa                                                                                                                          <  1:869402/54‑1 (MQ=255)
                                                                               aaaaaCTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCa                                                                                                                          >  2:869402/1‑54 (MQ=255)
                                                                               aaaaaCTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGgcg                                                                                                                      <  1:463155/58‑1 (MQ=255)
                                                                               aaaaaCTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGgcg                                                                                                                      >  2:463155/1‑58 (MQ=255)
                                                                                 aaaCTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGAcgtcgt                                                                                                   >  1:430390/1‑75 (MQ=255)
                                                                                 aaaCTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGAcgtcgt                                                                                                   <  2:430390/75‑1 (MQ=255)
                                                                                       aCAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGTCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATaa                             >  2:898987/1‑139 (MQ=255)
                                                                                             aCGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGa                                                                       <  1:303515/91‑1 (MQ=255)
                                                                                             aCGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGa                                                                       >  2:303515/1‑91 (MQ=255)
                                                                                               ggTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGgcg                                                                                                                      >  1:1173828/1‑42 (MQ=255)
                                                                                               ggTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGgcg                                                                                                                      <  2:1173828/42‑1 (MQ=255)
                                                                                                   atcatGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGc                                                                          >  2:907150/1‑82 (MQ=255)
                                                                                                   atcatGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGc                                                                          <  1:907150/82‑1 (MQ=255)
                                                                                                     catGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGcc                                                                               >  1:219460/1‑75 (MQ=255)
                                                                                                     catGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGcc                                                                               <  2:219460/75‑1 (MQ=255)
                                                                                                           gtgGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGAcgtcgt                                                                                                   >  1:1033986/1‑49 (MQ=255)
                                                                                                           gtgGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGAcgtcgt                                                                                                   <  2:1033986/49‑1 (MQ=255)
                                                                                                           gtgGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGAcgtcg                                                                                                    <  1:28491/48‑1 (MQ=255)
                                                                                                           gtgGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGAcgtcg                                                                                                    >  2:28491/1‑48 (MQ=255)
                                                                                                           gtgGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCaa                                                                                     >  1:370457/1‑63 (MQ=255)
                                                                                                           gtgGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCaa                                                                                     <  2:370457/63‑1 (MQ=255)
                                                                                                           gtgGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGa                                                                       >  2:1313895/1‑77 (MQ=255)
                                                                                                           gtgGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGa                                                                       <  1:1313895/77‑1 (MQ=255)
                                                                                                             ggCTGTCAATGAAATTATAACCCAGGCGCTTTCCATGACGAcgtcgt                                                                                                   >  1:856661/1‑47 (MQ=39)
                                                                                                             ggCTGTCAATGAAATTATAACCCAGGCGCTTTCCAAAACGAcgtcgt                                                                                                   <  2:856661/47‑1 (MQ=39)
                                                                                                              gCTGTCAATGAAATTATAACCCAGGCGCTTTCCATa                                                                                                             <  1:678544/36‑1 (MQ=38)
                                                                                                              gCTGTCAATGAAATTATAACCCAGGCGCTTTCCATa                                                                                                             >  2:678544/1‑36 (MQ=255)
                                                                                                              gCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGt                                                                                               <  2:1312891/50‑1 (MQ=255)
                                                                                                              gCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGt                                                                                               >  1:1312891/1‑50 (MQ=255)
                                                                                                              gCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGa                                                                       <  1:47841/74‑1 (MQ=255)
                                                                                                              gCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGa                                                                       >  2:47841/1‑74 (MQ=255)
                                                                                                               cTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTc               <  1:844761/129‑1 (MQ=255)
                                                                                                               cTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTc               >  2:844761/1‑129 (MQ=255)
                                                                                                                  tCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCTGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTCCGGGTCTGgtgtg  <  1:898987/139‑1 (MQ=255)
                                                                                                                           |                                                                                                                                 
TGGTGATGGAACGATGAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAACTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTCCGGGTCTGGTGTG  >  NC_000913/1598494‑1598746

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