Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,600,081	C→T	W43* (TGG→TGA)	lsrK ←	autoinducer‑2 kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,600,081	0	C	T	96.5%	90.1 / NA	29	W43* (TGG→TGA)	lsrK	autoinducer‑2 kinase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (12/16);  minor base G (0/1);  total (12/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTAAGTTCGCTAACTTCGCGTGCCGCTCTGGCATCCACATTGGCGCAGGCCCAGATCGGGGCTCCTTCATTATTATATAAAACAATGCCTTCACGCATCGAACATGCCGAAACGGCAGCGATATACTCCGGGGCTATGCCGGCGTTGTGCAGCGCCTGGCGCATACACTCACACGCCAGTTGCCAGTTTTTGTTGAGATCAAATTCCATAGAACCAGGAACGTCCGGTACTGCCAGATGCCGCCACTCCGCCTGTCCCACTGCTATTTGATTGCCTTCCAGGTCGAATATCACAGCCCGAATACTTCCGGTGCCTGCATCCAGCGCCATCAGGTAGTACTTTGATTCTGAAAGGGTAAAGAGTCGAGCCATTAGATTATCCTCGGCTTATGCGTTCTGTCTGTTCCTCTATACGTTCTCCATCATTCCCGGTAATAAGGTCATGCAAATTTAACTACGTAAAATCGCCGCTGCTGTGTCCT  >  NC_000913/1599839‑1600319                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                               ttAAGTTCGCTAACTTCGCGTGCCGCTCTGGCATCCACAGTGGCGCAGGCCCAGATCGGGTCTCCTTCATTATTATATAAAACAATGCCTTCACGCATCGAACATGCCGAAACGGCAGCGATATACTCCGGGGCTATGCCGGCGTTGTGCAGCGCCTGGCGCATACACTCACACGCCAGTGGCCAGTTTTTGTTGAGATCAAATTCCATAGAACCAGGAACGTCCGGTACTGCCAGATGCCGTCACTCCg                                                                                                                                                                                                                                         <  2:518729/250‑1 (MQ=255)      ttCGCTAACTTCGCGTGCCGCTCTGGCATCCACATTGGCGCAGGCCCAGATCGGGGCTCCTTCATTATTATATAAAACAATGCCTTCACGCATCGAACATGCCGAAACGGCAGCGATATACTCCGGGGCTATGCCGGCGTTGTGCAGCGCCTGGCGCATACACTCACACGCCAGTTGCCAGTTTTTGTTGAGATCAAATTCCATAGAACCAGGAACGTCCGGTACTGCCAGATGCCGTCACTCCGCCTGTc                                                                                                                                                                                                                                   >  2:423131/1‑251 (MQ=255)                                          ggCGCAGGCCCAGATCGGGGCTCCTTCATTATTATATAAAACAATGCCTTCACGCATCGAACATGCCGAAACGGCAGCGATATACTCCGGGGCTATGCCGGCGGTGTGCAGCGCCTGTCGCATACACTCACACGCCAGTTGCCAGTTTTTGTTGAGATCAAATTCCATAGAACCAGGAACGTCCGGTACTGCCAGATGCCGTCACTCCGCCTGTCCCACTGCTATTTGATTGCCTTCCAGGTCGAATATca                                                                                                                                                                                               <  1:489835/251‑1 (MQ=255)                                                    cAGATCGGGGCTCCTTCATTATTATATAAAACAATGCCTTCACGCATCGAACATGCCGAAACGGCAGCGATATACTCCGGGGCTATGCCGGCGTTGTGCAGCGCCTGGCGCATACACTCACACGCCAGTTGCCAGTTTTTGTTGAGATCAAATTCCATAGAACCAGGAACGTCCGGTACTGCCAGATGCCGTCACTCCGCCTGTCCCACTGCTATTTGATTGCCTTCCAGGTCGAATATCACAGCCCGAAt                                                                                                                                                                                     <  1:152562/251‑1 (MQ=255)                                                                        attatATAAAACAAGGCCTTCACGCATCGAACATGCCGAAACGGCAGCTATATACTCCGGGGCGATGCCGGCGTTGTGCAGCGCCTTGCGCATACACTCACACGCCAGTTGCCAGTTTTTGTTGAGATCAAATTCCATAGAACCAGGAACGTCCGGTACTGCCAGATTCCGTCACTCCGCCTGTCCCACTGCTATTTGATTGCCTTCCAGGTCGAATATCACAGCCCGAATACTTCCGGTGCCTGCATcc                                                                                                                                                                  <  2:362554/250‑1 (MQ=255)                                                                               aaaaCAATGCCTTCACGCATCGAACATGCCGAAACGGCAGCGATATACTCCGGGGCTATGCCGGCGTTGTGCAGCGCCTGGCGCATACACTCACACGCCAGTTGCCAGTTGTTGTTGAGAGCAAAGTCCATAGAACCAGGAACGTCCGGTACTGCCAGATGCCGTCACTCCGCCTGTCCCACTGCTATTGGATTGCCGTCCAGGTCGAATATCACAGCCCGAATACTTCCGGTGCCTGCAGCCATCGCCAt                                                                                                                                                          <  2:513237/251‑1 (MQ=255)                                                                                 aaCAATGCCTTCACGCATCGAACATGCCGAAACGGCAGCGATATACTCCGGGGCTATGCCGGCGTTGTGCAGCGCCTGGCGCATACACTCACACGCCAGTTGCCAGTTTTTGTTGAGATCAAATTCCATAGAACCAGGAACGTCCGGTACTGCCAGATGCCGTCACTCCGCCTGTCCCACTGCTATTTGATTGCCTTCCAGGTCGAATATCACAGCCCGAATACTTCCGGTGCCTGCATCCAGCGCCATCa                                                                                                                                                        <  1:143779/251‑1 (MQ=255)                                                                                  aCAATGCCTTCACGCATCGAACATGCCGAAACGGCAGCGATATACTCCGGGGCTATGCCGGCGTTGTGCAGCGCCTGGCGCATACACTCACACGCCAGTTGCCAGTTTTTGTTGAGATCAAATTCCATAGAACCAGGAACGTCCGGTACTGCCAGATGCCGTCACTCCGCCTGTCCCACTGCTATTTGATTGCCTTCCAGGTCGAATATCACAGCCCGAATACTTCCGGTGCCTGCATCCAGCGCCATc                                                                                                                                                         >  2:417386/1‑249 (MQ=255)                                                                                       gCCTTCACGCATCGAACATGCCGAAACGGCAGCGATATACTCCGGGGCTATGCCGGCGTTGTGCAGCGCCTGGCGCATACACTCACACGCCAGTTGCCAGTTTTTGTTGAGATCAAATTCCATAGAACCAGGAACGTCCGGTACTGCCAGATGCCGTCACTCCGCCTGTCCCACTGCTATTTGATTGCCTTCCAGGTCGAATATCACAGCCCGAATACTTCCGGTGCCTGCATCCAGCGCCATCAGgtagt                                                                                                                                                  >  1:429536/1‑251 (MQ=255)                                                                                        ccTTAACGCATCGAACATGCCGAAACGGCAGCGATATACTCCGGGGCTATGCCGGCGTTGTGCAGCGCCTGGCGCATACACTCACACGCCAGTTGCCAGTTTTTGTTGAGATCAAATTCCATAGAACCAGGAACGTCCGGTACTGCCAGATGCCGTCACTCCGCCTGTCCCACTGCTATTTGATTGCCTTCCAGGTCGAATATCACAGCCCGAATACTTCCGGTGCCTGCATCCAGCGCCATCAGgtagt                                                                                                                                                  >  1:12666/1‑250 (MQ=255)                                                                                             aCGCATCGAACATGCCGAAACGGCAGCGATATACTCCGGGGCTATGCCGGCGTTGTGCAGCGCCTGGCGCATACACTCACACGCCAGTTGCCAGTTTTTGTTGAGATCAAATTCCATAGAACCAGGAACGTCCGGTACTGCCAGATGCCGTCACTCCGCCTGTCCCACTGCTATTTGATTGCCTTCCAGGTCGAATATCACAGCCCGAATACTTCCGGTGCCTGCATCCAGCGCCATCAGGTAGTACTTTg                                                                                                                                            <  1:107319/251‑1 (MQ=255)                                                                                                 aTCGAACCTGCAGAACCGGCAGCGATCTAATCCGGGGCTATGCCGGCGTTGTGCAGCGCCTGGCGCATACACTCACACGCCAGTTGCCAGTTTTTGTTGAGATCAAATTCCATAGAACCAGGAACGTCCGGTACTGCCAGATGCCGTCACTCCGCCTGTCCCACTGCTATTTGATTGCCTTCCCGGTCGAATATCACAGCCCGAATACTTCCGGTGCCTGCATCCAGCGCCATCAGGCAGTACTTTGATTc                                                                                                                                        >  2:228307/1‑251 (MQ=255)                                                                                                       cATGCCGAAACGGCAGCGATATACGCCGGGGCTATGCCGGCGTTGTGCAGCGCCTGGCGCATACACTCACACGCCAGTTGCCAGTTTTTGTTGAGATCAAATTCCATAGAACCAGGAACGTCCGGTACTGCCAGATGCCGTCACTCCGCCTGTCCCACTGCTATTTGATTGCCTTCCAGGTCGAATATCACAGCCCGAATACTTCCGGTGCCTGCATCCAGCGCCATCAGGTAGTACTTTGATTCTGAAAg                                                                                                                                  <  1:255404/251‑1 (MQ=255)                                                                                                             gaaccGGCAGCGATATACTCCGGGGCTATGCCGGCCTTTTGCAGCGCCTGGCGCATACACTCACACGCCAGTTGCCAGTTTTTGTTGAGATCAAATTCCATAGAACCAGGAACGTCCTGTACTGCCAGATGCCGTCACTCCGCCTGTCCCACTGCTATTTGATTGCCTTCCAGTTCGAATATCACAGCCCGAATACTTCCGGTGCCTGCATCCAGCGCCATCAGGTAGTACTTTGATTCTGAAAGGGTaaa                                                                                                                            <  1:88057/247‑1 (MQ=255)                                                                                                                     aGCGATATACTCCGGGGCTATGCCGGCGTTGTGCAGCGCCTGGCGCATACACTCACACGCCAGTTGCCAGTTTTTGTTGAGATCAAATTCCATAGAACCAGGAACGTCCGGTACTGCCAGATGCCGTCACTCCGCCTGTCCCACTGCTATTTGATTACCTTCCAGGTCGAATATCACAGCCCGAATACTTCCGGTGCCTGCATCCAGCGCCATCAGGTAGTACTTTGATTCTGAAAGGGTAAAGAGTCg                                                                                                                      >  1:118069/1‑249 (MQ=255)                                                                                                                      gCGATATACTCCGGGGCTAGGCCGGCGTTGTGCAGCGCCTGGCGCATACACTCACACGCCAGTTGCCAGTTTTTGTTGAGATCAAATTCCATAGAACCAGGAACGTCCGGTACTGCCAGATGCCGTCCCTCCGCCTGTCCCACTGCTATTTGATTGCCTTCCAGGTCGAATATCACAGCCCGAAGACTTCCGGTGCCTGCATCCAGCGCCATCAGGTAGTACTTTGATTCTGAAAGGGTAAATAGTCGAGc                                                                                                                   <  2:207948/251‑1 (MQ=255)                                                                                                                                       tATGCCGGCGTTGTGCAGCGCCTGGCGCATACACTCACACGCCAGTTGCCAGTTTTTGTTGAGATCAAATTCCATAGAACCAGGAACGTCCGGTACTGCCAGATGCCGTCACTCCGCCTGTCCCACTGCTATTTGATTGCCTTCCAGGTCGAATATCACAGCCCGAATACTTCCGGTGCCTGCATCCAGCGCCATCAGGTAGTACTTTGATTCTGAAAGGGTAAAGAGTCGAGCCATTAGATTATCCTCgg                                                                                                  >  1:487884/1‑251 (MQ=255)                                                                                                                                        aTGCCGGCGTTGTGCAGCGCCTGGCGCATACACTCACACGCCAGTTGCCAGTTTTTGTTGAGATCAAATTCCATAGAACCAGGAACGTCCGGTACTGCCAGATGCCGTCACTCCGCCTGTCCCACTGCTATTTGATTGCCTTCCAGGTCGAATATCACAGCCCGAATACTTCCGGTGCCTGCATCCAGCGCCATCAGGTAGTACTTTGATTCTGAAAGGGTAAAGAGTCGAGCCATTAGATTATCCTCgg                                                                                                  <  1:228307/250‑1 (MQ=255)                                                                                                                                            cGGCGTTGGGCAGCCCCTGGCCCATCCACTCCAACGCCAGTTGCCAGTTTTTGTTGATATCAAATTCCTTAGACCCAGGAACGTCCGGTACTGCCAGATGCCGTCACTCCGCCTGTCCCACTGCTATTTTATTGCCTTCCAGGTCGAATATCACAGCCCGAATACTTCCGGTGCCTGCATCCAGCGCCATCAGGTAGTACTTTGATTCTGAAAGGGAAAAGAGTCGAGCCATTAGATTATCCTCGGCTTAt                                                                                             <  2:447170/251‑1 (MQ=255)                                                                                                                                                            cTGGCGCATACACTCACACGCCAGTTGCCAGTTTTTGTTGAGATCAAATTCCATAGAACCAGGAACGTCCGGTACTGCCAGATGCCGTCACTCCGCCTGTCCCACTGCTATTTGATTGCCTTCCAGGTCGAATATCACAGCCCGAATACTTCCGGTGCCTGCATCCAGCGCCATCAGGTAGTACTTTGATTCTGAAAGGGTAAAGAGTCGAGCCATTAGATTATCCTCGGCTTATGCGTTCTGTCTGTTcc                                                                             >  2:395028/1‑251 (MQ=255)                                                                                                                                                                   aTACACTCACACGCCAGTTGCCAGTTTTTGTTGAGATCAAATTCCATAGAACCAGGAACGTCCGGTACTGCCAGATGCCGTCACTCCGCCTGTCCCACTGCTATTTGATTGCCTTCCAGGTCGAATATCACAGCCCGAATACTTCCGGTGCCTGCATCCAGCGCCATCAGGTAGTACTTTGATTCTGAAAGGGTAAAGAGTCGAGCCATTAGATTATCCTCGGCTTATGCGTTCTGTCTGTTCCTCTATAc                                                                      >  1:288040/1‑251 (MQ=255)                                                                                                                                                                      catctacaagccatttgccagttTTTGGTGAGATCAAATTCCATAGCACCAGGAACGTCCGTTACTGCCAGATGTCGTCACTCCGCCTGTCCCCCTGCTATTTGATTGCCTTCCAGATCGAATATCACAGCCCGAATATTTCCGTTGCCTGCATCCATCGCCATCATGTAGTACTTTGATTCTGAAAGGGTAAAGAGTCGAGCCATTAGATTATCCTCGGCTTATGCGTTCTGTCTGTTCCTCTATACGtt                                                                   <  1:20615/242‑1 (MQ=255)                                                                                                                                                                           acacGCCAGTTGCCAGTTTTTGTTGAGATCAAATTCCATAGAACCAGGAACGTCCGGTACTGCCAGATGCCGTCACTCCGCCTGTCCCACTGCTATTTGATTGCCTTCCAGGTCGAATATCACAGCCCGAATACTTCCGGTGCCTGCATCCAGCGCCATCAGGTAGTACTTTGATTCTGAAAGGGTAAAGAGTCGAGCCATTAGATTATCCTCGGCTTATGCGTTCTGTCTGTTCCTCTATACGTTCTcc                                                               >  2:511695/1‑250 (MQ=255)                                                                                                                                                                                   gttgcaagttTTTGTTGAAATAAAATTCCATAGAACCAGCACCTTTCGGTACTGCCAGATGCCTTCACTCCGCCTGCCCCCCTGCTATTTGACTGCCTCCCAGGTCGAATATCACAGACCGAATCCTTCCGGTTCCTGCATCCAGCGCCATCAGTTATTTCTTGGATTCTGAAAGGGTAAAGAGTCGAGCCTTTAGATTATCCTCGGCTTATGCGTTCTGTCTGTGCCTCTATTCGTTCACCATCATTccc                                                      <  1:85750/251‑1 (MQ=255)                                                                                                                                                                                                                  aGAACCAGGAACGTCCGGTACTGCCAGATGCCGTCACTCCGCCTGTCCCACTGCTATTTGATTGCCTTCCAGGTCGAATATCACAGCCCGAATACTTCCGGTGCCTGCATCCAGCGCCATCAGGTAGTACTTTGATTCTGAAAGGGTAAAGAGTCGAGCCATTAGATTATCCTCGGCTTATGCGTTCTGTCTGTTCCTCTATACGTTCTCCATCATTCCCGGTAATAAGGTCATGCAAATTTAACTACGt                        >  2:411741/1‑250 (MQ=255)                                                                                                                                                                                                                   gAACCAGGATCGTCCGGCTCTGCCAGATGCCGGCACTCCGCCTGTCCCCCTGCTATTTGATTGCCTTCCAGGTCGAATATCACAGCCCGAATACTTCCGGTGCATGCATCCAGCGCCATCAGGTAGTACTTTGATTCTGAAAGGGGAAAGAGTCGAGCCATTAGATTATCCTCGGCTTATGCGTTCTGTCTGTTCCTCTATACGTTCTCCATCATTCCCGGTAATAAGGTCATGCAAATTTAACTACGTaa                      <  1:328089/251‑1 (MQ=255)                                                                                                                                                                                                                         ggAACGTCCGGTACTGCCAGATGCCGTCACTCCGCCTGTCCCACTGCTATTTGATTACCTTCCAGTGCGAATATCACAGCCCGAATACTTCCGGTGCCTGCATCCAGCGCCATCAGGTAGTACTTTGATTCTGAAAGGGTAAAGAGTCGAGCCATTAGATTATCCTCGGCTTATGCGTTCTGTCTGTTCCTCTATACGTTCTCCATCATTCCCGGTAATAAGGTCATGCAAATTTAACTACGTAAAATcgc                <  2:118069/251‑1 (MQ=255)                                                                                                                                                                                                                              gTCCGGTACTGCCAGATGCCGTCACTCCGCCTGTCCCACTGCTATTTGATTGCCTTCCAGGTCGAATATCACAGCCCGAATACTTCCGGTGCCTGCATCCAGCGCCATCAGGTAGTACTTTGATTCTGAAAGGGTAAAGAGTCGAGCCATTAGATTATCCTCGGCTTATGCGTTCTGTCTGTTCCTCTATACGTTCTCCATCATTCCCGGTAATAAGGTCATGCAAATTTAACTACGTAAAATCGCCgctg           >  2:54171/1‑251 (MQ=255)                                                                                                                                                                                                                                        gCCAGATGCCGTCACTCCGCCTGTCCCACTGCTATTTGATTGCCTTCCAGGTCGAATATCACAGCCCGAATACTTCCGGTGCCTGCATCCAGCGCCATCAGGTAGTACTTTGATTCTGAAAGGGTAAAGAGTCGAGCCATTAGATTATCCTCGGCTTATGCGTTCTGTCTGTTCCTCTATACGTTCTCCATCATTCCCGGTAATAAGGTCATGCAAATTTAACTACGTAAAATCGCCGCTGCTGTGTCCt  <  1:411741/250‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                               TTAAGTTCGCTAACTTCGCGTGCCGCTCTGGCATCCACATTGGCGCAGGCCCAGATCGGGGCTCCTTCATTATTATATAAAACAATGCCTTCACGCATCGAACATGCCGAAACGGCAGCGATATACTCCGGGGCTATGCCGGCGTTGTGCAGCGCCTGGCGCATACACTCACACGCCAGTTGCCAGTTTTTGTTGAGATCAAATTCCATAGAACCAGGAACGTCCGGTACTGCCAGATGCCGCCACTCCGCCTGTCCCACTGCTATTTGATTGCCTTCCAGGTCGAATATCACAGCCCGAATACTTCCGGTGCCTGCATCCAGCGCCATCAGGTAGTACTTTGATTCTGAAAGGGTAAAGAGTCGAGCCATTAGATTATCCTCGGCTTATGCGTTCTGTCTGTTCCTCTATACGTTCTCCATCATTCCCGGTAATAAGGTCATGCAAATTTAACTACGTAAAATCGCCGCTGCTGTGTCCT  >  NC_000913/1599839‑1600319

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑