Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	966,040	A→G	D399G (GAT→GGT)	ycaQ →	winged helix DNA‑binding protein YcaQ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	966,040	0	A	G	100.0%	117.3 / NA	31	D399G (GAT→GGT)	ycaQ	winged helix DNA‑binding protein YcaQ
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (15/16);  total (15/16)

TTACAGGAAGGTATTAAACCAACGACAACGCTGCAAAAAGGGCTACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCTCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCATAACCCGATCCATCTGGAGGAACTATGGATCATCGTCTGCTTGAAATCATTG  >  NZ_CP009273/965913‑966156                                                                                                                                |                                                                                                                     ttACAGGAAGGTATTAAACCAACGACAACGCTGCAAAAAGGGCTACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCTCTTTACTGGTTGCCGCACCg                                                                                                           >  1:547679/1‑139 (MQ=255)    cAGGAAGGTATTAAACCAACGACAACGCTGCAAAAAGGGCTACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCTCTTTACTGGTTGCCGCACCGGCt                                                                                                        >  1:35340/1‑139 (MQ=255)    cAGGAAGGTATTAAACCAACGACAACGCTGCAAAAAGGGCTACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCTCTTTACTGGTTGCCGCACCGGCt                                                                                                        <  2:184159/139‑1 (MQ=255)        aaGGTATTAAACCAACGACAACGCTGCAAAAAGGGCTACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCTCTTTACTGGTTGCCGCACCGGCTggg                                                                                                     >  1:435289/1‑138 (MQ=255)                   ccAACGACAACGCTGCAAAAAGGGCTACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCTCTTTACTGGTTGCCGCACCGGCTGGGAAATAGACCCCg                                                                                         <  1:246035/139‑1 (MQ=255)                   ccAACGACAACGCTGCAAAAAGGGCTACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCTCTTTACTGGTTGCCGCACCGGCTGGGAAATAGACCCCg                                                                                         <  1:536223/139‑1 (MQ=255)                      aCGACAACGCTGCAAAAAGGGCTACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCTCTTTACTGGTTGCCGCACCgg                                                                                                          <  2:376955/119‑1 (MQ=255)                      aCGACAACGCTGCAAAAAGGGCTACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCTCTTTACTGGTTGCCGCACCgg                                                                                                          >  1:376955/1‑119 (MQ=255)                      aCGACAACGCTGCAAAAAGGGCTACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCTCTTTACTGGTTGCCGCACCGGCTGGGAAATAGACCCCGt                                                                                        >  2:399767/1‑137 (MQ=255)                      aCGACAACGCTGCAAAAAGGGCTACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCTCTTTACTGGTTGCCGCACCGGCTGGGAAATAGACCCCGt                                                                                        <  1:399767/137‑1 (MQ=255)                      aCGACAACGCTGCAAAAAGGGCTACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCTCTTTACTGGTTGCCGCACCGGCTGGGAAATAGACCCCGTCg                                                                                      >  2:255496/1‑139 (MQ=255)                                tGCAAAAAGGGCTACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCTCTTTACTGGTTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGAt                                                                            >  2:625445/1‑139 (MQ=255)                                tGCAAAAAGGGCTACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCTCTTTACTGGTTGCCGCACCGGCTGGGAAATAGA‑CCCGTCGCATAAATGAt                                                                            >  2:195127/1‑138 (MQ=255)                                  cAAAAAGGGCTACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCTCTTTACTGGTTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGAtt                                                                           <  1:255496/138‑1 (MQ=255)                                  cAAAAAGGGCTACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCTCTTTACTGGTTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGAt                                                                            <  1:59014/137‑1 (MQ=255)                                      aaGGGCTACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCTCTTTACTGGTTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGAt                                                                      <  1:515405/139‑1 (MQ=255)                                          gCTACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCTCTTTACTGGTTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGAt                                                                      >  1:257296/1‑135 (MQ=255)                                               gTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCTCTTTACTGGTTGCCGCACCGGCTGGGAAATAGa                                                                                              <  1:593326/106‑1 (MQ=255)                                               gTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCTCTTTACTGGTTGCCGCACCGGCTGGGAAATAGa                                                                                              >  2:593326/1‑106 (MQ=255)                                               gTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCTCTTTACTGGTTGCCGCACCGGCTGGGAAATAGA‑CCCGTCGCATAAATGATTATGATAAGATTTATg                                                            <  1:195127/139‑1 (MQ=255)                                                tCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCTCTTTACTGGTTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATg                                                            <  2:449867/139‑1 (MQ=255)                                                    gCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCTCTTTACTGGTTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCAtt                                                        <  1:625445/139‑1 (MQ=255)                                                              aTTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCTCTTTACTGGTTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCATAACCCGa                                              <  2:183377/139‑1 (MQ=255)                                                                        cTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCTCTTTACTGGTTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCATAACCCGATCCATCTgga                                    <  1:286618/139‑1 (MQ=255)                                                                                 ggCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCTCTTTACTGGTTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCATAACCCGATCCATCTGGAGGAACTATg                           >  2:596225/1‑139 (MQ=255)                                                                                                         ctgcCCGCAAGGCCTCTTTACTGGTTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCATAACCCGATCCATCTGGAGGAACTATGGATCATCGTCTGCTTGAAATCAtt   >  2:320786/1‑139 (MQ=255)                                                                                                          tgcCCGCAAGGCCTCTTTACTGGTTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGc                                                           <  2:532649/82‑1 (MQ=255)                                                                                                          tgcCCGCAAGGCCTCTTTACTGGTTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGc                                                           >  1:532649/1‑82 (MQ=255)                                                                                                          tgcCCGCAAGGCCTCTTTACTGGTTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCATAACCCGATCCATCTGGAGGAACTATGGATCATCGTCTGCTTGAAATCATtg  >  2:378782/1‑139 (MQ=255)                                                                                                            cccGCAAGGCCTCTTTACTGGTTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGAt                                                                      >  2:546367/1‑69 (MQ=255)                                                                                                            cccGCAAGGCCTCTTTACTGGTTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGAt                                                                      <  1:546367/69‑1 (MQ=255)                                                                                                                                |                                                                                                                     TTACAGGAAGGTATTAAACCAACGACAACGCTGCAAAAAGGGCTACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCTCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCATAACCCGATCCATCTGGAGGAACTATGGATCATCGTCTGCTTGAAATCATTG  >  NZ_CP009273/965913‑966156

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