Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,358,354	A→G	T88T (ACT→ACC)	cadA ←	lysine decarboxylase, acid‑inducible

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,358,354	0	A	G	100.0%	177.4 / NA	50	T88T (ACT→ACC)	cadA	lysine decarboxylase, acid‑inducible
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (30/20);  total (30/20)

AACAGTGCTTTAGTCAGCGGAGGCAGAATAGTGTTGATATATTCGTCAGTGGTCTGCTTGATCTTATTAGCAATATCTTCAGCAGCACCCAGCGCATATTCAAAGAAGCTAATCTGTAAACGCAGGTCATTCAGGCTTACATCGAGAGTGGAATACGTATTAGCGAACGCGTACAACGGCAGGTTCTCGTTCATTTTGCTAATTTCTTCGCACAGCTCGAGATTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGATCAGTTTTAATAAGTCG  >  NC_000913/4358208‑4358498                                                                                                                                                   |                                                                                                                                                 aaCAGTGCTTTAGTCAGCGGAGGCAGAATAGTGTTGATATATTCGTCAGTGGTCTGCTTGATCTTATTAGCAATATCTTCAGCAGCACCCAGCGCATATTCAAAGAAGCTAATCTGTAAACGCAGGTCATTCAGGCTTACATCGAGGGTgg                                                                                                                                              <  1:822329/151‑1 (MQ=255) aaCAGTGCTTTAGTCAGCGGAGGCAGAATAGTGTTGATATATTCGTCAGTGGTCTGCTTGATCTTATTAGCAATATCTTCAGCAGCACCCAGCGCATATTCAAAGAAGCTAATCTGTAAACGCAGGTCATTCAGGCTTACATCGAGGGTgg                                                                                                                                              >  1:789314/1‑151 (MQ=255)   cAGTGCTTTAGTCAGCGGAGGCAGAATAGTGTTGATATATTCGTCAGTGGTCTGCTTGATCTTATTAGCAATATCTTCAGCAGCACCCAGCGCATATTCAAAGAAGCTAATCTGTAAACGCAGGTCATTCAGGCTTACATCGAGGGTGGa                                                                                                                                             >  1:792846/1‑150 (MQ=255)        cTTTAGTCAGCGGAGGCAGAATAGTGTTGATATATTCGTCAGTGGTCTGCTTGATCTTATTAGCAATATCTTCAGCAGCACCCAGCGCATATTCAAAGAAGCTAATCTGTAAACGCAGGTCATTCAGGCTTACATCGAGGGTGGAATACGt                                                                                                                                       >  2:139352/1‑151 (MQ=255)            aGTCAGCGGAGGCAGAATAGTGTTGATATATTCGTCAGTGGTCTGCTTGATCTTATTAGCAATATCTTCAGCAGCACCCAGCGCATATTCAAAGAAGCTAATCTGTAAACGCAGGTCATTCAGGCTTACATCGAGGGTGGAATACGTAtt                                                                                                                                    >  2:802267/1‑150 (MQ=255)             gTCAGCGGAGGCAGAATAGTGTTGATATATTCGTCAGTGGTCTGCTTGAGCTTATTAGCAATATCTTCAGCAGCACCCAGCGCATATTCAAAGAAGCTAATCTGTAAACGCAGGTCATTCAGGCTTACATCGAGGGTGGAATACGTATTAg                                                                                                                                  <  1:668252/151‑1 (MQ=255)                aGCGGAGGCAGAATAGTGTTGATATATTCGTCAGTGGTCTGCTTGATCTTATTAGCAATATCTTCAGCAGCACCCAGCGCATATTCAAAGAAGCTAATCTGTAAACGCAGGTCATTCAGGCTTACATCGAGGGTGGAATACGTATTAGc                                                                                                                                 <  1:509287/149‑1 (MQ=255)                   ggAGGCAGAATAGTGTTGATATATTCGTCAGTGGTCTGCTTGATCTTATTAGCAATATCTTCAGCAGCACCCAGCGCATATTCAAAGAAGCTAATCTGTAAACGCAGGTCATTCAGGCTTACATCGAGGGTGGAATACGTATTAGCGAAcg                                                                                                                            >  1:454885/1‑151 (MQ=255)                       gCAGAATAGTGTTGATATATTCGTCAGTGGTCTGCTTGATCTTATTAGCAATATCTTCAGCAGCACCCAGCGCATATTCAAAGAAGCTAATCTGTAAACGCAGGTCATTCAGGCTTACATCGAGGGTGGAATACGTATTAGCGAACGCGt                                                                                                                         >  1:142402/1‑150 (MQ=255)                            aTAGTGTTGATATATTCGTCAGTGGTCTGCTTGATCTTATTAGCAATATCTTCAGCAGCACCCAGCGCATATTCAAAGAAGCTAATCTGTAAACGCAGGTCATTCAGGCTTACATCGAGGGTGGAATACGTATTAGCGAACGCGTACAACg                                                                                                                   >  1:728052/1‑151 (MQ=255)                              aGTGTTGATATATTCGTCAGTGGTCTGCTTGATCTTATTAGCAATATCTTCAGCAGCACCCAGCGCATATTCAAAGAAGCTAATCTGTAAACGCAGGTCATTCAGGCTTACATCGAGGGTGGAATACGTATTAGCGAACGCGTACAACGGc                                                                                                                 <  2:328207/151‑1 (MQ=255)                                   tGATATATTCGTCAGTGGTCTGCTTGATCTTATTAGCAATATCTTCAGCAGCACCCAGCGCATATTCAAAGAAGCTAATCTGTAAACGCAGGTCATTCAGGCTTACATCGAGGGTGGAATACGTATTAGCGAACGCGTACAACGGCAGGtt                                                                                                            >  2:222993/1‑151 (MQ=255)                                     atatatTCGTCAGTGGTCTGCTTGATCTTATTAGCAATATCTTCAGCAGCACCCAGCGCATATTCAAAGAAGCTAATCTGTAAACGCAGGTCATTCAGGCTTACATCGAGGGTGGAATACGTATTAGCGAACGCGTACAACGGCAGGTtct                                                                                                          >  2:365331/1‑151 (MQ=255)                                      tatatTCGTCAGTGGTCTGCTTGATCTTATTAGCAATATCTTCAGCAGCACCCAGCGCATATTCAAAGAAGCTAATCTGTAAACGCAGGTCATTCAGGCTTACATCGAGGGTGGAATACGTATTAGCGAACGCGTACAACGGCAGGTtctc                                                                                                         >  2:371900/1‑151 (MQ=255)                                         atTCGTCAGTGGTCTGCTTGATCTTATTAGCAATATCTTCAGCAGCACCCAGCGCATATTCAAAGAAGCTAATCTGTAAACGCAGGTCATTCAGGCTTACATCGAGGGTGGAATACGTATTAGCGAACGCGTACAACGGCAGGTTCTCGtt                                                                                                      <  1:619114/151‑1 (MQ=255)                                         atTCGTCAGTGGTCTGCTTGATCTTATTAGCAATATCTTCAGCAGCACCCAGCGCATATTCAAAGAAGCTAATCTGTAAACGCAGGTCATTCAGGCTTACATCGAGGGTGGAATACGTATTAGCGAACGCGTACAACGGCAGGTTCTCGtt                                                                                                      >  2:711201/1‑151 (MQ=255)                                             gTCAGTGGTCTGCTTGATCTTATTAGCAATATCTTCAGCAGCACCCAGCGCATATTCAAAGAAGCTAATCTGTAAACGCAGGTCATTCAGGCTTACATCGAGGGTGGAATACGTATTAGCGAACGCGTACAACGGCAGGTTCTCGTTCAtt                                                                                                  >  1:636161/1‑151 (MQ=255)                                               cAGTGGTCTGCTTGATCTTATTAGCAATATCTTCAGCAGCACCCAGCGCATATTCAAAGAAGCTAATCTGTAAACGCAGGTCATTCAGGCTTACATCGAGGGTGGAATACGTATTAGCGAACGCGTACAACGGCAGGTTCTCGTTCAtttt                                                                                                >  1:819630/1‑151 (MQ=255)                                                aGTGGTCTGCTTGATCTTATTAGCAATATCTTCAGCAGCACCCAGCGCATATTCAAAGAAGCTAATCTGTAAACGCAGGTCATTCAGGCTTACATCGAGGGTGGAATACGTATTAGCGAACGCGTACAACGGCAGGTTCTCGTTCATTTTg                                                                                               <  2:541146/151‑1 (MQ=255)                                                             aTCTTATTAGCAATATCTTCAGCAGCACCCAGCGCATATTCAAAGAAGCTAATCTGTAAACGCAGGTCATTCAGGCTTACATCGAGGGTGGAATACGTATTAGCGAACGCGTACAACGGCAGGTTCTCGTTCATTTTGCTAATTTCTTCGc                                                                                  >  2:222217/1‑151 (MQ=255)                                                                ttattaGCAATATCTTCAGCAGCACCCAGCGCATATTCAAAGAAGCTAATCTGTAAACGCAGGTCATTCAGGCTTACATCGAGGGTGGAATACGTATTAGCGAACGCGTACAACGGCAGGTTCTCGTTCATTTTGCTAATTTCTTCGcac                                                                                >  2:832154/1‑150 (MQ=255)                                                                 tattaGCAATATCTTCAGCAGCACCCAGCGCATATTCAAAGAAGCTAATCTGTAAACGCAGGTCATTCAGGCTTACATCGAGGGTGGAATACGTATTAGCGAACGCGTACAACGGCAGGTTCTCGTTCATTTTGCTAATTTCTTCg                                                                                   <  2:9157/146‑1 (MQ=255)                                                                 tattaGCAATATCTTCAGCAGCACCCAGCGCATATTCAAAGAAGCTAATCTGTAAACGCAGGTCATTCAGGCTTACATCGAGGGTGGAATACGTATTAGCGAACGCGTACAACGGCAGGTTCTCGTTCATTTTGCTAATTTCTTCGcac                                                                                >  2:430867/1‑149 (MQ=255)                                                                 tattaGCAATATCTTCAGCAGCACCCAGCGCATATTCAAAGAAGCTAATCTGTAAACGCAGGTCATTCAGGCTTACATCGAGGGTGGAATACGTATTAGCGAACGCGTACAACGGCAGGTTCTCGTTCATTTTGCTAATTTCTTCGcac                                                                                <  2:292928/149‑1 (MQ=255)                                                                  attaGCAATATCTTCAGCAGCACCCAGCGCATATTCAAAGAAGCTAATCTGTAAACGCAGGTCATTCAGGCTTACATCGAGGGTGGAATACGTATTAGCGAACGCGTACAACGGCAGGTTCTCGTTCATTTTGCTAATTTCTTCGCACAg                                                                              <  2:180957/150‑1 (MQ=255)                                                                         atatCTTCAGCAGCACCCAGCGCATATTCAAAGAAGCTAATCTGTAAACGCAGGTCATTCAGGCTTACATCGAGGGTGGAATACGTATTAGCGAACGCGTACAACGGCAGGTTCTCGTTCATTTTGCTAATTTCTTCGCACAGCTCg                                                                          >  1:166876/1‑147 (MQ=255)                                                                              ttCAGCAGCACCCAGCGCATATTCAAAGAAGCTAATCTGTAAACGCAGGTCATTCAGGCTTACATCGAGGGTGGAATACGTATTAGCGAACGCGTACAACGGCAGGTtctc                                                                                                         >  2:784280/1‑111 (MQ=255)                                                                              ttCAGCAGCACCCAGCGCATATTCAAAGAAGCTAATCTGTAAACGCAGGTCATTCAGGCTTACATCGAGGGTGGAATACGTATTAGCGAACGCGTACAACGGCAGGTtctc                                                                                                         <  1:784224/111‑1 (MQ=255)                                                                              ttCAGCAGCACCCAGCGCATATTCAAAGAAGCTAATCTGTAAACGCAGGTCATTCAGGCTTACATCGAGGGTGGAATACGTATTAGCGAACGCGTACAACGGCAGGTTCTCGTTCATTTTGCTAATTTCTTCGCACAGCTCGAGATTatat                                                                 >  1:108390/1‑151 (MQ=255)                                                                               tCAGCAGCACCCAGCGCATATTCAAAGAAGCTAATCTGTAAACGCAGGTCATTCAGGCTTACATCGAGGGTGGAATACGTATTAGCGAACGCGTACAACGGCAGGTTCTCGTTCATTTTGCTAATTTCTTCGCACAGCTCGAGATtat                                                                   <  1:693434/148‑1 (MQ=255)                                                                                cagcagCACCCAGCGCATATTCAAAGAAGCTAATCTGTAAACGCAGGTCATTCAGGCTTACATCGAGGGTGGAATACGTATTAGCGAACGCGTACAACGGCAGGTTCTCGTTCATTTTGCTAATTTCTTCGCACAGCTCGAGATTATAttt                                                               <  2:837209/151‑1 (MQ=255)                                                                                cagcagCACCCAGCGCATATTCAAAGAAGCTAATCTGTAAACGCAGGTCATTCAGGCTTACATCGAGGGTGGAATACGTATTAGCGAACGCGTACAACGGCAGGTTCTCGTTCATTTTGCTAATTTCTTCGCACAGCTCGAGATTATAtt                                                                <  1:542415/150‑1 (MQ=255)                                                                                 agcagcACCCAGCGCATATTCAAAGAAGCTAATCTGTAAACGCAGGTCATTCAGGCTTACATCGAGGGTGGAATACGTATTAGCGAACGCGTACAACGGCAGGTTCTCGTTCATTTTGCTAATTTCTTCGCACAGCTCGAGATTATAttt                                                               >  1:550078/1‑150 (MQ=255)                                                                                         ccAGCGCATATTCAAAGAAGCTAATCTGTAAACGCAGGTCATTCAGGCTTACATCGAGGGTGGAATACGTATTAGCGAACGCGTACAACGGCAGGTTCTCGTTCATTTTGCTAATTTCTTCGCACAGCTCGAGATTATATTTATCCCAGTc                                                      >  2:92809/1‑151 (MQ=255)                                                                                             cgcATATTCAAAGAAGCTAATCTGTAAACGCAGGTCATTCAGGCTTACATCGAGGGTGGAATACGTATTAGCGAACGCGTACAACGGCAGGTTCTCGTTCATTTTGCTAATTTCTTCGCACAGCTCGAGATTATATTTATCCCAGTCaaaa                                                  <  2:821279/151‑1 (MQ=255)                                                                                             cgcATATTCAAAGAAGCTAATCTGTAAACGCAGGTCATTCAGGCTTACATCGAGGGTGGAATACGTATTAGCGAACGCGTACAACGGCAGGTTCTCGTTCATTTTGCTAATTTCTTCGCACAGCTCGAGATTATATTTATCCCAGTCaaaa                                                  <  1:371876/151‑1 (MQ=255)                                                                                               cATATTCAAAGAAGCTAATCTGTAAACGCAGGTCATTCAGGCTTACATCGAGGGTGGAATACGTATTAGCGAACGCGTACAACGGCAGGTTCTCGTTCATTTTGCTAATTTCTTCGCACAGCTCGAGATTATATTTATCCCAGTCAAAAAt                                                >  1:395172/1‑151 (MQ=255)                                                                                                          aagCTAATCTGTAAACGCAGGTCATTCAGGCTTACATCGAGGGTGGAATACGTATTAGCGAACGCGTACAACGGCAGGTTCTCGTTCATTTTGCTAATTTCTTCGCACAGCTCGAGATTATATTTATCCCAGTCAAAAATAACGCCGcac                                      >  2:729248/1‑150 (MQ=255)                                                                                                            gCTAATCTGTAAACGCAGGTCATTCAGGCTTACATCGAGGGTGGAATACGTATTAGCGAACGCGTACAACGGCAGGTTCTCGTTCATTTTGCTAATTTCTTCGCACAGCTCGAGATTATATTTATCCCAGTCAAAAATAACGCCGCACAGa                                   <  1:274161/151‑1 (MQ=255)                                                                                                               aaTCTGTAAACGCAGGTCATTCAGGCTTACATCGAGGGTGGAATACGTATTAGCGAACGCGTACAACGGCAGGTTCTCGTTCATTTTGCTAATTTCTTCGCACAGCTCGAGATTATATTTATCCCAGTCAAAAATAACGCCGCACAGAcgc                                >  2:506205/1‑151 (MQ=255)                                                                                                               aaTCTGTAAACGCAGGTCATTCAGGCTTACATCGAGGGTGGAATACGTATTAGCGAACGCGTACAACGGCAGGTTCTCGTTCATTTTGCTAATTTCTTCGCACAGCTCGAGATTATATTTATCCCAGTCAAAAATAACGCCGCACAGAc                                  <  2:206233/149‑1 (MQ=255)                                                                                                                aTCTGTAAACGCAGGTCATTCAGGCTTACATCGAGGGTGGAATACGTATTAGCGAACGCGTACAACGGCAGGTTCTCGTTCATTTTGCTAATTTCTTCGCACAGCTCGAGATTATATTTATCCCAGTCAAAAATAACGCCGCACAGAcgc                                >  1:526154/1‑150 (MQ=255)                                                                                                                  cTGTAAACGCAGGTCATTCAGGCTTACATCGAGGGTGGAATACGTATTAGCGAACGCGTACAACGGCAGGTTCTCGTTCATTTTGCTAATTTCTTCGCACAGCTCGAGATTATATTTATCCCAGTCAAAAATAACGCCGCACAGACgcgc                              >  2:549726/1‑150 (MQ=255)                                                                                                                   tGTAAACGCAGGTCATTCAGGCTTACATCGAGGGTGGAATACGTATTAGCGAACGCGTACAACGGCAGGTTCTCGTTCATTTTGCTAATTTCTTCGCACAGCTCga                                                                         >  2:330171/1‑106 (MQ=255)                                                                                                                   tGTAAACGCAGGTCATTCAGGCTTACATCGAGGGTGGAATACGTATTAGCGAACGCGTACAACGGCAGGTTCTCGTTCATTTTGCTAATTTCTTCGCACAGCTCga                                                                         <  1:330148/106‑1 (MQ=255)                                                                                                                                 aTTCAGGCTTACATCGAGGGTGGAATACGTATTAGCGAACGCGTACAACGGCAGGTTCTCGTTCATTTTGCTAATTTCTTCGCACAGCTCGAGATTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGATCAGt              <  2:417209/151‑1 (MQ=255)                                                                                                                                 aTTCAGGCTTACATCGAGGGTGGAATACGTATTAGCGAACGCGTACAACGGCAGGTTCTCGTTCATTTTGCTAATTTCTTCGCACAGCTCGAGATTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGATCAGt              >  1:1271/1‑151 (MQ=255)                                                                                                                                    cAGGCTTACATCGAGGGTGGAATACGTATTAGCGAACGCGTACAACGGCAGGTTCTCGTTCATTTTGCTAATTTCTTCGCACATCTCGAGATTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGATCAGtttt           >  2:838314/1‑151 (MQ=255)                                                                                                                                            cATCGAGGGTGGAATACGTATTAGCGAACGCGTACAACGGCAGGTTCTCGTTCATTTTGCTAATTTCTTCGCACAGCTCGAGATTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGATCAGTTTTAATAAgt    <  2:677385/150‑1 (MQ=255)                                                                                                                                             aTCGAGGGTGGAATACGTATTAGCGAACGCGTACAACGGCAGGTTCTCGTTCATTTTGCTAATTTCTTCGCACAGCTCGAGATTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGATCAGTTTTAATAAgtcg  >  1:859129/1‑151 (MQ=255)                                                                                                                                                   |                                                                                                                                                 AACAGTGCTTTAGTCAGCGGAGGCAGAATAGTGTTGATATATTCGTCAGTGGTCTGCTTGATCTTATTAGCAATATCTTCAGCAGCACCCAGCGCATATTCAAAGAAGCTAATCTGTAAACGCAGGTCATTCAGGCTTACATCGAGAGTGGAATACGTATTAGCGAACGCGTACAACGGCAGGTTCTCGTTCATTTTGCTAATTTCTTCGCACAGCTCGAGATTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGATCAGTTTTAATAAGTCG  >  NC_000913/4358208‑4358498

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