Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,086,634	G→A	A143A (GCC→GCT)	plaP ←	putrescine:H(+) symporter PlaP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,086,634	0	G	A	95.8%	68.2 / ‑5.7	24	A143A (GCC→GCT)	plaP	putrescine:H(+) symporter PlaP
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (10/13);  total (11/13)
Fisher's exact test for biased strand distribution p-value = 4.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

TCCGACAGGTTGCTGATGCCGTCAAAGCCGGTAAAGGAGAAGCACAGGATTGTCGCCCCGGTAATCATCGGGATAACATGCGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAG‑CGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAGGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGCTAATGGATTTCTGGGCGTAAGTGTATGCAGAGCCAGCAGAAGGATAGCGGCGAACCAGCTTCCCGTAGCTCAGAGCCGTAAAC  >  NC_000913/2086389‑2086874                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                 tcCGACAGGTTGCTGATGCCGTCAAAGCCGGTAAAGGAGAAGCACAGGATTGTCGCCCCGGTAATCATCGGGATAACATGCGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAG‑CGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCg                                                                                                                                                                                                                                               <  2:93533/249‑1 (MQ=255) tcCGACAGGTTGCTGATGCCGTCAAAGCCGGTAAAGGAGAAGCACAGGATTGTCGCCCCGGTAATAATCGGGATAACATGCGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAG‑CGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGt                                                                                                                                                                                                                                             >  1:202883/1‑251 (MQ=255)                  gCCGTCAAAGCCGGTAAAGGAGAAGCACAGGATTGTCGCACCGGTAATCATCGGGATAACATGCGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAG‑CGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAgc                                                                                                                                                                                                                             <  2:271466/250‑1 (MQ=255)                         aaGCCGGTAAAGGAGAAGCACAGGATTGTCGCCCCGGTAATCATCGGGATAACATGCGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAG‑CGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCcaccac                                                                                                                                                                                                                     >  2:269922/1‑251 (MQ=255)                             cGGTAAAGGAGAAGCACAGGATTGTCGCCCCGGTAATCATCGGGAAAACATGCGCATCGCCAGACAGGAAGGGACGAGGGCTCGCCAGCGT‑CCAGGCCCCTTCACCTTCAAATACGCCATAAACACCCATCCCAAGATTAACCGCGATCAGCCCTACCTGCAACACGCGGATACCAGGGTTGAGGTTCGCTCCGGTTTTCAAACTCCGCGGGTTAAAGGCGGTCATAAGGGCACCCAGC‑CCACCACCAAc                                                                                                                                                                                                                 >  1:265502/1‑250 (MQ=255)                                      agaAGCACAGGATTGTCTCCCCGGTAATCATCGGGATAACATGCGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAG‑CGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACgat                                                                                                                                                                                                        >  2:164118/1‑251 (MQ=255)                                           aacaGGATTGGCGCACCGGTCATCATCGGGATAACATGCGCTTCGCCAGACCAGAATGGACGAGTGCTAGCCAGCGTACCAG‑CGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATAACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCCCCAGCGCCACCTCAATCATCCACgatgggat                                                                                                                                                                                                   <  2:228881/250‑1 (MQ=255)                                              aGGATTGTCGCCCCGGTAATCATCGGGATAACATGCGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAG‑CGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACgatgggatgg                                                                                                                                                                                                 >  2:199085/1‑250 (MQ=255)                                                                   atcGGGATAACATGCGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAG‑CGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATa                                                                                                                                                                           <  2:196771/251‑1 (MQ=255)                                                                                    aTCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAG‑CGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGaa                                                                                                                                                          <  1:81890/251‑1 (MQ=255)                                                                                           gaccagaATGGACGAGTGCTCGCCAGCGTACCAG‑CGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGGCACCAGCGCCACCACGAACATCCACGAGGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATTTTGat                                                                                                                                                   <  2:137519/251‑1 (MQ=255)                                                                                                    ggACGAGTGCTCGCCAGCGTACCAG‑CGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGgcgc                                                                                                                                          <  1:269922/251‑1 (MQ=255)                                                                                                                 ccAGCGTACCAG‑CGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCg                                                                                                                             >  1:145147/1‑251 (MQ=255)                                                                                                                        cccAG‑CGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCagag                                                                                                                      <  2:145147/250‑1 (MQ=255)                                                                                                                                    caccttcaAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTATTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAAccc                                                                                                           >  2:153276/1‑251 (MQ=255)                                                                                                                                         tcaAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACccacca                                                                                                       <  1:40375/250‑1 (MQ=255)                                                                                                                                         tcaAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACccacca                                                                                                       <  2:78712/250‑1 (MQ=255)                                                                                                                                                  ccATAAACAACCATGCCCAGAATCACAGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTAAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACACACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTCGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGCAGACCAACCCACCATAAAGCCGAc                                                                                             >  2:283097/1‑251 (MQ=255)                                                                                                                                                        acaacCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCgg                                                                                       <  2:265502/251‑1 (MQ=255)                                                                                                                                                                                    gCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGCTAATGGATTTCTGGGCGTAAGTGTAt                                                           >  2:243710/1‑251 (MQ=255)                                                                                                                                                                                         aCCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGCCCAACCCACCATAAAGCCGACAGTCGGGCTAATGGATTTCTGGTCGTAAGTGTATGCaga                                                      <  1:71543/251‑1 (MQ=255)                                                                                                                                                                                                                     aGTTCGCTACGGATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGGACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGCGAATGGATTTCTGGGCTGAAGGGTATGCAGAGCCAGCAGAAGGATAGCGGCGACCCAGc                          <  1:198175/251‑1 (MQ=255)                                                                                                                                                                                                                                gATTTCAGACTACGCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGCTAATGGATTTCTGGGCGTAAGTGTATGCAGAGCCAGCAGAAGGATAGCGGCGAACCAGCTTCCCGTAGCt               >  2:193590/1‑251 (MQ=255)                                                                                                                                                                                                                                             gCAGGTTAAAAGCGGTCATAAAGGCCACCAGCGCCACCACAAAAATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGCTAATGGATTTATGGGCGTACGAGAATGCAGAGCCAGCAGAAGGATAGCGGCGAACCAGCTTCCCGTAGCTCAGAGACGTAAAc  >  1:273500/1‑251 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                 TCCGACAGGTTGCTGATGCCGTCAAAGCCGGTAAAGGAGAAGCACAGGATTGTCGCCCCGGTAATCATCGGGATAACATGCGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAG‑CGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAGGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGACAGTCGGGCTAATGGATTTCTGGGCGTAAGTGTATGCAGAGCCAGCAGAAGGATAGCGGCGAACCAGCTTCCCGTAGCTCAGAGCCGTAAAC  >  NC_000913/2086389‑2086874

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