Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	510,138	T→C	12.5%	I52T (ATC→ACC)	trpD →	anthranilate phosphoribosyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	510,138	0	T	C	12.5%	76.5 / 3.6	32	I52T (ATC→ACC)	trpD	anthranilate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base T (12/16);  new base C (2/2);  total (14/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.00e-01

CCGTATCGTCGGGCAGTTGGACCTGACTACCGAAGAAATGCGCGAGGTCATGCGCCAGATCATGACCGGTCAGTGCAGCGAGGCGCAGATAGGCGCTTTCCTGATGGGCATGCGCATGAAAAGCGAAAGCATCGACGAAATCGTCGGTGCGGTGTCGGTGATGCGTGAGCTGGCCGAAAAGGTCGAGCTGCAAAGCCTCGATGGTGTCGTCGATATTGTCGGGACCGGTGGTGATGGCGCCAACATCTTCAACGTGTCCACCGCTTCGTCCTTC  >  NC_002947/510007‑510280                                                                                                                                    |                                                                                                                                               ccGTATCGTCGGGCAGTTGGACCTGACTACCGAAGAAATGCGCGAGGTCATGCGCCAGATCATGACCGGTCAGTGCAGCGCGGGCCCCAAACGGGCCTTCCTTAGTGGCATGCGCATGAAAAGCGAAAGCATCGACGAAATCGTcggtgc                                                                                                                              >  1:228654/1‑150 (MQ=255)    tATCGTCGGGCAGTTGGACCTGACTACCGAAGAAATGCGCGAGGTCATGCGCCAGATCATGACCGGTCAGTGCAGCGAGGCGCAGATAGGCGCTTTCCTGATGGGCATGCGCATGAAAAGCGAAAGCATcgac                                                                                                                                            >  2:94871/1‑133 (MQ=255)    tATCGTCGGGCAGTTGGACCTGACTACCGAAGAAATGCGCGAGGTCATGCGCCAGATCATGACCGGTCAGTGCAGCGAGGCGCAGATAGGCGCTTTCCTGATGGGCATGCGCATGAAAAGCGAAAGCATcgac                                                                                                                                            <  1:94871/133‑1 (MQ=255)    tATCGTCGGGCAGTTGGACCTGACTACCGAAGAAATGCGCGAGGTCATGCGCCAGATCATGACCGGTCAGTGCAGCGAGGCGCAGATAGGCGCTTTCCTGATGGGCATGCGCATGAAAAGCGAAAGCACCGACGATTTCGTcgatgcttt                                                                                                                           >  1:365659/1‑147 (MQ=255)     aTCGTCGGGCAGTTGGACCTGACTACCGAAGAAATGCGCGAGGTCATGCGCCAGATCATGACCGGTCAGTGCAGCGAGGCGCAGATAGGCGCTTTCCTGATGGGCATGCGCATGAAAAGCGAAAGCACCGACGATTTCGTcgatgctttc                                                                                                                          <  2:365659/150‑5 (MQ=255)         tcgGGCAGTTGGACCTGACTACCGAAGAAATGCGCGAGGTCATGCGCCAGATCATGACCGGTCAGTGCAGCGAGGCGCAGATAGGCGCTTTCCTGATGGGCATGCGCATGAAAAGCGAAAGCACCGACGATTTCGTcgatgcttt                                                                                                                           <  1:202042/145‑4 (MQ=255)         tcgGGCAGTTGGACCTGACTACCGAAGAAATGCGCGAGGTCATGCGCCAGATCATGACCGGTCAGTGCAGCGAGGCGCAGATAGGCGCTTTCCTGATGGGCATGCGCATGAAAAGCGAAAGCACCGACGATTTCGTcgatgcttt                                                                                                                           >  2:202042/1‑142 (MQ=255)           gggCAGTTGGACCTGACTAGAGAAGAAATGACAGAGGTCATGCTCCAGATCATGACCTGTCAGTGCAGCGAGGCGCAGATAGGCGCTTTCCTGATGGGCATGCGCATGAAAAGCGAAAGCATCGACGAAATCGTCGGTGCGGTGTCGGTg                                                                                                                    >  2:315866/1‑150 (MQ=255)           gggCAGTTGGACCTGACTACCGAAGAAATGCGCGAGGTCATGCGCCAGATCATGACCGGTCAGTGCAGCGAGGCGCAGATAGGCGCTTTCCTGATGGGCATGCGCATGAAAAGCGAAAGCATCGACGAAATCGTCGGTGCGGTGTCGGTg                                                                                                                    <  2:179288/150‑1 (MQ=255)             gCAGTTGGACCTGACTACCGAAGAAATGCGCGAGGTCATGCGCCAGATCATGAACCAGCAGTGCAGCGAGGCGCAGATAGGCGCTTTCCTGATGGGCATGCGCATGAAAAGCGAAAGCATCGACGAAATCGTCGGTGCGGTGTCGGTGAt                                                                                                                  <  1:469364/150‑1 (MQ=255)               aGTTGGACCTGACTACCGAAGAAATGCGCGAGGTCATGCGCCAGATCATGACCGGTCAGTGCAGCGAGGCGCAGATAGGCGCTTTCCTGATGGGCATGCGCATGAAAAGCGAAAGCATCGACGAAATCGTCGGTGCGGTGTCGGTGATGc                                                                                                                <  1:2455/150‑1 (MQ=255)                       cTGACTACCGAAGAAATGCGCGAGGTCATGCGCCAGATCATGACCGGTCAGTGCAGCGAGGCGCAGATAGGCGCTTTCCTGATGGGCATGCGCATGAAAAGCGAAAGCATCGACGAAATCGTCGGTGCGGTGTCGGTGATGCGt                                                                                                              <  2:512958/144‑1 (MQ=255)                       cTGACTACCGAAGAAATGCGCGAGGTCATGCGCCAGATCATGACCGGTCAGTGCAGCGAGGCGCAGATAGGCGCTTTCCTGATGGGCATGCGCATGAAAAGCGAAAGCATCGACGAAATCGTCGGTGCGGTGTCGGTGATGCGt                                                                                                              >  1:512958/1‑144 (MQ=255)                            tACCGAAGAAATGCGCGAGGTCATGCGCCAGATCATGACCGGTCAGTGCAGCGAGGCGCAGATAGGCGCTTTCCTGATGGGCATGCGCATGAAAAGCGAAAGCATCGACGAAATCGTCGGTGCGGTGTCGGTGATGCGTGAGCTGGCCGa                                                                                                   >  2:488989/1‑150 (MQ=255)                              ccGAAGAAATGCGCGAGGTCATGCGCCAGATCATGACCGGTCAGTGCAGCGAGGCGCAGATAGGCGCTTTCCTGATGGGCATGCGCATGAAAAGCGAAAGCATCGACGAAATCGTCGGTGCGGTGTCGGTGATGCGTGAGCTGGCCGaaa                                                                                                 >  2:176117/1‑150 (MQ=255)                                gaagaaATGCGCGAGGTCATGCGCCAGATCATGACCGGTCAGTGCAGCGAGGCGCAGATAGGCGCTTTCCTGATGGGCATGCGCATGAAAAGCGAAAGCATCGACGAAATCGTCGGTGCGGTGTCGGTGATGCGTGAGCTGGCCGaaaaa                                                                                               >  1:71844/1‑149 (MQ=255)                                 aagaaATGCGCGAGGTCATGCGCCAGATCATGACCGGTCAGTGCAGCGAGGCGCAGATAGGCGCTTTCCTGATGGGCATGCGCATGAAAAGCGAAAGCATCGACGAAATCGTCGGTGCGGTGTCGGTGATGCGTGAGCTGGCCGAAAAgg                                                                                              <  1:488989/150‑1 (MQ=255)                                          gcgAGGTCATGCGCCAGATCATGACCGGTCAGTGCAGCGAGGCGCAGATAGGCGCTTTCCTGATGGGCATGCGCATGAAAAGCGAAAGCATCGACGAAATCGTCGGTGCGGTGTCGGTGATGCGTGAGCTGGCCGAAAAGGTCGAGCTGc                                                                                     <  1:455975/150‑1 (MQ=255)                                                  aTGCGCCAGATCATGACCGGTCAGTGCAGCGAGGCGCAGATAGGCGCTTTCCTGATGGGCATGCGCATGAAAAGCGAAAGCATCGACGAAATCGTCGGTGCGGTGTCGGTGATGCGTGAg                                                                                                           >  1:424297/1‑120 (MQ=255)                                                  aTGCGCCAGATCATGACCGGTCAGTGCAGCGAGGCGCAGATAGGCGCTTTCCTGATGGGCATGCGCATGAAAAGCGAAAGCATCGACGAAATCGTCGGTGCGGTGTCGGTGATGCGTGAg                                                                                                           <  2:424297/120‑1 (MQ=255)                                                                gACCGGTCAGTGCAGCGAGGCGCAGATAGGCGCTTTCCTGATGGGCATGCGCATGAAAAGCGAAAGCATCGACGAAATCGTCGGTGCGGTGTCGGTGATGCGt                                                                                                              <  2:372415/103‑1 (MQ=255)                                                                gACCGGTCAGTGCAGCGAGGCGCAGATAGGCGCTTTCCTGATGGGCATGCGCATGAAAAGCGAAAGCATCGACGAAATCGTCGGTGCGGTGTCGGTGATGCGt                                                                                                              >  1:372415/1‑103 (MQ=255)                                                                 aCCGGTCAGTGCAGCGAGGCGCAGATAGGCGCTTTCCTGATGGGCATGCGCATGAAAAGCGAAAGCATCGACGAAATCGTCGGTGCGGTGTCGGTGATGCGTGAGCTGGCCGAAAAGGTCGAGCTGCAAAGCCTCGATGGTGTCGTCGat                                                              <  1:9998/150‑1 (MQ=255)                                                                        aGTGCAGCGAGGCGCAGATAGGCGCTTTCCTGATGGGCATGCGCATGAAAAGCGAAAGCATCGACGAAATCGTCGGTGCGGTGTCGGTGATGCGTGAGCTGGCCGAAAAGGTCGAGCTGCAAAGCCTCGATGGTGTCGTCGATATTGTCg                                                       <  2:228654/150‑1 (MQ=255)                                                                         gTGCAGCGAGGCGCAGATAGGCGCTTTCCTGATGGGCATGCGCATGAAAAGCGAAAGAATCGACGAAATCGTCGGTGCGGTGTCGGTGATGCGTGAGCTGGCCGAAAAGGTCGAGCTGCAAAGCCTCGATGGTGTCGTCGATATTGTCgg                                                      <  1:315866/150‑1 (MQ=255)                                                                                      ctgATAGGCGCTTTCCTGATGGGCATGCGCATGAAAAGCGAAAGCATCGACGAAATCGTCGGTGCGGTGTCGGTGATGCGTGAGCTGGCCGAAAAGGTCGAGCTGCAAAGCCTCGATGGTGTCGTCGATATTGTCGGGACCGGTGGTGAt                                         >  2:471632/3‑150 (MQ=255)                                                                                        gATAGGCGCTTTCCTGATGGGCATGCGCATGAAAAGCGAAAGCATCGACGAAATCGTCGGTGCGGTGTCGGTGATGCGTGAGCTGGCCGAAAAGGTCGAGCTGCAAAGCCTCGATGGTGTCGTCGATATTGTCGGGACCGGTGGTGATgg                                       <  1:471632/150‑1 (MQ=255)                                                                                            ggCGCTTTCCTGATGGGCATGCGCATGAAAAGCGAAAGCATCGACGAAATCGTCGGTGCGGTGTCGGTGATGCGTGAGCTGGCCGAAAAGGTCGAGCTGCAAAGCCTCGATGGTGTCGTCGATATTGTCGGGACCGGTGGTGATGGCGcc                                   <  2:267731/150‑1 (MQ=255)                                                                                                        aTGGGCATGCGCATGAAAAGCGAAAGCATCGACGAAATCGTCGGTGCGGTGTCGGTGATGCGTGAGCTGGCCGAAAAGGACGAGCTGCAAAACCTCGATCGTGTCGTCGATATTGTATGGACCGGTCGTGATGTCGCCTAAATCTTCAAc                       >  2:61793/1‑150 (MQ=255)                                                                                                                             gAAAGCATCGACGAAATCGTCGGTGCGGTGTCGGTGATGCGTGAGCTGGCCGAAAAGGTCGAGCTGCAAAGCCTCGATGGTGTCGTCGATATTGTCGGGACCGGTGGTGATGGCGCCAACATCTTCAACGTGTCCACCGcttcgtccttc  <  1:176117/150‑1 (MQ=255)                                                                                                                                  cATCGACGAAATCGTCGGTGCGGTGTCGGTGATGCGTGAGCTGGCCGAAAAGGTCGAGCTGCAAAGCCTCGATGGTGTCGTCGATATTGTCGGGACCGGTGGTGATGGCGCCAACATCTTCAACGTGTCCACCGcttcgtccttc  >  1:592367/1‑145 (MQ=255)                                                                                                                                  cATCGACGAAATCGTCGGTGCGGTGTCGGTGATGCGTGAGCTGGCCGAAAAGGTCGAGCTGCAAAGCCTCGATGGTGTCGTCGATATTGTCGGGACCGGTGGTGATGGCGCCAACATCTTCAACGTGTCCACCGcttcgtccttc  <  2:592367/145‑1 (MQ=255)                                                                                                                                    |                                                                                                                                               CCGTATCGTCGGGCAGTTGGACCTGACTACCGAAGAAATGCGCGAGGTCATGCGCCAGATCATGACCGGTCAGTGCAGCGAGGCGCAGATAGGCGCTTTCCTGATGGGCATGCGCATGAAAAGCGAAAGCATCGACGAAATCGTCGGTGCGGTGTCGGTGATGCGTGAGCTGGCCGAAAAGGTCGAGCTGCAAAGCCTCGATGGTGTCGTCGATATTGTCGGGACCGGTGGTGATGGCGCCAACATCTTCAACGTGTCCACCGCTTCGTCCTTC  >  NC_002947/510007‑510280

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑