Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	989,482	(T)5→4	coding (875/1203 nt)	pncB ←	nicotinate phosphoribosyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	989,478	0	T	.	100.0%	159.1 / NA	37	coding (879/1203 nt)	pncB	nicotinate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (27/10);  total (27/10)

GCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGTCGCCAGAGTCATGACGCAGGCCCTGATACCGACTAGCGAACTCGA  >  NC_000913/989400‑989545                                                                               |                                                                    gCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCtt                                                                    <  1:202668/80‑1 (MQ=255)   gTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCttttc                                                                 <  1:202669/81‑2 (MQ=255)   gTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCtttt                                                                  >  2:202718/1‑80 (MQ=255)       aTCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGC‑TTTTCAccc                                                            >  2:202708/1‑81 (MQ=255)          cAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGC‑TTTTCACCCCAt                                                         <  1:202670/81‑1 (MQ=255)           aGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGC‑TTTTCACCCCAtt                                                        >  2:202750/1‑81 (MQ=255)            gATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGC‑TTTTCACCCCATTc                                                       <  2:853763/81‑1 (MQ=255)               tGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGC‑TTTTCACCCCATTCAAc                                                    >  1:202698/1‑81 (MQ=255)                gTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGC‑TTTTCACCCCATTCAAcc                                                   >  2:202724/1‑81 (MQ=255)                      gaAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGC‑TTTTCACCCCATTCAACCGGGTCg                                             <  1:202671/81‑1 (MQ=255)                      gaAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGC‑TTTTCACCCCATTCAACCGGGTCg                                             >  2:202714/1‑81 (MQ=255)                       aaaaCCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGC‑TTTTCACCCCATTCAACCGGGTCGc                                            <  2:202672/81‑1 (MQ=255)                       aaaaCCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGC‑TTTTCACCCCATTCAACCGGGTCGc                                            >  2:202706/1‑81 (MQ=255)                       aaaaCCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGC‑TTTTCACCCCATTCAACCGGGTCGc                                            >  1:202720/1‑81 (MQ=255)                             aGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGC‑TTTTCACCCCATTCAACCGGGTCGCCAGAGt                                      >  1:202702/1‑81 (MQ=255)                              gCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGC‑TTTTCACCCCATTCAACCGGGTCGCCAGAGt                                      >  1:202692/1‑80 (MQ=255)                                gTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGC‑TTTTCACCCCATTCAACCGGGTCGCCAGAGTCAt                                   <  1:202673/81‑1 (MQ=255)                                  tttACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGC‑TTTTCACCCCATTCAACCGGGTCGCCAGAGTCATGa                                 >  1:202723/1‑81 (MQ=255)                                  tttACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGC‑TTTTCACCCCATTCAACCGGGTCGCCAGAGTCATGa                                 >  2:202696/1‑81 (MQ=255)                                        ctGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGC‑TTTTCACCCCATTCAACCGGGTCGCCAGAGTCATGACGCAgg                           >  1:202695/1‑81 (MQ=255)                                        ctGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGC‑TTTTCACCCCATTCAACCGGGTCGCCAGAGTCATGACGCAgg                           >  1:202701/1‑81 (MQ=255)                                         tgtgGATCAATTCCCAGCTTTTCATAATGTGCAATGGC‑TTTTCACCCCATTCAACCGGGTCGCCAGAGTCATGACGCAGGc                          >  2:202728/1‑81 (MQ=255)                                            ggATCAATTCCCAGCTTTTCATAATGTGCAATGGC‑TTTTCACCCCATTCAACCGGGTCGCCAGAGTCATGACGCAGGCCCt                       >  1:202704/1‑81 (MQ=255)                                             gATCAATTCCCAGCTTTTCATAATGTGCAATGGC‑TTTTCACCCCATTCAACCGGGTCGCCAGAGTCATGACGCAGGCCCt                       >  1:202697/1‑80 (MQ=255)                                             gATCAATTCCCAGCTTTTCATAATGTGCAATGGC‑TTTTCACCCCATTCAACCGGGTCGCCAGAGTCATGACGCAGGCCCt                       <  1:202674/80‑1 (MQ=255)                                             gATCAATTCCCAGCTTTTCATAATGTGCAATGGC‑TTTTCACCCCATTCAACCGGGTCGCCAGAGTCATGACGCAGGCCCt                       >  1:202700/1‑80 (MQ=255)                                              aTCAATTCCCAGCTTTTCATAATGTGCAATGGC‑TTTTCACCCCATTCAACCGGGTCGCCAGAGTCATGACGCAGGCCCTGa                     >  2:202707/1‑81 (MQ=255)                                               tCAATTCCCAGCTTTTCATAATGTGCAATGGC‑TTTTCACCCCATTCAACCGGGTCGCCAGAGTCATGACGCAGGCCCTGa                     >  1:202716/1‑80 (MQ=255)                                                cAATTCCCAGCTTTTCATAATGTGCAATGGC‑TTTTCACCCCATTCAACCGGGTCGCCAGAGTCATGACGCAGGCCCTGATa                   >  2:202712/1‑81 (MQ=255)                                                  aTTCCCAGCTTTTCATAATGTGCAATGGC‑TTTTCACCCCATTCAACCGGGTCGCCAGAGTCATGACGCAGGCCCTGATAcc                 <  1:202675/81‑1 (MQ=255)                                                  aTTCCCAGCTTTTCATAATGTGCAATGGC‑TTTTCACCCCATTCAACCGGGTCGCCAGAGTCATGACGCAGGCCCTGATAcc                 >  2:202717/1‑81 (MQ=255)                                                   ttCCCAGCTTTTCATAATGTGCAATGGC‑TTTTCACCCCATTCAACCGGGTCGCCAGAGTCATGACGCAGGCCCTGATAcc                 >  1:202722/1‑80 (MQ=255)                                                    tCCCAGCTTTTCATAATGTGCAATGGC‑TTTTCACCCCATTCAACCGGGTCGCCAGAGTCATGACGCAGGCCCTGATACCGa               >  1:202730/1‑81 (MQ=255)                                                     cccAGCTTTTCATAATGTGCAATGGC‑TTTTCACCCCATTCAACCGGGTCGCCAGAGTCATGACGCAGGCCCTGATACCGAc              >  2:202713/1‑81 (MQ=255)                                                     cccAGCTTTTCATAATGTGCAATGGC‑TTTTCACCCCATTCAACCGGGTCGCCAGAGTCATGACGCAGGCCCTGATACCGAc              >  2:202719/1‑81 (MQ=255)                                                      ccAGCTTTTCATAATGTGCAATGGC‑TTTTCACCCCATTCAACCGGGTCGCCAGAGTCATGACGCAGGCCCTGATACCGACt             <  2:202676/81‑1 (MQ=255)                                                        aGCTTTTCATAATGTGCAATGGC‑TTTTCACCCCATTCAACCGGGTCGCCAGAGTCATGACGCAGGCCCTGATACCGACTAg           <  1:202677/81‑1 (MQ=255)                                                         gCTTTTCATAATGTGCAATGGC‑TTTTCACCCCATTCAACCGGGTCGCCAGAGTCATGACGCAGGCCCTGATACCGACTAGc          <  1:202678/81‑1 (MQ=255)                                                               cATAATGTGCAATGGC‑TTTTCACCCCATTCAACCGGGTCGCCAGAGTCATGACGCAGGCCCTGATACCGACTAGCGAACTc    >  2:202703/1‑81 (MQ=255)                                                                  aaTGTGCAATGGC‑TTTTCACCCCATTCAACCGGGTCGCCAGAGTCATGACGCAGGCCCTGATACCGACTAGCGAACTCGa  >  1:202729/1‑80 (MQ=255)                                                                               |                                                                    GCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGTCGCCAGAGTCATGACGCAGGCCCTGATACCGACTAGCGAACTCGA  >  NC_000913/989400‑989545

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