Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,060,839	G→A	W131* (TGG→TGA)	yihL →	putative DNA‑binding transcriptional regulator YihL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,060,839	0	G	A	93.8%	90.8 / ‑5.8	32	W131* (TGG→TGA)	yihL	putative DNA‑binding transcriptional regulator YihL
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (13/17);  minor base T (1/1);  total (14/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-01

GCGTTGTTGAACCTCGAGGCAGAATCCAAAATCTTCCGTAAGGATCGTAAGGGGTGGTTCGTGACCCAGCCGCGATTTAATTACAGTCCGGAGCTGTCGGCGAGCTTTCAGCGGGCCGCAATTGAGCAAGGGCGAGAGCCTTCCTGGGGGTTTACCGAGAAAAACCGTACCAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGAAAAGTGCTACCAAAAAGAGACGGTAGTAAAAAAATTGATTTTCAAACCCGTCAGAATGCCGGGCGATATCAGCAAGTATCTTGGCGGTTCTGCGGGTATGCCAGCGATCTTAATTGAAAAGCATCGC  >  NC_000913/4060594‑4061076                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                              gCGTTGTTGAACCTCGAGGCAGAATCCAAAATCTTCCGTAAGGATCGTAAGGGGTGGTTCGTGACCCAGCCGCGATTTAATTACAGTCCGGAGCTGTCGGCGAGCTTTCAGCGGGCCGCAATTGAGCAAGGGCGAGAGCCTTCCTGGGGGTTTACCGAGAAAAACCGTACCAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCCTCAACTGAACTCTATCGCATCACCGGCTGAGGGgc                                                                                                                                                                                                                                          >  2:415737/1‑251 (MQ=255)           aacTCGAGGCAGAATACAAAATCTTACGTAAGGATCGTAAGGGGTGGTTCGTGACTCAGCCGCGATTTAATTACAGTCCGGAGCTGTCGGCGAGCTTTCAGCGGGCCGCAATTGAGCAAGGGCGAGAGCCTTCCTGGGGGTTTACCGAGAAAAACCGTACCAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGTCATCAACTGAACTCTATCGCATCACCGTCTGAGGGGCGCTGGAAGGa                                                                                                                                                                                                                                <  2:282008/249‑1 (MQ=255)                   gCAGAATCCAAAATCTTCCGTAAGGATCGTAAGGGGTGGTTCGTGACCCAGCCGCGATTTAATTACAGTACGGAGCTGTCGGCGAGCTTTCAGCGGGCCGCAATTGAGCAAGGGCGAGAGCCTTCCTGGGGGTTTACCGAGAAAAACCGTACCAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGAGGGGCGCTGGAAGGACATAAAGt                                                                                                                                                                                                                        >  1:275340/1‑251 (MQ=255)                       aaTCCAAAATCTTCCGTAAGGATCGTAAGGGGTGGTTCGTGACCCAGCCGCGATTTAATTACAGTCCGGAGCTGTCGGCGAGCTTTCAGCGGGCCGCAATTGAGCAAGGGCGAGAGCCTTCCTGGGGGTTTACCGAGAAAAACCGTACCAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGAGGGGCGCTGGAAGGACATAAAGTTTTc                                                                                                                                                                                                                    >  1:334007/1‑251 (MQ=255)                           cAAAATCTTCCGTAAGGATCGTAAGGGGTGGTTCGTGACCCAGCCGCGATTTAATTACAGTCCGGAGCTGTCGGCGAGCTTTCAGCGGGCCGCAATTGAGCAAGGGCGAGAGCCTTCCTGGGGGTTTACCGAGAAAAACCGTACCAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGAGGGGCGCTGGAAGGACATAAAGTTTTCTATc                                                                                                                                                                                                                <  1:183895/251‑1 (MQ=255)                                        aaGGATCTTAAGGGGTGGTTCGTGACCCAGCCGCGATTTAATTACAGTCCGGAGCGGTCGGCGAGCTTTCAGCGGGCCGCAATTGAGCAAGGGCGAGAGCCTTCCTGGGGGTTTACCGAGAAAAACCGTACCAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGAGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACatatat                                                                                                                                                                                                   <  1:3241/251‑1 (MQ=255)                                                 aaGGGGTGGTTCGTGACCCAGCCGCGATTTAATTACAGTCCGGAGCTGTCGGCGAGCTTTCAGCGGGCCGCAATTGAGCAAGGGCGAGAGCCTTCCTGGGGGTTTACCGAGAAAAACCGTACCAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTAGCGCATCACCGGCTGAGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGa                                                                                                                                                                                          <  1:388617/251‑1 (MQ=255)                                                                  ccAGCCGCGATTTAATTACAGTCCGGAGCTGTCGGCGAGCTTTCCGCGGGCCGCAATTGAGCAAGGGCGAGAGCCTTCCTGGGGGTTTACCGAGAAAAAACGTACCAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGAGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCtttt                                                                                                                                                                         >  2:303634/1‑250 (MQ=255)                                                                  ccAGCCGCGATTTAATTACAGTCCGGAGCTGTCGGCGAGCTTTCAGCGGGCCGCAATTGAGCAAGGGCGAGAGCCTTCCTGGGGGTTTACCGAGAAAAACCGTACCAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGAGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTa                                                                                                                                                                         <  1:415737/251‑1 (MQ=255)                                                                                     aGTCCGGAGCTGTCGGCGAGCTTTCAGCGGGCCGCAATTGAGCAAGGGCGAGAGCCTTCCTGGGGGTTTACCGAGAAAAACCGTACCAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGAGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAAcc                                                                                                                                                       >  2:101922/1‑250 (MQ=255)                                                                                       tCCGGAGCTGTCGGCGAGCTTTCAGCGGGCCGCAATTGAGCAAGGGCGAGAGCCTTCCTGGGGGTTTACCGAGAAAAACCGTACCAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGAGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTCTTGAACAACTTGAAAACCACt                                                                                                                                                    >  1:404082/1‑251 (MQ=255)                                                                                         cGGAGCTGTCGGCGAGCTTTCCGCGGGCCGCAATTGAGCAAGGGCGAGATCCTCTCTGGGGGTTTACCGCGAAAAACCGTACCAGCGATCTTCCCGAGACGGTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTTAGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGACCAACTTGAAAACCACTCa                                                                                                                                                  <  1:269898/251‑1 (MQ=255)                                                                                              cggTCGGCGAGCTTGCAGCGGGCCGCAATTGAGCAAGGGCGAGAGCCTTCCTGGGGGTTTACCCAGAAAAACCGTACCAGTGATATTCCCGAGACGCTCGCGCCACTGATTGCTGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGAGGGGCGTTGGAAGGTCATAAAGTTTTCTATCAAGAAACATATATTTATCCTGAAGTTGCTCCGGGCTTTATTGAACTACTTGAAAACCACTCATTTTc                                                                                                                                             <  1:24190/249‑1 (MQ=255)                                                                                                   ggCGAGCTTTCAGCGGGCCGCAATTGAGCAAGGGCGAGAGCCTTCCTGGGGGTTTACCGAGAAAAACCGTACCAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGAGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAg                                                                                                                                        <  1:303634/251‑1 (MQ=255)                                                                                                                        aaTTGAGCAAGGGCGAGAGCCTTCCTGGGGGTTTACCGAGAAAAACCGTACCAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGAGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACATCATGAAAACCACTCATTGTCTGCAGTCTGGGAAAAGTGCTACCaaa                                                                                                                   <  2:124347/251‑1 (MQ=255)                                                                                                                                aaGGGCGAGAGCCTTCCTGGGGGTTTACCTAGAAAAACCGGACCAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCTCATCACCGTCTGTGGTGCGCTTGAAGGACATAAAGTGGTCTATCACGAAACATATATTAATCCTTAAGTTGCTCCGGGCGTTATTGAACAACTTGAAAACCACGCATTTTCGGCAGTCTGGGAAAAGTGCTACCAAAAAGAGACg                                                                                                           <  2:412402/251‑1 (MQ=255)                                                                                                                                aaGGGCGAGAGCCTTCCTGGGGGTTTACCGAGAAAAACCGTACCAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGAGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGAAAAGTGCTACCAAAAAGAGACg                                                                                                           >  1:159990/1‑251 (MQ=255)                                                                                                                                       atagCCTTCCTGGGGGTTTACCGAGAAAAACCGTACCAGCGATATTCCCGAGACGCTCTCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGAGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGAAAAGTGCTTCCAAAAAGAGACGGTAGTaa                                                                                                    <  2:123074/249‑1 (MQ=255)                                                                                                                                         agCCTTCCTGGGGGTTTACCGAGAAAAACCGTACCAGCGATATTCCCGCGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGAGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATCAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGAAAAGTGCTACCAAAAAGAGACGGTAGTaaa                                                                                                   >  2:42948/1‑250 (MQ=255)                                                                                                                                                           ccGAGAAAAACCGTACCAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGAGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGAAAAGTGCTACCAAAAAGAGACGGTAGTAAAAAAATTGATTTTCAAAccc                                                                                <  1:5906/251‑1 (MQ=255)                                                                                                                                                               gaAAAACCGTACCAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGAGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGAAAAGTGCTACCAAAAAGAGACGGTAGTAAAAAAATTGATTTTCAAACCCGTCa                                                                            <  1:355445/251‑1 (MQ=255)                                                                                                                                                                    aCCGTAACAGCGATATTCCCGAGACGCTCGCGCCACTGATTGAAGTGACGCCATAAACTGAACTCTATCGCATCACCGGCTGTGGGGCGATGGAAGGACATAAAGTTTTCTAACACGAAACTTATATTAAACCTGAAGTTGCTCGGGGCTTTATGGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGAAAAGTGCTACCAAAACTAGACGGAAGTAAACAAATTGATTTTCAAACCCGTAAGAATg                                                                       >  2:78533/1‑251 (MQ=255)                                                                                                                                                                     ccGTACCAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGAGGGGCGCTGGAAGGACATAAAGATTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGAAAAGTGCTACCAAAAAGAGACGGTAGTCAAAAAATTGATTTTCAAACCCGCCAGAATGc                                                                      >  2:23524/1‑251 (MQ=255)                                                                                                                                                                           cAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGAGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGAAAAGTGCTACCAAAAAGAGACGGTAGTAAAAAAATTGATTTTCAAACCCGTCAGAATGCCGGGCg                                                                >  2:44526/1‑251 (MQ=255)                                                                                                                                                                                         agaCGCTCGCGCCAATGATTGCAGTGCCGCCATCAACTGAACTCTATCGCCTCACCGGCTGAGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGAAAAGTGCTACCAAAAAGAGACGGTAGTAAAAAAATTGATTTTCAAACCCGTCAGAATGCCGGGCGATATCAGCAAGTAt                                                  >  2:193332/1‑251 (MQ=255)                                                                                                                                                                                                   gcCACTGATTGCAGTGACGCCCTCAACTGAACTCTATCGCATCACCGGCTGAGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGAAAAGTGCTACCAAAAAGAGACGGTAGTAAAAAAATTGATTTTCAAACCCGTCAGAATGCCGGGCGATATCAGCAAGTATCTTGGCGGtt                                        >  1:136682/1‑251 (MQ=255)                                                                                                                                                                                                    ccACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGAGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGAAAAGTGCTACCAAAAAGAGACGGTAGTAAAAAAATTGATTTTCAAACCCGTCAGAATGCCGGGCGATATCAGCAAGTATCTTGGCGGTTc                                       <  1:24133/251‑1 (MQ=255)                                                                                                                                                                                                         gATTGCAGTGACGCCATAAACTGAACTCTATCGCATCACCGGCTGAGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATGAATCCTGAAGTTGCTCCGGGCTTTATTGAACAGCTTGAAAACCACTCATTTTCTGCAGTCTGGGAAAAGTGCTAGCAAAAAGAGACGGTAGTGAAAAAATTGATTTCCAAACCCGTGAGAATGCCGGGCGAGATCAGcttttcgcttggctgttctgcgg                                  <  2:26117/251‑23 (MQ=255)                                                                                                                                                                                                          aTTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGAGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGAAAAGTGCTACCAAAAAGAGACGGTAGTAAAAAAATTGATTTTCAAACCCGTCAGAATGCCGGGCGATATCAGCAAGTATCTTGGCGGTTCTGCggg                                 <  2:343769/251‑1 (MQ=255)                                                                                                                                                                                                          aTTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGAGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGAAAAGTGCTACCAAAAAGAGACGGTAGTAAAAAAATTGATTGTCAAACCCGTCAGAATGCCGGGCGATATCAGCAAGTATCTTGGCGGTTCTGCggg                                 <  2:38383/251‑1 (MQ=255)                                                                                                                                                                                                                      ccATCAACGGAACTCTATCGCATCTCCTGCTGATGGGCGCTGGAATGACATAAAGTTTTCGATAACGAAACATATATTAAGCCTGAATTTGCCCCGGGCTTTATTGAACAACTTGAAAACCACTGAGTTTCTTCAGTCTGTGAAAAGTGCTACCAAAAAGAGACTGTAGTAAAAAAGTTGATTGTCAATCCCGTCAGAATGCCGGGCGATATCAGCAATTATCGGGGCTGTTCTGTGGATATGCTAGCGAt                     <  2:356416/251‑1 (MQ=255)                                                                                                                                                                                                                                         gCATCACCGGCTGAGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGAAAAGTGCTACCAAAAAGAGACGGTAGTAAAAAAATTGATTTTCAAACCCGTCAGAATGCCGGGCGATATCAGCAAGTATCTTGGCGGTTCTGCGGGTATGCCAGCGACCTTAATTGAAAAGCATcgc  >  2:235483/1‑251 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                              GCGTTGTTGAACCTCGAGGCAGAATCCAAAATCTTCCGTAAGGATCGTAAGGGGTGGTTCGTGACCCAGCCGCGATTTAATTACAGTCCGGAGCTGTCGGCGAGCTTTCAGCGGGCCGCAATTGAGCAAGGGCGAGAGCCTTCCTGGGGGTTTACCGAGAAAAACCGTACCAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGAAAAGTGCTACCAAAAAGAGACGGTAGTAAAAAAATTGATTTTCAAACCCGTCAGAATGCCGGGCGATATCAGCAAGTATCTTGGCGGTTCTGCGGGTATGCCAGCGATCTTAATTGAAAAGCATCGC  >  NC_000913/4060594‑4061076

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