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	97.6%	142.7 / ‑3.3	42	coding (747/1332 nt)	argA	N‑acetylglutamate synthase
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base T (21/20);  total (22/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.41e-01

TTGGTTTTTGCTCTTCCCAGGGCGTCACTAATGACGACGGTGATATTGTCTCCGAACTTTTCCCTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGC‑TTT‑TTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAG  >  NC_000913/2949858‑2950122                                                                                                                                    |                                                                                                                                        ttGGTTTTTGCTCTTCCCAGGGCGTCACTAATGACGACGGTGATATTGTCTCCGAACTTTTCCCTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGC‑ttt‑tttgcg                                                                                                                               <  2:249564/140‑5 (MQ=255)      ttttGCTCTTCCCAGGGCGTCACTAATGACGACGGTGATATTGTCTCCGAACTTTTCCCTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTT‑TTGCGTGgcg                                                                                                                           <  1:413610/140‑1 (MQ=255)          gCTCTTCCCAGGGCGTCACTAATGACGACGGTGATATTGTCTCCGAACTTTTCCCTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTT‑TTGCGTGGCGCAGt                                                                                                                       >  2:39943/1‑140 (MQ=255)               tCCCAGGGCGTCACTAATGACGACGGTGATATTGTCTCCGAACTTTTCCCTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTT‑TTGCGTGGCGCAGTGAAAg                                                                                                                  <  1:452701/140‑1 (MQ=255)                cccAGGGCGTCACTAATGACGACGGTGATATTGTCTCCGAACTTTTCCCTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTT‑TTGCGTGGCGCAGTGAAAGc                                                                                                                 >  1:468031/1‑140 (MQ=255)                         tCACTAATGACGACGGTGATATTGTCTCCGAACTTTTCCCTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTT‑TTGCGTGGCGCAGTGAAAGCCTGCCGCAg                                                                                                        >  1:874973/1‑140 (MQ=255)                          cACTAATGACGACGGTGATATTGTCTCCGAACTTTTCCCTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTT‑TTGCGTGGCGCAGTGAAAGcc                                                                                                                >  2:509651/1‑131 (MQ=255)                          cACTAATGACGACGGTGATATTGTCTCCGAACTTTTCCCTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTT‑TTGCGTGGCGCAGTGAAAGcc                                                                                                                <  1:509651/131‑1 (MQ=255)                          cACTAATGACGACGGTGATATTGTCTCCGAACTTTTCCCTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTT‑TTGCGTGGCGCAGTGAAAGCCTGCCGCAgc                                                                                                       >  2:63150/1‑140 (MQ=255)                          cACTAATGACGACGGTGATATTGTCTCCGAACTTTTCCCTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTT‑TTGCGTGGCGCAGTGAAAGCCTGCCGCAgc                                                                                                       >  2:843842/1‑140 (MQ=255)                          cACTAATGACGACGGTGATATTGTCTCCGAACTTTTCCCTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTT‑TTGCGTGGCGCAGTGAAAGCCTGCCGCAgc                                                                                                       >  2:185812/1‑140 (MQ=255)                          cACTAATGACGACGGTGATATTGTCTCCGAACTTTTCCCTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTT‑TTGCGTGGCGCAGTGAAAGCCTGCCGCAgc                                                                                                       >  1:486904/1‑140 (MQ=255)                          cACTAATGACGACGGTGATATTGTCTCCGAACTTTTCCCTAACCAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTT‑TTGCGTGGCGCAGTGAACGCCTGCCGCAgc                                                                                                       >  2:347629/1‑140 (MQ=255)                           aCTAATGACGACGGTGATATTGTCTCCGAACTTTTCCCTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTT‑TTGCGTGGCGCAGTGAAAGCCTGCCGCAgcg                                                                                                      <  2:731581/140‑1 (MQ=255)                           aCTAATGACGACGGTGATATTGTCTCCGAACTTTTCCCTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTT‑TTGCGTGGCGCAGTGAAAGCCTGCCGCAgcg                                                                                                      <  1:39943/140‑1 (MQ=255)                                    gacgGTGATATTGTCTCCGAACTTTTCCCTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTT‑TTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTc                                                                                             >  1:528294/1‑140 (MQ=255)                                         tGATATTGTCTCCGAACTTTTCCCTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTT‑TTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGt                                                                                        >  2:90435/1‑140 (MQ=255)                                         tGATATTGTCTCCGAACTTTTCCCTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTT‑TTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGt                                                                                        <  1:63150/140‑1 (MQ=255)                                              ttGTCTCCGAACTTTTCCCTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTT‑TTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCAttt                                                                                   <  1:185812/140‑1 (MQ=255)                                              ttGTCTCCGAACTTTTCCCTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTT‑TTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCAttt                                                                                   <  2:528294/140‑1 (MQ=255)                                                   tcCGAACTTTTCCCTAACGAAGCGCAAGCNCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTT‑TTGCGTGGCGCAGTNAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCa                                                                              <  2:486904/140‑1 (MQ=255)                                                    ccGAACTTTTCCCTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTT‑TTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAg                                                                             <  1:90435/140‑1 (MQ=255)                                                    ccGAACTTTTCCCTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTT‑TTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAg                                                                             <  2:874973/140‑1 (MQ=255)                                                    ccGAACTTTTCCCTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTT‑TTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAg                                                                             <  1:843842/140‑1 (MQ=255)                                                            ttCCCTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTT‑TTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAgg                                                                     >  2:38425/1‑140 (MQ=255)                                                                    cGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTT‑TTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGgc                                                             >  2:57925/1‑140 (MQ=255)                                                                      aagcgcaagcgcGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTT‑TTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGgcgc                                                           <  1:38425/140‑1 (MQ=255)                                                                        gcgcaagcgcGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTT‑TTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGgcgcgc                                                         <  2:929903/140‑1 (MQ=255)                                                                             agcgcGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTT‑TTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTg                                                    >  2:373430/1‑140 (MQ=255)                                                                                   ggTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTT‑TTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAagag                                              >  1:421166/1‑140 (MQ=255)                                                                                        aaGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTT‑TTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAgttgtt                                         >  1:789577/1‑140 (MQ=255)                                                                                              aGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTT‑TTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCAcg                                   >  2:101562/1‑140 (MQ=255)                                                                                               ggAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTT‑TTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCAcgc                                  >  2:293299/1‑140 (MQ=255)                                                                                                 aaGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTT‑TTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGa                                <  2:811491/140‑1 (MQ=255)                                                                                                      aaaGGCGATTACAACTCCGGTACGGTGCGC‑TTTCTTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTa                           >  2:796295/1‑140 (MQ=255)                                                                                                          gCGATTACAACTCCGGTACGGTGCGCTTTT‑TTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCgg                       <  2:421166/140‑1 (MQ=255)                                                                                                              ttACAACTCCGGTACGGTGCGCTTTT‑TTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACg                   <  2:789577/140‑1 (MQ=255)                                                                                                              ttACAACTCCGGTACGGTGCGCTTTT‑TTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACg                   <  1:293299/140‑1 (MQ=255)                                                                                                                            cGGTGCGCTTTT‑TTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGa     <  1:856405/140‑1 (MQ=255)                                                                                                                               tGCGCTTTT‑TTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATtgtg          <  2:749747/132‑1 (MQ=255)                                                                                                                               tGCGCTTTT‑TTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATtgtg          >  1:749747/1‑132 (MQ=255)                                                                                                                               tGCGCTTTT‑TTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAg  >  2:741820/1‑140 (MQ=255)                                                                                                                               tGCGCTTTT‑TTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAg  >  2:313929/1‑140 (MQ=255)                                                                                                                                    |                                                                                                                                        TTGGTTTTTGCTCTTCCCAGGGCGTCACTAATGACGACGGTGATATTGTCTCCGAACTTTTCCCTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGC‑TTT‑TTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAG  >  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: ‑