Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,234,632	A→G	20.0%	L29S (TTA→TCA)	nhaB ←	Na(+):H(+) antiporter NhaB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,234,632	0	A	G	20.0%	41.6 / 5.3	20	L29S (TTA→TCA)	nhaB	Na(+):H(+) antiporter NhaB
Reads supporting (aligned to +/- strand):  ref base A (10/6);  new base G (2/2);  total (12/8)
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.75e-01

CATAAAATAGATACCCGCCACCATAAACATCAGTAACAGCAAGACTTCAAGATTTGCCGCCACCTCTTCACGGACGTGTTCCGCGCTGGTCATGCCGATGAATACCGCTTCGATAGCCAACAGACCACCGGGGAGCAGCGGGTAGCATTTCAGGGCCATCGCCAGAGTGAAAATAAATTCCGCTACCAGCAACCAGCCCGCGACGAAAGGGCTGATGAGGAAAATTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGACTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCCCTTACCTTACAAATAATCAATGATGTTTTTATGTTTAAACGCAAAGCTTAACGGTCAGGCAGGAGTGAGGCAAGTCTTGATAGTCAAGGGGAAAGAGATGCGGAAAATGAAGCCTTGATCCCTTTTTCTTCTTTTTGTCTGCTATCAGCGTAGTTAGCCCTC  >  NC_000913/1234391‑1234889                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                  cATAAAATAGATACCCGCCACCATAAACATCAGTAACAGCAAGACTTCAAGATTTGCCGCCACCTCTTCACGGACGTGTTCCGCGCTGGTCATGCCGATGAATACCGCTTCGATAGCCAACAGACCACCGGGGAGCAGCGGGTAGCATTTCAGGGCCATCGCCAGAGTGAAAATAAATTCCGCTACCAGCAACCAGCCCGCGACGAAAGGGCTGATGAGGAAAATTAACGGGTTTACGATTAAGAAAatgatg                                                                                                                                                                                                                                                        >  1:254644/1‑253 (MQ=255) cATAAAATAGATACCCGCCACCATAAACATCAGTAACAGCAAGACTTCAAGATTTGCCGCCACCTCTTCACGGACGTGTTCCGCGCTGGTCATGCCGATGAATACCGCTTCGATAGCCAACAGACCACCGGGGAGCAGCGGGTAGCATTTCAGGGCCATCGCCAGAGTGAAAATAAATTCCGCTACCAGCAACCAGCCCGCGACGAAAGGGCTGATGAGGAAAATTAACGGGTTTACGATTAAGAAAAtgatga                                                                                                                                                                                                                                                       >  1:239500/1‑254 (MQ=255)                                 gTAACAGCAAGACTTCAAGATTTGCCGCCACCTCTTCACGGACGTGTTCCGCGCTGGTCATGCCGATGAATACCGCTTCGATAGCCAACAGACCACCGGGGAGCAGCGGGTAGCATTTCAGGGCCATCGCCAGAGTGAAAATAAATTCCGCTACCAGCAACCAGCCCGCGACGAAAGGGCTGATGAGGAAAATTAACGGGTTTACGATTAAGAAAATAAGAAGGGCGGGTTGTTACAAGTCGGGGAATGGGCCAAAAAAATTCCGCCATAgcgc                                                                                                                                                                                                   >  1:26694/1‑274 (MQ=255)                                  tAACAGCAAGACTTCAAGATTTGCCGCCACCTCTTCACGGACGTGTTCCGCGCTGGTCATGCCGATGAATACCGCTTCGATAGCCAACAGACCACCGGGGAGCAGCGGGTAGCATTTCAGGGCCATCGCCAGAGTGAAAATAAATTCCGCTACCAGCAACCAGCCCGCGACGAAAGGGCTGATGAGGAAAATTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACAAGTCGGGGGACTGGCCCAAAAAATTGCGCCATAg                                                                                                                                                                                                      >  1:144710/1‑270 (MQ=255)                                                                                                           gCTTCGATAGCCAACAGACCACCGGGGAGCAGCGGGTAGCATTTCAGGGCCATCGCCAGAGTGAAAATAAATTCCGCTACCAGCAACCAGCCCGCGACGAAAGGGCTGATGAGGAAAATTAACGGGGTTACGATTGAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGACTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCCCTTACCTTACAAATAATCAATGAATTTTTTATGTTTAAAAGCa                                                                                                                         >  1:36671/1‑274 (MQ=255)                                                                                                             ttCGATAGCCAACAGACCACCGGGGAGCAGCGGGTAGCATTTCAGGGCCATCGCCAGAGTGAAAATAAATTCCGCTACCAGCAACCAGCCCGCGACGAAAGGGCTGATGAGGAAAATTAACGGGTTTACGATTAAGAAAAtgatga                                                                                                                                                                                                                                                       >  2:105444/1‑146 (MQ=255)                                                                                                                gATAGCCAACAGACCACCGGGGAGCAGCGGGTAGCATTTCAGGGCCATCGCCAGAGTGAAAATAAATTCCGCTACCAGCAACCAGCCCGCGACGAAAGGGCTGATGAGGAAAATTAACGGGTTTACGATTAAGAAAATGATGAggg                                                                                                                                                                                                                                                    <  1:105444/146‑1 (MQ=255)                                                                                                                    caaaaggccccacgGGGGA‑CAGCGGTTGGACTTTCAGGGCCGACCCCAGAGGGAAAATAATATCCGCTACCAGCAACCAGCCCGCGACGAAATGGCTGATGAGGAAAATTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGACTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCCCTTACCTTACAAATAATCAATGATGTTTTTATGTTTAAACGCAAAGCTTAACg                                                                                                               <  2:143889/265‑1 (MQ=255)                                                                                                                             ccaccGGGGAGCAGCGGGTAGCATTTCAGGGCCATCGCCAGAGTGAAAATAAATTCCGCTACCAGCAACCAGCCCGCGACGAAAGGGCTGATGAGGAAAATTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGACTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCCCTTACCTTACAAATAATCAATGATGTTTTTATGTTTAAACGCAAAGCTTa                                                                                                                  >  1:143889/1‑263 (MQ=255)                                                                                                                                               tAGCATTTCAGGGCCACGGCCAGGGTGAAAATAAATTCCGCTACCAGCAACCAGCCCGGGACGAAAGGGCGGATGAGGAAAATTAACGGGTTTACGATTGAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGACTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCCCTTACCTTACAAATAATCAATGATGTTTTTATGTTTAAACGCAAAGCTTAACGGTCAGGCAGGAGTGAGGCAAGTCTTg                                                                                     <  2:36671/274‑1 (MQ=255)                                                                                                                                                           gCCATCGCCAGAGTGAAAATAAATTCCGCTACCAGCAACCAGCCCGCGACGAAAGGGCTGATGAGGAAAATTAACGGGTTTACGATTAAGAAAatgatg                                                                                                                                                                                                                                                        >  1:178018/1‑99 (MQ=255)                                                                                                                                                              aTCGCCAGAGTGAAAATAAATTCCGCTACCAGCAACCAGCCCGCGACGAAAGGGCTGATGAGGAAAATTAACGGGTTTACGATTAAGAAAATGATGAgg                                                                                                                                                                                                                                                     <  2:178018/99‑1 (MQ=255)                                                                                                                                                                cGCCAGAGTGAAAATAAATTCCGCTACCAGCAACCAGCCCGCGACGAAAGGGCTGATGAGGAAAATTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGACTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCCCTTACCTTACAAATAATCAATGATGTTTTTATGTTTAAACGCAAAGCTTaa                                                                                                                 >  2:263638/1‑229 (MQ=255)                                                                                                                                                                cGCCAGAGTGAAAATAAATTCCGCTACCAGCAACCAGCCCGCGACGAAAGGGCTGATGAGGAAAATTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGACTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCCCTTACCTTACAAATAATCAATGATGTTGTTATGTTTAAACGCAAAGCTTAACGGTCAGGCAGGAGTGAGGCAAGGCTTGATAGTCAAGGGGaagaa                                                                    >  1:15271/1‑271 (MQ=255)                                                                                                                                                                   cAGAGTGAAAATAAATTCCGCTACCAGCAACCAGCCCGCGACGAAAGGGCTGATGAGGAAAATTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGACTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCCCTTACCTTACAAATAATCAATGATGTTTTTATGTTTAAACGCAAAGCTTAACgg                                                                                                              <  1:263638/229‑1 (MQ=255)                                                                                                                                                                          aaaaTAAATTCCGATACCATCAACCCGCCCGGGACGAAAGGGCTGATGAGGAAAATTAACGGGTTTACGATTAAGAAAAGGATGAGGGCGAGTTTGTACCAGGCGGGGGACTGGCCAAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCCCTTACCTTACAaataa                                                                                                                                                    <  1:294365/184‑1 (MQ=255)                                                                                                                                                                                                     cccGCGACGAAAGGGCTGATGAGGAAAATTAACGGGTTTACGATTGAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGACTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCCCTTACCTTACAAATAATCAATGATGTTTTTATGTTTAAACGCAAAGCTTAACGGTcaggcag                                                                                                      >  1:290996/1‑203 (MQ=255)                                                                                                                                                                                                        gcgACGAAAGGGCTGATGAGGAAAATTAACGGGTTTACGATTGAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGACTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCCCTTACCTTACAAATAATCAATGATGTTTTTATGTTTAAACGCAAAGCTTAACGGTCAGGCAGGAg                                                                                                   <  2:290996/203‑1 (MQ=255)                                                                                                                                                                                                                                  tAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGACTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCCCTTACCTTACAAATAATCAATGATGTTTTTATGTTTAAACGCAAAGCTTAACGGTCAGGCAGGAGTGAGGCAAGTCTTGATAGTCAAGGGGAAAGAGATGCGGAAAATGAAACCTTGATCCCTTTTTCTTCTTTTTGTCTGTTATCCGCGTAGTTATccccc  >  2:283021/1‑272 (MQ=255)                                                                                                                                                                                                                                     cGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGACTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCCCTTACCTTACAAATAATCAATGATGTTTTTATGTTTAAACGCAAAGCTTAACGGTCAGGCAGGAGTGAGGCAAGTCTTGATAGTCAAGGGGAAAGAGATGCGGAAAATGAAGCCTTGATCCCTTTTTCTTCTTTTTGTCTGCTATCAGCGTAGTTAGccc    <  1:283021/269‑1 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                  CATAAAATAGATACCCGCCACCATAAACATCAGTAACAGCAAGACTTCAAGATTTGCCGCCACCTCTTCACGGACGTGTTCCGCGCTGGTCATGCCGATGAATACCGCTTCGATAGCCAACAGACCACCGGGGAGCAGCGGGTAGCATTTCAGGGCCATCGCCAGAGTGAAAATAAATTCCGCTACCAGCAACCAGCCCGCGACGAAAGGGCTGATGAGGAAAATTAACGGGTTTACGATTAAGAAAATGATGAGGGCGAGTTTGTACCAGTCGGGGGACTGGCCCAAAAAATTGCGCCATAGCGCGCGGCCCCAGGAGATCTCCATGATGGTTTCCCTTACCTTACAAATAATCAATGATGTTTTTATGTTTAAACGCAAAGCTTAACGGTCAGGCAGGAGTGAGGCAAGTCTTGATAGTCAAGGGGAAAGAGATGCGGAAAATGAAGCCTTGATCCCTTTTTCTTCTTTTTGTCTGCTATCAGCGTAGTTAGCCCTC  >  NC_000913/1234391‑1234889

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