Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,835,944	G→A	A144T (GCT→ACT)	ynjA →	carboxymuconolactone decarboxylase family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,835,944	0	G	A	100.0%	145.3 / NA	46	A144T (GCT→ACT)	ynjA	carboxymuconolactone decarboxylase family protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (22/24);  total (22/24)

GCTGACTGGCGGCAAAGCCCGCTCTTTAGCGATGAAGAACGGCTGGCGCTGGAGTACGCCGAAGCCGCAAGCGTAACGCCGCCAACGGTCGATGATGCCCTGCGTACCCGACTGGCTGCGCATTTTGACGCTCAGGCGCTCACCGAACTGACGGCATTGATCGGCCTGCAAAATCTGTCAGCCCGTTTTAATTCTGCCATGGACATTCCCGCTCAGGGGCTGTGCCGTATTCCTGAAAAACGTTCTTAAGGAGAGA  >  NC_000913/1835815‑1836070                                                                                                                                  |                                                                                                                               gCTGACTGGCGGCAAAGCCCGCTCTTTAGCGATGAAGAACGGCTGGCGCTGGAGTACGCCGAAGCCGCAAGCGTAACGCCGCCAACGGTCGATGATGCCCTGCGTACCCGACTGGCTGCGCATTTTGACACTCAGgcgc                                                                                                                       <  1:1172606/139‑1 (MQ=255)          cggcAAAGCCCGCTCTTTAGCGATGAAGAACGGCTGGCGCTGGAGTACGCCGAAGCCGCAAGCGTAACGCCGCCAACGGTCGATGATGCCCTGCGTACCCGACTGGCTGCGCATTTTGACACTCAgg                                                                                                                          <  2:281625/127‑1 (MQ=255)          cggcAAAGCCCGCTCTTTAGCGATGAAGAACGGCTGGCGCTGGAGTACGCCGAAGCCGCAAGCGTAACGCCGCCAACGGTCGATGATGCCCTGCGTACCCGACTGGCTGCGCATTTTGACACTCAgg                                                                                                                          >  1:281625/1‑127 (MQ=255)           ggcAAAGCCCGCTCTTTAGCGATGAAGAACGGCTGGCGCTGGAGTACGCCGAAGCCGCAAGCGTAACGCCGCCAACGGTCGATGATGCCCTGCGTACCCGACTGGCTGCGCATTTTGACACTCAGGCGCTCACCGTACt                                                                                                             >  2:1181629/1‑139 (MQ=255)           ggcAAAGCCCGCTCTTTAGCGATGAAGAACGGCTGGCGCTGGAGTACGCCGAAGCCGCAAGCGTAACGCCGCCAACGGTCGATGATGCCCTGCGTACCCGACTGGCTGCGCATTTTGACACTCAGGCGCTCACCGAACt                                                                                                             >  2:1602565/1‑139 (MQ=255)             cAAAGCCCGCTCTTTAGCGATGAAGAACGGCTGGCCCTGGAGTACGCCGAAGCCGCAAGCGTAACGCCGCCAACGGTCGATGATGCCCTGCGTACCCGACTGGCTGCGCATTTTGACACTCAGGCGCTCACCGAACTGa                                                                                                           <  2:431736/139‑1 (MQ=255)                       tctTTAGCGATGAAGAACGGCTGGCGCTGGAGTACGCCGAAGCCGCAAGCGTAACGCCGCCAACGGTCGATGATGCCCTGCGTACCCGACTGGCTGCGCATTTTGACACTCAGGCGCTCACCGAACTGACGGCATTGAt                                                                                                 >  2:101574/1‑139 (MQ=255)                        ctTTAGCGATGAAGAACGGCTGGCGCTGGAGTACGCCGAAGCCGCAAGCGTAACGCCGCCAACGGTCGATGATGCCCTGCGTACCCGACTGGCTGCGCATTTTGACACTCAGGCGCTCACCGAACTGACGGCATTGATc                                                                                                >  1:301230/1‑139 (MQ=255)                          ttAGCGATGAAGAACGGCTGGCGCTGGAGTACGCCGAAGCCGCAAGCGTAACGCCGCCAACGGTCGATGATGCCCTGCGTACCCGACTGGCTGCGCATTTTGACACTCAGGCGCTCACCGAACTGACGGCATTGATCgg                                                                                              <  1:1262265/139‑1 (MQ=255)                                aTGAAGAACGGCTGGCGCTGGAGTACGCCGAAGCCGCAAGCGTAACGCCGCCAACGGTCGATGATGCCCTGCGTACCCGACTGGCTGCGCATTTTGACACTCAGGCGCTCACCGAACTGa                                                                                                           <  1:215197/120‑1 (MQ=255)                                aTGAAGAACGGCTGGCGCTGGAGTACGCCGAAGCCGCAAGCGTAACGCCGCCAACGGTCGATGATGCCCTGCGTACCCGACTGGCTGCGCATTTTGACACTCAGGCGCTCACCGAACTGa                                                                                                           >  2:215197/1‑120 (MQ=255)                                 tGAAGAACGGCTGGCGCTGGAGTACGCCGAAGCCGCAAGCGTAACGCCGCCAACGGTCGATGATGCCCTGCGTACCCGACTGGCTGCGCATTTTGACACTCAGGCGCTCACCGAACTGACGGCATTGATCGGCCTGCaa                                                                                       >  1:1714369/1‑139 (MQ=255)                                   aagaaCGGCTGGCGCTGGAGTACGCCGAAGCCGCAAGCGTAACGCCGCCAACGGTCGATGATGCCCTGCGTACCCGACTGGCTGCGCATTTTGACACTCAGGCGCTCACCGAACTg                                                                                                            >  2:1720087/1‑116 (MQ=255)                                   aagaaCGGCTGGCGCTGGAGTACGCCGAAGCCGCAAGCGTAACGCCGCCAACGGTCGATGATGCCCTGCGTACCCGACTGGCTGCGCATTTTGACACTCAGGCGCTCACCGAACTg                                                                                                            <  1:1720087/116‑1 (MQ=255)                                   aagaaCGGCTGGCGCTGGAGTACGCCGAAGCCGCAAGCGTAACGCCGCCAACGGTCGATGATGCCCTGCGTACCCGACTGGCTGCGCATTTTGACACTCAGGCGCTCACCGAACTGACGGCATTGATCGGCCTGCaaaa                                                                                     >  2:1365036/1‑139 (MQ=255)                                       aCGGCTGGCGCTGGAGTACGCCGAAGCCGCAAGCGTAACGCCGCCAACGGTCGATGATGCCCTGCGTACCCGACTGGCTGCGCATTTTGACACTCAGGCGCTCACCGAACTGACGGCATTGATCGGCCTGCAAAATCTg                                                                                 >  1:446805/1‑139 (MQ=255)                                             ggCGCTGGAGTACGCCGAAGCCGCAAGCGTAACGCCGCCAACGGTCGATGATGCCCTGCGTACCCGACTGGCTGCGCATTTTGACACTCAGGCGTTCACCGAACTGACGGCATTGATCGGCCTGCAAAATCTGTCAGcc                                                                           <  2:323266/139‑1 (MQ=255)                                                gcTGGAGTACGCCGAAGCCGCAAGCGTAACGCCGCCAACGGTCGATGATGCCCTGCGTACCCGACTGGCTGCGCATTTTGACACt                                                                                                                              <  2:1133419/85‑1 (MQ=255)                                                     aGTACGCCGAAGCCGCAAGCGTAACGCCGCCAACGGTCGATGATGCCCTGCGTACCCGACTGGCTGCGCATTTTGACACTCAGGCGCTCACCGAACTGACGGCATTGATCGGCCTGCAAAATCTGTCAGCCCGTTTTaa                                                                   >  2:158700/1‑139 (MQ=255)                                                        aCGCCGAAGCCGCAAGCGTAACGCCGCCAACGGTCGATGATGCCCTGCGTACCCGACTGGCTGCGCATTTTGACACTCAGGCGCTCACCGAACTGACGGCATTGATCGGCCTGCAAAATCTGTCAGCCCGTTTTAATTc                                                                >  1:790595/1‑139 (MQ=255)                                                                 ccGCAAGCGTAACGCCGCCAACGGTCGATGATGCCCTGCGTACCCGACTGGCTGCGCATTTTGACACTCAGGCGCTCACCGAACTGACGGCATTGATCGGCCTGCAAAATCTGTCAGCCCGTTTTAATTCTGCCATGGa                                                       <  2:1729871/139‑1 (MQ=255)                                                                      aGCGTAACGCCGCCAACGGTCGATGATGCCCTGCGTACCCGACTGGCTGCGCATTTTGACACTCAGGCGCTCACCGAACTGACGGCATTGATCGGCCTGCAAAATCTGTCAGCCCGTTTTAATTCTGCCATGGACATTc                                                  <  1:995086/139‑1 (MQ=255)                                                                       gCGTAACGCCGCCAACGGTCGATGATGCCCTGCGTACCCGACTGGCTGCGCATTTTGACACTCAGGCGCTCACCGAACTGACGGCATTGATCGGCCTGCAAAATCTGTCAGCCCGTTTTAATTCTGCCATGGACATTcc                                                 <  1:1365036/139‑1 (MQ=255)                                                                                  ccAACGGTCGATGATGCCCTGCGTACCCGACTGGCTGCGCATTTTGACACTCAGGCGCTCACCGAACTGACGGCATTGATCGGCCTGCAAAATCTGTCAGcc                                                                           <  1:1445810/102‑1 (MQ=255)                                                                                  ccAACGGTCGATGATGCCCTGCGTACCCGACTGGCTGCGCATTTTGACACTCAGGCGCTCACCGAACTGACGGCATTGATCGGCCTGCAAAATCTGTCAGcc                                                                           >  2:1445810/1‑102 (MQ=255)                                                                                  ccAACGGTCGATGATGCCCTGCGTACCCGACTGGCTGCGCATTTTGACACTCAGGCGCTCACCGAACTGACGGCATTGATCGGCCTGCAAAATCTGTCAGCCCGTTTTAATTCTGCCATGGACATTCCCGCTCAGGGGc                                      >  2:921123/1‑139 (MQ=255)                                                                                    aaCGGTCGATGATGCCCTGCGTACCCGACTGGCTGCGCATTTTGACACTCAGGCGCTCACCGAACTGACGGCATTGATCGGCCTGCAAAATCTGTCAGCCCGTTTTAATTCTGCCATGGACATTCCCGCTCAGGGGCtg                                    <  2:790595/139‑1 (MQ=255)                                                                                      cGGTCGATGATGCCCTGCGTACCCGACTGGCTGCGCATTTTGACACTCAGGCGCTCACCGAACTGACGGCATTGATCGGCCTGCAAAATCTGTCAGCCCGTTTTAATTCTGCCATGGACATTCCCGCTCAGGGGCtgtg                                  <  1:921123/139‑1 (MQ=255)                                                                                              gatgCCCTGCGTACCCGACTGGCTGCGCATTTTGACACTCAGGCGCTCACCGAACTGACGGCATTGATCGGCCTGCAAAATCTGTCAGCCCGTTTTAATTCTGCCATGGACATTCCCGCTCAGGGGCTGTGCCGTGTTc                          <  1:1181629/139‑1 (MQ=255)                                                                                                 gCCCTGCGTACCCGACTGGCTGCGCATTTTGACACTCAGGCGCTCACCGAACTGACGGCATTGATCGGCCTGCAAAATCTGTCAGCCCGTTTTAATTCTGCCATGGACATTCCCGCTCAGGGGCTGTGCCGTATTCCTg                       >  1:1412169/1‑139 (MQ=255)                                                                                                   ccTGCGTACCCGACTGGCTGCGCATTTTGACACTCAGGCGCTCACCGAACTGACGGCATTGATCGGCCTGCAAAATCTGTCAGCCCGTTTTAATTCTGCCATGGACATTCCCGCTCAGGGGCTGTGCCGTATTCCTGaa                     >  2:1487611/1‑139 (MQ=255)                                                                                                      gCGTACCCGACTGGCTGCGCATTTTGACACTCAGGCGCTCACCGAACTGACGGCATTGATCGGCCTGCAAAATCTGTCAGCCCGTTTTAATTCTGCCATGGACATTCCCGCTCAGGGGCTGTGCCGTATTCCTGaaaaa                  <  2:446805/139‑1 (MQ=255)                                                                                                       cGTACCCGACTGGCTGCGCATTTTGACACTCAGGCGCTCACCGAACTGACGGCATTGATc                                                                                                >  1:41439/1‑60 (MQ=255)                                                                                                       cGTACCCGACTGGCTGCGCATTTTGACACTCAGGCGCTCACCGAACTGACGGCATTGATc                                                                                                <  2:41439/60‑1 (MQ=255)                                                                                                            ccGACTGGCTGCGCATTTTGACACTCAGGCGCTCACCGAACTGACGGCATTGATCGGCCTGCAAAATCTGTCAGCCCGTTTTAATTCTGCCATGGACATTCCCGCTCAGGGGCTGTGCCGTa                             <  1:869346/122‑1 (MQ=255)                                                                                                            ccGACTGGCTGCGCATTTTGACACTCAGGCGCTCACCGAACTGACGGCATTGATCGGCCTGCAAAATCTGTCAGCCCGTTTTAATTCTGCCATGGACATTCCCGCTCAGGGGCTGTGCCGTa                             >  2:869346/1‑122 (MQ=255)                                                                                                              gACTGGCTGCGCATTTTGACACTCAGGCGCTCACCGAACTGACGGCATTGATCGGCCTGCAAAATCTGTCAGCCCGTTTTAATTCTGCCATGGACATTCCCGCTCAGGGGCTGTGCCGTATTCCTGAAAAACGTTCTTa          <  1:1602565/139‑1 (MQ=255)                                                                                                              gACTGGCTGCGCATTTTGACACTCAGGCGCTCACCGAACTGACGGCATTGATCGGCCTGCAAAATCTGTCAGCCCGTTTTAATTCTGCCATGGACATTCCCGCTCAGGGGCTGTGCCGTATTCCTGAAAAACGTTCTTa          <  1:158700/139‑1 (MQ=255)                                                                                                               aCTGGCTGCGCATTTTGACACTCAGGCGCTCACCGAACTGACGGCATTGATCGGCCTGCAAAATCTGTCAGCCCGTTTTAATTCTGCCATGGACATTCCCGCTCAGGGGCTGTGCCGTATTCCTGAAAAACg                >  1:1468792/1‑132 (MQ=255)                                                                                                               aCTGGCTGCGCATTTTGACACTCAGGCGCTCACCGAACTGACGGCATTGATCGGCCTGCAAAATCTGTCAGCCCGTTTTAATTCTGCCATGGACATTCCCGCTCAGGGGCTGTGCCGTATTCCTGAAAAACg                <  2:1468792/132‑1 (MQ=255)                                                                                                                 tGGCTGCGCATTTTGACACTCAGGCGCTCACCGAACTGACGGCATTGATCGGCCTGCAAAATCTGTCAGCCCGTTTTAATTCTGCCATGGACAtt                                                   <  2:408728/95‑1 (MQ=255)                                                                                                                 tGGCTGCGCATTTTGACACTCAGGCGCTCACCGAACTGACGGCATTGATCGGCCTGCAAAATCTGTCAGCCCGTTTTAATTCTGCCATGGACAtt                                                   >  1:408728/1‑95 (MQ=255)                                                                                                                      gcgcATTTTGACACTCAGGCGCTCACCGAACTGACGGCATTGATCGGCCTGCAAAATCTGTCAGCCCGTTTTAATTCTGCCAt                                                          <  1:122120/83‑1 (MQ=255)                                                                                                                      gcgcATTTTGACACTCAGGCGCTCACCGAACTGACGGCATTGATCGGCCTGCAAAATCTGTCAGCCCGTTTTAATTCTGCCAt                                                          >  2:122120/1‑83 (MQ=255)                                                                                                                      gcgcATTTTGACACTCAGGCGCTCACCGAACTGACGGCATTGATCGGCCTGCAAAATCTGTCAGCCCGTTTTAATTCTGCCATGGACATTCCCGCTCAGGGGCTGTGCCGTATTCCTGAAAAACGTTCTTAAGgagaga  >  1:1456114/1‑139 (MQ=255)                                                                                                                      gcgcATTTTGACACTCAGGCGCTCACCGAACTGACGGCATTGATCGGCCTGCAAAATCTGTCAGCCCGTTTTAATTCTGCCATGGACATTCCCGCTCAGGGGCTGTGCCGTATTCCTGAAAAACGTTCTTAAGgagaga  <  1:943053/139‑1 (MQ=255)                                                                                                                                  |                                                                                                                               GCTGACTGGCGGCAAAGCCCGCTCTTTAGCGATGAAGAACGGCTGGCGCTGGAGTACGCCGAAGCCGCAAGCGTAACGCCGCCAACGGTCGATGATGCCCTGCGTACCCGACTGGCTGCGCATTTTGACGCTCAGGCGCTCACCGAACTGACGGCATTGATCGGCCTGCAAAATCTGTCAGCCCGTTTTAATTCTGCCATGGACATTCCCGCTCAGGGGCTGTGCCGTATTCCTGAAAAACGTTCTTAAGGAGAGA  >  NC_000913/1835815‑1836070

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