Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,949,993	(T)5→6	coding (752/1332 nt)	argA →	N‑acetylglutamate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,949,988	1	.	T	96.5%	194.4 / ‑0.5	57	coding (747/1332 nt)	argA	N‑acetylglutamate synthase
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base T (32/23);  total (33/24)
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.69e-01

TTGGTTTTTGCTCTTCCCAGGGCGTCACTAATGACGACGGTGATATTGTCTCCGAACTTTTCCCTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGC‑TTTTTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAG  >  NC_000913/2949858‑2950122                                                                                                                                    |                                                                                                                                       ttGGTTTTTGCTCTTCCCAGGGCGTCACTAATGACGACGGTGATATTGTCTCCGAACTTTTCCCTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGC‑ttttttgcg                                                                                                                               <  1:850073/140‑5 (MQ=255)               tCCCAGGGCGTCACTAATGACGACGGTGATATTGTCTCCGAACTTTTCCCTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTTTTGCGTGgcg                                                                                                                           <  2:590537/131‑1 (MQ=255)               tCCCAGGGCGTCACTAATGACGACGGTGATATTGTCTCCGAACTTTTCCCTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTTTTGCGTGgcg                                                                                                                           >  1:590537/1‑131 (MQ=255)               tCCCAGGGCGTCACTAATGACGACGGTGATATTGTCTCCGAACTTTTCCCTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTTTTGCGTGGCGCAGTGAAAg                                                                                                                  <  2:749127/140‑1 (MQ=255)                 ccAGGGCGTCACTAATGACGACGGTGATATTGTCTCCGAACTTTTCCCTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTTTTGCGTGGCGCAGTGAAAGcc                                                                                                                <  2:874616/140‑1 (MQ=255)                   aGGGCGTCACTAATGACGACGGTGATATTGTCTCCGAACTTTTCCCTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTTTTGCGTGGCGCAGTGAAAGCCTg                                                                                                              >  2:902509/1‑140 (MQ=255)                     ggCGTCACTAATGACGACGGTGATATTGTCTCCGAACTTTTCCCTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTTTTGCGTGGCGCAGTGAAAGCCTGcc                                                                                                            >  1:11430/1‑140 (MQ=255)                         tCACTAATGACGACGGTGATATTGTCTCCGAACTTTTCCCTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTTTTGCGTGGCGCAGTGAAAGcc                                                                                                                >  2:596558/1‑132 (MQ=255)                         tCACTAATGACGACGGTGATATTGTCTCCGAACTTTTCCCTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTTTTGCGTGGCGCAGTGAAAGcc                                                                                                                <  1:596558/132‑1 (MQ=255)                         tCACTAATGACGACGGTGATATTGTCTCCGAACTTTTCCCTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTTTTGCGTGGCGCAGTGAAAGCCTGCCGCa                                                                                                         >  1:554474/1‑139 (MQ=255)                         tCACTAATGACGACGGTGATATTGTCTCCGAACTTTTCCCTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTTTTGCGTGGCGCAGTGAAAGCCTGCCGCa                                                                                                         <  2:554474/139‑1 (MQ=255)                          cACTAATGACGACGGTGATATTGTCTCCGAACTTTTCCCTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTTTTGCGTGGCGCAGTGAAAGCCTGCCGCAgc                                                                                                       >  1:50715/1‑140 (MQ=255)                          cACTAATGACGACGGTGATATTGTCTCCGAACTTTTCCCTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTTTTGCGTGGCGCAGTGAAAGCCTGCCGCAgc                                                                                                       >  1:487008/1‑140 (MQ=255)                          cACTAATGACGACGGTGATATTGTCTCCGAACTTTTCCCTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTTTTGCGTGGCGCAGTGAAAGCCTGCCGCAgc                                                                                                       >  1:402063/1‑140 (MQ=255)                          cACTAATGACGACGGTGATATTGTCTCCGAACTTTTCCCTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTTTTGCGTGGCGCAGTGAAAGCCTGCCGCAgc                                                                                                       >  2:994951/1‑140 (MQ=255)                          cACTAATGACGACGGTGATATTGTCTCCGAACTTTTCCCTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTTTTGCGTGGCGCAGTGAAAGCCTGCCGCAgc                                                                                                       >  1:7583/1‑140 (MQ=255)                          cACTAATGACGACGGTGATATTGTCTCCGAACTTTTCCCTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTTTTGCGCGGCGCAGTGAAAGCCTGCCGCAgc                                                                                                       >  2:362342/1‑140 (MQ=255)                               aTGACGACGGTGATATTGTCTCCGAACTTTTCCCTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTTTTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGgcgt                                                                                                  >  2:31412/1‑140 (MQ=255)                               aTGACGACGGTGATATTGTCTCCGAACTTTTCCCTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTTTTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGgcgt                                                                                                  >  2:696157/1‑140 (MQ=255)                               aTGACGACGGTGATATTGTCTCCGAACTTTTCCCTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTTTTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGgcgt                                                                                                  >  2:1040215/1‑140 (MQ=255)                                     acgGTGATATTGTCTCCGAACTTTTCCCTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTTTTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCg                                                                                            <  1:1040215/140‑1 (MQ=255)                                     acgGTGATATTGTCTCCGAACTTTTCCCTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTTTTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCg                                                                                            <  2:11430/140‑1 (MQ=255)                                         tGATATTGTCTCCGAACTTTTCCCTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTTTTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGt                                                                                        >  2:58429/1‑140 (MQ=255)                                                 tctcCGAACTTTTCCCTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTTTTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAAt                                                                                >  2:937465/1‑140 (MQ=255)                                                 tctcCGAACTTTTCCCTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTTTTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAAt                                                                                >  2:622483/1‑140 (MQ=255)                                                  ctcCGAACTTTTCCCTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTTTTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATc                                                                               >  2:660431/1‑140 (MQ=255)                                                   tcCGAACTTTTCCCTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTTTTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCt                                                                                          <  2:520490/128‑1 (MQ=255)                                                   tcCGAACTTTTCCCTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTTTTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCt                                                                                          >  1:520490/1‑128 (MQ=255)                                                   tcCGAACTTTTCCCTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTTTTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCa                                                                              >  1:214846/1‑140 (MQ=255)                                                   tcCGAACTTTTCCCTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTTTTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCa                                                                              >  2:753174/1‑140 (MQ=255)                                                   tcCGAACTTTTCCCTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTTTTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCa                                                                              >  2:383914/1‑140 (MQ=255)                                                   tcCGAACTTTTCCCTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTTTTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCa                                                                              >  2:392882/1‑140 (MQ=255)                                                    ccGAACTTTTCCCTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTTTTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAg                                                                             <  2:402063/140‑1 (MQ=255)                                                           tttCCCTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTTTTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAg                                                                      >  1:993588/1‑140 (MQ=255)                                                                  aaCGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTTTTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATg                                                               <  1:397551/140‑1 (MQ=255)                                                                      aagcgcaagcgcGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTTTTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGgcgc                                                           <  2:821188/140‑1 (MQ=255)                                                                      aagcgcaagcgcGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTTTTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGgcgc                                                           <  2:821187/140‑1 (MQ=255)                                                                      aagcgcaagcgcGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTTTTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGgcgc                                                           >  1:722232/1‑140 (MQ=255)                                                                       agcgcaagcgcGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTTTTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGgcgcg                                                          <  1:622483/140‑1 (MQ=255)                                                                       agcgcaagcgcGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTTTTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGgcgcg                                                          <  1:937465/140‑1 (MQ=255)                                                                        gcgcaagcgcGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTTTTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGgcgcgc                                                         <  2:301346/140‑1 (MQ=255)                                                                        gcgcaagcgcGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTTTTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGgcgcgc                                                         <  1:902509/140‑1 (MQ=255)                                                                           caagcgcGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTTTTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGt                                                      <  1:58429/140‑1 (MQ=255)                                                                                 cgGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTTTTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAag                                                >  1:596267/1‑140 (MQ=255)                                                                                               ggAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTTTTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCAcgc                                  <  1:765676/140‑1 (MQ=255)                                                                                                 aaGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTTTTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGa                                <  1:31412/140‑1 (MQ=255)                                                                                                 aaGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTTTTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGa                                <  1:660431/140‑1 (MQ=255)                                                                                                     gaAAGGCGATTACAACTCCGGTACGGTGCGCTTTTTTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGt                            >  2:284441/1‑140 (MQ=255)                                                                                                       aaGGCGATTACAACTCCGGTACGGTGCGCTTTTTTGCGTGGCCCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTAt                          <  2:722232/140‑1 (MQ=255)                                                                                                           cGATTACAACTCCGGTACGGTGCGCTTTTTTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGt                      >  2:998046/1‑140 (MQ=255)                                                                                                                aCAACTCCGGTACGGTGCGCTTTTTTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCa                 <  1:392882/140‑1 (MQ=255)                                                                                                                aCAACTCCGGTACGGTGCGCTTTTTTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCa                 <  2:487008/140‑1 (MQ=255)                                                                                                                    cTCCGGTACGGTGCGC‑TTTTTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTa                           >  1:611777/1‑125 (MQ=255)                                                                                                                    cTCCGGTACGGTGCGC‑TTTTTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTa                           <  2:611777/125‑1 (MQ=255)                                                                                                                         gTACGGTGCGCTTTTTTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGAt        >  2:637485/1‑140 (MQ=255)                                                                                                                               tGCGCTTTTTTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAg  >  1:1005239/1‑140 (MQ=255)                                                                                                                               tGCGCTTTTTTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAg  >  2:456102/1‑140 (MQ=255)                                                                                                                               tGCGCTTTTTTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAg  >  2:493528/1‑140 (MQ=255)                                                                                                                                    |                                                                                                                                       TTGGTTTTTGCTCTTCCCAGGGCGTCACTAATGACGACGGTGATATTGTCTCCGAACTTTTCCCTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGC‑TTTTTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAG  >  NC_000913/2949858‑2950122

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