Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	2,057,664	T→C	G95G (GGA→GGG)	cobT ←	nicotinate‑nucleotide‑‑dimethylbenzimidazole phosphoribosyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	2,057,664	0	T	C	100.0%	137.9 / NA	39	G95G (GGA→GGG)	cobT	nicotinate‑nucleotide‑‑dimethylbenzimidazole phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (21/18);  total (21/18)

GGGATAGGCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAG‑‑CACACACACGCCGGTTGTTCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATCGGC  >  NZ_CP009273/2057578‑2057747                                                                                         |                                                                                    gggATAGGCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAG‑‑CACACACACGCCGGTTGTCCCa                                                                                  >  2:569363/1‑90 (MQ=255)      aGGCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAG‑‑CACACACACGCCGGTTGTCCCACGGGt                                                                             >  1:566651/1‑90 (MQ=255)       ggCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAG‑‑CACACACACGCCGGTTGTCCCACGGGTc                                                                            >  1:417686/1‑90 (MQ=255)       ggCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAG‑‑CACACACACGCCGGTTGTCCCACGGGTc                                                                            >  2:344790/1‑90 (MQ=255)         cTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAG‑‑CACACACACGCCGGTTGTCCCACGGGTCat                                                                          <  2:832848/90‑1 (MQ=255)         cTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAG‑‑CACACACACGCCGGTTGTCCCACGGGTCat                                                                          <  1:54828/90‑1 (MQ=255)           cagcagTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAG‑‑CACACACACGCCGGTTGTCCCACGGGTCatat                                                                        >  1:129499/1‑90 (MQ=255)           cagcagTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAG‑‑CACACACACGCCGGTTGTCCCACGGGTCatat                                                                        >  1:168981/1‑90 (MQ=255)             gcagTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAG‑‑CACACACACGCCGGTTGTCCCACGGGTCATAttt                                                                      >  1:755708/1‑90 (MQ=255)                 tATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAG‑‑CACACACACGCCGGGTGTCCCACGGGTCATATTTTCAg                                                                  >  2:214553/1‑90 (MQ=255)                       tACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAG‑‑CACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTa                                                            >  2:162628/1‑90 (MQ=255)                       tACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAG‑‑CACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTa                                                            >  2:540400/1‑90 (MQ=255)                       tACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAG‑‑CACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTa                                                            >  1:562532/1‑90 (MQ=255)                          cAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAG‑‑CACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATgg                                                         >  1:437908/1‑90 (MQ=255)                          cAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAG‑‑CACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATgg                                                         >  2:593146/1‑90 (MQ=255)                          cAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAG‑‑CACACACACGCCGGTTGTCCCACGGGGCATATTTTCAGCCTGTATgg                                                         >  2:886171/1‑90 (MQ=255)                            aCATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATgg                                                         >  1:285062/1‑90 (MQ=255)                            aCATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAG‑‑CACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCt                                                       >  2:810317/1‑90 (MQ=255)                            aCATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAG‑‑CACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCt                                                       >  2:45623/1‑90 (MQ=255)                                cAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAG‑‑CACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCa                                                   >  1:22549/1‑90 (MQ=255)                                  aTTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAG‑‑CACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCAct                                                 <  1:683449/90‑1 (MQ=255)                                        tGGACGTTAGCGCCCGCTTGTTCTGCCAG‑‑CACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCtttt                                           <  1:886171/90‑1 (MQ=255)                                         ggACGTTAGCGCCCGCTTGTTCTGCCAG‑‑CACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCAct                                                 >  1:11666/1‑83 (MQ=255)                                         ggACGTTAGCGCCCGCTTGTTCTGCCAG‑‑CACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCAct                                                 <  2:11666/83‑1 (MQ=255)                                              ttAGCGCCCGCTTGTTCTGCCAG‑‑CACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGaa                                     <  1:569363/90‑1 (MQ=255)                                                  cgcCCGCTTGTTCTGCCAG‑‑CACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAg                                 <  2:891974/90‑1 (MQ=255)                                                     ccGCTTGTTCTGCCAG‑‑CACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGa                              <  1:569723/90‑1 (MQ=255)                                                       gCTTGTTCTGCCAG‑‑CACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGAcc                            >  2:753430/1‑90 (MQ=255)                                                         ttGTTCTGCCAGCACACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGAcc                            <  2:285062/90‑1 (MQ=255)                                                          tGTTCTGCCAG‑‑CACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCt                         <  1:665596/90‑1 (MQ=255)                                                           gTTCTGCCAG‑‑CACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCtt                        <  1:110562/90‑1 (MQ=255)                                                                    g‑‑CACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGa               <  2:553082/90‑1 (MQ=255)                                                                         cacacacGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGc            <  1:450128/90‑1 (MQ=255)                                                                                gCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATc     <  1:881547/90‑1 (MQ=255)                                                                                gCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATc     <  1:358970/90‑1 (MQ=255)                                                                                   ggTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCtt                        <  1:682157/68‑1 (MQ=255)                                                                                   ggTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCtt                        >  2:682157/1‑68 (MQ=255)                                                                                   ggTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATCGGc  <  2:174463/90‑1 (MQ=255)                                                                                   ggTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATCGGc  <  1:810317/90‑1 (MQ=255)                                                                                         |                                                                                    GGGATAGGCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAG‑‑CACACACACGCCGGTTGTTCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATCGGC  >  NZ_CP009273/2057578‑2057747

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