Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,184,432	T→A	15.8%	I5F (ATT→TTT)	potB ←	spermidine preferential ABC transporter membrane subunit PotB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,184,432	0	T	A	15.8%	60.4 / 3.3	24	I5F (ATT→TTT)	potB	spermidine preferential ABC transporter membrane subunit PotB
Reads supporting (aligned to +/- strand):  ref base T (8/12);  new base A (2/2);  total (10/14)
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.83e-01

GCAACTTCGCCAGAAACCAGGCAAACGGGTAGCCCAGCACCAGGCAGGCGAGGGTGGCGATCAGCGCCATATTCAGCGAGTGCAATAGCACTTCAAAATAGAGCGGATCGAGCAGACGCGTGTAGTTATCCAGCGTAAAGACCATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAGGTTGGGCAGAAAGACAAACAACACAAGCCAACCGACAATAGTGACAATCACTACATTCTGGAACTTACTTGTGTTCTTCATCAGCCAGTACGACCTCCCAGCTTTCTACCCAATTAATGGCCATTTTTTGGTCGAGAGAGTGGTCAAAGTCAGGATCGTCTTCATTGAAGAATTCGCTGACCATCACCATCTTGCCATTTTCCAGTTCAACAACCGACTCCAGCGTCATGCCTTTGTAGTTACGCTCGCGAACGTAACCAATCAGCCCTTCAGCGTGGTTGTCGTCGTTAATCTCTTCAACACGTA  >  NC_000913/1184195‑1184699                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                            gCAACTTCGCCAGAAACCAGGCAAACGGGTAGCCCAGCACCAGGCAGGCGAGGGTGGCGATCAGCGCCATATTCAGCGAGTGCAATAGCACTTCAAAATAGAGCGGATCGAGCAGACGCGTGTAGTTATCCAGCGTAAAGACCATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAGGTTGGGCAGAAAGACAAACAACACAAGCCAACCGACAATAGTGACAATCa                                                                                                                                                                                                                                                                  <  1:166335/249‑1 (MQ=255)                                               ggcGAGGGTGGCGATCAGCGCCATATTCAGCGAGTGCAATAGCACTTCAAAATAGAGCGGATCGAGCAGACGCGTGTAGTTATCCAGCGTAAAGACCATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAGGTTGGGCAGAAAGACAAACAACACAAGCCAACCGACAATAGTGACAATCACTACATTCTg                                                                                                                                                                                                                                                        <  1:39108/213‑1 (MQ=255)                                               ggcGAGGGTGGCGATCAGCGCCATATTCAGCGAGTGCAATAGCACTTCAAAATAGAGCGGATCGAGCAGACGCGTGTAGTTATCCAGCGTAAAGACCATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAGGTTGGGCAGAAAGACAAACAACACAAGCCAACCGACAATAGTGACAATCACTACATTCTGGAacttactt                                                                                                                                                                                                                                              <  1:50853/223‑1 (MQ=255)                                                 cGAGGGTGGCGATCAGCGCCATATTCAGCGAGTGCAATAGCACTTCAAAATAGAGCGGATCGAGCAGACGCGTGTAGTTATCCAGCGTAAAGACCATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAGGTTGGGCAGAAAGACAAACAACACAAGCCAACCGACAATAGTGACAATCACTACATTCTGGAACTTCCTTTGGTTCTTCATCAGCCAGTACGACCCCACAGCTTTCCACCCAATTAATTGc                                                                                                                                                                                           >  2:150959/1‑272 (MQ=255)                                                    gggTGGCGATCAGCGCCATATTCAGCGAGTGCAATAGCACTTCAAAATAGAGCGGATCGAGCAGACGCGTGTAGTTATCCAGCGTAAAGACCATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAGGTTGGGCAGAAAGACAAACAACACAAGCCAACCGACAATAGTGACAATCACTACATTCTGGAACTTACTTGTGttc                                                                                                                                                                                                                                        >  2:73758/1‑224 (MQ=255)                                                     ggTGGCGATCAGCGCCATATTCAGCGAGTGCAATAGCACTTCAAAATAGAGCGGATCGAGCAGACGCGTGTAGTTATCCAGCGTAAAGACCATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAGGTTGGGCAGAAAGACAAACAACACAAGCCAACCGACAATAGTGACAATCACTACATTCTGGAACTTACTTGTGttct                                                                                                                                                                                                                                       <  1:73758/224‑1 (MQ=255)                                                             tCAGCGCCATATTCAGCGAGTGCAATAGCACTTCAAAATAGAGCGGATCGAGCAGACGCGTGTAGTTATCCAGCGTAAAGACCATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAGGTTGGGTAGAAAGACAAACAACACAAGCCAACCGACAATAGTGa                                                                                                                                                                                                                                                                        >  2:63324/1‑183 (MQ=255)                                                              cAGCGCCATATTCAGCGAGTGCAATAGCACTTCAAAATAGAGCGGATCGAGCAGACGCGTGTAGTTATCCAGCGTAAAGACCATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAGGTTGGGTAGAAAGACAAACAACACAAGCCAACCGACAATAGTGAc                                                                                                                                                                                                                                                                       <  1:63324/183‑1 (MQ=255)                                                                                gTGCAATAGCACTTCAAAATAGAGCGGATCGAGCAGACGCGTGTAGTTATCCAGCGTAAAGACCATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAGGTTGGGCAGAAAGACAAACAACACAAGCCAACCGACAATAGTGACAATCACTACAt                                                                                                                                                                                                                                                            >  1:109430/1‑176 (MQ=255)                                                                                 tGCAATAGCACTTCAAAATAGAGCGGATCGAGCAGACGCGTGTAGTTATCCAGCGTAAAGACCATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAGGTTGGGCAGAAAGACAAACAACACAAGCCAACCGACAATAGTGACAATCACTACAtt                                                                                                                                                                                                                                                           <  2:109430/176‑1 (MQ=255)                                                                                              cAAAATAGAGCGGATCGAGCAGACGCGTGTAGTTATCCAGCGTAAAGACCATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAGGTTGGGCAGAAAGACAAACAACACAAGCCAACCGACAATAg                                                                                                                                                                                                                                                                           >  2:153448/1‑147 (MQ=255)                                                                                               aaaaTAGAGCGGATCGAGCAGACGCGTGTAGTTATCCAGCGTAAAGACCATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAGGTTGGGCAGAAAGACAAACAACACAAGCCAACCGACAATAGt                                                                                                                                                                                                                                                                          <  1:153448/147‑1 (MQ=255)                                                                                                 aaTAGAGCGGATCGAGCAGACGCGTGTAGTTATCCAGCGTAAAGACCATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAGGTTGGGCAGAAAGACAAACAACACAAGCCAACCGACAAAAGTGACAATCACTACATTCTGGAACTTACTTGTGTTCTTCATCAGCCAGTACGACCTCCCAGCTTTCTACCCAATTAATGGCCATTTTTTGGTCgaga                                                                                                                                                                           >  2:95600/1‑240 (MQ=255)                                                                                                 aaTAGAGCGGATCGAGCAGACGCGTGTAGTTATCCAGCGTAAAGACCATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAGGTTGGGCAGAAAGACAAACAACACAAGCCAACCGACAAAAGTGACAATCACTACATTCTGGAACTTACTTGTGTTCTTCATCAGCCAGTACGACCTCCCAGCTTTCTACCCAATTAATGGCCATTTTTTGGTCgaga                                                                                                                                                                           <  1:95600/240‑1 (MQ=255)                                                                                                                        atgtAGTTATCCAGCGTAAAGACCATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAGGTTGGGCAGAAAGACAAACAACACAAGCCAACCGACAATAGTGACAATCACTACATTCTGGAACTTACTTGTGTTCTTCATCAGCCAGTACGACCTCCCAGc                                                                                                                                                                                                              <  2:65314/181‑1 (MQ=255)                                                                                                                            aGTTATCCAGCGTAAAGACCATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAGGTTGGGCAGAAAGACAAACAACACAAGCCAACCGACAATAGTGACAATCACTACATTCTGGAACTTACTTGTGTTCTTCATCAGCCAGTACGACCTccc                                                                                                                                                                                                                 >  1:9264/1‑175 (MQ=255)                                                                                                                             gTTATCCAGCGTAAAGACCATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAGGTTGGGCAGAAAGACAAACAACACAAGCCAACCGACAATAGTGACAATCACTACATTCTGGAACTTACTTGTGTTCTTCATCAGCCAGTACGACCTCCCa                                                                                                                                                                                                                <  2:9264/175‑1 (MQ=255)                                                                                                                                         aaaGACCATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAGGTTGGGCAGAAAGACAAACAACACAAGCCAACCGACAAAAGTGACAATCACTACATTCTGGAACTTACTTGTGTTCTTCATCAGCCAGTACGACCTCCCAGCTTTCTACCCAATTAATGGCCATTTTTTGGTCGAGAGAGTGGTCAAAGTCAGGATCGTCCTCATTGAATAATTTGCTAACCCACCCTAACTTGACATTTTcctcttta                                                                                                   >  2:153079/1‑266 (MQ=255)                                                                                                                                                         aaaCTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAGGTTGGGCAGAAAGACAAACAACACAAGCCAACCGACAATAGTGACAATCACTACATTCTGGAACTTACTTGTGTTCTTCATCAGCCAGTACGACCTCCCAGCTTTCTACCCAATTAATGGCCATTTTTTGGTCGAGAGAGTGGTCAAAGTCAGGATCGTCTTCATTGAAGAATTCGCTGACCATCACCATCTTGCCATTTTCCAGTTCAACAACCGACTCCAGCg                                                                                   <  1:150959/272‑1 (MQ=255)                                                                                                                                                          aaCGGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAGGTTGGGCAGAAAGACAAACAACACAAGCCAACCGACAAAAGTGACAATCACTACATTCTGGAACTTACTTGTGTTCTTCATCAGCCAGTACGACCTCCCAGCTTTCTACCCAATTAATGGCCATTTTTTGGTCGAGAGAGTGGTCAAAGTCAGGATCGTCTTCATTGAAGAATTCGCTGACCATCACCATCTTGCCATTTTCCAGTTCAACAACCGACTCCAGCGt                                                                                  <  1:153079/272‑1 (MQ=255)                                                                                                                                                                                                                 gaaGACAAACAACACAAGCCAACCGACAATAGTGACAATCACTACATTCTGGAACTTACTTGTGTTCTTCATCAGCCAGTACGACCTCCCAGCTTTCTACCCAATTAATGGCCATTTTTTGGTCGAGAGAGTGGTCAAAGTCAGGATCGTCTTCATTGAAGAATTCGCTGACCATCACCATCTTGCCATTTTCCAGTTCAACAACCGACTCCAGCGTCATGCCTTTGTAGTTACGCTCGCGAACGTAACCAATCAGCCCTTCAGCGTGGTTg                           <  1:127296/271‑1 (MQ=255)                                                                                                                                                                                                                    gACAAACAACACAAGCCAACCGACAATAGTGACAATCACTACATTCTGGAACTTACTTGTGTTCTTCATCAGCCAGTACGAc                                                                                                                                                                                                                      >  2:205948/1‑82 (MQ=255)                                                                                                                                                                                                                     acaaacaaCACAAGCCAACCGACAATAGTGACAATCACTACATTCTGGAACTTACTTGTGTTCTTCATCAGCCAGTACGAcc                                                                                                                                                                                                                     <  1:205948/82‑1 (MQ=255)                                                                                                                                                                                                                                          aCAATAGTGACAATCACTACATTCTGGAACTTACTTGTGTTCTTCATCAGCCAGTACGACCTCCCAGCTTTCTACCCAATTAATGGCCATTTTTTGGTCGAGAGAGTGGTCAAAGTCAGGATCGTCTTCATTGAAGAATTCGCTGACCATCACCATCTTGCCATTTTCCAGTTCAACAACCGACTCCAGCGTCATGCCTTTGTAGTTACGCTCGCGAACGTAACCAATCAGCCCTTCAGCGTGGTTGTCGTCGTTAATCTCTTCAACACGTa  >  1:125319/1‑272 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                            GCAACTTCGCCAGAAACCAGGCAAACGGGTAGCCCAGCACCAGGCAGGCGAGGGTGGCGATCAGCGCCATATTCAGCGAGTGCAATAGCACTTCAAAATAGAGCGGATCGAGCAGACGCGTGTAGTTATCCAGCGTAAAGACCATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAGGTTGGGCAGAAAGACAAACAACACAAGCCAACCGACAATAGTGACAATCACTACATTCTGGAACTTACTTGTGTTCTTCATCAGCCAGTACGACCTCCCAGCTTTCTACCCAATTAATGGCCATTTTTTGGTCGAGAGAGTGGTCAAAGTCAGGATCGTCTTCATTGAAGAATTCGCTGACCATCACCATCTTGCCATTTTCCAGTTCAACAACCGACTCCAGCGTCATGCCTTTGTAGTTACGCTCGCGAACGTAACCAATCAGCCCTTCAGCGTGGTTGTCGTCGTTAATCTCTTCAACACGTA  >  NC_000913/1184195‑1184699

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