Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,162,505	A→G	16.6%	I95V (ATT→GTT)	murI →	glutamate racemase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	1,162,505	0	A	G	16.6%	47.5 / 4.7	25	I95V (ATT→GTT)	murI	glutamate racemase
Reads supporting (aligned to +/- strand):  ref base A (9/11);  new base G (2/2);  total (11/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.66e-01

CAGTTGAAATCGCTCGTAAATTAATGGAATTTGATATAAAAATGCTCGTGATTGCTTGTAATACAGCAACTGCTGTAGCTTTAGAATATTTACAAAAGACCTTATCAATCCCAGTGATTGGCGTAATTGAACCAGGTGCTAGAACAGCAATAATGACTACTAGAAATCAAAATGTATTAGTACTAGGAACGGAAGGCACAATTAAATCTGAAGCATATCGTACGCATATTAAACGTATCAATCCACATGTAGAGGTAC  >  CP000730/1162380‑1162637                                                                                                                              |                                                                                                                                     cAGTTGAAATCGCTCGTAAATTAATGGAATTTGATATAAAAATGCTCGTGATTGCTTGTAATACAGCAACTGCTGTAGCTTTAGAATATTTACAAAAGACCTTATCAATCCCAGTGATTGGCGTAATTGAACAAGGTGCTa                                                                                                                       >  2:222118/1‑141 (MQ=255)  aGTTGAAAGCGCCCGTAAATTAATGGAATTTGATATAAAAATGCTCGTGATTGCTTGTAATAACGCAACTGCTGTAGCTTTAGAATATTTACAAAAGACCTTATCAATCCCAGTGATTGGCGTAATTGAACCAGGTGCTAg                                                                                                                      <  1:119958/141‑1 (MQ=255)         aTCGCTCGTAAATTAATGGAATTTGATATAAAAATGCTCGTGATTGCTTGTAATACAGCAACTGCTGTAGCTTTAGAATATTTACAAAAGACCTTATCAATCCCAGTGATTGGCGTAATTGACCCAGGTGCTAGAACAGCa                                                                                                               >  2:108406/1‑141 (MQ=255)         aTCGCTCGTAAATTAATGGAATTTGATATAAAAATGCTCGTGATTGCTTGTAATACAGCAACTGCTGTAGCTTTAGAATATTTACAAAAGACCTTATCAATCCCAGTGATTGGCGTAATTGAACCAGGTGCTAGAACAGCa                                                                                                               >  2:60120/1‑141 (MQ=255)                  aaaTTAATGGAATTTGATATAAAAATGCTCGTGATTGCTTGTAATACAGCAACTGCTGTAGCTTTAGAATATTTACAAAAGACCTTATCAATCCCAGTGATTGGCGTAGTTGAACCAGGTGCTAGAACAGCAATAATGact                                                                                                      >  2:105156/1‑141 (MQ=255)                  aaaTTAATGGAATTTGATATAAAAATGCTCGTGATTGCTTGTAATACAGCAACTGCTGTAGCTTTAGAATATTTACAAAAGACCTTATCAATCCCAGTGATTGGCGTAGTTGAACCAGGTGCTAGAACAGCAATAATGact                                                                                                      >  2:105155/1‑141 (MQ=255)                                 gATATAAAAATGCTCGTGATTGCTTGTAATACAGCAACTGCTGTAGCTTTAGAATATTTACAAAAGACCTTATCAATCCCAGTGATTGGCGTAATTGAACCAGGTGCTAGAACAGCAATAATGACTACTAGAAATCAAAAt                                                                                       >  1:256026/1‑141 (MQ=255)                                 gATATAAAAATGCTCGTGATTGCTTGTAATACAGCAACTGCTGTAGCTTTAGAATATTTACAAAAGACCTTATCAATCCCAGTGATTGGCGTAATTGAACCAGGTGCTAGAACAGCAATAATGACTACTAGAAATCAAAAt                                                                                       <  1:17658/141‑1 (MQ=255)                                    ataAAAATGCTCGTGATTGCTTGTAATACAGCAACTGCTGTAGCTTTAGAATATTTACAAAAGACCTTATCAATCCCAGTGATTGGCGTAATTGAACCAGGTGCTAGAACAGCAATAATGACTACTAGAAATCAAAATGTa                                                                                    >  1:54804/1‑141 (MQ=255)                                                     tGCTTGTAATACAGCAACTGCTGTAGCTTTAGAATATTTACAAAAGACCTTATCAATCCCAGTGATTGGCGTAATTGAACCAGGTGCTAGAACAGCAATAATGACTACTAGAAATCAAAATGTATTAGTACTAGGAACGGa                                                                   >  1:78600/1‑141 (MQ=255)                                                        tggTAAAACAGGAACCGCTGTAACGTTAGAATATTTACAAAAGGCCGTATCAAGCCCAGTGATTGGCGTAATTGAACCAGGGGCTAGAACTGCAATAATGACTACTAGAAATCAAAATGTATGAGTACTAGGAACGGAAgg                                                                <  2:87115/139‑1 (MQ=255)                                                                    aaCTGCTGTAGCTTTAGAATATTTACAAAAGACCTTATCAATCCCAGTGATTGGCGTAATTGAACCAGGTGCTAGAACAGCAATAATGACTACTAGAAATCAAAATGTATTAGTACTAGGAACGGAAGGCACAATTAAATc                                                    >  1:7051/1‑141 (MQ=255)                                                                              gCTTTAGAATATTTACAAAAGATCTCAACAATCCCAGTGAGTGGCGTAATTGAACCAGGTGCTAGAACAGCAATAATTACCACCAACATTCATGCTGTGTCATAATGACGAACGGAAGGCACAATTAAATCTGAAGCatat                                          <  1:265961/141‑1 (MQ=255)                                                                               cTTTAGAATATTTACAAAAGACCTTATCAATCCCAGTGATTGGCGTACGTGTACCAGGTGCTAGAACAGCAATAATGACTACTAGAAATCAAAATGTATTAGTACTAGGAACGGAAGGCACAATTAAATCTGAAGCATATc                                         <  1:222118/141‑1 (MQ=255)                                                                                 ttAGAATATTTACAAAAGACCTTATCAATCCCAGTGATTGGCGTAATTGAACCAGGTGCTAGAACAGCAATAATGACTACTAGAAATCAAAATGTATTAGTACTAGGAACGGAAGGCACAATTAAATCTGAAGCATATCGt                                       <  1:210433/141‑1 (MQ=255)                                                                                 ttAGAATATTTACAAAAGACCTTATCAATCCCAGTGATTGGCGTAATTGAACCAGGGGCTAGAACAGCAATAATGACTACTAGAAATCAAAATGTATTAGTACTAGGAACGGAAGGCACAATTAAATCTGAAGCATATCGt                                       <  1:46639/141‑1 (MQ=255)                                                                                    gAATAATTACAAAAGACCTAATCAATCCGAGTGATTGGCGTAATTGAACCAGGTGCTAGAACAGCAATAATGACTACTAGAAATCAAAATGTATTAGTACTAGGAACGGAAGGCACAATTAAATCTGAAGCATATCGTACg                                    <  1:187926/141‑1 (MQ=255)                                                                                     aaTATTTACAAAAGACCTTATCAATCCCAGTGATTGGCGTAGTTGAACCAGGTGCTAGAACAGCAATAATGACTACTAGAAATCAAAATTTATTAGTACTAGGAACGGAAGGCACAATTAAATCTGAAGCATATCGTACGc                                   <  1:105155/141‑1 (MQ=255)                                                                                     aaTATTTACAAAAGACCTTATCAATCCCAGTGATTGGCGTAGTTGAACCAGGTGCTAGAACAGCAATAATGACTACTAGAAATCAAAATGTATTAGTACTAGGAACGGAAGGCACAATTAAATCTGAAGCATATCGTACGc                                   <  1:105156/141‑1 (MQ=255)                                                                                       tatTTACAAACGACCTTATCAATCCCAGGGATTGGCGTAATTGAACCAGGTGCTAGAACAGCAATAATGACTACTAGAAATCAAAATGTATTAGTACTAGGAACGGAAGGCACAATTAAATCTGAAGCATATCGTACGCat                                 <  1:90826/141‑1 (MQ=255)                                                                                         tttACAAAAGACCTTATCAATCCCAGTGATTGGCGTAATTGAACCAGGTGCTAGAACAGCAATAATGACTACTAGAAATCACAATGTATTAGTACTAGGAACGGAAGGCACAATTAAATCTGAAGCATATCGTACGCatat                               <  1:184613/141‑1 (MQ=255)                                                                                          ttACAAAAGACCTTATCAATCCCAGTGATTGGCGTAATTGAACCAGGTGCTAGAACAGCAATAATGACTACTAGAAATCAAAATGTATTAGTACTAGGAACGGAAGGCACAATTAAATCTGAAGCATATCGTACGCATAtt                              <  1:108406/141‑1 (MQ=255)                                                                                                   aCCTTATCAATCCCAGTGATTGGCGTAATTGAACCAGGTGCTAGAACAGCAATAATGACTACTAAAAATCAAAATGTATTAGTACTAGGAACGGAAGGCACAATTAAATCTGAAGCATATCGTACGCATATTAAACGTATc                     >  1:218328/1‑141 (MQ=255)                                                                                                                  gTGATTGGCGTAATTGAACCAGGTGCTAGAACAGCAATAAGGACTACTAGAAATCAAAATGTATTAGTACTAGGAACGGAAGGCACAGTTAAATCTGAAGCATATCGTACGCATATTAAAGGCATCAACACACCTGTagag      >  2:51390/1‑141 (MQ=255)                                                                                                                      ttGGCGTTATTGAACCAGATGCCAGAACAGCAATAATAAATACAAGAAATCAAAATGTATTAGTACTAGGAACGGAAGGCACAATTAAATCTGAAGCATATCGTACGCATATTAAACGTATCAATCCACATGTAGAGGTAc  <  1:178825/141‑1 (MQ=255)                                                                                                                              |                                                                                                                                     CAGTTGAAATCGCTCGTAAATTAATGGAATTTGATATAAAAATGCTCGTGATTGCTTGTAATACAGCAACTGCTGTAGCTTTAGAATATTTACAAAAGACCTTATCAATCCCAGTGATTGGCGTAATTGAACCAGGTGCTAGAACAGCAATAATGACTACTAGAAATCAAAATGTATTAGTACTAGGAACGGAAGGCACAATTAAATCTGAAGCATATCGTACGCATATTAAACGTATCAATCCACATGTAGAGGTAC  >  CP000730/1162380‑1162637

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