Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,204,645	G→A	*442* (TAG→TAA)  N429N (AAC→AAT)	zraR → purD ←	fused DNA‑binding response regulator in two‑component regulatory system with ZraS: response regulator/sigma54 interaction protein phosphoribosylglycinamide synthetase phosphoribosylamine‑glycine ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,204,645	0	G	A	98.3%	183.3 / NA	58	*442* (TAG→TAA)  N429N (AAC→AAT)	zraR purD	fused DNA‑binding response regulator in two‑component regulatory system with ZraS: response regulator/sigma54 interaction protein phosphoribosylglycinamide synthetase phosphoribosylamine‑glycine ligase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (28/29);  minor base C (1/0);  total (29/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 = 9.99e-01

GGATATTCAGCCGTTGGTGGAAGTGGAAAAAGAGGTGATTCTGGCGGCGCTGGAGAAAACGGGCGGCAACAAAACCGAAGCCGCCCGTCAGTTAGGGATCACGCGCAAAACGCTATTGGCAAAACTGTCGCGTTAGTTCTGCTCGCGTTCGATAGCGCGCCAGCCGATATCTTTCCGGCAGAAGCAGTCGTCCCAGTGGATATCGGTCATTAAGGCATAGGCGCGTTTCTGCGCTTCTGCCACGGTATGACCCAGCGCGGTGACGCAC  >  NC_000913/4204510‑4204777                                                                                                                                        |                                                                                                                                     ggATATTCAGCCGTTGGTGGAAGTGGAAAAAGAGGTGATTCTGGCGGCGCTGGAGAAAACGGGCGGCAACAAAACCGAAGCCGCCCGTCAGTTAGGGATCACGCGCAAAACGCTATTGGCAAAACTGTCGCGTTAATTc                                                                                                                                   >  1:189219/1‑139 (MQ=255)      ttCAGCCGTTGGTGGAAGTGGAAAAAGAGGTGATTCTGGCGGCGCTGGAGAAAACGGGCGGCAACAAAACCGAAGCCGCCCGTCAGTTAGGGATCACGCGCAAAACGCTATTGGCAAAACTGTCGCGTTAATTCTGCTc                                                                                                                              <  2:1092935/139‑1 (MQ=255)        cAGCCGTTGGTGGAAGTGGAAAAAGAGGTGATTCTGGCGGCGCTGGAGAAAACGGGCGGCAACAAAACCGAAGCCGCCCGTCAGTTAGGGATCACGCGCAAAACGCTATTGGCAAAACTGTCGCGTTAATTCTGCTcgc                                                                                                                            >  2:513339/1‑139 (MQ=255)            cGTTGGTGGAAGTGGAAAAAGAGGTGATTCTGGCGGCGCTGGAGAAAACGGGCGGCAACAAAACCGAAGCCGCCCGTCAGTTAGGGATCACGCGCAAAACGCTATTGGCAAAACTGTCGCGTTAATTCTGCTCGCGTTc                                                                                                                        >  2:446887/1‑139 (MQ=255)              ttGGTGGAAGTGGAAAAAGAGGTGATTCTGGCGGCGCTGGAGAAAACGGGCGGCAACAAAACCGAAGCCGCCCGTCAGTTAGGGATCACGCGCAAAACGCTATTGGCAAAACTGTCGCGTTACTTCTGCTCGCGTTCGa                                                                                                                      >  1:888548/1‑139 (MQ=255)              ttGGTGGAAGTGGAAAAAGAGGTGATTCTGGCGGCGCTGGAGAAAACGGGCGGCAACAAAACCGAAGCCGCCCGTCAGTTAGGGATCACGCGCAAAACGCTATTGGCAAAACTGTCGCGTTAATTCTGCTCGCGTTCGa                                                                                                                      >  1:200891/1‑139 (MQ=255)              ttGGTGGAAGTGGAAAAAGAGGTGATTCTGGCGGCGCTGGAGAAAACGGGCGGCAACAAAACCGAAGCCGCCCGTCAGTTAGGGATCACGCGCAAAACGCTATTGGCAAAACTGTCGCGTTAATTCTGCTCGCGTTCGa                                                                                                                      >  1:371226/1‑139 (MQ=255)                              aaGAGGTGATTCTGGCGGCGCTGGAGAAAACGGGCGGCAACAAAACCGAAGCCGCCCGTCAGTTAGGGATCACGCGCAAAACGCTATTGGCAAAACTGTCGCGTTAATTCTGCTCGCGTTCGATAGCGCGCCAGCCGat                                                                                                      <  1:513339/139‑1 (MQ=255)                               agagGTGATTCTGGCGGCGCTGGAGAAAACGGGCGGCAACAAAACCGAAGCCGCCCGTCAGTTAGGGATCACGCGCAAAACGCTATTGGCAAAACTGTCGCGTTAATTCTGCTcgcg                                                                                                                           <  1:802334/117‑1 (MQ=255)                               agagGTGATTCTGGCGGCGCTGGAGAAAACGGGCGGCAACAAAACCGAAGCCGCCCGTCAGTTAGGGATCACGCGCAAAACGCTATTGGCAAAACTGTCGCGTTAATTCTGCTcgcg                                                                                                                           >  2:802334/1‑117 (MQ=255)                               agagGTGATTCTGGCGGCGCTGGAGAAAACGGGCGGCAACAAAACCGAAGCCGCCCGTCAGTTAGGGATCACGCGCAAAACGCTATTGGCAAAACTGTCGCGTTAATTCTGCTCGCGTTCGATAGCGCGCCAGCCGata                                                                                                     >  2:563018/1‑139 (MQ=255)                                   gTGATTCTGGCGGCGCTGGAGAAAACGGGCGGCAACAAAACCGAAGCCGCCCGTCAGTTAGGGATCACGCGCAAAACGCTATTGGCAAAACTGTCGCGTTAATTCTGCTCGCGTTCGATAgcg                                                                                                                 <  1:1046132/123‑1 (MQ=255)                                   gTGATTCTGGCGGCGCTGGAGAAAACGGGCGGCAACAAAACCGAAGCCGCCCGTCAGTTAGGGATCACGCGCAAAACGCTATTGGCAAAACTGTCGCGTTAATTCTGCTCGCGTTCGATAgcg                                                                                                                 >  2:1046132/1‑123 (MQ=255)                                           ggcggcGCTGGAGAAAACGGGCGGCAACAAAACCGAAGCCGCCCGTCAGTTAGGGATCACGCGCAAAACGCTATTGGCAAAACTGTCGCGTTAATTCTGCTCGCGTTCGATAGCGCGCCAGCCGATATCTTTCCGGCAg                                                                                         >  1:872870/1‑139 (MQ=255)                                                 gcTGGAGAAAACGGGCGGCAACAAAACCGAAGCCGCCCGTCAGTTAGGGATCACGCGCAAAACGCTATTGGCAAAACTGTCGCGTTAATTCTGCTCGCGTTCGATAGCGCGCCAGCCGATATCTTTCCGGCAGAAGCAg                                                                                   >  1:5962/1‑139 (MQ=255)                                                  cTGGAGAAAACGGGCGGCAACAAAACCGAAGCCGCCCGTCAGTTAGGGATCACGCGCAAAACGCTATTGGCAAAACTGTCGCGTTAATTCTGCTCGCGTTCGATAGCGCGCCAGCCGATATCTTTCCGGCAGAAGCAgt                                                                                  <  1:481166/139‑1 (MQ=255)                                                  cTGGAGAAAACGGGCGGCAACAAAACCGAAGCCGCCCGTCAGTTAGGGATCACGCGCAAAACGCTATTGGCAAAACTGTCGCGTTAATTCTGCTCGCGTTCGATAGCGCGCCAGCCGATATCTTTCCGGCAGAAGCAgt                                                                                  <  2:388755/139‑1 (MQ=255)                                                  cTGGAGAAAACGGGCGGCAACAAAACCGAAGCCGCCCGTCAGTTAGGGATCACGCGCAAAACGCTATTGGCAAAACTGTCGCGTTAATTCTGCTCGCGTTCGATAGCGCGCCAGCCGATATCTTTCCGGCAGAAGCAgt                                                                                  <  2:388756/139‑1 (MQ=255)                                                        aaaaCGGGCGGCAACAAAACCGAAGCCGCCCGTCAGTTAGGGATCACGCGCAAAACGCTATTGGCAAAACTGTCGCGTTAATTCTGCTCGCGTTCGATAGCGCGCCAGCCGATATCTTTCCGGCAGAAGCAGTCGTccc                                                                            <  2:189219/139‑1 (MQ=255)                                                        aaaaCGGGCGGCAACAAAACCGAAGCCGCCCGTCAGTTAGGGATCACGCGCAAAACGCTATTGGCAAAACTGTCGCGTTAATTCTGCTCGCGTTCGATAGCGCGCCAGCCGATATCTTTCCGGCAGAAGCAGTCGTccc                                                                            <  2:872870/139‑1 (MQ=255)                                                         aaaCGGGCGGCAACAAAACCGAAGCCGCCCGTCAGTTAGGGATCACGCGCAAAACGCTATTGGCAAAACTGTCGCGTTAATTCTGCTCGCGTTCGATAGCGCGCCAGCCGATATCTTTCCGGCAGAAGCAGTCGTCCCa                                                                           <  1:446887/139‑1 (MQ=255)                                                               gcggcAACAAAACCGAAGCCGCCCGTCAGTTAGGGATCACGCGCAAAACGCTATTGGCAAAACTGTCGCGTTAATTCTGCTCGCGTTCGATAGCGCGCCAGCCGATATCTTTCCGGCAGAAGCAGTCGTCCCAGTGGat                                                                     >  1:666730/1‑139 (MQ=255)                                                                cggcAACAAAACCGAAGCCGCCCGTCAGTTAGGGATCACGCGCAAAACGCTATTGGCAAAACTGTCGCGTTAATTCTGCTCGCGTTCGATAGCGCGCCAGCCGATATCTTTCCGGCAGAAGCAGTCGTCCCAGTGGata                                                                    <  2:666730/139‑1 (MQ=255)                                                                 ggcAACAAAACCGAAGCCGCCCGTCAGTTAGGGATCACGCGCAAAACGCTATTGGCAAAACTGTCGCGTTAATTCTGCTCGCGTTCGATAGCGCGCCAGCCGATATCTTTCCGGCAGAAGCAGTCGTCCCAGTGGatat                                                                   <  1:156161/139‑1 (MQ=255)                                                                 ggcAACAAAACCGAAGCCGCCCGTCAGTTAGGGATCACGCGCAAAACGCTATTGGCAAAACTGTCGCGTTAATTCTGCTCGCGTTCGATAGCGCGCCAGCCGATATCTTTCCGGCAGAAGCAGTCGTCCCAGTGGatat                                                                   >  2:156161/1‑139 (MQ=255)                                                                                   gccCGTCAGTTAGGGATCACGCGCAAAACGCTATTGGCAAAACTGTCGCGTTAATTCTGCTCGCGTTCGATAGCGCGCCAGCCGATATCTTTCCGGCAGAAGCAGTCGTCCCAGTGGATATCGGTCATTAAGGCAt                                                    >  1:63041/1‑136 (MQ=255)                                                                                   gccCGTCAGTTAGGGATCACGCGCAAAACGCTATTGGCAAAACTGTCGCGTTAATTCTGCTCGCGTTCGATAGCGCGCCAGCCGATATCTTTCCGGCAGAAGCAGTCGTCCCAGTGGATATCGGTCATTAAGGCAt                                                    <  2:63041/136‑1 (MQ=255)                                                                                    cccGTCAGTTAGGGATCACGCGCAAAACGCTATTGGCAAAACTGTCGCGTTAATTCTGCTCGCGTTCGATAGCGCGCCAGCCGATATCTTTCCGGCAGAAGCAGTCGTCCCAGTGGATATCGGTCATTAAGGCATAGgc                                                >  2:615186/1‑139 (MQ=255)                                                                                     ccGTCAGTTAGGGATCACGCGCAAAACGCTATTGGCAAAACTGTCGCGTTAATTCTGCTCGCGTTCGATAGCGCGCCAGCCGATATCTTTCCGGCAGAAGCAGTCGTCCCAGTGGATATCGGTCATTAAGGCATAGgcg                                               <  1:615186/139‑1 (MQ=255)                                                                                                 gATCACGCGCAAAACGCTATTGGCAAAACTGTCGCGTTAATTCTGCTCGCGTTCGATAGCGCGCCAGCCGATATCTTTCCGGCAGAAGCAGTCGTCCCAGTGGATATCTGTCATTAAGGCATAGGCGCGTTTCTGCGCt                                   >  2:354020/1‑139 (MQ=255)                                                                                                   tCACGCGCAAAACGCTATTGGCAAAACTGTCGCGTTAATTCTGCTCGCGTTCGATAGCGCGCCAGCCGATATCTTTCCGGCAGAAGCAGTCGTCCCAGTGGATATCGGTCATTAAGGCATAGGCGCGTTTCTGCGCTTc                                 >  2:297437/1‑139 (MQ=255)                                                                                                     aCGCGCAAAACGCTATTGGCAAAACTGTCGCGTTAATTCTGCTCGCGTTCGATAGCGCGCCAGCCGATATCTTTCCGGCAGAAGCa                                                                                    <  1:169116/86‑1 (MQ=255)                                                                                                     aCGCGCAAAACGCTATTGGCAAAACTGTCGCGTTAATTCTGCTCGCGTTCGATAGCGCGCCAGCCGATATCTTTCCGGCAGAAGCa                                                                                    >  2:169116/1‑86 (MQ=255)                                                                                                      cgcgCAAAACGCTATTGGCAAAACTGTCGCGTTAATTCTGCt                                                                                                                               >  1:1166067/1‑42 (MQ=255)                                                                                                      cgcgCAAAACGCTATTGGCAAAACTGTCGCGTTAATTCTGCt                                                                                                                               <  2:1166067/42‑1 (MQ=255)                                                                                                      cgcgCAAAACGCTATTGGCAAAACTGTCGCGTTAATTCTGCTCGCGTTCGATAGCGCGCCAGCCGATATCTTTCCGGCAGAAGCAGTCGTCCCAGTGGATATCGGTCATTa                                                          >  1:196632/1‑111 (MQ=255)                                                                                                      cgcgCAAAACGCTATTGGCAAAACTGTCGCGTTAATTCTGCTCGCGTTCGATAGCGCGCCAGCCGATATCTTTCCGGCAGAAGCAGTCGTCCCAGTGGATATCGGTCATTa                                                          <  2:196632/111‑1 (MQ=255)                                                                                                      cgcgCAAAACGCTATTGGCAAAACTGTCGCGTTAATTCTGCTCGCGTTCGATAGCGCGCCAGCCGATATCTTTCCGGCAGAAGCAGTCGTCCCAGTGGATATCGGTCATTAAGGCATAGGCGCGTTTCTGCGCTTCTGc                              >  1:1292793/1‑139 (MQ=255)                                                                                                        cgcAAAACGCTATTGGCAAAACTGTCGCGTTAATTCTGCTCGCGTTCGATAGCGCGCCAGCCGATATCTTTCCGGCAGAAGCAGTCGTCCCAGTGGATATCGGTCATTAAGGCATAGGCGCGTTTCTGCGCTTCTGCCa                            <  2:888548/139‑1 (MQ=255)                                                                                                        cgcAAAACGCTATTGGCAAAACTGTCGCGTTAATTCTGCTCGCGTTCGATAGCGCGCCAGCCGATATCTTTCCGGCAGAAGCAGTCGTCCCAGTGGATATCGGTCATTAAGGCATAGGCGCGTTTCTGCGCTTCTGCCa                            <  2:1292793/139‑1 (MQ=255)                                                                                                        cgcAAAACGCTATTGGCAAAACTGTCGCGTTAATTCTGCTCGCGTTCGATAGCGCGCCAGCCGATATCTTTCCGGCAGAAGCAGTCGTCCCAGTGGATATCGGTCATTAAGGCATAGGCGCGTTTCTGCGCTTCTGCCa                            <  1:563018/139‑1 (MQ=255)                                                                                                         gcAAAACGCTATTGGCAAAACTGTCGCGTTAATTCTGCTCGCGTTCGATAGCGCGCCAGCCGATATCTTTCCGGCAGAAGc                                                                                     >  1:609740/1‑81 (MQ=255)                                                                                                         gcAAAACGCTATTGGCAAAACTGTCGCGTTAATTCTGCTCGCGTTCGATAGCGCGCCAGCCGATATCTTTCCGGCAGAAGc                                                                                     <  2:609740/81‑1 (MQ=255)                                                                                                          cAAAACGCTATTGGCAAAACTGTCGCGTTAATTCTGCTCGCGTTCGATAGCGCGCCAGCCGat                                                                                                      >  2:379386/1‑63 (MQ=255)                                                                                                          cAAAACGCTATTGGCAAAACTGTCGCGTTAATTCTGCTCGCGTTCGATAGCGCGCCAGCCGat                                                                                                      <  1:379386/63‑1 (MQ=255)                                                                                                          cAAAACGCTATTGGCAAAACTGTCGCGTTAATTCTGCTCGCGTTCGATAGCGCGCCAGCCGATATCTTTCCGGCAGAAGCAGTCGTCCCAGTGGATATCGGTCATTAAGGCATAGGCGCGTTTCTGCGCtt                                  >  1:476495/1‑131 (MQ=255)                                                                                                          cAAAACGCTATTGGCAAAACTGTCGCGTTAATTCTGCTCGCGTTCGATAGCGCGCCAGCCGATATCTTTCCGGCAGAAGCAGTCGTCCCAGTGGATATCGGTCATTAAGGCATAGGCGCGTTTCTGCGCtt                                  <  2:476495/131‑1 (MQ=255)                                                                                                                   aTTGGCAAAACTGTCGCGTTAATTCTGCTCGCGTTCGATAGCGCGCCAGCCGATATCTTTCCGGCAGAAGCAGTCGTCCCAGTGGATATCGGTCATTAAGGCATAGGCGCGTTTCTGCGCTTCTGc                              >  2:122572/1‑126 (MQ=255)                                                                                                                   aTTGGCAAAACTGTCGCGTTAATTCTGCTCGCGTTCGATAGCGCGCCAGCCGATATCTTTCCGGCAGAAGCAGTCGTCCCAGTGGATATCGGTCATTAAGGCATAGGCGCGTTTCTGCGCTTCTGc                              <  1:122572/126‑1 (MQ=255)                                                                                                                       gCAAAACTGTCGCGTTAATTCTGCTCGCGTTCGATAGCGCGCCAGCCGATATCTTTCCGGCAGAAGCAGTCGTCCCAGTGGATATCGGTCATTAAGGCATAGGCGCGTTTCTGCGCTTCTGCCACGGTATGACCCAgcg             <  2:5962/139‑1 (MQ=255)                                                                                                                            aCTGTCGCGTTAATTCTGCTCGCGTTCGATAGCGCGCCAGCCGATATCTTTCCGGCAGAAGCAGTCGTCCCAGTGGATATCGGTCATTAAGGCATAGGCGCGTTTCTGCGCTTCTGCCACGGTATGACCCAGCGCGGTg        >  1:691102/1‑139 (MQ=255)                                                                                                                            aCTGTCGCGTTAATTCTGCTCGCGTTCGATAGCGCGCCAGCCGATATCTTTCCGGCAGAAGCAGTCGTCCCAGTGGATATCGGTCATTAAGGCATAGGCGCGTTTCTGCGCTTCTGCCACGGTATGACCCAGCGCGGTg        >  1:363505/1‑139 (MQ=255)                                                                                                                              tGTCGCGTTAATTCTGCTCGCGTTCGATAGCGCGCCAGCCGATATCTTTCCGGCAGAAGCAGTCGTCCCAGTGGATATCGGTCATTAAGGCATAGGCGCGTTTCTGCGCTTCTGc                              >  2:1024351/1‑115 (MQ=255)                                                                                                                              tGTCGCGTTAATTCTGCTCGCGTTCGATAGCGCGCCAGCCGATATCTTTCCGGCAGAAGCAGTCGTCCCAGTGGATATCGGTCATTAAGGCATAGGCGCGTTTCTGCGCTTCTGc                              <  1:1024351/115‑1 (MQ=255)                                                                                                                                tgggtttAATTCTGCTCGCTTTCAATAGCGCGCCAGCCGATATTTTCCGGGCAGAAGCAGTCGTCCCGGTGGATATCGGTCATTAAGGCATAGGCGCGTTTCTGCGCTTCTGCCACGGTATGACCCAGCGCGGTGACGc    <  2:691102/134‑1 (MQ=255)                                                                                                                                tCGCGTTAATTCTGCTCGCGTTCGATAGCGCGCCAGCCGATATCTTTCCGGCAGAAGCAg                                                                                   <  1:727211/60‑1 (MQ=255)                                                                                                                                tCGCGTTAATTCTGCTCGCGTTCGATAGCGCGCCAGCCGATATCTTTCCGGCAGAAGCAg                                                                                   >  2:727211/1‑60 (MQ=255)                                                                                                                                  gcgTTAATTCTGCTCGCGTTCGATAGCGCGCTAGCCGATATCTTTCCGGCAGAAGCAGTCGTCCCAGTGGATATCGGTCATTAAGGCATAGGCGCGTTTCTGCGCTTCTGCCACGGTATGACCCAGCGCGGTGACGcac  <  2:837613/139‑1 (MQ=255)                                                                                                                                  gcgTTAATTCTGCTCGCGTTCGATAGCGCGCCAGCCGATATCTTTCCGGCAGAAGCAGTCGTCCCAGTGGATATCTGTCATTAAGGCATAGGCGCGTTTCTGCGCTTCTGCCACGGTATGACCCAGCGCGGTGACGcac  <  1:354020/139‑1 (MQ=255)                                                                                                                                        |                                                                                                                                     GGATATTCAGCCGTTGGTGGAAGTGGAAAAAGAGGTGATTCTGGCGGCGCTGGAGAAAACGGGCGGCAACAAAACCGAAGCCGCCCGTCAGTTAGGGATCACGCGCAAAACGCTATTGGCAAAACTGTCGCGTTAGTTCTGCTCGCGTTCGATAGCGCGCCAGCCGATATCTTTCCGGCAGAAGCAGTCGTCCCAGTGGATATCGGTCATTAAGGCATAGGCGCGTTTCTGCGCTTCTGCCACGGTATGACCCAGCGCGGTGACGCAC  >  NC_000913/4204510‑4204777

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