Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,419,414	C→T	G15S (GGT→AGT)	ydaG ←	uncharacterized protein YdaG

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,419,414	0	C	T	96.9%	91.9 / NA	33	G15S (GGT→AGT)	ydaG	uncharacterized protein YdaG
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (17/14);  total (17/16)
Fisher's exact test for biased strand distribution p-value = 4.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGTATACGTCATGGCAGTGAACGTTCCATCCTGGTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAGCGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGCAACAATCCATCAAATACTGAACTACTTCGTAATGCACCATCTTATCTCTCCCCTTAACGCCGGGTGGCGGAACTAACTGCTGCACTGCAAAATTTGAATCCCGCCGTCATGTTCATACGCCTCGGGCTGGCTACTTAACCCCTTACCACTGCCTGGTAACTCGAAGTATTGCCCGGCGTTCTGTGGGGCGGGGTGGGTGGTATGCTGGAACTATAGGTAATGCCTAATTGATTGTC  >  NC_000913/1419184‑1419652                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                               cGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGTATACGTCATGGCAGTGAACGTTCCATCCTGGTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAGCGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACTGCAACAATCCATCAAATACt                                                                                                                                                                                                                            <  1:261512/251‑1 (MQ=255)        cACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGTATACGTCATGGCAGTGAACGTTCCATCCTGGTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAGCGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACTGCAACAATCCATCAAATACTGAactac                                                                                                                                                                                                                     >  1:423599/1‑251 (MQ=255)              tttCGGTTTTAAATGTTTTGCTTTTGGTATACGTCATGGCAGTGAACGTTCCATCCTGGTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAGCGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACTGCAACAATCCATCAAATACTGAACTACTTCGt                                                                                                                                                                                                                >  2:256062/1‑250 (MQ=255)                tCGGTTTTAAATGTTTTGCTTTTGGTATACGTCATGGCAGTGAACGTTCCATCCTGGTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAGCGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACTGCAACAATCCATCAAATACTGAACTACTTCGTaa                                                                                                                                                                                                              <  1:369586/250‑1 (MQ=255)                         aaTGTTTTGCTTTTGGTATACGTCATGGCAGTGAACGTTCCATCCTGGTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAGCGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACTGCAACAATCCATAAAATACTGAACTACTTCGTAATGCACCATCt                                                                                                                                                                                                    >  2:218157/1‑251 (MQ=255)                            gTTTTGCTTTTGGTATACGTCATGGCAGTGAACGTTCCATCCTGGTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAGCGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACTGCAACAATCCATCAAATACTGAACTACTTCGTAATGCACCATCTTAt                                                                                                                                                                                                 >  1:42208/1‑251 (MQ=255)                                      tGGTATACGTCATGGCAGTGAACGTTCCATCCTGGTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAGCGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACTGCAACAATCCATCAAATACTGAACTACTTCGTAATGCACCATCTTATCTCTCCCCtt                                                                                                                                                                                       >  2:220244/1‑251 (MQ=255)                                             cGTCATGGCAGTGAACGTTCCATCCTGGTTTGGGAACACGTCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAGCGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACTGCAACAATCCATCAAATACTGAACTACTTCGTAATGCACCATCTTATCTCTCCCCTTAACGCCg                                                                                                                                                                                <  1:537338/251‑1 (MQ=255)                                                cATGGCAGTGAACGTTCCATCCTGTTTGGGGAACACGCCGCCCACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCAGTTTGCGAATCTCACATCTTGTTGCTACGTATAGCGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACTGCAACAATCCATCAAATACTGAACTACTTCGTAATGCACCATCTTATCTCTCCCCTTAACGCCggg                                                                                                                                                                              <  2:42208/250‑1 (MQ=255)                                                   ggcggTGAACGTTCCATGCTGGTTGGGGTACACGCCGAACACCAGGGATTCGTTGTTGCCGAGGTCGAGTTTTTGCATTGTGCGAATCTCACATCTGGTTGCTACGTATAGCGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTTGTTGTAAGTGATACTGCAACAGTCCATCAAATACTGAACTACTTCGTAATGCTCCATCTTATCTCTCCCCTTAACGCAGGGTgg                                                                                                                                                                           <  2:546278/246‑1 (MQ=255)                                                    gCAGTGAACGTTCCATCCGTGTGGGGGAACGCGCCTCACCTCAGGGATGCGTTGTTGCCGAGGTCGATTTTTTGGAATTTGCGAATCTCACATCTTGTTGCTACGTATAGCGCCTTCTGCCTGCCAGAGATCCCATTCGTTGCTGCGTAAAGTCTGCACAGCCTGGTTGTAAGTGATACTGCAACAATCCATCAAATACGGAACTACTTCGTAATGCACCCTCTTATCTCTACCCTTAGCGCCGGGTgtcg                                                                                                                                                                         <  2:212894/251‑4 (MQ=255)                                                       gTGAACGTTCCATCCTGGTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAGCGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACTGCAACAATCCATCAAATACTGAACTACTTCGTAATGCACCATCTTATCTCTCCCCTTAACGCCGGGTGGCGGa                                                                                                                                                                       >  1:395769/1‑250 (MQ=255)                                                                       ggTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAGCGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACTGCAACAATCCATCAAATACTGAACTACTTCGTAATGCACCATCTTATCTCTCCCCTTAACGCCGGGTGGCGGAACTAACTGctgcactgc                                                                                                                                                      <  1:256062/251‑1 (MQ=255)                                                                         ttGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAGCGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACTGCAACAATCCATCAAATACTGAACTACTTCGTAATGCACCATCTTATCTCTCCCCTTAACGCCGGGTGGCGGAACTAACTGctgcactgca                                                                                                                                                     >  1:16457/1‑250 (MQ=255)                                                                                  acGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAGCGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACTGCAACAATCCATCAAATACTGAACTACTTCGTAATGCACCATCTTATCTCTCCCCTTAACGCCGGGTGGCGGAACTAACTGCTGCACTGCAAAATTTGAAt                                                                                                                                           >  1:144499/1‑251 (MQ=255)                                                                                          cacCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAGCGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACTGCAACAATCCATCAAATACTGAACTACTTCGTAATGCACCATCTTATCTCTCCCCTTAACGCCGGGTGGCGGAACTAACTGCTGCACTGCAAAATTTGAATCCCGCCGt                                                                                                                                   >  1:414881/1‑251 (MQ=255)                                                                                                   ttCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAACCTCACATCTTGTTGCTCCGTATAGCGACTTCTGCGTGCCAGAGAACCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTGGTAAGTGATACTGCAACAATCCATCAAATACTGAACTTCTTCGTAATGCCCCATCTTATCTCTCCCCTTAACGCCGGGTGGCGGAACTAACTGCTGCCCTGCAAATCTTGAATCCCGCCGTCATGTTCa                                                                                                                           <  2:467030/250‑1 (MQ=255)                                                                                                                  gTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAGCGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACTGCAACAATCCATCAAATACTGAACTACTTCGTAATGCACCATCTTATCTCTCCCCTTAACGCCGGGTGGCGGAACTAACTGCTGCACTGCAAAATTTGAATCCCGCCGTCATGTTCATACGCCTCGggctgg                                                                                                            >  1:366830/1‑250 (MQ=255)                                                                                                                   ccGATTTTTTGCATTTTGCGAATCTCTCATCTTGTTGCTTCGTATAGCGACTTCTGCCTGCCAGAGATCCCATTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATCCTGCAACAATCCACCAAATACGGAACTTCTTCGTAATGCACCATCTTATCTCTCCCCTTAACGCCGGGTGGCGGAACTAACTGCTGCACTGCAAAATTTGAATCCCGCCGTCATGTTCATACGCCTCGggctggct                                                                                                          <  2:172599/250‑1 (MQ=255)                                                                                                                                              caTCTTGTTGCTACGTATAGCGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACTGCAACAATCCATCAAATACTGAACTACTTCGTAATGCACCCTCTTATCTCTCCCCTTAACGCCGGGTGGCGGAACTAACTGCTGCACTGCTAAATTTGAATCACGCCGTCATGTTCATACGCATCGTGCTGGCTACTTAACCCCTTACCACTGCCTGGTa                                                                                >  2:151345/1‑250 (MQ=255)                                                                                                                                                          acggataGCGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACTGCAACAATCCATCAAATACTGAACTACTTCGTAATGCACCATCTTATCTCTCCCCTTAACGCCGGGTGGCGGAACTAACTGCTGCACTGCAAAATTTGAATCCCGCCGTCATGTTCATACGCCTCGGGCTGGCTACTTAACCCCTTACCACTGCCTGGTAACTCGACGTAttt                                                                   >  2:531834/5‑250 (MQ=255)                                                                                                                                                            gTATAGCGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACTGCAACAATCCATCAAATACTGAACTACTTCGTAATGCACCAGCTTATCTCTCCCCTTAACGCCGGGTGGCGGAACTAACTGCTGCACTGCAAAATTTGAATCCCGCCGTCATGTTCATACGCCTCGGGCTGGCTACTTAACCCCTTACCACTGCCTGGTAACTCGAAGTATTGcc                                                                 <  1:348416/251‑1 (MQ=255)                                                                                                                                                              ataGCGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACTGCAACAATCCATCAAATACTGAACTACTTCGTAATGCACCATCTTATCTCTCCCCTTAACGCCGGGTGGCGGAACTAACTGCTGCACTGCAAAATTTGAATCCCGCCGTCATGTTCATACGCCTCGGGCTGGCTACTTAACCCCTTACCACTGCCTGGTAACTCGAAGTATTGccc                                                                <  1:220244/250‑1 (MQ=255)                                                                                                                                                              ataGCGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACTGCAACAATCCATCAAATACTGAACTACTTCGTAATGCACCATCTTATCTCTCCCCTTAACGCCGGGTGGCGGAACTAACTGCTGCACTGCAAAATTTGAATCCCGCCGTCATGTTAATACGCCTCGGGCTGGCTACTTAACCCCTTACCACTGCCTGGTAACTCGAAGAATTGCCCg                                                               >  2:451916/1‑251 (MQ=255)                                                                                                                                                                       cctgcctgcCAGGGAACCCAGTCGTCGCTGCTTAAAGCCTGCACAGCCTGGTTGTAAGTGTTACCGCAACAACCCATCAAATACTGACGTGCTTCGTAATGCACCAGCTTCTCTCCCCCCGTAACGCCGGGTGGCGGACCTACCTGCTGCTCTGCAAAATTTGAATCCCGCCGTCATGTTCATACGCCTCGGGCTGGCTACTTAACCCCTTACCACTGCCTGGTAACTCGAAGTATTGCCCGGCGTTCtgt                                                      <  2:144499/250‑1 (MQ=255)                                                                                                                                                                        ctgcctgcCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACTGCAACAATCCATCAAATACTGAACTACTTCGTAATGCACCATCTTATCTCTCCCCTTAACGCCGGGTGGCGGAACTAACTGCTGCACTGCAAAATTTGAATCCCGCCGTCATGTTCATACGCCTCGGGCTGGCTACTTAACCCCTTACCACTGCCTGGTAACTCGAAGTATTGcc                                                                 >  2:159236/1‑239 (MQ=255)                                                                                                                                                                        ctgcctgcCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACTGCAACAATCCATCAAATACTGAACTACTTCGTAATGCACCATCTTATCTCTCCCCTTAACGCCGGGTGGCGGAACTAACTGCTGCACTGCAAAATTTGAATCCCGCCGTCATGTTCATACGCCTCGGGCTGGCTACTTAACCCCTTACCACTGCCTGGTAACTCGAAGTATTGcc                                                                 <  1:159236/239‑1 (MQ=255)                                                                                                                                                                         tgcctgccAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACTGCAACAATCCATCAAATACTGAACTACTTCGTAATGCACCATCTTATCTCTCCCCTTAACGCCGGGTGGCGGAACTAACTGCTGCACTGCAAAATTTGAATCCCGCCGTCATGTTCATACGCCTCGGGCTGGCTACTTAACCCCTTACCACTGCCTGGTAACTCGAAGTATTGCCCGGCGTTCTGTgg                                                    >  1:57845/1‑251 (MQ=255)                                                                                                                                                                                        ccAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACTGCAACAATCCATCAAATACTGAACTACTTCGTAATGCACCATCTTATCTCTCCCCTTAACGCCGGGTGGCGGAACTAACTGCTGCACTGCAAAATTTGAATCCCGCCGTCATGTTCATACGCCTCGGGCTGGCTACTTAACCCCTTACCACTGCCTGGTAACTCGAAGTATTGCCCGGCGTTCTGTGGGGCGGGGTGggtggt                                     >  1:351637/1‑251 (MQ=255)                                                                                                                                                                                        ccAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACTGCAACAATCCATCAAATACTGAACTACTTCGTAATGCACCATCTTATCTCTCCCCTTAACGCCGGGTGGCGGAACTAACTGCTGCACTGCAAAATTTGAATCCCGCCGTCATGTTCATACGCCTCGGGCTGGCTACTTAACCCCTTACCACTGCCTGGTAACTCGAAGTATTGCCCGGCGTTCTGTGGGGCGGGGTGggtggt                                     <  1:218157/251‑1 (MQ=255)                                                                                                                                                                                             cGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACTGCAACAATCCATCAAATACTGAACTACTTCGTAATGCACCATCTTATCTCTCCCCTTAACGCCGGGTGGCGGAACTAACTGCTGCACTGCAAAATTTGAATCCCGCCGTCATGTTCATACGCCTCGGGCTGGCTACTTAACCCCTTACCACTGCCTGGAAACTCGAAGTATTGCCCGGCGTTCTGTGGGGCGGGGGGGGTGGTATGCt                                >  2:285178/1‑251 (MQ=255)                                                                                                                                                                                                    gtttAAAGCCTGCCCAGCTGGGTTGTAGGTGATTCGTCAACAATCCACAAAACGCGAAACTTCTTCGTAATGCACTATGTGATTCCTCCCCTTAACGCCGGGTGGCGTAACTAACTGCGGCACTGCAAATTTTGAACCCCGCCGTCATGTTCATACGCCTCGGGCTGGCTACTTTACCCCTTACCACTGCCTGGTAACTCGAAGTATTGCCCGGCGTTCTGTGGGGCGGGGTGGGTGGTATGCTGGAACta                         <  2:395769/248‑1 (MQ=255)                                                                                                                                                                                                                            gTAAGTGATACTGCAACAACCCATCAAATACTGAACTACTTCGTAATGCACCATCTTATCTCCCCCCTTAACGCCGGGTGGCGGAACTAACTGCTGCACTGCAAAATTTGAATCCCGCCGTCATGTTCATACGCCTCGGGCTGGCTACTTAACCCCTTACCACTGCCTGGTAACTCGAAGTATTGCCCGGCGTTCTGTGGGGCGGGGTGGGTGGTATGCTGGAACTATAGGTAATGCCTAATTGATTGTc  <  1:490579/250‑1 (MQ=255)                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                               CGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGTATACGTCATGGCAGTGAACGTTCCATCCTGGTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAGCGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGCAACAATCCATCAAATACTGAACTACTTCGTAATGCACCATCTTATCTCTCCCCTTAACGCCGGGTGGCGGAACTAACTGCTGCACTGCAAAATTTGAATCCCGCCGTCATGTTCATACGCCTCGGGCTGGCTACTTAACCCCTTACCACTGCCTGGTAACTCGAAGTATTGCCCGGCGTTCTGTGGGGCGGGGTGGGTGGTATGCTGGAACTATAGGTAATGCCTAATTGATTGTC  >  NC_000913/1419184‑1419652

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