Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,606,332	T→G	L78R (CTG→CGG)	lsrF →	3‑hydroxy‑2,4‑pentadione 5‑phosphate thiolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,606,332	0	T	G	100.0%	213.9 / NA	55	L78R (CTG→CGG)	lsrF	3‑hydroxy‑2,4‑pentadione 5‑phosphate thiolase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (27/28);  total (27/28)

GAATGCAGTCACGCTTATCGCGGATATTTAATCCGAAAACGGGTAAAACCGTGATGCTGGCTTTTGACCATGGTTATTTTCAGGGACCGACTACCGGACTTGAACGCATTGATATAAATATCGCCCCGCTGTTTGAACATGCCGATGTATTAATGTGTACGCGCGGCATTTTGCGCAGCGTAGTTCCCCCTGCGACCAATAGGCCGGTGGTACTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTG  >  NC_000913/1606203‑1606468                                                                                                                                  |                                                                                                                                         gAATGCAGTCACGCTTATCGCGGATATTTAATCCGAAAACGGGTAAAACCGTGATGCTGGCTTTTGACCATGGTTATTTTCAGGGACCGACTACCGGACTTGACCGCATTGATATAAATATCGCCCCGCGGTTTGAACAt                                                                                                                                <  2:873687/140‑1 (MQ=255)   aTGCAGTCACGCTTATCGCGGATATTTAATCCGAAAACGGGTAAAACCGTGATGCTGGCTTTTGACCATGGTTATTTTCAGGGACCGACTACCGGACTTGAACGCATTGATATAAATATCGCCCCGCGGTTTGAACATGc                                                                                                                              <  2:1046259/140‑1 (MQ=255)   aTGCAGTCACGCTTATCGCGGATATTTAATCCGAAAACGGGTAAAACCGTGATGCTGGCTTTTGACCATGGTTATTTTCAGGGACCGACTACCGGACTTGAACGCATTGATATAAATATCGCCCCGCGGTTTGAACATGc                                                                                                                              <  1:828510/140‑1 (MQ=255)   aTGCAGTCACGCTTATCGCGGATATTTAATCCGAAAACGGGTAAAACCGTGATGCTGGCTTTTGACCATGGTTATTTTCAGGGACCGACTACCGGACTTGAACGCATTGATATAAATATCGCCCCGCGGTTTGAACATGc                                                                                                                              <  1:563510/140‑1 (MQ=255)   aTGCAGTCACGCTTATCGCGGATATTTAATCCGAAAACGGGTAAAACCGTGATGCTGGCTTTTGACCATGGTTATTTTCAGGGACCGACTACCGGACTTGAACGCATTGATATAAATATCGCCCCGCGGTTTGAACATGc                                                                                                                              <  2:136315/140‑1 (MQ=255)   aTGCAGTCACGCTTATCGCGGATATTTAATCCGAAAACGGGTAAAACCGTGATGCTGGCTTTTGACCATGGTTATTTTCAGGGACCGACTACCGGACTTGAACGCATTAATATAAATATCGCCCCGCGGTTTGAACATGc                                                                                                                              <  1:525068/140‑1 (MQ=255)    tGCAGTCACGCTTATCGCGGATATTTAATCCGAAAACGGGTAAAACCGTGATGCTGGCTTTTGACCATGGTTATTTTCAGGGACCGACTACCGGACTTGAACGCATTGATATAAATATCGCCCCGCGGTTTGAACa                                                                                                                                 <  2:367036/136‑1 (MQ=255)    tGCAGTCACGCTTATCGCGGATATTTAATCCGAAAACGGGTAAAACCGTGATGCTGGCTTTTGACCATGGTTATTTTCAGGGACCGACTACCGGACTTGAACGCATTGATATAAATATCGCCCCGCGGTTTGAACa                                                                                                                                 >  1:367036/1‑136 (MQ=255)    tGCAGTCACGCTTATCGCGGATATTTAATCCGAAAACGGGTAAAACCGTGATGCTGGCTTTTGACCATGGTTATTTTCAGGGACCGACTACCGGACTTGAACGCATTGATATAAATATCGCCCCGCGGTTTGAACATGcc                                                                                                                             >  2:1053648/1‑140 (MQ=255)    tGCAGTCACGCTTATCGCGGATATTTAATCCGAAAACGGGTAAAACCGTGATGCTGGCTTTTGACCATGGTTATTTTCAGGGACCGACTACCGGACTTGAACGCATTGATATAAATATCGCCCCGCGGTTTGAACATGcc                                                                                                                             >  2:193725/1‑140 (MQ=255)    tGCAGTCACGCTTATCGCGGATATTTAATCCGAAAACGGGTAAAACCGTGATGCTGGCTTTTGACCATGGTTATTTTCAGGGACCGACTACCGGACTTGAACGCATTGATATAAATATCGCCCCGCGGTTTGAACATGcc                                                                                                                             >  2:7948/1‑140 (MQ=255)     gCAGTCACGCTTATCGCGGATATTTAATCCGAAAACGGGTAAAACCGTGATGCTGGCTTTTGACCATGGTTATTTTCAGGGACCGACTACCGGACTTGAACGCATTGATATAAATATCGCCCCGCGGTTTGAACATGCCg                                                                                                                            <  1:193725/140‑1 (MQ=255)     gCAGTCACGCTTATCGCGGATATTTAATCCGAAAACGGGTAAAACCGTGATGCTGGCTTTTGACCATGGTTATTTTCAGGGACCGACTACCGGACTTGAACGCATTGATATAAATATCGCCCCGCGGTTTGAACATGCCg                                                                                                                            <  1:1054521/140‑1 (MQ=255)             gCTTATCGCGGATATTTAATCCGAAAACGGGTAAAACCGTGATGCTGGCTTTTGACCATGGTTATTTTCAGGGACCGACTACCGGACTTGAACGCATTGATATAAATATCGCCCCGCGGTTTGAACATGCCGATGTATTa                                                                                                                    <  1:862009/140‑1 (MQ=255)                 aTCGCGGATATTTAATCCGAAAACGGGTAAAACCGTGATGCTGGCTTTTGACCATGGTTATTTTCAGGGACCGACTACCGGACTTGAACGCATTGATATAAATATCGCCCCGCGGTTTGAACATGc                                                                                                                              >  2:120440/1‑126 (MQ=255)                 aTCGCGGATATTTAATCCGAAAACGGGTAAAACCGTGATGCTGGCTTTTGACCATGGTTATTTTCAGGGACCGACTACCGGACTTGAACGCATTGATATAAATATCGCCCCGCGGTTTGAACATGc                                                                                                                              <  1:120440/126‑1 (MQ=255)                     cgGATATTTAATCCGAAAACGGGTAAAACCGTGATGCTGGCTTTTGACCATGGTTATTTTCAGGGACCGACTACCGGACTTGAACGCATTGATATAAATATCGCCCCGCGGTTTGAACATGCCGATGTATTAATGTGTAc                                                                                                            <  2:250775/140‑1 (MQ=255)                         tatTTAATCCGAAAACGGGTAAAACCGTGATGCTGGCTTTTGACCATGGTTATTTTCAGGGACCGACTACCGGACTTGAACGCATTGATATAAATATCGCCCCGCGGTTTGAACATGCCGATGTATTAATGTGTAcgcgc                                                                                                        >  1:783582/1‑140 (MQ=255)                                    aaaaCGGGTAAAACCGTGATGCTGGCTTTTGACCATGGTTATTTTCAGGGACCGACTACCGGACTTGAACGCATTGATATAAATATCGCCCCGCGGTTTGAACATGCCGATGTATTAATGTGTACGCGCGGCATTTTgcg                                                                                             >  1:650086/1‑140 (MQ=255)                                     aaaCGGGTAAAACCGTGATGCTGGCTTTTGACCATGGTTATTTTCAGGGACCGACTACCGGACTTGAACGCATTGATATAAATATCGCCCCGCGGTTTGAACATGCCGATGTATTAATGTGTACGCGCGGCATTTTgcgc                                                                                            >  1:383300/1‑140 (MQ=255)                                        cGGGTAAAACCGTGATGCTGGCTTTTGACCATGGTTATTTTCAGGGACCGACTACCGGACTTGAACGCATTGATATAAATATCGCCCCGCGGTTTGAACATGCCGATGTATTAATGTGTACGCGCGGCATTTTGCGCAGc                                                                                         <  1:737990/140‑1 (MQ=255)                                        cGGGTAAAACCGTGATGCTGGCTTTTGACCATGGTTATTTTCAGGGACCGACTACCGGACTTGAACGCATTGATATAAATATCGCCCCGCGGTTTGAACATGCCGATGTATTAATGTGTACGCGCGGCATTTTGCGCAGc                                                                                         <  2:135736/140‑1 (MQ=255)                                        cGGGTAAAACCGTGATGCTGGCTTTTGACCATGGTTATTTTCAGGGACCGACTACCGGACTTGAACGCATTGATATAAATATCGCCCCGCGGTTTGAACATGCCGATGTATTAATGTGTACGCGCGGCATTTTGCGCAGc                                                                                         <  1:1053648/140‑1 (MQ=255)                                           gTAAAACCGTGATGCTGGCTTTTGACCATGGTTATTTTCAGGGACCGACTACCGGACTTGAACGCATTGATATAAATATCGCCCCGCGGTTTGAACATGCCGATGTATTAATGTGTACGCGCGGCATTTTGCGCAGCGTa                                                                                      >  1:413183/1‑140 (MQ=255)                                             aaaaCCGTGATGCTGGCTTTTGACCATGGTTATTTTCAGGGACCGACTACCGGACTTGAACGCATTGATATAAATATCGCCCCGCGGTTTGAACATGCCGATGTATTAATGTGTACGCGCGGCATTTTGCGCAGCGTAGt                                                                                    <  2:1239238/140‑1 (MQ=255)                                                   gTGATGCTGGCTTTTGACCATGGTTATTTTCAGGGACCGACTACCGGACTTGAACGCATTGATATAAATATCGCCCCGCGGTTTGAACATGCCGATGTATTAATGTGTACGCGCGGCATTTTGCGCa                                                                                           >  1:480860/1‑127 (MQ=255)                                                   gTGATGCTGGCTTTTGACCATGGTTATTTTCAGGGACCGACTACCGGACTTGAACGCATTGATATAAATATCGCCCCGCGGTTTGAACATGCCGATGTATTAATGTGTACGCGCGGCATTTTGCGCa                                                                                           <  2:480860/127‑1 (MQ=255)                                                     gATGCTGGCTTTTGACCATGGTTATTTTCAGGGACCGACTACCGGACTTGAACGCATTGATATAAATATCGCCCCGCGGTTTGAACATGCCGATGTATTAATGTGTACGCGCGGCATTTTGCGCAGCGTAGTTCCCCCTg                                                                            >  2:193783/1‑140 (MQ=255)                                                            gCTTTTGACCATGGTTATTTTCAGGGACCGACTACCGGACTTGAACGCATTGATATAAATATCGCCCCGCGGTTTGAACATGCCGATGTATTAATGTGTACGCGCGGCATTTTGCGCAGCGTAGTTCCCCCTGCGACCaa                                                                     <  2:1271145/140‑1 (MQ=255)                                                                 tGACCATGGTTATTTTCAGGGACCGACTACCGGACTTGAACGCATTGATATAAATATCGCCCCGCGGTTTGAACATGCCGATGTATTAATGTGTACGCGCGGCATTTTGCGCAGCGTAGTTCCCCCTGCGACCAATAGGc                                                                >  1:454503/1‑140 (MQ=255)                                                                 tGACCATGGTTATTTTCAGGGACCGACTACCGGACTTGAACGCATTGATATAAATATCGCCCCGCGGTTTGAACATGCCGATGTATTAATGTGTACGCGCGGCATTTTGCGCAGCGTAGTTCCCCCTGCGACCAATAGGc                                                                <  1:553902/140‑1 (MQ=255)                                                                      aTGGTTATTTTCAGGGACCGACTACCGGACTTGAACGCATTGATATAAATATCGCCCCGCGGTTTGAACATGCCGATGTATTAATGTGTACGCGCGGCATTTTGCGCAGCGTAGTTCCCCCTGCGACCAATAGGCCggtg                                                           >  2:139783/1‑140 (MQ=255)                                                                                   gggACCGACTACCGGACTTGAACGCATTGATATAAATATCGCCCCGCGGTTTGAACATGCCGATGTATTAATGTGTACGCGCGGCATTTTGCGCAGCGTAGTTCCCCCTGCGACCAATAGGCCGGTGGTACTGCGGGCGt                                              >  1:1046219/1‑140 (MQ=255)                                                                                         gACTACCGGACTTGAACGCATTGATATAAATATCGCCCCGCGGTTTGAACATGCCGATGTATTAATGTGTACGCGCGGCATTTTGCGCAGCGTAGTTCCCCCTGCGACCAATAGGCCGGTGGTACTGCGGGCGTCAGGTg                                        >  2:1041733/1‑140 (MQ=255)                                                                                          aCTACCGGACTTGAACGCATTGATATAAATATCGCCCCGCGGTTTGAACATGCCGATGTATTAATGTGTACGCGCGGCATTTTGCGCAGCGTAGTTCCCCCTGCGACCAATAGGCCGGTGGTACTGCGGGCGTCAGGTGc                                       <  2:393413/140‑1 (MQ=255)                                                                                          aCTACCGGACTTGAACGCATTGATATAAATATCGCCCCGCGGTTTGAACATGCCGATGTATTAATGTGTACGCGCGGCATTTTGCGCAGCGTAGTTCCCCCTGCGACCAATAGGCCGGTGGTACTGCGGGCGTCAGGTGc                                       <  2:383300/140‑1 (MQ=255)                                                                                          aCTACCGGACTTGAACGCATTGATATAAATATCGCCCCGCGGTTTGAACATGCCGATGTATTAATGTGTACGCGCGGCATTTTGCGCAGCGTAGTTCCCCCTGCGACCAATAGGCCGGTGGTACTGCGGGCGTCAGGTGc                                       <  2:1290066/140‑1 (MQ=255)                                                                                                    ttGAACGCATTGATATAAATATCGCCCCGCGGTTTGAACATGCCGATGTATTAATGTGTACGCGCGGCATTTTGCGCAGCGTAGTTCCCCCTGCGACCAATAGGCCGGTGGTACTGCGGGCGTCAGGTGCGAACTCTAtt                             >  2:1144261/1‑140 (MQ=255)                                                                                                      gAACGCATTGATATAAATATCGCCCCGCGGTTTGAACATGCCGATGTATTAATGTGTACGCGCGGCATTTTGCGCAGCGTAGTTCCCCCTGCGACCAATAGGCCGGTGGTACTGCGGGCGTCAGGTGCGAACTCTATTCt                           >  1:261641/1‑140 (MQ=255)                                                                                                        aCGCATTGATATAAATATCGCCCCGCGGTTTGAACATGCCGATGTATTAATGTGTACGCGCGGCATTTTGCGCAGCGTAGTTCCCCCTGCGACCAATAGGCCGGTGGTACTGCGGGCGTCAGGTGCGAACTCTATTCTgg                         >  1:325461/1‑140 (MQ=255)                                                                                                        aCGCATTGATATAAATATCGCCCCGCGGTTTGAACATGCCGATGTATTAATGTGTACGCGCGGCATTTTGCGCAGCGTAGTTCCCCCTGCGACCAATAGGCCGGTGGTACTGCGGGCGTCAGGTGCGAACTCTATTCTgg                         <  2:661720/140‑1 (MQ=255)                                                                                                         cGCATTGATATAAATATCGCCCCGCGGTTTGAACATGCCGATGTATTAATGTGTACGCGCGGCATTTTGCGCAGCGTAGTTCCCCCTGCGACCAATAGGCCGGTGGTACTGCGGGCGTCAGGTGCGAACTCTATTCTGGc                        >  2:253072/1‑140 (MQ=255)                                                                                                             ttGATATAAATATCGCCCCGCGGTTTGAACATGCCGATGTATTAATGTGTACGCGCGGCATTTTGCGCAGCGTAGTTCCCCCTGCGACCAATAGGCCGGTGGTACTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGaa                    <  2:454503/140‑1 (MQ=255)                                                                                                             ttGATATAAATATCGCCCCGCGGTTTGAACATGCCGATGAATTAATGTGTACGCGCGGCATTTTGCGCAGCGTAGTTCCCCCTGCGACCAATAGGCCGGTGGTACTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGaa                    >  1:1200585/1‑140 (MQ=255)                                                                                                                atatAAATATCGCCCCGCGGTTTGAACATGCCGATGTATTAATGTGTACGCGCGGCATTTTGCGCAGCGTAGTTCCCCCTGCGACCAATAGGCCGGTGGTACTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGAATTa                 <  1:1041733/140‑1 (MQ=255)                                                                                                                 tataAATATCGCCCCGCGGTTTGAACATGCCGATGTATTAATGTGTACGCGCGGCATTTTGCGCAGCGTAGTTCCCCCTGCGACCAATAGGCCGGTGGTACTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGAATTaa                <  1:7948/140‑1 (MQ=255)                                                                                                                      atatCGCCCCGTGGTTTGAACATGCCGATGTATTAATGTGTACGCGCGGCATTTTGCGCAGCGTAGTTCCCCCTGCGACCAATAGGCCGGTGGTACTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAAt           >  2:278833/1‑140 (MQ=255)                                                                                                                          cGCCCCGCGGTTTGAACATGCCGATGTATTAATGTGTACGCGCGGCATTTTGCGCAGCGTAGTTCCCCCTGCGACCAATAGGCCGGTGGTACTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGaa                    <  1:525940/127‑1 (MQ=255)                                                                                                                          cGCCCCGCGGTTTGAACATGCCGATGTATTAATGTGTACGCGCGGCATTTTGCGCAGCGTAGTTCCCCCTGCGACCAATAGGCCGGTGGTACTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGaa                    >  2:525940/1‑127 (MQ=255)                                                                                                                          cGCCCCGCGGTTTGAACATGCCGATGTATTAATGTGTACGCGCGGCATTTTGCGCAGCGTAGTTCCCCCTGCGACCAATAGGCCGGTGGTACTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAAt           >  2:1158695/1‑136 (MQ=255)                                                                                                                          cGCCCCGCGGTTTGAACATGCCGATGTATTAATGTGTACGCGCGGCATTTTGCGCAGCGTAGTTCCCCCTGCGACCAATAGGCCGGTGGTACTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAAt           <  1:1158695/136‑1 (MQ=255)                                                                                                                          cGCCCCGCGGTTTGAACATGCCGATGTATTAATGTGTACGCGCGGCATTTTGCGCAGCGTAGTTCCCCCTGCGACCAATAGGCCGGTGGTACTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGaaa       >  2:895193/1‑139 (MQ=255)                                                                                                                          cGCCCCGCGGTTTGAACATGCCGATGTATTAATGTGTACGCGCGGCATTTTGCGCAGCGTAGTTCCCCCTGCGACCAATAGGCCGGTGGTACTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAg       >  1:1334958/1‑140 (MQ=255)                                                                                                                          cGCCCCGCGGTTTGAACATGCCGATGTATTAATGTGTACGCGCGGCATTTTGCGCAGCGTAGTTCCCCCTGCGACCAATAGGCCGGTGGTACTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAg       >  1:423488/1‑140 (MQ=255)                                                                                                                            ccccGCGGTTTGAACATGCCGATGTATTAATGTGTACGCGCGGCATTTTGCGCAGCGTAGTTCCCCCTGCGACCAATAGGCCGGTGGTACTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGcc     >  2:509216/1‑140 (MQ=255)                                                                                                                               cgcgGTTTGAACATGCCGATGTATTAATGTGTACGCGCGGCATTTTGCGCAGCGTAGTTCCCCCTGCGACCAATAGGCCGGTGGTACTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCg    <  1:777265/138‑1 (MQ=255)                                                                                                                               cgcgGTTTGAACATGCCGATGTATTAATGTGTACGCGCGGCATTTTGCGCAGCGTAGTTCCCCCTGCGACCAATAGGCCGGTGGTACTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTg  <  1:253072/140‑1 (MQ=255)                                                                                                                               cgcgGTTTGAACATGCCGATGTATTAATGTGTACGCGCGGCATTTTGCGCAGCGTAGTTCCCCCTGCGACCAATAGGCCGGTGGTACTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTg  <  1:790912/140‑1 (MQ=255)                                                                                                                                  |                                                                                                                                         GAATGCAGTCACGCTTATCGCGGATATTTAATCCGAAAACGGGTAAAACCGTGATGCTGGCTTTTGACCATGGTTATTTTCAGGGACCGACTACCGGACTTGAACGCATTGATATAAATATCGCCCCGCTGTTTGAACATGCCGATGTATTAATGTGTACGCGCGGCATTTTGCGCAGCGTAGTTCCCCCTGCGACCAATAGGCCGGTGGTACTGCGGGCGTCAGGTGCGAACTCTATTCTGGCGGAATTAAGTAATGAAGCCGTG  >  NC_000913/1606203‑1606468

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