Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,539,947	C→G	A882P (GCG→CCG)	narZ ←	nitrate reductase Z subunit alpha

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,539,947	0	C	G	96.5%	88.5 / ‑6.2	29	A882P (GCG→CCG)	narZ	nitrate reductase Z subunit alpha
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base G (22/6);  minor base T (1/0);  total (23/6)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGGTGTCAATCAGCGGACGCCCCTGTGCCGGGCCATCACGCTTGGTGTAATTGAGTTTGCCGAGGAAATCGATTTCATCCTGGGTATTCCACGAAATCCCTTTACCGCCATTGCC‑AAGTTTGTCCATCAATGGCCCGAGCGAGGTAAAGCGTTCATACGTAGCAGGGTAGTCGCGCTCCACCGCCACAATATTCGGCGCGG‑TTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCGGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCTTTGGCGATACCTTTGTAGATTTCCCAGTCGCTGCGTGATTCCCACGCCGGATCGACCGCCGCAGAAAGCGGATGAATAAACGGATGCATA  >  NC_000913/1539717‑1540188                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   cagtgccagAATCACTTCTGAAGCGTCAATGGCGGTGTCAATCCGCGGGCGCCCCTGTGCCGGGCCCTAACACTTGGGGTAATTGAGTTTGCCGCGGAAATCGATTTCATCCCGGGTATTCCCCGACATCCCTTTCCCGCCCTTGCCCAAGTT‑GTCCATCAATTGCCCGCGCGAGGGAAAGCGTTCATACGTAGCAGGGTAGTCTCGCTCCACCGCCCCAATATTCGGCGGGGGTTTTGctgggatccgt                                                                                                                                                                                                                                 >  2:313198/1‑240 (MQ=255)       cagAATCACTTCTGAAGCGTCAATGGCGGTGTCAATCAGCGGACGCCCCTGTGCCGGGCCATCACGCTTGGTGTAATTGAGTTTGCCGAGGAAATCGATTTCATCCTGGGTATTCCACGAAATCCCTTTACCGCCATTGCC‑AAGTTTGTCCATCAATGGCCCGAGCGAGGTAAAGCGTTCATACGTAGCAGGGTAGTCGCGCTCCACCGCCACAATATTCGGCGGGG‑GTTTGCCTGGGATCAGATCGCAtt                                                                                                                                                                                                                         >  1:10608/1‑251 (MQ=255)             cACTTCTGAAGCGTCAATGGCGGGGTCCATCAGCGGACGCCCCTGTGCCGGGCCATCACGCATGGGGTAATTGAGTTTGCCGAGGAAATCGATTTCATACTGGGGATTTCCCCAAAACCCTTTACCGCCATTTCC‑AAGGTAGTCCATCAATGGCCCCAGCGAGGTAAAACGTTCATACGTAGCAGGGTGGTCGCGCCCCCCCGCCACAATATTCGGCGGGG‑TTTTGcccgggattcgaatggattctgctt                                                                                                                                                                                                                   >  1:63698/1‑228 (MQ=255)              aCTTCTGAAGCGTCAATGGCGGTGTCAATCAGCGGACGCCCCTGTGCCGGGCCATCACGCTTGGTGTAATTGAGTTTGCCGAGGAAATCGATTTCATCCTGGGGATTCCACGAAATCCCTTTACCGCCATTGCC‑AAGTTTGTCCATCAATGGCCCGAGCGAGGGAAAGCGTTCCTACGTAGCAGGGTAGTCGCGCTCCCCCGCCACAATATTCGGCGGGG‑TTTTGCCTGGCATCAGATCGCCATCGCCttt                                                                                                                                                                                                                  >  1:234637/1‑251 (MQ=255)                               ggCGGTGTCAATCAGCGGACGCCCCTGTGCCGGGCCATCACGCTTGGTGTAATTGAGTTTGCCGAGGGAATCGATTTCATCCTGGGTATTCCACGAAATCCCTTTACCGCCATTGCC‑AAGTTTGTCCCTCAATGGCTCGAGCGAGGTAAAGCGTTCATACGTAGCAGTGTAGTCGCGCTCCACCGCCACAATATTCGGCGTGG‑TTTTGCCTGGAATCAGATCGCATTCCCCTTTGCGCCAgacaggccctc                                                                                                                                                                                                 >  2:7312/1‑241 (MQ=255)                                            aGCGGCCGCCCCTGTGCCGGGCCATCACGCTTGGTGTAATTGCGTTTGCCGAGGAAATCGATTTCATCCTGGGTATTCCACGAAATCCCTTTACCGCCATTGCC‑AAGTTTGTCCATCAATGGCCCGAGCGAGGTAAAGCGTTCATACGTAGCAGGGTAGTCGCGCTCCACCGCCACAATATTCGGCGGGG‑TTTTGCCTGGAATCAGAGCGCCATCGCCCTTGCGCCAGTCGCGCCCTTCACACGGCtgtg                                                                                                                                                                                     >  2:37066/1‑250 (MQ=255)                                                    ccccTGTGCCGGGCCATCACGCTTGGTGTAATTGAGTTTGCCGAGGAAATCGATTTCATCCTGGGTATTCCACGAAATCCCTTTACCGCCATTGCC‑AAGTTTGTCCATCAATGGCCCGAGCGCGGTAAAGCGTTCATACGTAGCAGGGTAGTCGCGCTCCACCGCCACAATATTCTGCGGGG‑TTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGGGCCCTTCCCACGGCTGTGAGAGCTCTg                                                                                                                                                                            >  2:88315/1‑251 (MQ=255)                                                                   atcccGCTTGGTGTAATTGAGTTTGCCGAGGAAATCGATTTCATCCTGGGTATTCCACGAAATCCCTTTACCGCCATTGCC‑AAGTTTGTCCATCAATGGCCCGAGCGAGGTAAAGCGTTCATACGTAGCAGGGTAGTCGCGCTCCACCGCCACAATATTCGGCGGGG‑TTTTGCCTGGAATCAGCTCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCGGAGAGTCATga                                                                                                                                                              >  2:66923/5‑249 (MQ=255)                                                                                       gTTTGCCGAGGAAATCGATTTCATCCTGGGTATTCCACGAAATCCCTTTACCGCCATTGCC‑AAGTTTGTCCATCAATGGCCCGAGCGAGGTAAAGCGTTCATACGTAGCAGGGTAGTCGCGCTCCACCGCCACAATATTCGGCGGGG‑TTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCGGAGAGTCATGCAGTAGGGGTTGTAATACCACg                                                                                                                                         >  1:267488/1‑251 (MQ=255)                                                                                                       gATTTCATCCTGGGTATTCCACGAAATCCCTTTACCGCCATTGCC‑AAGTTTGTCCATCAATGGCCCGAGCGAGGTAAAGCGTTCATACGTAGCAGGGTAGTCGCGCTCCACCGCCACAATATTCGGCGGGG‑TTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCGGAGAGGCATGCAGTAGGGGTTGTAATAACACGTCGGTTTCTTTGCCa                                                                                                                          >  1:57812/1‑250 (MQ=255)                                                                                                             aTCCTGGGTATTCCACGAAATCCCTTTACCGCCATTGCC‑AAGTTTGTCCATCAATGGCCCGAGCGAGGTAAAGCGTTCATACGTAGCAGGGTAGTCGCGCTCCACCGCCACAATATTCGGCGGGG‑TTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCGGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGAc                                                                                                                   >  1:174165/1‑251 (MQ=255)                                                                                                                      aTTCCACGAAATCCCTTTACCGCCATTGCC‑AAGTTTGTCCATCAATGGCCCGAGCGAGGTAAAGCGTTCATACGTAGCAGGGTAGTCGCGCTCCACCGCCACAATATTCGGCGGGG‑TTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCGGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGcac                                                                                                           >  1:298639/1‑250 (MQ=255)                                                                                                                             gAAATCCCTTTACCGCCATTGCC‑AAGTTTGTCCATCAATGGCCCGAGCGAGGTAAAGCGTTCATACGTAGCAGGGTAGTCGCGCTCCACCGCCACAATATTCGGCGGGG‑tttt                                                                                                                                                                                                                                             >  1:162937/1‑113 (MQ=255)                                                                                                                               actcccttccccGCCATTGCA‑AAGTTTTCCCATAAAGGGCCGGAGAGGGGAAAAGCTTTCATCCGTGGCAGGGTTGTCGCGCCCCACCGCCACAATATTCGGCGGGG‑TTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCGGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGaaa                                                                                                 <  1:116523/241‑1 (MQ=255)                                                                                                                                           gCCATTGCC‑AAGTTTGTCCATCAATGGCCCGAGCGAGGTAAAGCGTTCATACGTAGCAGGGTAGTCGCGCTCCACCGCCACAATATTCGGCGGGG‑TTTTGCCTGGAATCAGATCGCATTCGACTTTGCGCCAGTCGAGCACTTCACACGGGTGTGAGAGCTCTGCCGGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCTTTGGCgg                                                                                     >  1:123185/1‑250 (MQ=255)                                                                                                                                                   c‑AAGTTTGTCCATCAATGGCCCGAGCGAGGTAAAGCGTTCATACGTAGCAGGGTAGTCGCGCGCCACCGCCACAATATTCGGCGGGG‑TTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGCGCTCTGCCGGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGCACATGCTTTGGCGactcctttt                                                                             >  2:220801/1‑243 (MQ=255)                                                                                                                                                    aaattttgccaaaaaagggCCCGAGGGAGGAAAAGCGTTCATACGTAGCAGGGTAGTCGCCCCCCTCCGCCACAATATTCGGCGGGG‑TTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCGGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCTTTGGCGATACCTTTg                                                                             <  1:86934/234‑1 (MQ=255)                                                                                                                                                                  aaTGGCCCGAGCGAGGTAAAGCGTTCATACGTAGCAGGGTAGTCGCGCTCCACCGCCACAATATTCGGCGGGG‑TTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCGGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGCAAATGCTTTGGCGATACCTTTGTAGATTTCCCAGTc                                                               >  1:58908/1‑251 (MQ=255)                                                                                                                                                                        ccGAGCGAGGTAAAGCGTTCATACGTAGCAGGGTAGTCGCGCTCCACCGCCACAATATTCGGCGGGG‑TTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCGGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCTTTGGCGATACCTTTGTAGATTTCCCAGTCGCTGCg                                                         >  1:266303/1‑251 (MQ=255)                                                                                                                                                                        ccGAGCGAGGTAAAGCGTTCATACGTAGCAGGGTAGTCGCGCTCCACCGCCACAATATTCGGCGGGG‑TTTTGCATGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCGGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGCTGACCCACGCACACTTGCGAAAATGCTTTGGCGATACCTTTGTAGATTTCCCAGTCGCTGCg                                                         >  1:370473/1‑251 (MQ=255)                                                                                                                                                                          gagagaggaaaagtgtcatctgtgaGCAGGGTAGTCTCTCTCCTCCGCCACAATATTCGGCGGGG‑TTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCGGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCTTTGGCGATACCTTTGTAGATTTCCCAGTCGCTGCGTg                                                       <  1:81435/227‑1 (MQ=255)                                                                                                                                                                              gagGTAAAGCGTTCATACGTAGCAGGGTAGTCGCGCTCCACCGCCACAATATTCGGCGGGG‑TTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCGGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCCCACTTGCGAAAATGCTTTGGCGATACCTTTGTAGATTTCCCAGTCGCTGCGTGGTTc                                                   >  2:154143/1‑251 (MQ=255)                                                                                                                                                                              gagGTAAAGCGTTCATACGTAGCAGGGTAGTCGCGCTCCACCGCCACAATATTCGGCGGGG‑TTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCGGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCTTTGGCGATACCTTTGTAGATTTCCCAGTCGCTGCGTGATTc                                                   >  1:261752/1‑251 (MQ=255)                                                                                                                                                                                     aGCGTTCATACGTAGCAGGGTAGTCGCGCTCCACCGCCACAATATTCGGCGGGG‑TTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCGGAGAGTCATGAAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCTTTGGCGATACCTTTGTAGATTTCCCAGTCGCCGCGTGATTCCCACGcc                                            >  2:108191/1‑251 (MQ=255)                                                                                                                                                                                                 aaggaggttgtccgccccccgccGCCACAATATTCGGCGGGG‑TTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCGGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCTTTGGCGATACCTTTGTAGATTTCCCAGTCGCTGCGTGATTCCCACGCCGGATCGAccgcc                                <  1:37066/230‑1 (MQ=255)                                                                                                                                                                                                     aGGGTAGTCGCGCTCCACCGCCACAATATTCGGCGGGG‑TTTTGACTGGAATCAGATAGCCTACGCCTTTGCGACAGTCGAGCACTTCACACGGCTGAGAGAACTCTGACGGAGAGTCATGCAGTAAGGGGTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCTTTGGCGATACCTTTGTAGATTTCCCAGTCGCTGCGGGATTCCCACACCGGATCGACcgccgcgg                            >  1:366353/1‑249 (MQ=255)                                                                                                                                                                                                            cccccccaccccGCCACAATATCCGGGGGGG‑GTTTGCCTGGAATCAGATCGCATTCTCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGATAGCTCTTCCGGAGCGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGCCCCACGCCCCCTTGCGAAAATGCTTTGGCGATACCTTTGTAGATTTCCCAGTCGCTGCGTGATTCCCACGCCGGATCGACCGCCGCAGAAAGCgg                     <  2:131925/241‑1 (MQ=255)                                                                                                                                                                                                                 cTCCACCGCCACAATATTCGGCGGGG‑TTTTGCCTGGAATCAGATCGCCTTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCGGAGAGTCATGCAGAAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCTTTGGCGATACCTTTTTAGATTTCCCAGTCGCTGCTTGCTTCCCACGCCGGATCGACCGCCGCCGAAAGCGGGTGa                 >  2:133010/1‑250 (MQ=255)                                                                                                                                                                                                                                atTCGGCGGGG‑TTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCGGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCTTGGGCGATACCTTTGTAGATTTCCCAGTCGCTGCGTGATTCCCACGCCGGATCGACCGCCGCAGAAAGCGGATGAATAAACGGATGCata  <  1:218878/250‑1 (MQ=255)                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   CAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGGTGTCAATCAGCGGACGCCCCTGTGCCGGGCCATCACGCTTGGTGTAATTGAGTTTGCCGAGGAAATCGATTTCATCCTGGGTATTCCACGAAATCCCTTTACCGCCATTGCC‑AAGTTTGTCCATCAATGGCCCGAGCGAGGTAAAGCGTTCATACGTAGCAGGGTAGTCGCGCTCCACCGCCACAATATTCGGCGCGG‑TTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCGGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCTTTGGCGATACCTTTGTAGATTTCCCAGTCGCTGCGTGATTCCCACGCCGGATCGACCGCCGCAGAAAGCGGATGAATAAACGGATGCATA  >  NC_000913/1539717‑1540188

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