Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,028,546	G→A	A156V (GCC→GTC)	ygfS ←	putative electron transport protein YgfS

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,028,546	0	G	A	96.4%	72.7 / NA	28	A156V (GCC→GTC)	ygfS	putative electron transport protein YgfS
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (18/9);  minor base C (0/1);  total (18/10)
Fisher's exact test for biased strand distribution p-value = 3.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCTGGCCTGGTGTGTGCCGTGTTCATCGTTCTGACTTGTAATATCGTTACCGGTATTATGCTGGGCTTTGTGACACTGGTCGTAGGCCGCGTCTTTGCACGCGAATGGCAAAAGCTGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGGCAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGTGATCATTCCAAACGGGCACGCACTGACGCAACTCTGACAGCCAATACAACGGGCAGAATTGGTCTGCACCACCT  >  NC_000913/3028305‑3028783                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                              ggctggcCTGGTGTGTGCCGTGTTCATCGTTCTGACTTGTAATATCGTTACCGGTATTATGCTGGGCTTTGTGACACTGGTCGTAGGCCGCGTCTTTGCACGCGAATGGCAAAAGCTGATTATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGCTTCTCACCGCTGACAACAATAc                                                                                                                                                                                                                                      >  1:37250/1‑251 (MQ=255)                       ttCATCGTTCTGACTTTTAATATCGTTCCCGGTAGTATGCTGGGCTTTGTGACACGGGTCGTAGGCCTCGTCTTTGCACGCGAATGGCAAAAGCTGAATATTGGTACGGTGATCAGTACTGCCGCACTGGTCGCAGTTTACGCGGGTGGGTTGTCAATCTAATATATTCTCCGCGCCTTGCTGCGCGTATATTTCGTGCAGAGAGGATTCTTTCCGCTGACAACAATCCGCCGCTggtgaacttccctga                                                                                                                                                                                                                 <  2:386750/250‑14 (MQ=255)                                       gTAATATCGTTACCGGTATTATGCTGAGCTTTGTGACACTTGCCGTAGGCCGCGTCTTTGAACGCGAATGGCGAACGCTGACAATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGCTTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGCTTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGc                                                                                                                                                                                                >  2:247408/1‑251 (MQ=255)                                         aaTATCGTTACCGGGATTATGCTGGGCTTTGTGACACTGGTCGTAGGCCGCGTCTTTGCACGCGAATGGCAAAAGCTGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCa                                                                                                                                                                                               >  2:315746/1‑250 (MQ=255)                                                                              tGGTCGTAGGCCAAGTCTTTCCACCAGAATGGCAAAAGCTGTATATTGGTACGGTGATCATTACTGCCGAAATGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCACGATTCAACgaagcct                                                                                                                                                         >  2:396541/1‑245 (MQ=255)                                                                                         cgcgTCTTTGCACGCGAATGGCAAAAGCTGAATATTTGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTc                                                                                                                                              <  1:7204/251‑1 (MQ=255)                                                                                              cTTTGCACGCGACTGGCAAAAGCTGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCt                                                                                                                                         >  2:324431/1‑251 (MQ=255)                                                                                                      gcgAATCGCCAAAGCTGAATATTGGTGCGGTGATCATGACTGCCGCCCTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGCATCTCCGCGCCTTTCAGCGCGGTGATTTCTTTCAGAGAGGATTCTCACCGCTGACCACAATACGCTGCTGGAGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCCAGATTCAACGCCGGCTGGCCCCTCTTCCCGCTGTTcac                                                                                                                                   >  2:173414/1‑249 (MQ=255)                                                                                                                   gCTGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCCTTTc                                                                                                                     >  2:182863/1‑250 (MQ=255)                                                                                                                              gTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTtgttgt                                                                                                         <  1:362295/251‑1 (MQ=255)                                                                                                                               tACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTc                                                                                                        >  1:177525/1‑251 (MQ=255)                                                                                                                                                           gCATTTTACGCAGGTGGATGGGCAATCTACTAGATTCTACGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACAGCTGACAACAATACGCTGCTGGCGCACTCGCCTGGGTTCTCTTTCAATGAGCAACTGCAGCGCCTGGGTTCGGCAATATTCAACGAACGCTGGCCCCTCTACCCGCTGTTCACAGCGACCGCATTTAAAAAAGTGTTTTCGGCTATCCCCCCGCAACGACTGaagg                                                                              >  2:285028/1‑247 (MQ=255)                                                                                                                                                                 tACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGTGATCa                                                                       <  1:263010/250‑1 (MQ=255)                                                                                                                                                                  aCGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGTGATCAtt                                                                     >  1:211907/1‑251 (MQ=255)                                                                                                                                                                    gcgGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGTGATCATTcc                                                                   >  1:395639/1‑251 (MQ=255)                                                                                                                                                                          ggtTGGACAATCTAATAGATTCTCCGCGCATTTCAGCGCGGAGATTTCTTTCAGAGAGGCTTCTCACCGCTGACACCAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGCATGGTGATCATTCCACACgg                                                             >  2:169939/1‑251 (MQ=255)                                                                                                                                                                                       aaTAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGTGATCATTCCAAACGGGCACGCACTg                                                   <  1:19547/248‑1 (MQ=255)                                                                                                                                                                                             ttctcagagcCTTTAAGAGCGGAGATTTCTTTCAGAGAGGATTATCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGAGATCATTCCAAACGGGCACGCACTGACGCAActc                                          >  1:114957/9‑251 (MQ=255)                                                                                                                                                                                                 ccGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGTGATCATTCCAAACGGGCACGCACTGACGCAACTCTg                                        <  1:225142/249‑1 (MQ=255)                                                                                                                                                                                                         ttCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCACGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGTGATCATTCCAAACGGGCACGCACTAACGCAACTCTGACAGCCAAt                               >  1:26037/1‑250 (MQ=255)                                                                                                                                                                                                          tCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACATCAATAAGCTGCTGGCGGACTCGCCCGAGTTCTCTTTCAGTAAGCACCTGCAGCACCTGTGTTGAGCACGATTCAACTCAGACTGGACCCTCTTCCCGCTGTTCACAGCGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGCCAACGACTAACTGGTGATCATTCCACACGGGCAAGCCCTGACGCAAACCTGCCAGCCAATa                              >  1:400967/1‑250 (MQ=255)                                                                                                                                                                                                          tCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGCTTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGTGATCATTCCAAACGGGCACGCACTGACGCAACTCTCACAGCCAATAc                             >  1:229926/1‑251 (MQ=255)                                                                                                                                                                                                          tCAGCGCGGAGATTTCTTGCAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCTGACTCGCCTGAGTTCTCTGTCAGTAAGCAACTGCAGCTCCTGTGTTGGGCAAGATTCAACTCAGGCTGTCCCCGCTTCCCGCTGTTCACAGAGAGCGCATTTCACAATTTGTTGACGGGTATCCCCCGGCAACGACTGAATGGTTATCATTCCAAACGGGCACGCACTGACGCAACTCTGACAGTCAATa                              <  2:201777/250‑1 (MQ=255)                                                                                                                                                                                                           cAGCGCTGAGATTTCTTTCATAGAGGATTATCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCGCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTTGCCCCGTTTCCCGCTCTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGTGATCATTCCAAACGGGCTCGCCCTGAGGCAATTCGGACAGCCAATAc                             <  2:5322/250‑1 (MQ=255)                                                                                                                                                                                                             gcgcGGAGATTTCTTTCTGAGAGGATTCTCACCGCTGCCAACAATACGCTTCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTTTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGTGATCATTCCAAACGGGCACTCACTGACTCAACTCTGACAGCCAATACAAc                          <  2:15359/251‑1 (MQ=255)                                                                                                                                                                                                                       ttctttcAGAGAGGATTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTAAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCATCCCGCTGGTCACAGCGATCGCATTTCACAATTTGTTGTCGGGAATCCCCCGGCAACGACTGAATGGGGATCATTCCAAACGAGCACGCACTGACGCAACTCTGACAGCCAATACAACGGGCAGAAtt                >  2:273730/1‑251 (MQ=255)                                                                                                                                                                                                                                     ggtctcACCGCTGACAACAATCCGCTGCTGGCGGACTCGCCTGAGTTCTCTTGCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATCCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGTGATCATTCCAAACGGGCACGCACGGACGCAACTCTGACAGCCAATACAACGGGCAGAATTGGTCTGCACCACCt  <  2:356126/249‑1 (MQ=255)                                                                                                                                                                                                                                     aTTCTCACCGCTGACAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGTGATCATTCCAAACGGGCACGCACTGACGCAACTCTGACAGCCAATACAACGGGCAGAATTGGTCTGCACCACCt  >  1:100814/1‑251 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                              GGCTGGCCTGGTGTGTGCCGTGTTCATCGTTCTGACTTGTAATATCGTTACCGGTATTATGCTGGGCTTTGTGACACTGGTCGTAGGCCGCGTCTTTGCACGCGAATGGCAAAAGCTGAATATTGGTACGGTGATCATTACTGCCGCACTGGTCGCATTTTACGCGGGTGGTTGGGCAATCTAATAGATTCTCCGCGCCTTTCAGCGCGGAGATTTCTTTCAGAGAGGATTCTCACCGCTGGCAACAATACGCTGCTGGCGGACTCGCCTGAGTTCTCTTTCAGTAAGCAACTGCAGCGCCTGTGTTGGGCAAGATTCAACGCAGGCTGGCCCCTCTTCCCGCTGTTCACAGAGATCGCATTTCACAATTTGTTGTCGGGTATCCCCCGGCAACGACTGAATGGTGATCATTCCAAACGGGCACGCACTGACGCAACTCTGACAGCCAATACAACGGGCAGAATTGGTCTGCACCACCT  >  NC_000913/3028305‑3028783

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