Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	989,500	T→C	D286G (GAC→GGC)	pncB ←	nicotinate phosphoribosyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	989,500	0	T	C	100.0%	185.5 / NA	52	D286G (GAC→GGC)	pncB	nicotinate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (26/26);  total (26/26)

GAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGTCGCCAGAGTCATGACGCAGGCCCTGATACCGACTAGCGAACTCGACACCGAAATCACGCAGGAAAGCATCCATAGTGATGCAGTCGGTTAATGCAATGCCAAGTTGGTCGGGATACTCTTCCAGCCAGGCAGCA  >  NC_000913/989373‑989634                                                                                                                                |                                                                                                                                       gaAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGCCGCCAGAGTCAt                                                                                                                            <  1:258844/140‑1 (MQ=255)     tGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGCCGCCAGAGTCATGACg                                                                                                                        <  2:771042/140‑1 (MQ=255)       gCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGCCGCCAGAGTCATGACGCa                                                                                                                      <  1:910758/140‑1 (MQ=255)          gTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGCCGCCAGAGTCATGACGCAGGc                                                                                                                   >  1:86680/1‑140 (MQ=255)           tataGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGCCGCCAGAGTCATGACGCAGGcc                                                                                                                  <  2:600780/140‑1 (MQ=255)                cTCAACCGCTTTGCCTAAATCCAGATTGTCAGAGAAAAACAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCCGGCCGCCAGAGCCATGACGCAGCCCCTGAt                                                                                                             >  2:522347/1‑140 (MQ=255)                     ccGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGCCGCCAGAGTCATGACg                                                                                                                        >  1:205486/1‑124 (MQ=255)                     ccGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGCCGCCAGAGTCATGACg                                                                                                                        <  2:205486/124‑1 (MQ=255)                       gCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGCCGCCAGAGTCATGACGCAGGCCCTGAt                                                                                                             <  1:618596/133‑1 (MQ=255)                       gCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGCCGCCAGAGTCATGACGCAGGCCCTGAt                                                                                                             >  2:618596/1‑133 (MQ=255)                        cTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAA‑TCCCATCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCCGGCCGCCAGAGTCATGACGCAGGCCCTGATACCGACTAg                                                                                                    >  2:547450/1‑140 (MQ=255)                         tttGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTATGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGCCGCCAGAGTCATGACGCAGGCCCTGATACCGACTAg                                                                                                    >  1:407286/1‑140 (MQ=255)                             cGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGCCGCCAGAGTCATGACGCAGGCCCTGATACCGACTAGCGaa                                                                                                <  2:565584/140‑1 (MQ=255)                                aaaTCCAGATTGTCAGAGAACCCCAGCGTTTTCCTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCATCCGGGCCGTCAGAGTCATGACGCAGGCCCTGATCCCGACTAGCGAACTc                                                                                             <  2:956228/140‑1 (MQ=255)                                       gATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGCCGCCAGAGTCATGACGCAGGCCCTGATACCGACTAGCGAACTCGACACCg                                                                                      >  1:956228/1‑140 (MQ=255)                                          tGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAA‑TCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGCCGCCAGAGTCATGACGCAGGCCCTGATACCGACTAGCGAACTCGACACCGAAAt                                                                                  <  1:547450/140‑1 (MQ=255)                                                 gaAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGCCGCCAGAGTCATGACGCAGGCCCTGATACCGACTAGCGAACTCGACACCGAAATCACGCa                                                                            <  1:235877/140‑1 (MQ=255)                                                  aaaaCCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGCCGCCAGAGTCATGACGCAGGCCCTGATACCGACTAGCGAACTCGACACCGAAATCACGCAg                                                                           >  1:865240/1‑140 (MQ=255)                                                    aaCCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGCCGCCAGAGTCATGACGCAGGCCCTGATACCGACTAGCGAACTCGACACCGAAATCACGCAGGa                                                                         >  2:568505/1‑140 (MQ=255)                                                     aCCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGCCGCCAGAGTCATGACGCAGGCCCTGATACCGACTAGCGAACTCGACACCGAAATCACGCAGGaa                                                                        <  2:314088/140‑1 (MQ=255)                                                      ccAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGCCGCCAGAGTCATGACGCAGGCCCTGATACCGACTAGCGAACTCGACACCGAAATCACGCAGGaaa                                                                       >  2:675638/1‑140 (MQ=255)                                                      ccAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGCCGCCAGAGTCATGACGCAGGCCCTGATACCGACTAGCGAACTCGACACCGAAATCACGCAGGaaa                                                                       <  2:283136/140‑1 (MQ=255)                                                      ccAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGCCACCAGAGTCATGACGCAGGCCCTGATACCGACTAGCGAACTCGACACCGAAATCACGCAGGaaa                                                                       >  2:260287/1‑140 (MQ=255)                                                          cGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGCCGCCAGAGTCATGACGCAGGCCCTGATACCGACTAGCGAACTCGACACCGAAATCACGCAGGAAAGCAt                                                                   <  2:228333/140‑1 (MQ=255)                                                          cGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGCCGCCAGAGTCATGACGCAGGCCCTGATACCGACTAGCGAACTCGACACCGAAATCACGCAGGAAAGCAt                                                                   >  1:228333/1‑140 (MQ=255)                                                           gTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGCCGCCAGAGTCATGACGCAGGCCCTGATACCGACTAGCGAACTCGACACCGAAATCACGCAGGAAAGCATc                                                                  <  2:440287/140‑1 (MQ=255)                                                               tACTATGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGCCGCCAGAGTCATGACGCAGGCCCTGATACCGACTAGCGAACTCGACACCGAAATCACGCAGGAAAGCATCCATa                                                              <  2:407286/140‑1 (MQ=255)                                                                        gATCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGCCGCCAGAGTCATGACGCAGGCCCTGATACCGACTAGCGAACTCGACACCGAAATCACGCAGGAAAGCATCCATAGTGATGCAg                                                     >  1:80403/1‑140 (MQ=255)                                                                        gATCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGCCGCCAGAGTCATGACGCAGGCCCTGATACCGACTAGCGAACTCGACACCGAAATCACGCAGGAAAGCATCCATAGTGATGCAg                                                     <  1:675638/140‑1 (MQ=255)                                                                         aTCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGCCGCCAGAGTCATGACGCAGGCCCTGATACCGACTAGCGAACTCGACACCGAAATCACGCAGGAAAGCATCCATAGTGATGCAGt                                                    >  1:152731/1‑140 (MQ=255)                                                                          tCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGCCGCCAGAGTCATGACGCAGGCCCTGATACCGACTAGCGAACTCGACACCGAAATCACGCAGGAAAGCATCCATAGTGATGCAGTc                                                   <  1:892469/140‑1 (MQ=255)                                                                              ttCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGCCGCCAGAGTCATGACGCAGGCCCTGATACCGACTAGCGAACTCGACACCGAAATCACGCAGGAAAGCATCCATAGTGATGCAGTCGGtt                                               >  2:145/1‑140 (MQ=255)                                                                                cccAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGCCGCCAGAGTCATGACGCAGGCCCTGATACCGACTAGCGAACTCGACACCGAAATCACGCAGGAAGGCATCCATAGTGATGCAGTCGGTTaa                                             >  1:673872/1‑140 (MQ=255)                                                                                cccAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGCCGCCAGAGTCATGACGCAGGCCCTGATACCGACTAGCGAACTCGACACCGAAATCACGCAGGAAAGCATCCATAGTGATg                                                        >  1:816396/1‑129 (MQ=255)                                                                                cccAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGCCGCCAGAGTCATGACGCAGGCCCTGATACCGACTAGCGAACTCGACACCGAAATCACGCAGGAAAGCATCCATAGTGATg                                                        <  2:816396/129‑1 (MQ=255)                                                                                cccAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGCCGCCAGAGTCATGACGCAGGCCCTGATACCGACTAGCGAACTCGACACCGAAATCACGCAGGAAAGCATCCATAGTGATGCAGTCGGTTaa                                             >  2:102178/1‑140 (MQ=255)                                                                                  cAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGCCGCCAGAGTCATGACGCAGGCCCTGATACCGACTAGCGAACTCGACACCGAAATCACGCAGGAAAGCATCCATAGTGATGCAGTCGGTTaatg                                           <  1:522347/140‑1 (MQ=255)                                                                                   aGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGCCGCCAGAGTCATGACGCAGGCCCTGATACCGACTAGCGAACTCGACACCGAAATCACGCAGGAAAGCATCCATAGTGATGCAGTCGGTTaatgc                                          >  2:80022/1‑140 (MQ=255)                                                                                   aGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGCCGCCAGAGTCATGACGCAGGCCCTGATACCGACTAGCGAACTCGACACCGAAATCACGCAGGAAAGCATCCATAGTGATGCAGTCGGTTaatgc                                          >  1:137759/1‑140 (MQ=255)                                                                                   aGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGCCGCCAGAGTCATGACGCAGGCCCTGATACCGACTAGCGAACTCGACACCGAAATCACGCAGGAAAGCATCCATAGTGATGCAGTCGGTTaatgc                                          >  1:562056/1‑140 (MQ=255)                                                                                            tAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGCCGCCAGAGTCATGACGCAGGCCCTGATACCGACTAGCGAACTCGACACCGAAATCACGCAGGAAAGCATCCATAGTGATGCAGTCGGTTAATGCAATGCCAAg                                 <  2:152731/140‑1 (MQ=255)                                                                                            tAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGCCGCCAGAGTCATGACGCAGGCCCTGATACCGACTAGCGAACTCGACACCGAAATCACGCAGGAAAGCATCCATAGTGATGCAGTCGGTTAATGCAATGCCAAg                                 <  2:80403/140‑1 (MQ=255)                                                                                               tgtgCAATGGCTTTTTCACCCCATTCAACCGGGCCGCCAGAGTCATGACGCAGGCCCTGATACCGACTAGCGAACTCGACACCGAAATCACGCAGGAAAGCATCCATAGTGATGCAGTCGGTTAATGCAATGCCAAGt                                <  2:357686/138‑1 (MQ=255)                                                                                               tgtgCAATGGCTTTTTCACCCCATTCAACCGGGCCGCCAGAGTCATGACGCAGGCCCTGATACCGACTAGCGAACTCGACACCGAAATCACGCAGGAAAGCATCCATAGTGATGCAGTCGGTTAATGCAATGCCAAGt                                >  1:357686/1‑138 (MQ=255)                                                                                               tgtgCAATGGCTTTTTCACCCCATTCAACCGGGCCGCCAGAGTCATGACGCAGGCCCTGATACCGACTAGCGAACTCGACACCGAAATCACGCAGGAAAGCATCCATAGTGATGCAGTCGGTTAATGCAATGCCAAGTTg                              >  2:126944/1‑140 (MQ=255)                                                                                                 tgCAATGGCTTTTTCACCCCATTCAACCGGGCCGCCAGAGTCATGACGCAGGCCCTGATACCGACTAGCGAACTCGACACCGAAATCACGCAGGAAAGCATCCATAGTGATGCAGTCTGTTAATGCAATGCCAAGTTGGt                            >  1:681269/1‑140 (MQ=255)                                                                                                   cAATGGCTTTTTCACCCCATTCAACCGGGCCGCCAGAGTCATGACGCAGGCCCTGATACCGACTAGCGAACTCGACACCGATATCACGCAGGAAAGCATCCATAGTGATGCAGTCGGTTAATGCAATGCCAAGTTGGTCg                          >  2:842909/1‑140 (MQ=255)                                                                                                      tGGCTTTTTCACCCCATTCAACCGGGCCGCCAGAGTCATGACGCAGGCCCTGATACCGACTAGCGAACTCGACACCGAAATCACGCAGGAAAGCATCCATAGTGATGCAGTCGGTTAATGCAATGCCAAGTTGGTCGGGa                       <  2:681269/140‑1 (MQ=255)                                                                                                        gCTTTTTCACCCCATTCAACCGGGCCGCCAGAGTCATGACGCAGGCCCTGATACCGACTAGCGAACTCGACACCGAAATCACGCAGGAAAGCATCCATAGTGATGCAGTCGGTTAATGCAATGCCAAGTTGGTCGGGATa                     <  1:145/140‑1 (MQ=255)                                                                                                              tCACCCCATTCAACCGGGCCGCCAGAGTCATGACGCAGGCCCTGATACCGACTAGCGAACTCGACACCGAAATCACGCAGGAAAGCATCCATAGTGATGCAGTCGGTTAATGCAATGCCAAGTTGGTCGGGATACTCTTc               <  2:865240/140‑1 (MQ=255)                                                                                                                  cccATTCAACCGGGCCGCCAGAGTCATGACGCAGGCCCTGATACCGACTAGCGAACTCGACACCGAAATCACGCAGGAAAGCATCCATAGTGATGCAGTCGGTTAATGCAATGCCAAGTTGGTCGGGATACTCTTccagc           <  1:126944/140‑1 (MQ=255)                                                                                                                          aCCGGGCCGCCAGAGTCATGACGCAGGCCCTGATACCGACTAGCGAACTCGACACCGAAATCACGCAGGAAAGCATCCATAGTGATGCAGTCGGTTAATGCAATGCCAAGTTGGTCGGGATa                     >  2:180386/1‑122 (MQ=255)                                                                                                                          aCCGGGCCGCCAGAGTCATGACGCAGGCCCTGATACCGACTAGCGAACTCGACACCGAAATCACGCAGGAAAGCATCCATAGTGATGCAGTCGGTTAATGCAATGCCAAGTTGGTCGGGATa                     <  1:180386/122‑1 (MQ=255)                                                                                                                          aCCGGGCCGCCAGAGTCATGACGCAGGCCCTGATACCGACTAGCGAACTCGACACCGAAATCACGCAGGAAAGCATCCATAGTGATGCAGTCGGTTAATGCAATGCCAAGTTGGTCGGGATACTCTTCCAGCCAGgcagc   >  1:958301/1‑140 (MQ=255)                                                                                                                           ccGGGCCGCCAGAGTCATGACGCAGGCCCTGATACCGACTAGCGAACTCGACACCGAAATCACGCAGGAAAGCATCCATAGTGATGCAGTCGGTTAATGCAATGCCAAGTTGGTCGGGATACTCTTCCAGCCAGgcagca  >  2:11658/1‑140 (MQ=255)                                                                                                                                |                                                                                                                                       GAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGTCGCCAGAGTCATGACGCAGGCCCTGATACCGACTAGCGAACTCGACACCGAAATCACGCAGGAAAGCATCCATAGTGATGCAGTCGGTTAATGCAATGCCAAGTTGGTCGGGATACTCTTCCAGCCAGGCAGCA  >  NC_000913/989373‑989634

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