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_000913	1,598,617	0	C	T	100.0%	154.8 / NA	48	*531* (TAG→TAA)	lsrK	autoinducer‑2 (AI‑2) kinase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (26/22);  total (26/22)

TGAATATCGTACTGGTGATGGAACGATGAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAACTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTCCGGGTCTGGTGTG  >  NC_000913/1598482‑1598746                                                                                                                                        |                                                                                                                                  tGAATATCGTACTGGTGATGGAACGATGAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATtat                                                                                                                                <  2:587584/139‑1 (MQ=255)                  aTGGAACGATGAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAgg                                                                                                                        <  2:368376/130‑1 (MQ=255)                  aTGGAACGATGAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAgg                                                                                                                        >  1:368376/1‑130 (MQ=255)                              aTCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGAcgtcgt                                                                                                   <  2:1462413/139‑1 (MQ=255)                              aTCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGAcgtcgt                                                                                                   <  1:366749/139‑1 (MQ=255)                                                    aGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGcccc                                                                             >  2:689748/1‑139 (MQ=255)                                                                    agGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAacgac                                                                                                        <  1:54225/96‑1 (MQ=255)                                                                    agGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAacgac                                                                                                        >  2:54225/1‑96 (MQ=255)                                                                       gAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAAc                                                                                    >  1:690026/1‑113 (MQ=255)                                                                       gAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAAc                                                                                    <  2:690026/113‑1 (MQ=255)                                                                        aaaTAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAActgc                                                         <  2:273606/139‑1 (MQ=255)                                                                                  aaCTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAAccc                                                                                                                           <  2:1112623/63‑1 (MQ=255)                                                                                  aaCTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAAccc                                                                                                                           >  1:1112623/1‑63 (MQ=255)                                                                                      ataCTCAAAAACTAACAGCCACGGTCATCATGATGTGGTTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATcgc                                           >  2:1300658/1‑139 (MQ=255)                                                                                       taCTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATcgcg                                          >  1:415016/1‑139 (MQ=255)                                                                                         cTCAAAAACTAACAGCCATGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTg                                        >  1:1466265/1‑139 (MQ=255)                                                                                         cTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGa                                                                       >  2:1139782/1‑108 (MQ=255)                                                                                         cTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGa                                                                       <  1:1139782/108‑1 (MQ=255)                                                                                            aaaaaCTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGa                                                                       >  2:861789/1‑105 (MQ=255)                                                                                            aaaaaCTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGa                                                                       <  1:861789/105‑1 (MQ=255)                                                                                                 cTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGa                                                                       <  2:284598/100‑1 (MQ=255)                                                                                                 cTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGa                                                                       >  1:284598/1‑100 (MQ=255)                                                                                                    aCAGCCACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGCCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATaa                             >  2:731724/1‑139 (MQ=255)                                                                                                        ccACGGTCATCATGATGTGGCTGTCAATGAAATTATAAccc                                                                                                                           >  1:17228/1‑41 (MQ=255)                                                                                                        ccACGGTCATCATGATGTGGCTGTCAATGAAATTATAAccc                                                                                                                           <  2:17228/41‑1 (MQ=255)                                                                                                        ccACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTc                                                                                              >  2:795581/1‑70 (MQ=255)                                                                                                        ccACGGTCATCATGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTc                                                                                              <  1:795581/70‑1 (MQ=255)                                                                                                                  catGATGTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTc               <  2:1466265/139‑1 (MQ=255)                                                                                                                  catGATGTGGCTGTCAATGAAATTATAACCCAGCCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTc               <  1:731724/139‑1 (MQ=255)                                                                                                                    tgatgTGGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTCCg             >  1:82924/1‑139 (MQ=255)                                                                                                                       tgtgGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTc               >  1:1086405/1‑134 (MQ=255)                                                                                                                       tgtgGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTCCGGGt          >  2:550708/1‑139 (MQ=255)                                                                                                                       tgtgGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACGGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTc               <  2:1086405/134‑1 (MQ=255)                                                                                                                        gtgGCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTCCGGGTc         >  2:553144/1‑139 (MQ=255)                                                                                                                           gCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTAt                                              <  1:736492/99‑1 (MQ=255)                                                                                                                           gCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTAt                                              >  2:736492/1‑99 (MQ=255)                                                                                                                           gCTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTCCGGGTCTgg      >  2:587454/1‑139 (MQ=255)                                                                                                                            cTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCttt                 >  2:761301/1‑127 (MQ=255)                                                                                                                            cTGTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCttt                 <  1:761301/127‑1 (MQ=255)                                                                                                                              gTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGt                                                                                               >  2:1460030/1‑47 (MQ=255)                                                                                                                              gTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGt                                                                                               <  1:1460030/47‑1 (MQ=255)                                                                                                                              gTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATc                                                                     >  1:1275558/1‑73 (MQ=255)                                                                                                                              gTCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATc                                                                     <  2:1275558/73‑1 (MQ=255)                                                                                                                               tCAATGAAATTATAACCCAGGCGCTTTCCATAACGAcgtcgt                                                                                                   <  1:483954/42‑1 (MQ=255)                                                                                                                               tCAATGAAATTATAACCCAGGCGCTTTCCATAACGAcgtcgt                                                                                                   >  2:483954/1‑42 (MQ=255)                                                                                                                               tCAATGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTCCGGGTCTGgtgtg  >  1:1532979/1‑139 (MQ=255)                                                                                                                                 aaTGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTc               <  1:1470695/124‑1 (MQ=255)                                                                                                                                 aaTGAAATTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTc               >  2:1470695/1‑124 (MQ=255)                                                                                                                                        |                                                                                                                                  TGAATATCGTACTGGTGATGGAACGATGAATCTGCAAGCTGGCTTTTAACAAGCCAGCTCTAAAAGAAGGGAAATAAGAATAACTATACTCAAAAACTAACAGCCACGGTCATCATGATGTGGCTGTCAATGAAACTATAACCCAGGCGCTTTCCATAACGACGTCGTCAGTCCATGATCAACCAGCCCCAGCTGATCCTGATAAACTGCCTGCCACTTATCGCGTGAATCCTGATAAAGTTCATGCTTTTCCGGGTCTGGTGTG  >  NC_000913/1598482‑1598746

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: ‑