Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,337,394	A→G	S522G (AGC→GGC)	acnA →	aconitate hydratase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,337,394	0	A	G	97.0%	99.8 / ‑6.2	33	S522G (AGC→GGC)	acnA	aconitate hydratase 1
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (12/20);  total (12/21)
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.98e-01

CCCAAGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAAGCGAT‑TTAACCGTCGGTGCGGTGCTGT‑CCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGCACAAGAAATTGCCCGTGCGGTAGA  >  NC_000913/1337150‑1337626                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                           cccAAGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGAt‑t                                                                                                                                                                                                                                     <  1:317910/251‑1 (MQ=255)                   cGCAGTCTTGTTGGCTAAAATAGCCGTAACTGTGGGCTTCAAACGGCAACCTCGGGGCACAGCGGCGCGTGCACGGGGTTCGAAAGTTGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGGATCTCGACGAACTGAGGTTTAACCTTGTGGTATACGTTTGTACCACCTGTATTGGTAACGCTGGGCCGCTGCCCGATCCTATCGAAACGGCACGCAAAAAAGGCGAT‑TCACCCGTcgtagctgtgc                                                                                                                                                                                                                   <  2:136025/251‑10 (MQ=255)                          ttGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCACGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGCCGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGAT‑TTAACCGTCGATGCGGTGCTGT‑CCgg                                                                                                                                                                                                           >  1:99855/1‑251 (MQ=255)                                         gCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGAT‑TTAACCGTCGGTGCGGTGC‑GTGCCGGCAACCGGAACTTTg                                                                                                                                                                                             >  2:326626/1‑250 (MQ=255)                                                      gCCTCAAGCGGTAACCATGGCTCAAAGCGCCGCTGGCACCGGGTTCGAAAGGCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGAT‑TTAACCGTCGGTGCGGTGCTGT‑CCGGCACCCGTAACTTTGAAGGCCTTatccat                                                                                                                                                                               <  1:346953/251‑1 (MQ=255)                                                                 cAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGAT‑TTATCCGTCGGTGCGGTGCTGT‑CCGGCAACCGTATCTTTGAAGGCCGTATCCATCCGCTGGTTaa                                                                                                                                                                    <  1:330171/251‑1 (MQ=255)                                                                              aaGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGAT‑TTAACCGTCGGTGCGGTGCTGT‑CCGGCAACCGTAACTTTGAAGGCCGTATCCATACGCTGGTTAAAACTAActggctg                                                                                                                                                       <  1:3105/251‑1 (MQ=255)                                                                                           aCCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGAT‑TTAACCGTCGGTGCGGTGCTGT‑CCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTcgccgccgc                                                                                                                                          >  1:212913/1‑251 (MQ=255)                                                                                                gttagAAAGTCGTTTCTGATTAACGTGCAAAAGCGAAACTGACACCGTATCTCGACGAACGGGGGTTTAACCTTGGGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGAT‑TTAACCGTCGGTGCGGTGCTGT‑CCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCtggtg                                                                                                                                     <  2:332167/247‑1 (MQ=255)                                                                                                        gTCGTTTCTGATTATCTGGCAAAAGAGAAACTGACATCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGAT‑TTAACCGTCGGTGCGGTGCTGT‑CCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTa                                                                                                                             <  1:156753/251‑1 (MQ=255)                                                                                                                tGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGAT‑TTAACCGTCGGTGCGGTGCTGT‑CCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTgg                                                                                                                     >  1:47838/1‑251 (MQ=255)                                                                                                                 gATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAAATGGGGTTTAACCTTGTGGGCTACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAACAAAGGCGAT‑TTACCCGGCGGTGCGGTACTGT‑CCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTCCCTGTGCGCTGGc                                                                                                                    >  1:131570/1‑251 (MQ=255)                                                                                                                  aTTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATTCGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGAT‑TTAACCGTCGGTGCGGTGCTGT‑CCGGCAACCGTAACTTTGAAGGCCGTATCCTTCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCg                                                                                                                   <  1:289333/251‑1 (MQ=255)                                                                                                                    tATCTGGCAAAAGCGAAACTGACACCGTATCGCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGAT‑TTAACCGTCGGTTCGGTGCGGT‑CCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTATCGGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCgg                                                                                                                  <  1:413843/250‑1 (MQ=255)                                                                                                                                  gAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGa                                                                                                                                                                                                                                        <  1:186740/120‑1 (MQ=255)                                                                                                                                  gAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGa                                                                                                                                                                                                                                        >  2:186740/1‑120 (MQ=255)                                                                                                                                     aCTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGAT‑TTAACCGTCGGTGCGGTGCTGT‑CCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAAcc                                                                                                >  1:361675/1‑251 (MQ=255)                                                                                                                                       tGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGAT‑TTAACCGTCGGTGCGGTGCTGT‑CCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTg                                                                                              >  1:199310/1‑251 (MQ=255)                                                                                                                                         acacCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGAT‑TTAACCGTCGGTGCGGTGCTGT‑CCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGc                                                                                            <  1:326626/251‑1 (MQ=255)                                                                                                                                         acacCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGAT‑TTAACCGTCGGTGCGGTGCTGT‑CCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGc                                                                                            <  1:293984/251‑1 (MQ=255)                                                                                                                                         acacCGAATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCCCCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGAT‑TTAACCGTCGGTGCGGTGCTGT‑CCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTCTGCGCTGGCGGGAGATATGTATATCACCCTGGc                                                                                            >  2:182779/1‑251 (MQ=255)                                                                                                                                            ccGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGAT‑TTAACCGTCGGTGCGGTGCTGT‑CCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTc                                                                                         <  1:322020/251‑1 (MQ=255)                                                                                                                                                       cgaACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGAT‑TTAACCGTCGGTGCGGTGCTGT‑CCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCg                                                                              <  1:41925/251‑1 (MQ=255)                                                                                                                                                                gTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGAT‑TTAACCGTCGGTGCGGTGCTGT‑CCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGGGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATc                                                                     <  1:201626/251‑1 (MQ=255)                                                                                                                                                                                       gTACAACCTGTATTGGTAACTCTGGGCCGCGGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGAT‑TTAACCGTCGGTGCGGTGCTGT‑CCGGCAACCGTAACTTTGAAGGCCGTATTCAGCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCt                                               <  1:272171/250‑1 (MQ=255)                                                                                                                                                                                         tgccccTGTATTTGTAACTCTTGGCGGTCGCCCGCCCCTATCGAAACTGCAA‑CAAAAAAAGTGGTATTAACCGTCGGTGCGGTGCTGT‑CCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTGAAAACTAACTTGCTGGCCTTGCCGCCGCTGGTGGTTGTCTATGCGATGGCGGGAAATATGAATATCAACCGGGCGTCGGAGCCGATCGGGCATGATCTCaaatgcttgccgtttgatctga                                             <  2:192966/246‑20 (MQ=255)                                                                                                                                                                                          ccaccTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCTAAACGGCAATCAAAAAAGGCGAT‑TTAACCGTCGGTGCGGTGCTGT‑CCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGaaa                                           <  1:314063/251‑1 (MQ=255)                                                                                                                                                                                                  aTTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGAT‑TTAACCGTCGGTGCGGTGCTGT‑CCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCt                                    <  1:271600/250‑1 (MQ=255)                                                                                                                                                                                                       tAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGAT‑TTAACCGTCGGTGCGGTGCTGT‑CCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTCAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGTTCCCAAAGGCGATCCGGTTTATCTGAAAGATATCCGGCCAt                              >  1:267408/1‑251 (MQ=255)                                                                                                                                                                                                        aatgctGGGCCGTTTCCCGATCTTACTGAAACGGAAATCAAAAAAGGCGAT‑TTAACCGTCGGTGCGGTGCTGT‑CCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACGGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCAt                              <  2:310031/246‑1 (MQ=255)                                                                                                                                                                                                                         cGATCCTATCGAAACGGCAATCAAAAAAGGCGAT‑TTAACCGTCGGTGCGGTGCTGT‑CCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGCACAAGAAATTGccc            >  1:244836/1‑251 (MQ=255)                                                                                                                                                                                                                                  cGAAAAGGCAATCAAAAAAGGCGAT‑TTAACCGTCGGTGCGGTGCTGT‑CCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGGTGCCTATGCGCTGGTGGGAAATATGAATATCAACCTTGCTTCTGATCCTATCGGCCATGATCGCAAAGGCGATTCGGTTTATCTGAAAGATATCTGGCCATCGGCACAAGAAATTGCCCGTGCGGTa    <  1:440886/250‑1 (MQ=255)                                                                                                                                                                                                                                    aaaCGGCAATCAAAAAAGGCGAT‑TTAACCGTCGGTGCGGTGCTGT‑CCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGCACAAGAAATTGCCCGTGCGGTAgc  >  1:19907/1‑249 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                           CCCAAGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAAGCGAT‑TTAACCGTCGGTGCGGTGCTGT‑CCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGCACAAGAAATTGCCCGTGCGGTAGA  >  NC_000913/1337150‑1337626

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