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	86.7%	73.1 / 6.0	31	W43* (TGG→TGA)	lsrK	autoinducer‑2 kinase
Reads supporting (aligned to +/- strand):  ref base C (1/3);  new base T (17/9);  total (19/12)
Fisher's exact test for biased strand distribution p-value = 2.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.77e-01

TCTTTAAGTTCGCTAACTTCGCGTGCCGCTCTGGCATCCACATTGGCGCAGGCCCAGATCGGGGCTCCTTCATTATTATATAAAACAATGCCTTCACGCATCGAACATGCCGAAACGGCAGCGATATACTCCGGGGCTATGCCGGCGTTGTGCAGCGCCTGGCGCATACACTCACACGCCAGTTGCCAGTTTTTGTTGAGATCAAATTCCATAGAACCAGGAACGTCCGGTACTGCCAGATGCCGCCACTCCGCCTGTCCCACTGCTATTTGATTGCCTTCCAGGTCGAATATCACAGCCCGAATACTTCCGGTGCCTGCATCCAGCGCCATCAGGTAGTACTTTGATTCTGAAAGGGTAAAGAGTCGAGCCATTAGATTATCCTCGGCTTATGCGTTCTGTCTGTTCCTCTATACGTTCTCCATCATTCCCGGTAATAAGGTCATGCAAATTTAACTACGTAAAATCGCCGCTGCTGTGTCCT  >  NC_000913/1599836‑1600319                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                               tCTTTAAGTGCTCTAACTTCGCGTGCCGCTCTGGCATCCACATTGGCGCAGGCCCAGATCGGGGCTCCTTCATTATTATATAAAACAATGCCTTCACGCATCGAACATGCCGAAACGGCAGCGATATACTCCGGGGCTATGCCGGCGTTGTGCAGCGCCTGGCGCATACACTCACACGCCAGTTGCCAGTTTTTGTTGAGATCAAATTCCATAGAACCAGGAACGTCCGGTACTGCCAGATGCCGTCACt                                                                                                                                                                                                                                            <  2:255069/250‑1 (MQ=255)                  ttCGCGTGCCGCTCTGGCATCCACATTAGCGCAGGCCCAGATCGCGGCTCCTTCATTATTATATAAAACAATGCCTTCACGCATCGAACCTGCCGAAACGGCAGCGATATACTCCGGGGCTATGCCGGCGTAGTGCAGCCCCTGGCGCCTACACTCACACGCCAGTTGCCAGTTTTTGTTGAGATCCAATTCCATAGAACAAGGAACGTCCGGTACTGACAGATGACGTCACTCCGCCTGTCCCACTggt                                                                                                                                                                                                                           >  2:41361/1‑248 (MQ=255)                      cgTGCCGCTCTGGCATCCACATTGGCGCAGGCCCAGATCGGGGCTCCTTCATTATTATATAAAACAATGCCTTCACGCATCGAACATGCCGAAACGGCAGCGATATACTCCGGGGCTATGCCGGCGTTGTGCAGCGCCTGGCGCATACACTCACACGCCAGTTGCCAGTTTTTGTTGAGATCAAATTCCATAGAACCAGGAACGTCCGGTACTGCCAGATGCCGCCACTCCGCCTGTCCCACTGCTATTTg                                                                                                                                                                                                                      <  1:222725/251‑1 (MQ=255)                                                gcAGGCCCAGATCGGGGCTCCTTCATTATTATATAAAACAATGCCTTCACGCATCGAACATGCCGAAACGGCAGCGATATACTCCGGGGCTATGCCGGCGTTGTGCAGCGCCTGGCGCATACACTCACACGCCAGTTGCCAGTTTTTGTTGAGATCAAATTCCATAGAACCAGGAACGTCCGGTACTGCCAGATGCCGTCACTCCGCCTGTCCCACTGCTATTTGATTGCCTTCCAGGTCGAATATCACAg                                                                                                                                                                                            <  2:121779/251‑1 (MQ=255)                                                         gATCGGGGCTCCTTCATTATTATATAAAACAATGCCTTCACGCATCGAACATGCCGAAACGGCAGCGATATACTCCGGGGCTATGCCGGCGTTGTGCAGCGCCTGGCGCATACACTCACACGCCAGTTGCCAGTTTTTGTTGAGATCAAATTCCATAGAACCAGGAACGTCCGGTACTGCCAGATGCCGTCACTCCGCCTGTCCCACTGCTATTTGATTGCCTTCCAGGTCGAATATCACAGCCCGAATAc                                                                                                                                                                                   >  1:241074/1‑251 (MQ=255)                                                             ggggCTCCTTCATTATTATATAAAACAATGCCTTCACGCATCGAACATGCCGAAACGGCAGCGATATACTCCGGGGCTATGCCGGCGTTGTGCAGCGCCTGGCGCATACACTCACACGCCAGTTGCCAGTTTTTGTTGAGATCAAATTCCATAGAACCAGGAACGTCCGGTACTGCCAGATGCCGCCACTCCGCCTGTCCCACTGCTATTTGATTGCCTTCCAGGTCGAATATCACAGCCCGAATACTTcc                                                                                                                                                                               <  1:242328/251‑1 (MQ=255)                                                                               tataAAACAATGCCTTCACGCATCGAACATGCCGAAACGGCAGCGATATACTCCGGGGCTATGCCGGCGTTGTGCAGCGCCTGGCGCATACACTCACACGCCAGTTGCCAGTTTTTGTTGAGATCAAATTCCATAGAACCAGGAACGTCCGGTACTGCCAGATGCCGTCACTCCGCCTGTCCCACTGCTATTTGATTGCCTTCCAGGTCGAATATCACAGCCCGAATACTTCCGGGGCCTGCATCCAgcgc                                                                                                                                                             >  2:114939/1‑251 (MQ=255)                                                                                        aTGCCTTCACGCATCGAACATGCCGAAACGGCAGCGATATACTCCGGGGCTATGCCGGCGTTGTGCAGCGCCTGGCGCATACACTCACACGCCAGTTGCCAGTTTTTGTTGAGATCAAATTCCATAGAACCAGGAACGTCCGGTACTGCCAGATGCCGTCACTCCGCCTGTCCCACTGCTATTTGATTGCCTTCCAGGTCGAATATCACAGCCCGAATACTTCCGGTGCCTGCATCCAGCGCCATCAGgt                                                                                                                                                     >  1:288281/1‑250 (MQ=255)                                                                                                  gCATCGAACATGCCGAAACGGCAGCGATATACTCCGGGGCTATGCCGGCGTTGTGCAGCGCCTGGCGCATACACTCACACGCCAGTTGCCAGTTTTTGTTGAGATCAAATTCCATAGAACCAGGAACGTCCGGTACTGCCAGATGCCGTCACTCCGCCTGTCCCACTGCTATTTGATTGCCTTCCAGGTCGAATATCACAGCCCGAATACTTCCGGTGCCTGCATCCAGCGCCATCAGGTAGTACTTTGAt                                                                                                                                          >  1:303995/1‑251 (MQ=255)                                                                                                        aaCATGCCGAAGCGGCCGCGATATACTCCTGGGTGATGCCGGCGTTGTTCAGCGCCTGGGGCGTACACTCACACGCCAGTTGCCAGGTTTTGTTGAGATCAAAGTCCATAGAACCAGGAACGTCCGGTACTGCCAGATGCCGTCACTCCGCCTGGCCCACTGCTATTTGATTGCCTTCCAGGTCGAATATCACAGCCCGAATACTTCCGGTTCCTGCATCCAGCGCCATCAGGTAGTACTTTGATGCTGaa                                                                                                                                    <  2:236302/251‑1 (MQ=255)                                                                                                              ccGAAACGGCAGCGATATACTCCGGGGCTATGCCGGCGTTGTGCAGCGCCTGGCGCATACACTCACACGCCAGTTGCCAGTTTTTGTTGAGCTCAAATTCCATAGAACCAGGAACGTCCGGTACTGCCAGATGCCGTCACTCCGCCTGTCCCACTGCTATTTGATTGCCTTCCAGGTCGAATATCACAGCCCGAATACTTCCGGTGCCTGCATCCAGCGCCATCAGGTCGTACTTTGATTCTGAAAGGGt                                                                                                                               >  1:443861/1‑250 (MQ=255)                                                                                                                  aaCGGCAGCGATATACTCCGGGGCTATGCCGGCGTTGTGCAGCGCCTGGCGCATACACTCACACGCCAGTTGCCAGTTTTTGTTGAGATCAAATTCCATAGAACCAGGAACGTCCGGTACTGCCAGATGCCGTCACTCCGCCTGTCCCACTGCTATTTGATTGCCTTCCAGGTCGAATATCACAGCCCGAATACTTCCGGTGCCTGCATCCAGCGCCATCAGGTAGTACTTTGATTCTGAAAGGGTAAaga                                                                                                                          <  1:306653/251‑1 (MQ=255)                                                                                                                    cGGCAGCGATATACTCCGGGGCTATGCCGGCGTTGTGCAGCGCCTGGCGCATACCCTCACACGCCAGTTGCCAGTTTTTGTTGAGATCAAATTCCATAGAACCAGGAACGTCCGGTACTGCCAGATGCCGTCACTCCGCCTGTCCCACTGCTATTTGATTGCCTTCCAGGTCGAATATCACAGCCCGAATACTTCCGGTGCCTGCATCCAGCGCCATCAGGTAGTACTTTGATTCTGAAAGGGTAAAGAGt                                                                                                                        <  1:389743/251‑1 (MQ=255)                                                                                                                           gATATACTCCGGGGCTATGCCGGCGTTGTGCAGCGCCTGGCGCATACACTCACACGCCAGTTGCCAGTTTTTGTTGAGATCAAATTCCATAGAACCAGGAACGTCCGGTACTGCCAGATGCCGTCACTCCGCCTGTCCCACTGCTATTTGATTGCCTTCCAGGTCGAATATCACAGCCCGAATACTTCCGGTGCCTGCATCCAGCGCCATCAGGTAGTACTTTGATTCTGAAAGGGTAAAGAGTCGAGcc                                                                                                                  >  1:120654/1‑250 (MQ=255)                                                                                                                           gATATACTCCGGGGCTATGCCGGCGTTGTGCAGCGCCTGGCGCATACACTCACACGCCAGTTGCCAGTTTTTGTTGAGATCAAATTCCATAGAACCAGGAACGTCCGGTACTGCCAGATGCCGCCACTCCGCCTGTCCCACTGCTATTTGATTGCCTTCCAGGTCGAATATCACAGCCCGAATACTTCCGGTGCCTGCATCCAGCGCCATCAGGTAGTACTTTGATTCTGAAAGGGTAAAGAGTCGAGcc                                                                                                                  >  2:137953/1‑250 (MQ=255)                                                                                                                                              ccGGCGTTGTGCAGCGCCTGGCGCATACACTCACACGCCAGTTGCCAGTTTTTGTTGAGATCAAATTCCAAAGAACCAGGAACGTCCGGTACTGCCAGATGCCGCCACTCCGCCTGTCCCACTGCTATTTGATTGCCTTCCAGGTCGAATATCACAGCCCGAATACTTCCGGTGCCTGCATCCAGCGCCATCAGGTAGTACTTTGATTCTGAAAGGGTAAAGAGTCGAGCCATTAGATTATCCTCGGCTTa                                                                                              <  1:172439/251‑1 (MQ=255)                                                                                                                                                              ccTGGCGCATACACTCACAAGCCAGTTTCCAGTTTTTGTTAAGCTCAAATTCCATAGAACCAGGAACGTCCGGTACTGCCAGATGCCGTCACTCCGCCTGTCCCACTGCTATTTGATTGCCTTCCCGGCCGAATATCACAGCCCGAATACTTCCGGTGCCTGCATCCAGCGCCCTCAGGAATTACTTTGATTCTGAAAGGTTAAAGAGACGAGCCATTAGATTATCCTCGGCTTATGCGTTCTGTCTGTTc                                                                              >  2:343445/1‑251 (MQ=255)                                                                                                                                                                  gcgcATACCCTCACACGCCAGTTGCCAGTTTTTGTTGAGATCAAATTCCATAGAACCAGGAACGTCCGGTACTGCCAGATGCCGTCACTCCGCCTGTCCCACTGCTATTTGATTGCCTTCCAGGTCGAATATCACAGCCCGAATACTTCCGGTGCCTGCATCCAGCGCCATCAGGTAGTACTTTGATTCTGAAAGGGTAAAGAGTCGAGCCATTAGATTATCCTCGGCTTATGCGTTCTGTCTGTTCctct                                                                          <  1:247999/251‑1 (MQ=255)                                                                                                                                                                      aTACACTCACACGCCAGTTGCCAGTTTTTGTTGAGATCAAATTCCATAGAACCAGGAACGTCCGGTACTGCCAGATGCCGTCACTCCGCCTGTCCCACTGCTATTTGATTGCCTTCCAGGTCGAATATCACAGCCCGAATACTTCCGGGGCCTGCATCCAGCGCCCTCAGGTAGTACTTTGATTCTGAACGGGGAAAGAGTCGAGCCATTAGATTATCCTCGGCTTATGCGTTCTGTCTGTTCCTCTATAc                                                                      >  2:477540/1‑251 (MQ=255)                                                                                                                                                                         cacTCACACGCCAGTTGCCAGTTTTTGTTGAGATCAAATTCCATAGAACCAGGAACGTCCGGTACTGCCAGATGCCGTCACTCCGCCTGTCCCACTGCTATTTGATTGCCTTCCAGGTCGAATATCACAGCCCGAATACTTCCGGTGCCTGCATCCAGCGCCATCAGGTAGTACTTTGATTCTGAAAGGGTAAAGAGTCGAGCCATTAGATTATCCTCGGCTTATGCGTTCTGTCTGTTCCTCTATACGtt                                                                   <  1:278502/251‑1 (MQ=255)                                                                                                                                                                          acTCACACGCCAGTTGCCAGTTTTTGTTGAGATCAAATTCCATAGAACCAGGAACGTCCGGTACTGCCAGATGCCGTCACTCCGCCTGTCCCACTGCTATTTGATTGCCTTCCAGGTCGAATATCACAGCCCGAATACTTCCGGTGCCTGCATCCAGCGCCATCAGGTAGTACTTTGATTc                                                                                                                                        <  1:442822/181‑1 (MQ=255)                                                                                                                                                                          acTCACACGCCAGTTGCCAGTTTTTGTTGAGATCAAATTCCATAGAAACAGGAACGTCCGGTACTGCCAGATGCCGTCACTCCGCCTGTCCCACTGCTATTTGATTGCCTTCCAGGTCGAATATCACAGCCCGAATACTTCCGGTGCCTGCATCCAGCGCCATCAGGTAGTACTTTGATTc                                                                                                                                        >  2:442822/1‑181 (MQ=255)                                                                                                                                                                             cacaCGCCAGTTGCCAGTTTTTGTTGAGATCAAATTCCATAGAACCAGGAACGTCCGGTACTGCCAGATGCCGTCACTCCGCCTGTCCCACTGCTATTTGATTGCCTTCCAGGTCGAATATCACAGCCCGAATACTTCCGGTGCCTGCATCCAGCGCCATCAGGTAGTACTTTGATTCTGAAAGGGTAAAGAGTCGAGCCATTAGATTATCCTCGGCTTATGCGTTCTGTCTGTTCCTCTATACGTTCTcc                                                               >  2:13183/1‑251 (MQ=255)                                                                                                                                                                                   ccagttgccagttTTTGTTGAGATCAAATTCCATAGAACCAGGAACGTCCGGTACTGCCAGATGCCGTCACTCCGCCTGTCCCACTGCTATTTGATTGCCTTCCAGGTCGAATATCACAGCCCGAATACTTCCGGTGCCTGCATCCAGCGCCATCAGGGAGTACTTTGATTCTGAAAGGGTAAAGAGTCGAGCCATTAGATTATCCTCGGCTTATGCGTTCTGTCTGTTCCTCTATACGTTCTCCATCAtt                                                         >  1:247689/1‑251 (MQ=255)                                                                                                                                                                                              tttttGTTGAGATCAAATTCCATAGAACCAGGAACGTCCGGTACTGCCAGATGCCGTCACTCCGCCTGTCCCACTGCTATTTGATTGCCTTCCAGGTCGAATATCACAGCCCGAATACTTCCGGTGCCTGCATCCAGCGCCATCAGGTAGTACTTTGATTCTGAAAGGGTAAAGAGTCGAGCCATTAGATTATCCTCGGCTTATGCGTTCTGTCTGTTCCTCTATACGTTCTCCATCATTCCCGGtaata                                               >  2:404073/1‑250 (MQ=255)                                                                                                                                                                                                 ttgttgAGATCAAATTCCATAGAACCAGGAACGTCCGGTACTGCCAGATGCCGTCACTCCGCCTGTCCCACTGCTATTTGATTGCCTTCCAGGTCGAATATCACAGCCCGAATACTTCCGGTGCCTGCATCCAGCGCCATCAGGTAGTACTTTGATTCTGAAAGGGTAAAGAGTCGAGCCATTAGATTATCCTCGGCTTATGCGTTCTGTCTGTTCCTCTATACGTTCTCCATCATTCCCGGTAATAAGGt                                           >  1:371099/1‑251 (MQ=255)                                                                                                                                                                                                           caaatccaaTAGACCCAGGACCGGCCGGTACGGCCAGATGCCGTCCCTCCGCCTGTCCCACTGCTATTTGATTGCCTTCCAGGTCGAATATCACAGCCCGAATACTTCCGGTGCCTGCATCCAGCGCCATCAGGTAGTACTTTGATTCTGAAAGGGTAAAGAGTCGAGCCATTAGATTATCCTCGGCTTATGCGTTCTGTCTGTTCCTCTATACGTTCTCCATCATTCCCGGTAATAAGGTCATGCAAAt                                  <  2:398336/242‑1 (MQ=255)                                                                                                                                                                                                                  cATAGAACCAGGAACGTCCGGTACTGCCAGATGCCGTCACTCCGCCTGTCCCACTGCTATTTGATTGCCTTCCAGGTCGAATATCACAGCCCGAATACTTCCGGTGCCTGCATCCAGCGCCATCAGGTAGTACTTTGATTCTGAAAGGGTAAAGAGTCGAGCCATTAGATTATCCTCGGCTTATGCGTTCTGTCTGTTCCTCTATACGTTCTCCATCATTCCCGGTAATAAGGTCATGCAAATTTAACTAc                          >  1:72235/1‑251 (MQ=255)                                                                                                                                                                                                                          cAGGAACGTCCGGTACTGCCAGATGCCGGCACTCCGCCTGTCCCACTGCTATTTGATTGCCTTCCAGGTCGAATATCACAGCCCGAATACTTCCGGTGCCTGCATCCAGCGCCATCAGGTAGTACTTTGATTCTGCAAGGGTAAAGAGTCGAGCCATTAGATTATCCTCGGCTTATGCGTTCTGTCTGTTCCTCTATACGTTCTCCATCATTCCCGGTAATAAGGTCATGCAAATTTAACTACGTAAAATc                  >  2:38435/1‑251 (MQ=255)                                                                                                                                                                                                                           aGGAACGTCCGGTACTGCCAGATGCCGTCACTCCGCCTGTCCCACTGCTATTTGATTGACTTCCAGGTCGAATATCACAGCCCGAATACTTCCGGTGCCTGCATCCAGCGCCATCAGGTAGTACTTTGATTCTGAAAGGGTAAAGAGTCGAGCCATTAGATTATCCTCGGCTTATGCGTTCTGTCTGTTCCTCTATACGTTCTCCATCATTCCCGGTAATAAGGTCATGCAAATTTAACTACGTAAAATcg                 >  1:263150/1‑251 (MQ=255)                                                                                                                                                                                                                                          tGCCAGATGCCGTCACTCCGCCTGTCCCCCTGCTATTTGATTGCCTTACAGGTCGAATATCACAGCCCGAATACTTCCGGTGCCTGCATCCAGCGCCATCAGGTAGTACTTTGATTCTGAAAGGGAAAAGAGTCGAGCCATTAGATTATCCTCGGCTTATGCGTTCTGTCTGTTCCTCTATACGTTCTCCATCATTCCCGGTAATAAGGTCATGCAAATTTAACTACGTAAAATCGCCGCTGCTGTGTCCt  >  1:29864/1‑251 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                               TCTTTAAGTTCGCTAACTTCGCGTGCCGCTCTGGCATCCACATTGGCGCAGGCCCAGATCGGGGCTCCTTCATTATTATATAAAACAATGCCTTCACGCATCGAACATGCCGAAACGGCAGCGATATACTCCGGGGCTATGCCGGCGTTGTGCAGCGCCTGGCGCATACACTCACACGCCAGTTGCCAGTTTTTGTTGAGATCAAATTCCATAGAACCAGGAACGTCCGGTACTGCCAGATGCCGCCACTCCGCCTGTCCCACTGCTATTTGATTGCCTTCCAGGTCGAATATCACAGCCCGAATACTTCCGGTGCCTGCATCCAGCGCCATCAGGTAGTACTTTGATTCTGAAAGGGTAAAGAGTCGAGCCATTAGATTATCCTCGGCTTATGCGTTCTGTCTGTTCCTCTATACGTTCTCCATCATTCCCGGTAATAAGGTCATGCAAATTTAACTACGTAAAATCGCCGCTGCTGTGTCCT  >  NC_000913/1599836‑1600319

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