Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,462,065	T→G	G391G (GGT→GGG)	paaJ →	3‑oxoadipyl‑CoA/3‑oxo‑5,6‑dehydrosuberyl‑CoA thiolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,462,065	0	T	G	100.0%	97.6 / NA	28	G391G (GGT→GGG)	paaJ	3‑oxoadipyl‑CoA/3‑oxo‑5,6‑dehydrosuberyl‑CoA thiolase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (9/19);  total (9/19)

GAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTACCCTACAATGATAACCAATACAAAGCTTGACCCG  >  NC_000913/1461975‑1462151                                                                                           |                                                                                       gAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGGCAGGGCATc                                                                               <  2:1271896/100‑1 (MQ=255)          ccGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGGCAGGGCATCGCCATGAtt                                                                      >  1:371234/1‑100 (MQ=255)           cGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGGCAGGGCATCGCCATGATTc                                                                     <  2:596104/100‑1 (MQ=255)                ggCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGGCAGGGCATCGCCATGATTCTGGAg                                                                <  1:236219/100‑1 (MQ=255)                  cACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGGCAGGGCATCGCCATGATTCTGGAGCg                                                              <  2:74937/100‑1 (MQ=255)                    cTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGGCAGGGCATCGCCATGATTCTGGAGCgtg                                                            <  1:96316/100‑1 (MQ=255)                        cTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGGCAGGGCATCGCCATGATTCTGGAGCGTGTTTg                                                        >  2:707383/1‑100 (MQ=255)                        cTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGGCAGGGCATCGCCATGATTCTGGAGCGTGTTTg                                                        <  2:843242/100‑1 (MQ=255)                          gccagccaTGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGGCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAg                                                      >  2:1257754/1‑100 (MQ=255)                            cagccaTGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGGCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCa                                                    <  2:371231/100‑1 (MQ=255)                             agccaTGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGGCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCat                                                   <  2:213567/100‑1 (MQ=255)                                     gCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGGCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAAcc                                            <  1:450963/99‑1 (MQ=255)                                       cttccgatcttAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGGCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGc                                         <  2:847032/90‑1 (MQ=255)                                        gCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGGCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCg                                        >  1:276880/1‑100 (MQ=255)                                              gCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGGCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTAcc                                  <  1:546471/100‑1 (MQ=255)                                                 tAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGGCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGcagatcggaag                               >  1:847038/1‑90 (MQ=255)                                                 tAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGGCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTACCCTa                               >  1:1149426/1‑100 (MQ=255)                                                 tAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGGCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTACCCTa                               <  2:1004571/100‑1 (MQ=255)                                                      gTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGGCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTACCCTACAATg                          <  1:771862/100‑1 (MQ=255)                                                      gTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGGCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTACCCTACAATg                          <  1:771856/100‑1 (MQ=255)                                                       tCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGGCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTACCCTACAATGa                         <  2:876371/100‑1 (MQ=255)                                                           tACGCATTGTGCACCATGTGCATCGGTGTCGGGCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTACCCTACAATGATAAc                     <  2:201970/100‑1 (MQ=255)                                                            aCGCATTGTGCACCATGTGCATCGGTGTCGGGCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTACCCTACAATGATAAcc                    <  2:1211681/100‑1 (MQ=255)                                                               cATTGTGCACCATGTGCATCGGTGTCGGGCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTACCCTACAATGATAACCAAt                 <  2:52239/100‑1 (MQ=255)                                                                   gtgCACCATGTGCATCGGTGTCGGGCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTACCCTACAATGATAACCAATACaa             >  1:1216090/1‑100 (MQ=255)                                                                     gCACCATGTGCATCGGTGTCGGGCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTACCCTACAATGATAACCAATACAAAg           <  1:300355/100‑1 (MQ=255)                                                                             tgCATCGGTGTCGGGCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTACCCTACAATGATAACCAATACAAAGCTTGAccc   >  1:1206778/1‑100 (MQ=255)                                                                              gCATCGGTGTCGGGCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTACCCTACAATGATAACCAATACAAAGCTTGACCCg  >  1:948884/1‑100 (MQ=255)                                                                                           |                                                                                       GAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTACCCTACAATGATAACCAATACAAAGCTTGACCCG  >  NC_000913/1461975‑1462151

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑