Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	136,457	(T)8→7	coding (8/867 nt)	speE ←	spermidine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	136,450	0	T	.	96.9%	126.6 / ‑2.2	32	coding (15/867 nt)	speE	spermidine synthase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base . (17/14);  minor base C (1/0);  total (18/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 = 7.11e-01

CCAGCGCCATTACGCGACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGAT‑ACCT‑CCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTA  >  NC_000913/136317‑136581                                                                                                                                      |                                                                                                                                      ccAGCGCCATGACGCGACCAAATGCGGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAGAACGTTATCTACCGCAAAGTACTGCCCAAACGGGTCGTGTAGCGTTTCATGCCACTGtttttttcgg                                                                                                                              <  2:625019/143‑4 (MQ=255)       ccATTCCTCGACGAATTGCCCCGTTCTCGAACATGAGCAACTCCTCGGCATCGGCCTGTCCATGATACAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGCCGTGTAGCGTTTCATGCCACTG‑TTTTTTTCGGCCATgg                                                                                                                       <  2:813963/143‑1 (MQ=255)                  ccAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTG‑TTTTTTTCGGCCATGGGTTGAT‑Acct‑c                                                                                                          <  2:760726/143‑1 (MQ=255)                     aaTGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGCTACCGGTCGTTAGCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTG‑TTTTTTTCGGCCATGGGTTGAT‑ACCT‑CCtt                                                                                                       <  2:105150/143‑1 (MQ=255)                     aaTGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTGCATGCCGCTG‑TTTTTTTCGGCCATGGGTTGAT‑Acat‑ctct                                                                                                       >  2:95755/1‑140 (MQ=255)                           gCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTG‑TTTTTTTCGGCCATGGGTTGAT‑ACCT‑CCTTTGTTaa                                                                                                 >  1:706678/1‑143 (MQ=255)                              ttCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTG‑TTTTTTTCGGGCATGGGTTGAT‑ACGT‑CCTTTGttgaaag                                                                                              >  1:480450/1‑138 (MQ=255)                              ttCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAGAACGTTATCTACCGCAAAGTACTGCCCAAAATGGTCGGGTAGCGTTTCATGCCACTG‑TTTTTTTCGGCCATGGGTTGAT‑ACAT‑GCGTTGTTGacac                                                                                              >  1:161180/1‑143 (MQ=255)                                ctcAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTG‑TTTTTTTCGGCCATGGGTGGAT‑ACCT‑CCTTTGTAAACAccc                                                                                            >  2:853453/1‑143 (MQ=255)                                  ccaaaaTGATCAGTTCCTTGCGATCGGGCGTTGCATGATACATAACGTTATCTACCGCAGAGTACTGCCCACGCTGGTCGTGTAGCGTTCCATGCCACTG‑TTTTTTTCGGCCATGGGTTGAT‑ACCT‑CCTTTGTTAACACCCGt                                                                                          <  2:660408/141‑1 (MQ=255)                                  cAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTG‑TTTTTTTCGGCCATGGGATGAT‑ACCT‑CCGTTGTTAACACCCGt                                                                                          >  1:151075/1‑143 (MQ=255)                                                  cTGGTGATCGGTCTTTTCATGATACAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTG‑TTTTTTTCGGCCATGGGTTGAT‑ACCT‑CCTTTGTTAACACCCGTAAAAAAAGGGCGCAAc                                                                          <  1:191618/143‑1 (MQ=255)                                                         ccGGTCTCTCCATCATACAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTG‑TTTTTTTCGGCCATGGGTTGAT‑ACCT‑CCTTTGTTAACTCCCGTAAAAAAAGGGCGCAACATAATAg                                                                   <  1:119413/142‑1 (MQ=255)                                                           ggTCTGTTCATGATGCTGTACGTTATCTACCGCGAAGTACTGCGCAAACTGGTCGTGTAGCGTTTCACGCCGCTG‑TTTTTTTCGGCCATGGGTTGAT‑ACCT‑CCTTTGTTAACACCCGTTAAAAAAGGGCGCAACATAATAGCt                                                                 <  2:661384/143‑1 (MQ=255)                                                             tCTTTTCATGATACAGAACGGTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTG‑TTTTTTTCGGCCATGGGTTGAT‑ACCT‑CCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTaa                                                               <  1:403227/143‑1 (MQ=255)                                                              cTTTTCATGATACAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTG‑TTTTTTTCGGCCATGGGTTGAT‑ACCT‑CCTTTGTTAACACCCGTGAAACAAGGGCGCAACATCATAGCTAAc                                                              >  2:207259/1‑143 (MQ=255)                                                               ttttCATGATACAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTG‑TTTTTTTCGGCCATGGGTTGAT‑ACCT‑CCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACa                                                             >  1:606673/1‑143 (MQ=255)                                                                  tCATGATAAAGAACGTTATCTACCACACAGTAATGCTCAAACTGGTCGTATAGCGTTTCATGACACTG‑TTTTTTTCAGCCATGGGTTGAT‑ACCT‑CCTTTGTTAACACCAGTAAAAAAAGGGCGCAACATAATAGCTAACATTg                                                          <  2:95311/143‑1 (MQ=255)                                                                  ccATGATACAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTG‑TTTTTTTCGGCCATGGGTTGAT‑ACCT‑CCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTg                                                          <  1:78328/142‑1 (MQ=255)                                                                     tGATACAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTG‑TTTTTTTCGGCCATGGGTTGAT‑ACCT‑CCTTTGTTAACACCCGTAAACAAAGGGCGCAACATAATAGCAAACATTGAcc                                                       >  1:653513/1‑143 (MQ=255)                                                                         aCAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTG‑TTTTTTTCGGCCATGGGTTGGTAAGATGACCAAGTT‑‑CTCGCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCgg                                                   >  1:207896/1‑143 (MQ=255)                                                                         aCAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTG‑TTTTTTTCGACCATGGGTGGAT‑ACCT‑CCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCgg                                                   >  1:818896/1‑143 (MQ=255)                                                                             aaCGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTG‑TTTTTTTCGGCCATGGGTTGAT‑ACCT‑CCTTTGTTAACACCCGT‑AAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCa                                              >  1:796258/1‑143 (MQ=255)                                                                                           ctacgcaCTGCCCAAACTCGTCGTGTAGCGTTTCATGCCACCG‑CTCTTTTCGGCCATGGGTTGAT‑ACCT‑CCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTc                                 <  2:872115/137‑1 (MQ=255)                                                                                           cAAACTCCTGCCCAAACTCGTCGTGTAGCGTTTCATGCCACTG‑TTTTTTTCGGCCATGGGTTGAT‑ACCT‑CCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTc                                 <  1:366097/143‑1 (MQ=255)                                                                                            aaaCTACTGCCCCAACTGGTCGTGTAGCGCTTCATGCCACTG‑TTTTTTTCGGCCATGGGTTGAT‑ACCT‑CCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCt                                <  1:248528/143‑1 (MQ=255)                                                                                             aaGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGCTTTTTTT‑TTTTTTTTTTTGAT‑ACCT‑ACTTTGTTAACACCCGTAAATAAAGGGCGCGAGATAATAGCTAATATTGACCGCGGATGCACGGTCAACATTCCTa                               >  2:369372/1‑143 (MQ=255)                                                                                                     cccAAACTGGTCGTGTAGCGTTTCATGCCACTG‑TTTTTTTCGGCCATGGGTTGTT‑ACCT‑CCTTTGTTCACATCCGTTAAAAGTGGGCGCATCGTAATAGCTAACATTGACCGCGGATGCACGGTCGATATTTCTAGCAAAAgg                       >  2:760865/1‑143 (MQ=255)                                                                                                             ggTCGTGTAGCGTTTCATGCCACTG‑TTTTTTTCGGCCATGGGTTGAT‑ACCT‑CCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGCCCGCGGATGCACGGTCAATATTTCTTCACAAAGGTATCAggg               >  2:466446/1‑143 (MQ=255)                                                                                                                      gCGTTTCATGCCACTG‑TTTTTTTCGGCCATGGGTTGAT‑ACCT‑CCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACTGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATGATTTgg      >  2:534174/1‑142 (MQ=255)                                                                                                                      gCGTTTCATGCCACTG‑TTTTTTTCGGCCATGGGTTGAT‑ACCT‑CCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGa      >  1:595912/1‑143 (MQ=255)                                                                                                                      gCGTTTCATGCCACTG‑TTTTTGTCGGCCATGGGTTGAT‑ACCG‑CCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGa      <  1:357524/143‑1 (MQ=255)                                                                                                                          ttCATGCCACTG‑TTTTTTTCGGCCATGGGTTGAT‑ACCT‑CCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCGATATTTCTACTAAAAGGTATCAGGCCTTATTTGACGTa  >  2:287068/1‑143 (MQ=255)                                                                                                                                      |                                                                                                                                      CCAGCGCCATTACGCGACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGAT‑ACCT‑CCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTA  >  NC_000913/136317‑136581

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