Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	2,057,587	T→C	E121G (GAG→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,587	0	T	C	98.3%	214.4 / ‑3.9	61	E121G (GAG→GGG)	cobT	nicotinate‑nucleotide‑‑dimethylbenzimidazole phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (30/29);  total (32/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.42e-01

TTTTTTGCCAGCTCCTGCGTATAACATATGACGTCCAAAAGCAACTTTTCAGCCTGACGGCGACTCATTGCCGGAGCTGAAGCAATATTGCCGCTACCTCGTGCGACACGCATGTTGATAAGCCCGGGGATAGGCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTTCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGAC  >  NZ_CP009273/2057452‑2057720                                                                                                                                        |                                                                                                                                      ttttttGCCAGCTCCTGCGTATAACATATGACGTCCAAAAGCAACTTTTCAGCCTGACGGCGACTCATTGCCGGAGCTGAAGCAATATTGCCGCTACCTCGTGCGACACGCATGTTGATAAGCCCGGGGATAGGCCcag                                                                                                                                    >  2:1272754/1‑139 (MQ=255) ttttttGCCAGCTCCTGCGTATAACATATGACGTCCAAAAGCAACTTTTCAGCCTGACGGCGACTCATTGCCGGAGCTGAAGCAATATTGCCGCTACCTCGTGCGACACGCATGTTGATAAGCCCGGGGATAGGCCcag                                                                                                                                    <  2:905115/139‑1 (MQ=255)    tttGCCAGCTCCTGCGTATAACATATGACGTCCAAAAGCAACTTTTCAGCCTGACGGCGACTCATTGCCGGAGCTGAAGCAATATTGCCGCTACCTCGTGCGACACGCATGTTGATAAGCCCGGGGATAGGTCcagca                                                                                                                                  <  2:319429/138‑1 (MQ=255)    tttGCCAGCTCCTGCGTATAACATATGACGTCCAAAAGCAACTTTTCAGCCTGACGGCGACTCATTGCCGGAGCTGAAGCAATATTGCCGCTACCTCGTGCGACACGCATGTTGATAAGCCCGGGGATAGGCCcagca                                                                                                                                  <  1:351607/138‑1 (MQ=255)           gCTCCTGCGTATAACATATGACGTCCAAAAGCAACTTTTCAGCCTGACGGCGACTCATTGCCGGAGCTGAAGCAATATTGCCGCTACCTCGTGCGACACGCATGTTGATAAGCCCGGGGATAGGCCCAGCAGTATCaa                                                                                                                           >  2:404737/1‑138 (MQ=255)               cTGCGTATAACATATGACGTCCAAAAGCAACTTTTCAGCCTGACGGCGACTCATTGCCGGAGCTGAAGCAATATTGCCGCTACCTCGTGCGACACGCATGTTGATAAGCCCGGGGATAGGCCCAGCAGTATCAATAcc                                                                                                                       <  2:580331/138‑1 (MQ=255)               cTGCGTATAACATATGACGTCCAAAAGCAACTTTTCAGCCTGACGGCGACTCATTGCCGGAGCTGAAGCAATATTGCCGCTACCTCGTGCGACACGCATGTTGATAAGCCCGGGGATAGGCCCAGCAGTATCAATAcc                                                                                                                       >  2:847558/1‑138 (MQ=255)               cTGCGTATAACATATGACGTCCAAAAGCAACTTTTCAGCCTGACGGCGACTCATTGCCGGAGCTGAAGCAATATTGCCGCTACCTCGTGCGACACGCATGTTGATAAGCCCGGGGATAGGCCCAGCAGTATCAATAcc                                                                                                                       <  2:986126/138‑1 (MQ=255)                  cGTATAACATATGACGTCCAAAAGCAACTTTTCAGCCTGACGGCGACTCATTGCCGGAGCTGAAGCAATATTGCCGCTACCTCGTGCGACACGCATGTTGATAAGCCCGGGGATAGGCCCAGCAGTATCAATACCaa                                                                                                                     <  1:1272754/137‑1 (MQ=255)                       aaCATATGACGTCCAAAAGCAACTTTTCAGCCTGACGGCGACTCATTGCCGGAGCTGAAGCAATATTGCCGCTACCTCGTGCGACACGCATGTTGATAAGCCCGGGGATAGGCCcagc                                                                                                                                   <  2:562205/118‑1 (MQ=255)                       aaCATATGACGTCCAAAAGCAACTTTTCAGCCTGACGGCGACTCATTGCCGGAGCTGAAGCAATATTGCCGCTACCTCGTGCGACACGCATGTTGATAAGCCCGGGGATAGGCCcagc                                                                                                                                   >  1:562205/1‑118 (MQ=255)                       aaCATATGACGTCCAAAAGCAACTTTTCAGCCTGACGGCGACTCATTGCCGGAGCTGAAGCAATATTGCCGCTACCTCGTGCGACACGCATGTTGATAAGCCCGGGGATAGGCCCAGCAGTATCAATACCAACATCAAt                                                                                                              >  2:708107/1‑139 (MQ=255)                            atGACGTCCAAAAGCAACTTTTCAGCCTGACGGCGACTCATTGCCGGAGCTGAAGCAATATTGCCGCTACCTCGTGCGACACGCATGTTGATAAGCCCGGGGATAGGCCCAGCAGTATCAATACCAACATCAATTACGt                                                                                                         >  1:666886/1‑139 (MQ=255)                             tGACGTCCAAAAGCAACTTTTCAGCCTGACGGCGACTCATTGCCGGAGCTGAAGCAATATTGCCGCTACCTCGTGCGACACGCATGTTGATAAGCCCGGGGATAGGCCCAGCAGTATCAATACCAACATCAATTACGTg                                                                                                        <  2:736077/139‑1 (MQ=255)                                 gTCCAAAAGCAACTTTTCAGCCTGACGGCGACTCATTGCCGGAGCTGAAGCAATATTGCCGCTACCTCGTGCGACACGCATGTTGATAAGCCCGGGGATAGGCCCAGCAGTATCAATACCAACATCAATTACGTGGACg                                                                                                    >  2:881353/1‑139 (MQ=255)                                   ccAAAAGCAACTTTTCAGCCTGACGGCGACTCATTGCCGGAGCTGAAGCAATATTGCCGCTACCTCGTGCGACACGCATGTTGATAAGCCCGGGGATAGGCCCAGCAGTATCAATACCAACATCAATTACGTGGACGt                                                                                                   >  1:692855/1‑138 (MQ=255)                                    cAAAAGTAACTTTTCAGCTTGACGGCGACTCATTGCCGTAGCTGAAGCAATATTGCCGCTACCTCGTGCGACACGCATGTTGATAAGCCCGGGGATAGGCCCAGCAGTACCAATACCAACATGAATTACGTGGACGt                                                                                                   >  2:315251/1‑137 (MQ=255)                                    cAAAAGCAACTTTTCAGCCTGACGGCGACTCATTGCCGGAGCTGAAGCAATATTGCCGCTACCTCGTGCGACACGCATGTTGATAAGCCCGGGGATAGGCCCAGCAGTATCAATACCAACATCAATTACGTGGACGt                                                                                                   >  1:87714/1‑137 (MQ=255)                                    cAAAAGCAACTTTTCAGCCTGACGGCGACTCATTGCCGGAGCTGAAGCAATATTGCCGCTACCTCGTGCGACACGCATGTTGATAAGCCCGGGGATAGGCCCAGCAGTATCAATACCAACATCAATTACGTGGACGt                                                                                                   >  2:34276/1‑137 (MQ=255)                                                               aCTCATTGCCGGAGCTGAAGCAATATTGCCGCTACCTCGTGCGACACGCATGTTGATAAGCCCGGGGATAGGCCCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGcacaca                                                                       >  2:980723/1‑138 (MQ=255)                                                                    ttGCCGGAGCTGAAGCAATATTGCCGCTACCTCGTGCGACACGCATGTTGATAAGCCCGGGGATAGGCCCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGcc                                                                 <  2:480601/139‑1 (MQ=255)                                                                           aGCTGAAGCAATATTGCCGCTACCTCGTGCGACACGCATGTTGATAAGCCCGGGGATAGGCCCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGgttgtt                                                          >  2:291778/1‑139 (MQ=255)                                                                           aGCTGAAGCAATATTGCCGCTACCTCGTGCGACACGCATGTTGATAAGCCCGGGGAT‑GGCCCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTTc                                                         >  1:864560/1‑139 (MQ=255)                                                                             cTGAAGCAATATTGCCGCTACCTCGTGCGACACGCATGTTGATAAGCCCGGGGATAGGCCCAGCAGTATCAATACCAACATCaa                                                                                                               <  2:349421/84‑1 (MQ=255)                                                                             cTGAAGCAATATTGCCGCTACCTCGTGCGACACGCATGTTGATAAGCCCGGGGATAGGCCCAGCAGTATCAATACCAACATCaa                                                                                                               >  1:349421/1‑84 (MQ=255)                                                                             cTGAAGCAATATTGCCGCTACCTCGTGCGACACGCATGTTGATAAGCCCGGGGATAGGCCCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTTcc                                                        <  1:34276/139‑1 (MQ=255)                                                                             cTGAAGCAATATTGCCGCTACCTCGTGCGACACGCATGTTGATAAGCCCGGGGATAGGCCCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTTcc                                                        <  2:818226/139‑1 (MQ=255)                                                                             cTGAAGCAATATTGCCGCTACCTCGTGCGACACGCATGTTGATAAGCCCGGGGATAGGCCCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTTcc                                                        <  1:89140/139‑1 (MQ=255)                                                                               gAAGCAATATTGCCGCTACCTCGTGCGACACGCATGTTGATAAGCCCGGGGATAGGCCCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTTCCAc                                                      >  1:312340/1‑139 (MQ=255)                                                                                aaGCAATATTGCCGCTACCTCGTGCGACACGCATGTTGATAAGCCCGGGGAT‑GGCCCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTTCCAc                                                      <  2:864560/137‑1 (MQ=255)                                                                                aaGCAATATTGCCGCTACCTCGTGCGACACGCATGGTGATAAGCCCGGGGATAGGCCCAGCAG‑‑GCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTCCTGCCAGCACACACACGCCGGTTGTTCCACgg                                                    <  2:1172468/138‑1 (MQ=255)                                                                                  gCAATATTGCCGCTACCTCGTGCGACACGCATGTTGATAAGCCCGGGGATAGGCCCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTTCCACggg                                                   >  1:841429/1‑139 (MQ=255)                                                                                     atatTGCCGCTACCTCGTGCGACACGCATGTTGATAAGCCCGGGGATAGGCCCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTTCCACGGGTCa                                                <  2:1095908/139‑1 (MQ=255)                                                                                     atatTGCCGCTACCTCGTGCGACACGCATGTTGATAAGCCCGGGGATAGGCCCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTTCCACGGGTCa                                                >  2:217050/1‑139 (MQ=255)                                                                                     atatTGCCGCTACCTCGTGCGACACGCATGTTGATAAGCCCGGGGATAGGCCCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTTCCACGGGTCa                                                <  1:881353/139‑1 (MQ=255)                                                                                          gCCGCTACCTCGTGCGACACGCATGTTGATAAGCCCGGGGATAGGCCcag                                                                                                                                    >  2:1008648/1‑50 (MQ=255)                                                                                          gCCGCTACCTCGTGCGACACGCATGTTGATAAGCCCGGGGATAGGCCcag                                                                                                                                    <  1:1008648/50‑1 (MQ=255)                                                                                          gCCGCTACCTCGTGCGACACGCATGTTGATAAGCCCGGGGATAGGCCCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTTCCACGGGTCATAttt                                           >  2:668284/1‑139 (MQ=255)                                                                                           ccGCTACCTCGTGCGACACGCATGTTGATAAGCCCGGGGATAGGCCCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGgt                                                              <  2:963502/119‑1 (MQ=255)                                                                                           ccGCTACCTCGTGCGACACGCATGTTGATAAGCCCGGGGATAGGCCCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGgt                                                              >  1:963502/1‑119 (MQ=255)                                                                                           ccGCTACCTCGTGCGACACGCATGTTGATAAGCCCGGGGATAGGCCCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTTCCACGGGTCATAtttt                                          <  2:749132/139‑1 (MQ=255)                                                                                             gCTACCTCGTGCGACACGCATGTTGATAAGCCCGGGGATAGGCCCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTTCCACGGGTCATAtttt                                          >  1:176111/1‑137 (MQ=255)                                                                                             gCTACCTCGTGCGACACGCATGTTGATAAGCCCGGGGATAGGCCCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTTCCACGGGTCATATTTTCa                                        <  1:291778/139‑1 (MQ=255)                                                                                             gCTACCTCGTGCGACACGCATGTTGATAAGCCCGGGGATAGGCCCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTTCCACGGGTCATATTTTCa                                        <  2:666886/139‑1 (MQ=255)                                                                                              cTACCTCGTGCGACACGCATGTTGATAAGCCCGGGGATAGGCCCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTTCCACGGGTCATATTTTCAg                                       >  2:1209396/1‑139 (MQ=255)                                                                                              cTACCTCGTGCGACACGCATGTTGATAAGCCCGGGGATAGGCCCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTTCCACGGGTCATATTTTCAg                                       >  1:306059/1‑139 (MQ=255)                                                                                                 ccTCGTGCGACACGCATGTTGATAAGCCCGGGGATAGGCCCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTTCCACGGGTCATATTTTCAGc                                      <  1:217050/137‑1 (MQ=255)                                                                                                     gTGCGACACGCATGTTGATAAGCCCGGGGATAGGCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTTCCACGGGTCATATTTTCAGCCTGTAt                                >  1:295224/1‑139 (MQ=255)                                                                                                      tGCGACACGCATGTTGATAAGCCCGGGGATAGGCCCAGCAGTATCAATACCAACATCAATTACGTGTACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTTCCACGGGTCATATTTTCAGc                                      <  1:278073/132‑1 (MQ=255)                                                                                                      tGCGACACGCATGTTGATAAGCCCGGGGATAGGCCCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAg                                                                             <  2:997001/93‑1 (MQ=255)                                                                                                      tGCGACACGCATGTTGATAAGCCCGGGGATAGGCCCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAg                                                                             >  1:997001/1‑93 (MQ=255)                                                                                                      tGCGACACGCATGTTGATAAGCCCGGGGATAGGCCCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTTCCACGGGTCATATTTTCAGCCTGTATg                               >  2:1112487/1‑139 (MQ=255)                                                                                                      tGCGACACGCATGTTGATAAGCCCGGGGATAGGCCCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTTCCACGGGTCATATTTTCAGCCTGTATg                               <  1:980723/139‑1 (MQ=255)                                                                                                      tGCGACACGCATGTTGATAAGCCCGGGGATAGGCACAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTTCCACGGGTCATATTTTCAGc                                      >  2:278073/1‑132 (MQ=255)                                                                                                           cacGCATGTTGATAAGCCCGGGGATAGGCCCAGCAGGATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTTCCACGGGTCATATTTTCAGCCTGTATGGCTGt                          >  1:934110/1‑139 (MQ=255)                                                                                                             cGCATGTTGATAAGCCCGGGGATAGGCCCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTTCCACGGGTCATATTTTCAGCCTGTATGGCTGTCa                        >  2:367380/1‑139 (MQ=255)                                                                                                                          gCCCGGGGATAGGCCCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTTCCACGGGTCATAttt                                           >  2:462863/1‑107 (MQ=255)                                                                                                                          gCCCGGGGATAGGCCCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTTCCACGGGTCATAttt                                           <  1:462863/107‑1 (MQ=255)                                                                                                                                ggATAGGCCCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTTCCACGGGTCATATTTTCa                                        <  2:336839/104‑1 (MQ=255)                                                                                                                                ggATAGGCCCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTTCCACGGGTCATATTTTCa                                        >  1:336839/1‑104 (MQ=255)                                                                                                                                   tAGGCCCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTTCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGAc  <  1:1209396/139‑1 (MQ=255)                                                                                                                                        |                                                                                                                                      TTTTTTGCCAGCTCCTGCGTATAACATATGACGTCCAAAAGCAACTTTTCAGCCTGACGGCGACTCATTGCCGGAGCTGAAGCAATATTGCCGCTACCTCGTGCGACACGCATGTTGATAAGCCCGGGGATAGGCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTTCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGAC  >  NZ_CP009273/2057452‑2057720

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