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%	243.3 / NA	75	*187* (TAG→TAA)	speG	spermidine N(1)‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (39/36);  total (39/36)

TTGAAGGTGAATTGATGCACGAGTTCTTTATTAATGGTCAATATCGTAATGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAGCATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTCACAGCTCTGTTTTA  >  NC_000913/1656608‑1656874                                                                                                                                         |                                                                                                                                   ttGAAGGTGAATTGATGCACGAGTTCTTTATTAATGGTCAATATCGTAATGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACa                                                                                                                                  <  2:162675/139‑1 (MQ=255)                 gCACGAGTTCTTTATTAATGGTCAATATCGTAATGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAatcgatc                                                                                                                  >  2:248306/1‑139 (MQ=255)                 gCACGAGTTCTTTATTAATGGTCAATATCGTAATGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAatcgatc                                                                                                                  >  2:136604/1‑139 (MQ=255)                  cACGAGTTCTTTATTAATGGTCAATATCGTAATGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAatcgatcg                                                                                                                 >  2:1348353/1‑139 (MQ=255)                               ttaATGGTCAATATCGTAATGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGgt                                                                                                    >  2:1458309/1‑139 (MQ=255)                                                aaTGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGAcgg                                                                                   <  1:622373/139‑1 (MQ=255)                                                aaTGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGAcgg                                                                                   >  2:622373/1‑139 (MQ=255)                                                 aTGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCATAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGAcggc                                                                                  <  1:654000/139‑1 (MQ=255)                                                 aTGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGAcggc                                                                                  <  1:136604/139‑1 (MQ=255)                                                    ccATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCggg                                                                               >  1:1591413/1‑139 (MQ=255)                                                      aTTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAACACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGtt                                                                             >  1:1334319/1‑139 (MQ=255)                                                       ttCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGAcg                                                                                    >  2:571298/1‑131 (MQ=255)                                                       ttCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGAcg                                                                                    <  1:571298/131‑1 (MQ=255)                                                              gtgtATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCt                                                                     <  1:248306/139‑1 (MQ=255)                                                                gtATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGa                                                                                      <  2:545743/120‑1 (MQ=255)                                                                gtATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGa                                                                                      >  1:545743/1‑120 (MQ=255)                                                                     ttCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATcaa                                                              <  2:1334319/139‑1 (MQ=255)                                                                      tCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACAt                                                                                                                                 >  1:709310/1‑71 (MQ=255)                                                                      tCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACAt                                                                                                                                 <  2:709310/71‑1 (MQ=255)                                                                        cAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGa                                                                                                        <  1:203191/94‑1 (MQ=255)                                                                        cAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGa                                                                                                        >  2:203191/1‑94 (MQ=255)                                                                        cAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGAcggc                                                                                  <  1:220445/116‑1 (MQ=255)                                                                        cAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGAcggc                                                                                  >  2:220445/1‑116 (MQ=255)                                                                         aGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGa                                                                                                        >  2:544473/1‑93 (MQ=255)                                                                         aGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGa                                                                                                        >  2:239990/1‑93 (MQ=255)                                                                         aGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGa                                                                                                        <  1:239990/93‑1 (MQ=255)                                                                         aGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGa                                                                                                        <  1:544473/93‑1 (MQ=255)                                                                           cATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATg                                                                                                      <  2:381274/93‑1 (MQ=255)                                                                           cATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATg                                                                                                      >  1:381274/1‑93 (MQ=255)                                                                           cATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGa                                                        >  1:479422/1‑139 (MQ=255)                                                                                   tCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGa                                                <  2:1591413/139‑1 (MQ=255)                                                                                           agCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACa                                                                                                                                  <  2:896223/49‑1 (MQ=255)                                                                                           agCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACa                                                                                                                                  <  2:1556158/49‑1 (MQ=255)                                                                                           agCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACa                                                                                                                                  >  1:896223/1‑49 (MQ=255)                                                                                           agCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACa                                                                                                                                  >  1:1556158/1‑49 (MQ=255)                                                                                           agCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTg                                                                                                         <  2:912464/74‑1 (MQ=255)                                                                                           agCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGa                                                                                                        <  1:1384777/75‑1 (MQ=255)                                                                                           agCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGa                                                                                                        >  2:1384777/1‑75 (MQ=255)                                                                                           agCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTtcatc                                                                >  2:255048/1‑115 (MQ=255)                                                                                           agCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTtcatc                                                                <  1:255048/115‑1 (MQ=255)                                                                                           agCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCAc                                        >  2:764461/1‑139 (MQ=255)                                                                                           agCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCAc                                        >  1:1525135/1‑139 (MQ=255)                                                                                           agCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGAACGTATTTTTg                                                                                                         >  1:912464/1‑74 (MQ=255)                                                                                             cacaAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTg                                                                                                         <  2:299882/72‑1 (MQ=255)                                                                                             cacaAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTg                                                                                                         >  1:299882/1‑72 (MQ=255)                                                                                             cacaAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGa                                                                                                        >  1:1029789/1‑73 (MQ=255)                                                                                             cacaAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGa                                                                                                        <  2:1029789/73‑1 (MQ=255)                                                                                               caAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTTGAcggc                                                                                  <  2:1185157/93‑1 (MQ=255)                                                                                               caAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTTGAcggc                                                                                  >  1:1185157/1‑93 (MQ=255)                                                                                                 aaaCACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTAGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGAtt                                  >  1:20432/1‑139 (MQ=255)                                                                                                       cGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGttt                                                                            <  1:1117939/91‑1 (MQ=255)                                                                                                       cGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGttt                                                                            >  2:1117939/1‑91 (MQ=255)                                                                                                        gggTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGa                           <  1:764461/139‑1 (MQ=255)                                                                                                          gTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGTATTAGCGTGACt                         <  1:1458309/139‑1 (MQ=255)                                                                                                            cAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGa                           <  1:1236481/135‑1 (MQ=255)                                                                                                            cAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGa                           >  2:1236481/1‑135 (MQ=255)                                                                                                               aCTCTCCTGAAGCCGACCGCACAATAACATTAATAATAatcgatc                                                                                                                  >  2:1498864/1‑45 (MQ=255)                                                                                                               aCTCTCCTGAAGCCGACCGCACAATAACATTAATAATAatcgatc                                                                                                                  <  1:1498864/45‑1 (MQ=255)                                                                                                                   tcCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTc                >  2:1329732/1‑139 (MQ=255)                                                                                                                    ccTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGa                                                                                                        <  2:476583/50‑1 (MQ=255)                                                                                                                    ccTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGa                                                                                                        >  1:476583/1‑50 (MQ=255)                                                                                                                    ccTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACt                                                                        >  2:995569/1‑82 (MQ=255)                                                                                                                    ccTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACt                                                                        <  1:995569/82‑1 (MQ=255)                                                                                                                             gACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCggg                                                                               >  2:1449419/1‑66 (MQ=255)                                                                                                                             gACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCggg                                                                               <  1:1449419/66‑1 (MQ=255)                                                                                                                             gACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGt                             <  2:1038675/116‑1 (MQ=255)                                                                                                                             gACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGt                             >  1:1038675/1‑116 (MQ=255)                                                                                                                             gACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGa                           <  2:1869046/118‑1 (MQ=255)                                                                                                                             gACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGa                           >  1:1869046/1‑118 (MQ=255)                                                                                                                             gACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTCACAGCTCTGt      >  1:1324745/1‑139 (MQ=255)                                                                                                                             gACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCATCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTCACAGCTCTGt      >  1:834601/1‑139 (MQ=255)                                                                                                                             gACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCATCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTCACAGCTCTGt      <  2:834601/139‑1 (MQ=255)                                                                                                                               ccGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGa                                    >  2:23272/1‑107 (MQ=255)                                                                                                                               ccGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGa                                    <  1:23272/107‑1 (MQ=255)                                                                                                                                 gCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGATTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTCACAGCTCTGTTTTa  <  1:1348353/139‑1 (MQ=255)                                                                                                                                         |                                                                                                                                   TTGAAGGTGAATTGATGCACGAGTTCTTTATTAATGGTCAATATCGTAATGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAGCATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTCACAGCTCTGTTTTA  >  NC_000913/1656608‑1656874

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: ‑