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

CTTCCCAGGGCGTCACTAATGACGACGGTGATATTGTCTCCGAACTTTTCCCTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGC‑T‑TT‑TTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGA  >  NC_000913/2949870‑2950115                                                                                                                        |                                                                                                                                  ctTCCCAGGGGGTCACTAATGACGACCGTGATATTGTCTCCTAACTTTTCCCTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTT‑TT‑TTGCGTGgcg                                                                                                                    <  2:273737/133‑1 (MQ=255) ctTCCCAGGGCGTCACTAATGACGACGGTGATATTGTCTCCGAACTTTTCCCTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTT‑TT‑TTGCGTGgcg                                                                                                                    <  2:627285/133‑1 (MQ=255)   tCCCAGGGCGTCACTAATGACGACGGTGATATTGTCTCCGAACTTTTCCCTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTT‑TT‑TTGCGTGGCGCa                                                                                                                  >  1:460402/1‑133 (MQ=255)            gTCACTAATGACGACGGTGATATTGTCTCCGAACTTTTCCCTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTT‑TT‑TTGCGTGGCGCAGTGAAAGcc                                                                                                         >  2:63378/1‑133 (MQ=255)            gTCACTAATGACGACCGTTATATTGTCTCCCAACTTTTCCCTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTT‑TT‑TTGCGTGGCGCAGTGAAAGcc                                                                                                         <  1:205589/133‑1 (MQ=255)              cACTAATGACGACGGTGATATTGTCTCCGAACTTTTCCCTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTACGCTT‑TT‑TTGCGTGGCGCAGTAAAAGCCTg                                                                                                       >  1:292187/1‑133 (MQ=255)               aCTAATGACGACGGTGATATTGTCTCCGAACTTTTCCCTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTT‑TT‑TTGCGTGGCGCAGTGAAAGCCTGc                                                                                                      >  2:481494/1‑133 (MQ=255)                   aTGACGACGGTGATATTGTCTCCGAACTTTTCCCTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGAACGGTGCGG‑T‑TTATTGAGTGGCGCAATGAAGGCCTGCCGCa                                                                                                  >  2:205589/1‑133 (MQ=255)                    tGACGACGGTGATATTGTCTCCGAACTTTTCCCTAACGAAGCACAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTT‑TT‑TTGCGTGGCGCAGTGAAAGCCTGCCGCAg                                                                                                 >  2:633934/1‑133 (MQ=255)                        gacgGTGATATTGTCTCCGAACTTTTCCCTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTT‑TT‑TTGCGTGGCGCAGTGAAAGCCTGCCGCAgcggc                                                                                             >  1:402625/1‑133 (MQ=255)                             tGATATTGGCTCCGAGCTTTTCACTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTT‑TT‑TTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGgcgtgcg                                                                                        <  1:311396/133‑1 (MQ=255)                                     tctcCGAACTTTTCCCTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTT‑TT‑TTGCTTGGCGCAGTGAATGACAACCGCAGCGGCGTTGGTCGCTggc                                                                                >  2:311396/1‑131 (MQ=255)                                     tctcCGAACTTTTCCCTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTT‑TT‑TTGCGTGGCGCAGTGTAAGCCTGGCGCAGCAGCGTGCGTCGCTGTc                                                                                >  1:180609/1‑133 (MQ=255)                                     tctcCGAACTTTTCCCTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTT‑TT‑TTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTc                                                                                >  1:584453/1‑133 (MQ=255)                                      ctcCGAACTTTTCCCTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTT‑TT‑TTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCa                                                                               >  1:324278/1‑133 (MQ=255)                                      ctcCGAACTTTTCCCTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTT‑TT‑TTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCa                                                                               >  1:621153/1‑133 (MQ=255)                                         cGAACTTTTCCCTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTT‑TT‑TTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCAttt                                                                            >  2:125970/1‑133 (MQ=255)                                         cGAACTTTTCCCTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTT‑TT‑TTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCAttt                                                                            >  2:423775/1‑133 (MQ=255)                                         cGAACTTTTCCCTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTT‑TT‑TTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCAttt                                                                            >  1:500544/1‑133 (MQ=255)                                                   ccTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTT‑TT‑TTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTAt                                                                  <  2:460402/133‑1 (MQ=255)                                                       aCGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTT‑TT‑TTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAgg                                                              >  1:104550/1‑133 (MQ=255)                                                        cGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTT‑TT‑TTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGa                                                             >  1:90521/1‑133 (MQ=255)                                                          aagcgcaagcgcGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTT‑TT‑TTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAg                                                           >  2:523634/1‑133 (MQ=255)                                                                 agcgcGGGTACAAGCCCAGGAAGAAAAACGCGCTTACAACTCCGGTACGGTGCGCTT‑TT‑TTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGgcgc                                                    <  2:167463/133‑1 (MQ=255)                                                                   cgcgGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTT‑TT‑TTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGgcgcgc                                                  <  2:324278/133‑1 (MQ=255)                                                                   cgcgGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTT‑TT‑TTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGgcgcgc                                                  <  1:125970/133‑1 (MQ=255)                                                                   cgcgGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTT‑TT‑TTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGgcgcgc                                                  <  1:423775/133‑1 (MQ=255)                                                                   cgcgGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTT‑TT‑TTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGgcg                                                     >  2:140070/1‑130 (MQ=255)                                                                   cgcgGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTT‑TT‑TTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGgcg                                                     <  1:140070/130‑1 (MQ=255)                                                                     cgGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTT‑TT‑TTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAGGATGGCGCGCTg                                                >  1:31244/1‑133 (MQ=255)                                                                       ggTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTT‑TT‑TTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGtt                                              >  2:57142/1‑133 (MQ=255)                                                                         tAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTT‑TT‑TTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGc                                            <  2:90521/133‑1 (MQ=255)                                                                         tAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTT‑TT‑TTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGc                                            <  2:621153/133‑1 (MQ=255)                                                                         tAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTT‑TT‑TTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGc                                            <  1:481494/133‑1 (MQ=255)                                                                          aGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTT‑TT‑TTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCa                                           <  2:500544/133‑1 (MQ=255)                                                                                          aaaGGCGATTTCAACGCCGGTACGCTGCACTT‑TT‑TTGCGTGACGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCAcgc                           <  2:180609/133‑1 (MQ=255)                                                                                          aaaGGCGATTACAACTCCGGTACGGTGCGCTT‑TT‑TTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATTAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCAcgc                           <  2:584453/133‑1 (MQ=255)                                                                                          aaaGGCGATTACAACTCCGGTACGGTGCGCTT‑TT‑TTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCAcgc                           >  2:496841/1‑133 (MQ=255)                                                                                          aaaGGCGATTACAACTCCGGTACGGTGCGCTT‑TT‑TTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCAcgc                           <  1:523634/133‑1 (MQ=255)                                                                                           aaGGCGATGACAACTACGTTACAGTGCGC‑TCTT‑TTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCAcgcg                          <  1:231165/133‑1 (MQ=255)                                                                                            aGGCGATTACAACTCCGGTACGGTGCGCTT‑TT‑TTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGa                         <  1:63378/133‑1 (MQ=255)                                                                                                     cAACTCCGGTACGGTGCGCTT‑TT‑TTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCgg                <  1:633934/133‑1 (MQ=255)                                                                                                        cTCCGGTACGGTGCGCTT‑TT‑TTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCACGAAGATGGGGCGCTGTTCCAACAGTCGTTCTCACACGACTGTATCGGTAc             >  2:245421/1‑133 (MQ=255)                                                                                                         tCCGGTACGGTGCGCTT‑TT‑TTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACg            <  1:496841/133‑1 (MQ=255)                                                                                                          ccGGTACGGTGCGCTT‑TT‑TTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGc           <  2:402625/133‑1 (MQ=255)                                                                                                                   tGCGCTT‑TT‑TTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGa  >  1:410609/1‑133 (MQ=255)                                                                                                                        |                                                                                                                                  CTTCCCAGGGCGTCACTAATGACGACGGTGATATTGTCTCCGAACTTTTCCCTAACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGC‑T‑TT‑TTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGA  >  NC_000913/2949870‑2950115

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