Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,851,894	(C)5→6	coding (8/642 nt)	pncA →	nicotinamidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,851,889	1	.	C	100.0%	189.1 / NA	53	M1M (ATG→ATC) †	pncA	nicotinamidase
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (31/22);  total (31/22)

TGTAGAAGCCACGCTAACCAAACTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCCAAAACCTGCGCGGCGAACTGACGCCGGATGATTAAGGAGACTGTAATG‑CCCCCTCGCGCCCTGTTACTGGTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGC  >  NC_000913/1851762‑1852023                                                                                                                                 |                                                                                                                                       tGTAGAAGCCACGCTAACCTAACTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCCAAAACCTGCGCGGCGAACTGACGCCGGATGATTAAGGAGACTGTAATGCCCCCCTCgcgc                                                                                                                             <  2:448457/140‑1 (MQ=255)  gTAGAAGCCACGCTAACCAAACTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCCAAAACCTGCGCGGCGAACTGACGCCGGATGATTAAGGAGACTGTAATGCCCCCCTCGCGcc                                                                                                                            <  1:408702/140‑1 (MQ=255)  gTAGAAGCCACGCTAACCAAACTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCCAAAACCTGCGCGGCGAACTGACGCCGGATGATTAAGGAGACTGTAATGCCCCCCTCGCGcc                                                                                                                            <  2:855066/140‑1 (MQ=255)    aGAAGCCACGCTAACCAAACTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCCAAAACCTGCGCGGCGAACTGACGCCGGATGATTAAGGAGACTGTAATGCCCCCCTcgcg                                                                                                                              <  2:467713/136‑1 (MQ=255)    aGAAGCCACGCTAACCAAACTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCCAAAACCTGCGCGGCGAACTGACGCCGGATGATTAAGGAGACTGTAATGCCCCCCTcgcg                                                                                                                              >  1:467713/1‑136 (MQ=255)    aGAAGCCACGCTAACCAAACTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCCAAAACCTGCGCGGCGAACTGACGCCGGATGATTAAGGAGACTGTAATG‑cccccctcg                                                                                                                               <  2:877259/134‑5 (MQ=255)    aGAAGCCACGCTAACCAAACTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCCAAAACCTGCGCGGCGAACTGACGCCGGATGATTAAGGAGACTGTAATG‑cccccctcg                                                                                                                               >  1:877259/1‑130 (MQ=255)     gAAGCCACGCTAACCAAACTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCCAAAACCTGCGCGGCGAACTGACGCCGGATGATTAAGGAGACTGTAATGCCCCCCTCGCGCCCTg                                                                                                                         >  2:343067/1‑140 (MQ=255)      aaGCCACGCTAACCAAACTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCCAAAACCTGCGCGGCGAACTGACGCCGGATGATTAAGGAGACTGTAATGCCCCCCTCGCGCCCTGt                                                                                                                        >  2:1036579/1‑140 (MQ=255)        gCCACGCTAACCAAACTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCCAAAACCTGCGCGGCGAACTGACGCCGGATGATTAAGGAGACTGTAATGCCCCCCTcgcg                                                                                                                              <  1:292231/132‑1 (MQ=255)        gCCACGCTAACCAAACTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCCAAAACCTGCGCGGCGAACTGACGCCGGATGATTAAGGAGACTGTAATGCCCCCCTcgcg                                                                                                                              >  2:292231/1‑132 (MQ=255)        gCCACGCTAACCAAACTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCCAAAACCTGCGCGGCGAACTGACGCCGGATGATTAAGGAGACTGTAATGCCCCCCTcgcg                                                                                                                              >  2:759094/1‑132 (MQ=255)        gCCACGCTAACCAAACTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCCAAAACCTGCGCGGCGAACTGACGCCGGATGATTAAGGAGACTGTAATGCCCCCCTcgcg                                                                                                                              <  1:759094/132‑1 (MQ=255)        gCCACGCTAACCAAACTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCCAAAACCTGCGCGGCGAACTGACGCCGGATGATTAAGGAGACTGTAATGCCCCCCTCGCGCCCTGTTa                                                                                                                      >  2:339131/1‑140 (MQ=255)           aCGCTAACCAAACTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCCAAAACCTGCGCGGCGAACTGACGCCGGATGATTAAGGAGACTGTAATGCCCCCCTcgcg                                                                                                                              >  2:913730/1‑129 (MQ=255)           aCGCTAACCAAACTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCCAAAACCTGCGCGGCGAACTGACGCCGGATGATTAAGGAGACTGTAATGCCCCCCTcgcg                                                                                                                              >  1:9748/1‑129 (MQ=255)           aCGCTAACCAAACTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCCAAAACCTGCGCGGCGAACTGACGCCGGATGATTAAGGAGACTGTAATGCCCCCCTcgcg                                                                                                                              <  1:913730/129‑1 (MQ=255)           aCGCTAACCAAACTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCCAAAACCTGCGCGGCGAACTGACGCCGGATGATTAAGGAGACTGTAATGCCCCCCTcgcg                                                                                                                              <  2:9748/129‑1 (MQ=255)           aCGCTAACCAAACTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCCAAAACCTGCGCGGCGAACTGACGCCGGATGATTAAGGAGACTGTAATGCCCCCCTCGCGCCCTGTTACTg                                                                                                                   >  1:646725/1‑140 (MQ=255)                    aaaCTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCCAAAACCTGCGCGGCGAACTGACGCCGGATGATTAAGGAGACTGTAATGCCCCCCTCGCGCCCTGTTACTGGTCGATTTa                                                                                                          <  2:646725/140‑1 (MQ=255)                       cTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCCAAAACCTGCGCGGCGAACTGACGCCGGATGATTAAGGAGACTGTAATGCCCCCCTCGCGCCCTGTTACTGGTCGATTTACaa                                                                                                       <  2:498093/140‑1 (MQ=255)                       cTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCCAAAACCTGCGCGGCGAACTGACGCCGGATGATTAAGGAGACTGTAATGCCCCCCTCGCGCCCTGTTACTGGTCGATTTACaa                                                                                                       <  1:862384/140‑1 (MQ=255)                                cTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCCAAAACCTGCGCGGCGAACTGACGCCGGATGATTAAGGAGACTGTAATGCCCCCCTCGCGCCCTGTTACTGGTCGATTTACAAAATGATTTc                                                                                              <  1:1036579/140‑1 (MQ=255)                                 tGCTCAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCCAAAACCTGCGCGGCGAACTGACGCCGGATGATTAAGGAGACTGTAATGCCCCCCTCGCGCCCTGTTACTGGTCGATTTACAAAATGATTTCt                                                                                             >  1:466318/1‑140 (MQ=255)                                  gCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCCAAAACCTGCGCGGCGAACTGACGCCGGATGATTAAGGAGACTGTAATGCCCCCCTCGCGCCCTGTTACTGGTCGATTTACAAAATGATTTCtg                                                                                            >  1:942798/1‑140 (MQ=255)                                    ttAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCCAAAACCTGCGCGGCGAACTGACGCCGGATGATTAAGGAGACTGTAATGCCCCCCTCGCGCCCTGTTACTGGTCGATTTACAAAATGATTTCtgtg                                                                                          <  1:339131/140‑1 (MQ=255)                                             aaCTGGATACTGAAACCATTCGCAAGGCCATGAGCCAAAACCTGCGCGGCGAACTGACGCCGGATGATTAAGGAGACTGTAATGCCCCCCTCGCGCCCTGTTACTGGTCGATTTACAAAATGATTTCTGTGCTGGTGgcg                                                                                 >  2:843941/1‑140 (MQ=255)                                                       tGAAACCATTCGCAAGGCCATGAGCCAAAACCTGCGCGGCGAACTGACGCCGGATGATTAAGGAGACTGTAATGCCCCCCTCGCGCCCTGTTACTGGTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGt                                                                       >  2:309124/1‑140 (MQ=255)                                                            ccATTCGCAAGGCCATGAGCCAAAACCTGCGCGGCGAACTGACGCCGGATGATTAAGGAGACTGTAATGCCCCCCTCGCGCCCTGTTACTGGTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCgg                                                                  >  1:377898/1‑140 (MQ=255)                                                              aTTCGCAAGGCCATGAGCCAAAACCTGCGCGGCGAACTGACGCCGGATGATTAAGGAGACTGTAATGCCCCCCTCGCGCCCTGTTACTGGTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGaa                                                                <  2:785493/140‑1 (MQ=255)                                                                   cAAGGCCATGAGCCAAAACCTGCGCGGCGAACTGACGCCGGATGATTAAGGAGACTGTAATGCCCCCCTCGCGCCCTGTTACTGGTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGa                                                           <  2:942798/140‑1 (MQ=255)                                                                    aaGGCCATGAGCCAAAACCTGCGCGGCGAACTGACGCCGGATGATTAAGGAGACTGTAATGCCCCCCTCGCGCCCTGTTACTGGTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGAc                                                          >  2:335815/1‑140 (MQ=255)                                                                         cATGAGCCAAAACCTGCGCGGCGAACTGACGCCGGATGATTAAGGAGACTGTAATGCCCCCCTCGCGCCCTGTTACTGGTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTAc                                                     >  1:95395/1‑140 (MQ=255)                                                                          aTGAGCCAAAACCTGCGCGGCGAACTGACGCCGGATGATTAAGGAGACTGTAATGCCCCCCTCGCGCCCTGTTACTGGTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACg                                                    <  2:377898/140‑1 (MQ=255)                                                                          aTGAGCCAAAACCTGCGCGGCGAACTGACGCCGGATGATTAAGGAGACTGTAATGCCCCCCTCGCGCCCTGTTACTGGTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACg                                                    >  1:737153/1‑140 (MQ=255)                                                                            gAGCCAAAACCTGCGCGGCGAACTGACGCCGGATGATTAAGGAGACTGTAATGCCCCCCTCGCGCCCTGTTACTGGTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGt                                                  >  1:634792/1‑140 (MQ=255)                                                                              gCCAAAACCTGCGCGGCGAACTGACGCCGGATGATTAAGGAGACTGTAATGCCCCCCTCGCGCCCTGTTACTGGTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTgg                                                <  2:634792/140‑1 (MQ=255)                                                                                  aaaCCTGCGCGGCGAACTGACGCCGGATGATTAAGGAGACTGTAATGCCCCCCTCGCGCCCTGTTACTGGTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAgg                                                              >  1:80819/1‑122 (MQ=255)                                                                                  aaaCCTGCGCGGCGAACTGACGCCGGATGATTAAGGAGACTGTAATGCCCCCCTCGCGCCCTGTTACTGGTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAgg                                                              <  2:80819/122‑1 (MQ=255)                                                                                  aaaCCTGCGCGGCGAACTGACGCCGGATGATTAAGGAGACTGTAATGCCCCCCTCGCGCCCTGTTACTGGTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGt                                            >  1:847111/1‑140 (MQ=255)                                                                                         cgcgGCGAACTGACGCCGGATGATTAAGGAGACTGTAATGCCCCCCTCGCGCCCTGTTACTGGTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAAc                                     >  1:72807/1‑140 (MQ=255)                                                                                         cgcgGCGAACTGACGCCGGATGATTAAGGAGACTGTAATGCCCCCCTCGCGCCCTGTTACTGGTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAAc                                     >  2:32035/1‑140 (MQ=255)                                                                                          gcggcgAACTGACGCCGGATGATTAAGGAGACTGTAATGCCCCCCTCGCGCCCTGTTACTGGTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAAcc                                    <  2:95395/140‑1 (MQ=255)                                                                                             gcgAACTGACGCCGGATGATTAAGGAGACTGTAATGCCCCCCTCGCGCCCTGTTACTGGTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGcc                                 <  1:32035/140‑1 (MQ=255)                                                                                             gcgAACTGACGCCGGATGATTAAGGAGACTGTAATGCCCCCCTCGCGCCCTGTTACTGGTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGc                                  <  1:483003/139‑1 (MQ=255)                                                                                             gcgAACTGACGCCGGATGATTAAGGAGACTGTAATGCCCCCCTCGCGCCCTGTTACTGGTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGc                                  >  2:483003/1‑139 (MQ=255)                                                                                                                   aGGAGACTGTAATGCCCCCCTCGCGCCCTGTTTCTGGTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTcgcg           >  2:778845/1‑140 (MQ=255)                                                                                                                        aCTGTAATGCCCCCCTCGCGCCCTGTTACTGGTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTg      >  1:70734/1‑140 (MQ=255)                                                                                                                          tGTAATGCCCCCCTCGCGCCCTGTTACTGGTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGaa    >  1:1078266/1‑140 (MQ=255)                                                                                                                          tGTAATGCCCCCCTCGCGCCCTGTTACTGGTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGaa    >  1:309407/1‑140 (MQ=255)                                                                                                                          tGTAATGCCCCCCTCGCGCCCTGTTACTGGTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGaa    >  1:538111/1‑140 (MQ=255)                                                                                                                          tGTAATGCCCCCCTCGCGCCCTGTTACTGGTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGaa    >  1:416827/1‑140 (MQ=255)                                                                                                                          tGTAATGCCCCCCTCGCGCCCTGTTACTGGTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGaa    >  1:98825/1‑140 (MQ=255)                                                                                                                           gTAATGCCCCCCTCGCGCCCTGTTACTGGTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAg   <  1:343067/140‑1 (MQ=255)                                                                                                                            tAATGCCCCCCTCGCGCCCTGTTACTGGTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGc  >  1:148884/1‑140 (MQ=255)                                                                                                                                 |                                                                                                                                       TGTAGAAGCCACGCTAACCAAACTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCCAAAACCTGCGCGGCGAACTGACGCCGGATGATTAAGGAGACTGTAATG‑CCCCCTCGCGCCCTGTTACTGGTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGC  >  NC_000913/1851762‑1852023

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