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	98.6%	218.8 / NA	70	*187* (TAG→TAA)	speG	spermidine N(1)‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (35/34);  minor base T (0/1);  total (35/35)
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

TTGAAGGTGAATTGATGCACGAGTTCTTTATTAATGGTCAATATCGTAATGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAGCATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTCACAGCTCTGTTTTAC  >  NC_000913/1656608‑1656875                                                                                                                                         |                                                                                                                                    ttGAAGGTGAATTGATGCACGAGTTCTTTATTAATGGTCAATATCGTAATGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACa                                                                                                                                   <  2:973235/139‑1 (MQ=255)   gAAGGTGAATTGATGCACGAGTTCTTTATTAATGGTCAATATCGTAATGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACAtt                                                                                                                                 <  2:152192/139‑1 (MQ=255)     aGGTGAATTGATGCACGAGTTCTTTATTAATGGTCAATATCGTAATGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATtaa                                                                                                                               <  2:450133/139‑1 (MQ=255)                 gCACGAGTTCTTTATTAATGGTCAATATCGTAATGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTaataataat                                                                                                                        <  1:1176823/134‑1 (MQ=255)                 gCACGAGTTCTTTATTAATGGTCAATATCGTAATGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTaataataat                                                                                                                        >  2:1176823/1‑134 (MQ=255)                  cACGAGTTCTTTATTAATGGTCAATATCGTAATGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAatcgatcg                                                                                                                  >  2:287560/1‑139 (MQ=255)                           tttATTAATGGTCAATATCGTAATGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGa                                                                                                         <  1:1382976/139‑1 (MQ=255)                           tttATTAATGGTCAATATCGTAATGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGa                                                                                                         <  2:554409/139‑1 (MQ=255)                               ttaATGGTCAATATCGTAATGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGa                                                                                                         <  1:561834/135‑1 (MQ=255)                               ttaATGGTCAATATCGTAATGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGa                                                                                                         >  2:561834/1‑135 (MQ=255)                                taATGGTCAATATCGTAATGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACa                                                                                                                                   <  2:128085/108‑1 (MQ=255)                                     gTCAATATCGTAATGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGa                                                                                                         <  2:235024/129‑1 (MQ=255)                                     gTCAATATCGTAATGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGa                                                                                                         <  1:759278/129‑1 (MQ=255)                                     gTCAATATCGTAATGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGa                                                                                                         >  1:235024/1‑129 (MQ=255)                                     gTCAATATCGTAATGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGa                                                                                                         >  2:759278/1‑129 (MQ=255)                                         atatCGTAATGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGt                                                                                           >  1:571558/1‑139 (MQ=255)                                               tAATGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGAcg                                                                                     >  1:749734/1‑139 (MQ=255)                                                 aTGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGAcggc                                                                                   <  2:1133410/139‑1 (MQ=255)                                                  tGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGa                                                                                                         >  2:113876/1‑116 (MQ=255)                                                  tGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGa                                                                                                         <  1:113876/116‑1 (MQ=255)                                                  tGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGAcggcg                                                                                  >  2:1325046/1‑139 (MQ=255)                                                       ttCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGa                                                                                                         >  1:353264/1‑111 (MQ=255)                                                       ttCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGa                                                                                                         <  2:353264/111‑1 (MQ=255)                                                                  atatTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCACCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGCCGGCGGGTTTTACTCCTtctt                                                                  <  2:860473/139‑3 (MQ=255)                                                                    atTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAatcgatc                                                                                                                   <  1:553519/88‑1 (MQ=255)                                                                    atTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAatcgatc                                                                                                                   >  2:553519/1‑88 (MQ=255)                                                                       ccAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATcaaca                                                             <  1:287560/139‑1 (MQ=255)                                                                        cAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACa                                                                                                                                   >  2:817646/1‑68 (MQ=255)                                                                         aGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAAt                                                           >  2:175304/1‑139 (MQ=255)                                                                         aGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAAt                                                           >  2:1383111/1‑139 (MQ=255)                                                                          gCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATcaaca                                                             <  1:513392/136‑1 (MQ=255)                                                                          gCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATcaaca                                                             >  2:513392/1‑136 (MQ=255)                                                                           cATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGAcggc                                                                                   <  1:1196298/113‑1 (MQ=255)                                                                           cATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGAcggc                                                                                   >  2:1196298/1‑113 (MQ=255)                                                                           cATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGa                                                         >  1:1369106/1‑139 (MQ=255)                                                                           cATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGa                                                         >  1:303442/1‑139 (MQ=255)                                                                           cATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGa                                                         >  2:1249353/1‑139 (MQ=255)                                                                           cATCAGTACCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGa                                                         >  2:569548/1‑139 (MQ=255)                                                                                gTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAATCCGTTCGAcg                                                                                     >  2:928181/1‑106 (MQ=255)                                                                                gTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATcaa                                                               <  1:1405035/128‑1 (MQ=255)                                                                                gTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATcaa                                                               >  2:1405035/1‑128 (MQ=255)                                                                                gTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCAGTTCGAcg                                                                                     <  1:928181/106‑1 (MQ=255)                                                                                 tATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATATCATTAATAATAATCGTTCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCt                                                   <  2:303442/139‑1 (MQ=255)                                                                                   tCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGa                                                 <  2:749734/139‑1 (MQ=255)                                                                                   tCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGa                                                 <  2:571558/139‑1 (MQ=255)                                                                                   ccTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGa                                                 <  1:569548/138‑1 (MQ=255)                                                                                     tGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACa                                               <  1:1249353/139‑1 (MQ=255)                                                                                      ggCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACa                                                                                                                                   >  2:1192111/1‑54 (MQ=255)                                                                                      ggCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACa                                                                                                                                   <  1:1192111/54‑1 (MQ=255)                                                                                      ggCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGa                                                                                                         >  1:1277623/1‑80 (MQ=255)                                                                                      ggCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGa                                                                                                         <  2:1277623/80‑1 (MQ=255)                                                                                           agCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCAc                                         >  1:1330995/1‑139 (MQ=255)                                                                                           agCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCAc                                         >  1:469477/1‑139 (MQ=255)                                                                                           agCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCAc                                         >  2:773161/1‑139 (MQ=255)                                                                                             cacaAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCg                                       >  2:471835/1‑139 (MQ=255)                                                                                                       cGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTg                             >  2:778231/1‑139 (MQ=255)                                                                                                        gggTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGa                            >  1:892890/1‑139 (MQ=255)                                                                                                                        aaGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGa                            <  1:1031208/123‑1 (MQ=255)                                                                                                                        aaGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGa                            >  2:1031208/1‑123 (MQ=255)                                                                                                                          gCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACt                                                                         <  2:479603/76‑1 (MQ=255)                                                                                                                          gCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACt                                                                         >  1:479603/1‑76 (MQ=255)                                                                                                                           ccGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTtcat                                                                  <  2:248374/82‑1 (MQ=255)                                                                                                                           ccGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTtcat                                                                  >  1:248374/1‑82 (MQ=255)                                                                                                                             gACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTtcat                                                                  <  2:75899/80‑1 (MQ=255)                                                                                                                             gACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTtcat                                                                  >  1:75899/1‑80 (MQ=255)                                                                                                                             gACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTCACAGCTCTGt       <  1:778231/139‑1 (MQ=255)                                                                                                                              aCCGCACAATAACATTAATAATATTCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGGTTAGCGTGACTGTCATATTCACAGCTCTGtt      <  2:1369106/139‑1 (MQ=255)                                                                                                                                 gCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTCACAGCTCTGTTTTa   <  2:892890/139‑1 (MQ=255)                                                                                                                                  cacaATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTCACAGCTCTGTTTTAc  <  2:1032647/139‑1 (MQ=255)                                                                                                                                  cacaATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTCACAGCTCTGTTTTAc  >  1:1032647/1‑139 (MQ=255)                                                                                                                                         |                                                                                                                                    TTGAAGGTGAATTGATGCACGAGTTCTTTATTAATGGTCAATATCGTAATGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAGCATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTCACAGCTCTGTTTTAC  >  NC_000913/1656608‑1656875

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