Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	349,456	A→G	T259A (ACG→GCG)	prpB →	2‑methylisocitrate lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	349,456	0	A	G	100.0%	132.8 / NA	37	T259A (ACG→GCG)	prpB	2‑methylisocitrate lyase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (15/22);  total (15/22)

ACATCACCGAATTTGGTGCCACGCCGCTGTTTACCACCGACGAATTACGCAGCGCCCATGTCGCAATGGCGCTGTACCCACTTTCAGCGTTCCGCGCCATGAACCGCGCCGCTGAACATGTCTACAACGTCCTGCGCCAGGAAGGCACGCAGAAAAGCGTCATCGACACCATGCAGACCCGCAACGAGCTGTACGAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGC  >  NC_000913/349310‑349600                                                                                                                                                   |                                                                                                                                                 aCATCACCGAATTTGGTGCCACGCCGCTGTTTACCACCGACGAATTACGCAGCGCCCATGTCGCAATGGCGCTGTACCCACTTTCAGCGTTCCGCGCCATGAACCGCGCCGCTGAACATGTCTACAACGTCCTGCGCCAGGAAGGCGCGCa                                                                                                                                              <  2:492823/151‑1 (MQ=255)       ccGAATTTGGTGCCACGCCGCTGTTTACCACCGACGAATTACGCAGCGCCCATGTCGCAATGGCGCTGTACCCACTTTCAGCGTTCCGCGCCATGAACCGCGCCGCTGAACATGTCTACAACGTCCTGCGCCAGGAAGGCGCGCAGaaa                                                                                                                                          >  2:164466/1‑149 (MQ=255)        cGAATTTGGTGCCACGCCGCTGTTTACCACCGACGAATTACGCAGCGCCCATGTCGCAATGGCGCTGTACCCACTTTCAGCGTTCCGCGCCATGAACCGCGCCGCTGAACATGTCTACAACGTCCTGCGCCAGGAAGGCGCGCAGAAAAGc                                                                                                                                       >  1:8277/1‑151 (MQ=255)          aaTTTGGTGCCACGCCGCTGTTTACCACCGACGAATTACGCAGCGCCCATGTCGCAATGGCGCTGTACCCACTTTCAGCGTTCCGCGCCATGAACCGCGCCGCTGAACATGTCTACAACGTCCTGCGCCAGGAAGGCGCGCAGAAAAGCGt                                                                                                                                     <  1:666200/151‑1 (MQ=255)                     aCGCCGCTGTTTACCACCGACGAATTACGCAGCGCCCATGTCGCAATGGCGCTGTACCCACTTTCAGCGTTCCGCGCCATGAACCGCGCCGCTGAACATGTCTACAACGTCCTGCGCCAGGAAGGCGCGCAGAAAAGCGTCATCGACACCa                                                                                                                          <  1:473504/151‑1 (MQ=255)                               ttACCACCGACGAATTACGCAGCGCCCATGTCGCAATGGCGCTGTACCCACTTTCAGCGTTCCGCGCCATGAACCGCGCCGCTGAACATGTCTACAACGTCCTGCGCCAGGAAGGCGCGCAGAAAAGCGTCATCGACACCATGCAGAccc                                                                                                                 <  2:592228/150‑1 (MQ=255)                                  ccaccGACGAATTACGCAGCGCCCATGTCGCAATGGCGCTGTACCCACTTTCAGCGTTCCGCGCCATGAACCGCGCCGCTGAACATGTCTACAACGTCCTGCGCCAGGAAGGCGCGCAGAAAAGCGTCATCGACACCATGCAGACCCGCaa                                                                                                             <  2:22999/151‑1 (MQ=255)                                  ccaccGACGAATTACGCAGCGCCCATGTCGCAATGGCGCTGTACCCACTTTCAGCGTTCCGCGCCATGAACCGCGCCGCTGAACATGTCTACAACGTCCTGCGCCAGGAAGGCGCGCAGAAAAGCGTCATCGACACCATGCAGACCCGCa                                                                                                              >  1:527670/1‑150 (MQ=255)                                   caccGACGAATTACGCAGCGCCCATGTCGCAATGGCGCTGTACCCACTTTCAGCGTTCCGCGCCATGAACCGCGCCGCTGAACATGTCTACAACGTCCTGCGCCAGGAAGGCGCGCAGAAAAGCGTCATCGACACCATGCAGACCCGCAAc                                                                                                            <  2:100259/151‑1 (MQ=255)                                   caccGACGAATTACGCAGCGCCCATGTCGCAATGGCGCTGTACCCACTTTCAGCGTTCCGCGCCATGAACCGCGCCGCTGAACATGTCTACAACGTCCTGCGCCAGGAAGGCGCGCAGAAAAGCGTCATCGACACCATGCAGACCCGCAAc                                                                                                            >  1:483478/1‑151 (MQ=255)                                       gacgaATTACGCAGCGCCCATGTCGCAATGGCGCTGTACCCACTTTCAGCGTTCCGCGCCATGAACCGCGCCGCTGAACATGTCTACAACGTCCTGCGCCAGGAAGGCGCGCAGAAAAGCGTCATCGACACCATGCAGACCCGCAACGAGc                                                                                                        <  1:106581/151‑1 (MQ=255)                                           aaTTACGCAGCGCCCATGTCGCAATGGCGCTGTACCCACTTTCAGCGTTCCGCGCCATGAACCGCGCCGCTGAACATGTCTACAACGTCCTGCGCCAGGAAGGCGCGCAGAAAAGCGTCATCGACACCATGCAGACCCGCAACGAGCTGTa                                                                                                    <  1:549338/151‑1 (MQ=255)                                                cGCAGCGCCCATGTCGCAATGGCGCTGTACCCACTTTCAGCGTTCCGCGCCATGAACCGCGCCGCTGAACATGTCTACAACGTCCTGCGCCAGGAAGGCGCGCAGAAAAGCGTCATCGACACCATGCAGACCCGCAACGAGCTGTACGaa                                                                                                >  1:103789/1‑150 (MQ=255)                                                     cgcCCATGTCGCAATGGCGCTGTACCCACTTTCAGCGTTCCGCGCCATGAACCGCGCCGCTGAACATGTCTACAACGTCCTGCGCCAGGAAGGCGCGCAGAAAAGCGTCATCGACACCATGCAGACCCGCAACGAGCTGTACGAAAGCATc                                                                                          >  2:199170/1‑151 (MQ=255)                                                       cccATGTCGCAATGGCGCTGTACCCACTTTCAGCGTTCCGCGCCATGAACCGCGCCGCTGAACATGTCTACAACGTCCTGCGCCAGGAAGGCGCGCAGAAAAGCGTCATCGACACCATGCAGACCCGCAACGAGCTGTACGAAAGCATCaa                                                                                        <  1:479364/151‑1 (MQ=255)                                                        ccATGTCGCAATGGCGCTGTACCCACTTTCAGCGTTCCGCGCCATGAACCGCGCCGCTGAACATGTCTACAACGTCCTGCGCCAGGAAGGCGCGCAGAAAAGCGTCATCGACACCATGCAGACCCGCAACGAGCTGTACGAAAGCATCAac                                                                                       >  2:489500/1‑151 (MQ=255)                                                          aTGTCGCAATGGCGCTGTACCCACTTTCAGCGTTCCGCGCCATGAACCGCGCCGCTGAACATGTCTACAACGTCCTGCGCCAGGAAGGCGCGCAGAAAAGCGTCATCGACACCATGCAGACCCGCAACGAGCTGTACGAAAGCATCAacta                                                                                     <  1:118571/151‑1 (MQ=255)                                                               gCAATGGCGCTGTACCCACTTTCAGCGTTCCGCGCCATGAACCGCGCCGCTGAACATGTCTACAACGTCCTGCGCCAGGAAGGCGCGCAGAAAAGCGTCATCGACACCATGCAGACCCGCAACGAGCTGTACGAAAGCATCAACTACTAcc                                                                                >  1:116118/1‑151 (MQ=255)                                                                 aaTGGCGCTGTACCCACTTTCAGCGTTCCGCGCCATGAACCGCGCCGCTGAACATGTCTACAACGTCCTGCGCCAGGAAGGCGCGCAGAAAAGCGTCATCGACACCATGCAGACCCGCAACGAGCTGTACGAAAGCATCAACTACTACCAg                                                                              <  2:300469/151‑1 (MQ=255)                                                                      cgcTGTACCCACTTTCAGCGTTCCGCGCCATGAACCGCGCCGCTGAACATGTCTACAACGTCCTGCGCCAGGAAGGCGCGCAGAAAAGCGTCATCGACACCATGCAGACCCGCAACGAGCTGTACGAAAGCATCAACTACTACCAGTACg                                                                          >  1:583819/1‑150 (MQ=255)                                                                       gcTGTACCCACTTTCAGCGTTCCGCGCCATGAACCGCGCCGCTGAACATGTCTACAACGTCCTGCGCCAGGAAGGCGCGCAGAAAAGCGTCATCGACACCATGCAGACCCGCAACGAGCTGTACGAAAGCATCAACTACTACCAGTACGaa                                                                        <  1:667811/151‑1 (MQ=255)                                                                              ccACTTTCAGCGTTCCGCGCCATGAACCGCGCCGCTGAACATGTCTACAACGTCCTGCGCCAGGAAGGCGCGCAGAAAAGCGTCATCGACACCATGCAGACCCGCAACGAGCTGTACGAAAGCATCAACTACTACCAGTACGAAGAGAAGc                                                                 >  2:248705/1‑151 (MQ=255)                                                                                       gCGTTCCGCGCCATGAACCGCGCCGCTGAACATGTCTACAACGTCCTGCGCCAGGAAGGCGCGCAGAAAAGCGTCATCGACACCATGCAGACCCGCAACGAGCTGTACGAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAAcc                                                        <  1:501273/151‑1 (MQ=255)                                                                                              gcgcCATGAACCGCGCCGCTGAACATGTCTACAACGTCCTGCGCCAGGAAGGCGCGCAGAAAAGCGTCATCGACACCATGCAGACCCGCAACGAGCTGTACGAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGc                                                 >  2:126958/1‑151 (MQ=255)                                                                                                 ccATGAACCGCGCCGCTGAACATGTCTACAACGTCCTGCGCCAGGAAGGCGCGCAGAAAAGCGTCATCGACACCATGCAGACCCGCAACGAGCTGTACGAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCg                                              <  2:152340/151‑1 (MQ=255)                                                                                                         cgcgCCGCTGAACATGTCTACAACGTCCTGCGCCAGGAAGGCGCGCAGAAAAGCGTCATCGACACCATGCAGACCCGCAACGAGCTGTACGAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAgg                                      <  1:457807/151‑1 (MQ=255)                                                                                                           cgccgcTGAACATGTCTACAACGTCCTGCGCCAGGAAGGCGCGCAGAAAAGCGTCATCGACACCATGCAGACCCGCAACGAGCTGTACGAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTg                                    <  1:390563/151‑1 (MQ=255)                                                                                                              cgcTGAACATGTCTACAACGTCCTGCGCCAGGAAGGCGCGCAGAAAAGCGTCATCGACACCATGCAGACCCGCAACGAGCTGTACGAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGaa                                  >  2:235194/1‑150 (MQ=255)                                                                                                                    aCATGTCTACAACGTCCTGCGCCAGGAAGGCGCGCAGAAAAGCGTCATCGACACCATGCAGACCCGCAACGAGCTGTACGAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATaaaaa                           <  2:257752/151‑1 (MQ=255)                                                                                                                    aCATGTCTACAACGTCCTGCGCCAGGAAGGCGCGCAGAAAAGCGTCATCGACACCATGCAGACCCGCAACGAGCTGTACGAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATaaaa                            >  2:195100/1‑150 (MQ=255)                                                                                                                        gTCTACAACGTCCTGCGCCAGGAAGGCGCGCAGAAAAGCGTCATCGACACCATGCAGACCCGCAACGAGCTGTACGAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGc                       <  2:109106/151‑1 (MQ=255)                                                                                                                            aCAACGTCCTGCGCCAGGAAGGCGCGCAGAAAAGCGTCATCGACACCATGCAGACCCGCAACGAGCTGTACGAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCg                    <  2:441443/150‑1 (MQ=255)                                                                                                                                 gTCCTGCGCCAGGAAGGCGCGCAGAAAAGCGTCATCGACACCATGCAGACCCGCAACGAGCTGTACGAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGAtt              <  2:10068/151‑1 (MQ=255)                                                                                                                                        gcCAGGAAGGCGCGCAGAAAAGCGTCATCGACACCATGCAGACCCGCAACGAGCTGTACGAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCg       <  1:127142/151‑1 (MQ=255)                                                                                                                                        gcCAGGAAGGCGCGCAGAAAAGCGTCATCGACACCATGCAGACCCGCAACGAGCTGTACGAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCg       >  2:127145/1‑151 (MQ=255)                                                                                                                                         ccAGGAAGGCGCGCAGAAAAGCGTCATCGACACCATGCAGACCCGCAACGAGCTGTACGAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCg       >  1:300956/1‑150 (MQ=255)                                                                                                                                             gAAGGCGCGCAGAAAAGCGTCATCGACACCATGCAGACCCGCAACGAGCTGTACGAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGc  <  2:661389/151‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                 ACATCACCGAATTTGGTGCCACGCCGCTGTTTACCACCGACGAATTACGCAGCGCCCATGTCGCAATGGCGCTGTACCCACTTTCAGCGTTCCGCGCCATGAACCGCGCCGCTGAACATGTCTACAACGTCCTGCGCCAGGAAGGCACGCAGAAAAGCGTCATCGACACCATGCAGACCCGCAACGAGCTGTACGAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGC  >  NC_000913/349310‑349600

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