Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,656,744	G→A	*187* (TAG→TAA)	speG →	spermidine N(1)‑acetyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,656,744	0	G	A	100.0%	205.4 / NA	63	*187* (TAG→TAA)	speG	spermidine N(1)‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (32/31);  total (32/31)

ATTGATGCACGAGTTCTTTATTAATGGTCAATATCGTAATGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAGCATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCC‑TTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTCACAGCTCTGTT  >  NC_000913/1656618‑1656871                                                                                                                               |                                                                                                                                 aTTGATGCACGAGTTCTTTATTAATGGTCAATATCGTAATGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATtaataataa                                                                                                                      <  1:768429/139‑1 (MQ=255)     aTGCACGAGTTCTTTATTAATGGTCAATATCGTAATGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAatcga                                                                                                                  >  1:1015321/1‑139 (MQ=255)       gCACGAGTTCTTTATTAATGGTCAATATCGTAATGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAatcgatc                                                                                                                <  1:964887/139‑1 (MQ=255)                cTTTATTAATGGTCAATATCGTAATGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTg                                                                                                       <  1:588382/139‑1 (MQ=255)                       aaTGGTCAATATCGTAATGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGgtgt                                                                                                >  2:198597/1‑139 (MQ=255)                           gTCAATATCGTAATGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACa                                                                                                                                >  1:160094/1‑103 (MQ=255)                           gTCAATATCGTAATGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACa                                                                                                                                <  2:160094/103‑1 (MQ=255)                           gTCAATATCGTAATGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGa                                                                                                      <  2:347815/129‑1 (MQ=255)                           gTCAATATCGTAATGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGa                                                                                                      >  1:347815/1‑129 (MQ=255)                           gTCAATATCGTAATGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAAc                                                                                            >  1:228833/1‑139 (MQ=255)                                tatCGTAATGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATtaa                                                                                                                            <  2:39263/102‑1 (MQ=255)                                tatCGTAATGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATtaa                                                                                                                            >  1:39263/1‑102 (MQ=255)                                   cGTAATGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACAt                                                                                                                               >  2:425672/1‑96 (MQ=255)                                   cGTAATGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACAt                                                                                                                               <  1:425672/96‑1 (MQ=255)                                   cGTAATGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGa                                                                                                      <  2:439136/121‑1 (MQ=255)                                   cGTAATGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGa                                                                                                      >  1:439136/1‑121 (MQ=255)                                     tAATGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGAcg                                                                                  <  2:228833/139‑1 (MQ=255)                                      aaTGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACa                                                                                                                                >  1:499475/1‑92 (MQ=255)                                        tGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGa                                                                                                      >  2:919713/1‑116 (MQ=255)                                        tGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGa                                                                                                      <  1:919713/116‑1 (MQ=255)                                              tCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGtttt                                                                         >  1:64618/1‑139 (MQ=255)                                                        atatTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCC‑Ttcat                                                              <  1:329086/139‑1 (MQ=255)                                                         tatTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAatcga                                                                                                                  >  1:289383/1‑87 (MQ=255)                                                         tatTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAatcga                                                                                                                  <  2:289383/87‑1 (MQ=255)                                                            tCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCC‑TTCATcaac                                                          >  2:63753/1‑139 (MQ=255)                                                               aGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGa                                                                                                      >  1:474332/1‑93 (MQ=255)                                                               aGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGa                                                                                                      <  2:474332/93‑1 (MQ=255)                                                               aGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCC‑TTCATCAACAAt                                                       >  2:564128/1‑139 (MQ=255)                                                               aGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCCCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGa                                                                                                      <  2:30811/93‑1 (MQ=255)                                                               aGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCCCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGa                                                                                                      >  1:30811/1‑93 (MQ=255)                                                                 cATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACa                                                                                                                                <  2:87733/65‑1 (MQ=255)                                                                 cATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGa                                                                                                      >  2:279981/1‑91 (MQ=255)                                                                 cATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGa                                                                                                      <  1:279981/91‑1 (MQ=255)                                                                 cATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCC‑Ttcat                                                              <  1:129017/130‑1 (MQ=255)                                                                 cATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCC‑Ttcat                                                              >  2:129017/1‑130 (MQ=255)                                                                      gTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACt                                                                      <  1:808860/118‑1 (MQ=255)                                                                      gTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACt                                                                      >  2:808860/1‑118 (MQ=255)                                                                                 agCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTg                                                                                                       <  2:16660/74‑1 (MQ=255)                                                                                 agCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTg                                                                                                       >  1:16660/1‑74 (MQ=255)                                                                                 agCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGa                                                                                                      <  1:491996/75‑1 (MQ=255)                                                                                 agCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGa                                                                                                      >  2:561635/1‑75 (MQ=255)                                                                                 agCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGa                                                                                                      >  2:491996/1‑75 (MQ=255)                                                                                 agCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGa                                                                                                      <  1:561635/75‑1 (MQ=255)                                                                                 agCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACt                                                                      <  1:268244/107‑1 (MQ=255)                                                                                 agCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACt                                                                      >  2:268244/1‑107 (MQ=255)                                                                                 agCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCC‑TTCATCAACAATGATTAGCTGACAGTCCAc                                     >  1:573806/1‑139 (MQ=255)                                                                                             cGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCC‑TTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTg                         >  1:1115809/1‑139 (MQ=255)                                                                                              gggTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCC‑TTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGa                        <  2:573806/139‑1 (MQ=255)                                                                                              gggTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCC‑TTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGa                        <  1:63753/139‑1 (MQ=255)                                                                                                   aGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCC‑TTCATCACCAATGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTc                   <  2:207322/139‑1 (MQ=255)                                                                                                         tcCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCC‑TTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTc             >  1:494264/1‑139 (MQ=255)                                                                                                            tGAGGCCGCCCGAAAATTACCATTAATAATAATCATTCTTATTTTGGATGGGGT‑AACCCGTTCGACGGCGGGTTTTACTCC‑TTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTCACAg         <  2:177415/139‑1 (MQ=255)                                                                                                               aGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCC‑TTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTCACAGct       <  2:1115809/139‑1 (MQ=255)                                                                                                               aGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCC‑TTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTCACAGct       <  1:173359/139‑1 (MQ=255)                                                                                                                gCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCC‑TTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTCACAGctc      >  1:661895/1‑139 (MQ=255)                                                                                                                 ccGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGa                                                                                                      <  2:946685/43‑1 (MQ=255)                                                                                                                 ccGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGa                                                                                                      >  1:946685/1‑43 (MQ=255)                                                                                                                   gACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGAcggc                                                                                >  2:418671/1‑63 (MQ=255)                                                                                                                   gACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGAcggc                                                                                <  1:418671/63‑1 (MQ=255)                                                                                                                   gACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCC‑TTCATCAACAATGATTAGCTGACAg                                          >  1:936700/1‑100 (MQ=255)                                                                                                                   gACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCC‑TTCATCAACAATGATTAGCTGACAg                                          <  2:936700/100‑1 (MQ=255)                                                                                                                   gACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCC‑TTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCg                           >  2:6157/1‑115 (MQ=255)                                                                                                                   gACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCC‑TTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCg                           <  1:6157/115‑1 (MQ=255)                                                                                                                   gACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGT‑AACCCCTTTGACGGCGGGTTTTACTCCTTTAATAAACAA‑GATTAGCTTACAGTCCACCGGGTTAGCCTGCACGTCATATTTACAGCCCTGtt  >  1:177415/1‑139 (MQ=255)                                                                                                                               |                                                                                                                                 ATTGATGCACGAGTTCTTTATTAATGGTCAATATCGTAATGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAGCATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCC‑TTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTCACAGCTCTGTT  >  NC_000913/1656618‑1656871

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