Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,606,917	T→C	V273A (GTT→GCT)	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,917	0	T	C	100.0%	142.5 / NA	44	V273A (GTT→GCT)	lsrF	3‑hydroxy‑2,4‑pentadione 5‑phosphate thiolase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (16/28);  total (16/28)

GCGCGAGGCGCTGGAAATGTGCTGGCAGGCTATCGATCAGGGCGCTTCTGGTGTGGATATGGGGCGTAATATTTTCCAGTCTGACCATCCGGTGGCGATGATGAAAGCCGTACAGGCGGTGGTTCACCATAACGAAACGGCTGATCGGGCATATGAACTCTATCTGAGTGAAAAACAGTAACTGCGGATCTAAGGAGAAGAATTATGCACGTCACACTGGTTGAAATTAACGTTCATGAAGACAAGGTTGACGA  >  NC_000913/1606795‑1607048                                                                                                                           |                                                                                                                                    gtgcgAGGCGCTGGAAATGTGCTGGCAGGCTATCGATCAGGGCGCTTCTGGTGTGGATATGGGGCGTAATATTTTCCAGTCTGACCATCCGGTGGCGATGATGAAAGCCGTACAGGCGGTGGCTCACCATAACGAAACgg                                                                                                                    <  2:799391/138‑1 (MQ=255) gtgcgAGGCGCTGGAAATGTGCTGGCAGGCTATCGATCAGGGCGCTTCTGGTGTGGATATGGGGCGTAATATTTTCCAGTCTGACCATCCGGTGGCGATGATGAAAGCCGTACAGGCGGTGGCTCACCATAACGAAACgg                                                                                                                    <  2:242140/138‑1 (MQ=255)     gAGGCGCTGGAAATGTGCTGGCAGGCTATCGATCAGGGCGCTTCTGGTGTGGATATGGGGCGTAATATTTTCCAGTCTGACCATCCGGTGGCGATGATGAAAGCCGTACAGGCGGTGGCTCACCATAACGAAACGGCTGa                                                                                                                <  1:635867/140‑1 (MQ=255)      aGGCGCTGGAAATGTGCTGGCAGGCTATCGATCAGGGCGCTTCTGGTGTGGATATGGGGCGTAATATTTTCCAGTCTGACCATCCGGTGGCGATGATGAAAGCCGTACAGGCGGTGGCTCACCATAACGAAACGGCTGAt                                                                                                               <  1:167708/140‑1 (MQ=255)      aGGCGCTGGAAATGTGCTGGCAGGCTATCGATCAGGGCGCTTCTGGTGTGGATATGGGGCGTAATATTTTCCAGTCTGACCATCCGGTGGCGATGATGAAAGCCGTACAGGCGGTGGCTCACCATAACGAAACGGCTGAt                                                                                                               <  1:788449/140‑1 (MQ=255)        gcgcTGGAAATGTGCTGGCAGGCTATCGATCAGGGCGCTTCTGGTGTGGATATGGGGCGTAATATTTTCCAGTCTGACCATCCGGTGGCGATGATGAAAGCCGTACAGGCGGTGGCTCACCATAACGAAACGGCTGATCg                                                                                                             <  2:8684/140‑1 (MQ=255)          gcTGGAAATGTGCTGGCAGGCTATCGATCAGGGCGCTTCTGGTGTGGATATGGGGCGTAATATTTTCCAGTCTGACCATCCGGTGGCGATGATGAAAGCCGTACAGGCGGTGGCTCACCATAACGAAACGGCTGATCggg                                                                                                           <  2:20121/140‑1 (MQ=255)            tGGAAATGTGCTGGCAGGCTATCGATCAGGGCGCTTCTGGTGTGGATATGGGGCGTAATATTTTCCAGTCTGACCATCCGGTGGCGATGATGAAAGCCGTACAGGCGGTGGCTCACCATAACGAAACGGCTGATCGGGCa                                                                                                         >  2:925672/1‑140 (MQ=255)              gAAATGTGCTGGCAGGCTATCGATCAGGGCGCTTCTGGTGTGGATATGGGGCGTAATATTTTCCAGTCTGACCATCCGGTGGCGATGATGAAAGCCGTACAGGCGGTGGCTCACCATAACGAAACGGCTGATCGGGCata                                                                                                       <  1:571581/140‑1 (MQ=255)                aaTGTGCTGGCAGGCTATCGATCAGGGCGCTTCTGGTGTGGATATGGGGCGTAATATTTTCCAGTCTGACCATCCGGTGGCGATGATGAAAGCCGTACAGGCGGTGGCTCACCATAACGAAACGGCTGATCGGGCATATg                                                                                                     <  1:925672/140‑1 (MQ=255)                  tgtgCTGGCAGGCTATCGATCAGGGCGCTTCTGGTGTGGATATGGGGCGTAATATTTTCCAGTCTGACCATCCGGTGGCGATGATGAAAGCCGTACAGGCGGTGGCTCACCATAACGAAACGGCTGATCGGGCATATGaa                                                                                                   <  1:87456/140‑1 (MQ=255)                      cTGGCAGGCTATCGATCAGGGCGCTTCTGGTGTGGATATGGGGCGTAATATTTTCCAGTCTGACCATCTGGTGGCGATGATGAAAGCCGTACAGGCGGTGGCTCACCATAACGAAACGGCTGATCGGGCATATGAActct                                                                                               >  2:235143/1‑140 (MQ=255)                        ggCAGGCTATCGATCAGGGCGCTTCTGGTGTGGATATGGGGCGTAATATTTTCCAGTCTGACCATCCGGTGGCGATGATGAAAGCCGTACAGGCGGTGGCTCACCATAACGAAACGGCTGATCGGGCATATGAACTCTAt                                                                                             <  2:688586/140‑1 (MQ=255)                          cAGGCTATCGATCAGGGCGCTTCTGGTGTGGATATGGGGCGTAATATTTTCCAGTCTGACCATCCGGTGGCGATGATGAAAGCCGTACAGGCGGTGGCTCACCATAACGAAACGGCTGATCGGGCATATGAACTCTATCt                                                                                           <  1:206281/140‑1 (MQ=255)                              cTATCGATCAGGGCGCTTCTGGTGTGGATATGGGGCGTAATATTTTCCAGTCTGACCATCCGGTGGCGATGATGAAAGCCGTACAGGCGGTGGCTCACCATAACGAAACGGCTGATCGGGCATATGAACTCTATCTGAGt                                                                                       >  1:578/1‑140 (MQ=255)                                aTCGATCAGGGCGCTTCTGGTGTGGATATGGGGCGTAATATTTTCCAGTCTGACCATCCGGTGGCGATGATGAAAGCCGTACAGGCGGTGGCTCACCATAACGAAACGGCTGATCGGGCATATGAACTCTATCTGAGTGa                                                                                     <  2:731245/140‑1 (MQ=255)                                aTCGATCAGGGCGCTTCTGGTGTGGATATGGGGCGTAATATTTTCCAGTCTGACCATCCGGTGGCGATGATGAAAGCCGTACAGGCGGTGGCTCACCATAACGAAACGGCTGATCGGGCATATGAACTCTATCTGAGTGa                                                                                     <  2:803208/140‑1 (MQ=255)                                       aGGGCGCTTCTGGTGTGGATATGGGGCGTAATATTTTCCAGTCTGACCATCCGGTGGCGATGATGAAAGCCGTACAGGCGGTGGCTCACCATAACGAAACGGCTGATCGGGCATATGAACTCTATCTGAGTGAAAAACAg                                                                              >  1:667711/1‑140 (MQ=255)                                        gggCGCTTCTGGTGTGGATATGGGGCGTAATATTTTCCAGTCTGACCATCCGGTGGCGATGATGAAAGCCGTACAGGCGGTGGCTCACCATAACGAAACGGCTGATCGGGCATATGAACTCTATCTGAGTGAAAAACAGt                                                                             <  2:667711/140‑1 (MQ=255)                                           cgcTTCTGGTGTGGATATGGGGCGTAATATTTTCCAGTCTGACCATCCGGTGGCGATGATGAAAGCCGTACAGGCGGTGGCTCACCATAACGAAACGGCTGATCGGGCATATGAACTCTATCTGAGTGAAAAACAGTAAc                                                                          >  2:159035/1‑140 (MQ=255)                                               tCTGGTGTGGATATGGGGCGTAATATTTTCCAGTCTGACCATCCGGTGGCGATGATGAAAGCCGTACAGGCGGTGGCTCACCATAACGAAACGGCTGATCGGGCATATGAACTCTATCTGAGTGAAAAACAGTAACTGCg                                                                      >  1:30186/1‑140 (MQ=255)                                                       ggATATGGGGCGTAATATTTTCCAGTCTGACCATCCGGTGGCGATGATGAAAGCCGTACAGGCGGTGGCTCACCATAACGAAACGGCTGATCGGGCATATGAACTCTATCTGAGTGAAAAACAGTAACTGCGGATCTAAg                                                              <  2:202529/140‑1 (MQ=255)                                                        gATATGGGGCGTAATATTTTCCAGTCTGACCATCCGGTGGCGATGATGAAAGCCGTACAGGCGGTGGCTCACCATAACGAAACGGCTGATCGGGCATATGAACTCTATCTGAGTGAAAAACAGTAACTGCGGATCTAAgg                                                             >  1:57746/1‑140 (MQ=255)                                                            tGGGGCGTAATATTTTCCAGTCTGACCATCCGGTGGCGATGATGAAAGCCGTACAGGCGGTGGCTCACCATAACGAAACGGCTGATCGGGCATATGAACTCTATCTGAGTGAAAAACAGTAACTGCGGATCTAAGGagaa                                                         <  2:183632/140‑1 (MQ=255)                                                            tGGGGCGTAATATTTTCCAGTCTGACCATCCGGTGGCGATGATGAAAGCCGTACAGGCGGTGGCTCACCATAACGAAACGGCTGATCGGGCATATGAACTCTATCTGAGTGAAAAACAGTAACTGCGGATCTAAGGagaa                                                         <  1:648065/140‑1 (MQ=255)                                                             ggggCGTAATATTTTCCAGTCTGACCATCCGGTGGCGATGATGAAAGCCGTACAGGCGGTGGCTCACCATAACGAAACGGCTGATCGGGCATATGAACTCTATCTGAGTGAAAAACAGTAACTGCGGATCTAAGGagaag                                                        >  1:635484/1‑140 (MQ=255)                                                                  gTAATATTTTCCAGTCTGACCATCCGGTGGCGATGATGAAAGCCGTACAGGCGGTGGCTCACCATAACGAAACGGCTGATCGGGCATATGAACTCTATCTTAGTCAAAAACAGTAACAGCGGTTCTAAGGAGAAGAATTa                                                   >  2:416227/1‑140 (MQ=255)                                                                  gTAATATTTTCCAGTCTGACCATCCGGTGGCGATGATGAAAGCCGTACAGGCGGTGGCTCACCATAACGAAACGGCTGATCGGGCATATGAACTCTATCTGAGTGAAAAACAGTAACTGCGGATCTAAGGAGAAGAATTa                                                   >  1:583701/1‑140 (MQ=255)                                                                   tAATATTTTCCAGTCTGACCATCCGGTGGCGATGATGAAAGCCGTACAGGCGGTGGCCCACCATAACGAAACGGCTGGTCGGGCATATGAACCCTATCTGAGTGAAAAACAGGAACCGCGGATCTAAGGAGGAGAAatat                                                  >  2:648065/1‑140 (MQ=255)                                                                           tCCAGTCTGACCATCCGGTGGCGATGATGAAAGCCGTACAGGCGGTGGCTCACCATAACGAAACGGCTGATCGGGCATATGAACTCTATCTGAGTGAAAAACAGTAACTGCGGATCTAAGGAGAAGAATTATGCACGTca                                          <  2:907067/140‑1 (MQ=255)                                                                            ccAGTCTGACCATCCGGTGGCGATGATGAAAGCCGTACAGGCGGTGGCTCACCATAACGAAACGGCTGATCGGGCATATGAACTCTATCTGAGTGAAAAACAGTAACTGCGGATCTAAGGAGAAGAATTATGCACGTcac                                         <  2:57746/140‑1 (MQ=255)                                                                             cAGTCTGACCATCCGGTGGCGATGATGAAAGCCGTACAGGCGGTGGCTCACCATAACGAAACGGCTGATCGGGCATATGAACTCTATCTGAGTGAAAAACAGTAACTGCGGATCTAAGGAGAAGAATTATGCACGTcaca                                        <  2:578/140‑1 (MQ=255)                                                                              aGTCTGACCATCCGGCGGCGATCATCAAAACCGTACAGGCGGTGGCTCACCATAACGAAACGGCTGATCGGGCATATGAACTCTATCTGAGTGAAAAACAGTAACTGCGGATCTAAGGAGAAGAATTATGCACGTCacac                                       <  1:470085/140‑1 (MQ=255)                                                                                tCTGACCATCCGGTGGCGATGATGAAAGCCGTACAGGCGGTGGCTCACCCTAACGAAACGGCCGATCGGGCATATTAAATCTATCTGAGGGAAAAACAGTAACTACGGATCTAAAGAGAAAAATTATTGACGTCACACTg                                     >  2:995330/1‑140 (MQ=255)                                                                                tCTGACCATCCGGTGGCGATGATGAAAGCCGTACAGGCGGTGGCTCACCATAACGAAACGGCTGATCGGGCATATGAACTCTATCTGAGTGAAAAACAGTAACTGCGGATCTAAGGAGAAGAATTATGCACGTCACACTg                                     >  1:789269/1‑140 (MQ=255)                                                                                tCTGACCATCCGGTGGCGATGATGAAAGCCGTACAGGCGGTGGCTCACCATAACGAAACGGCTGATCGGGCATATGAACTCTATCTGAGTGAAAAACAGTAACTGCGGATCTAAGGAGAAGAATTATGCACGTCACACTg                                     >  1:800286/1‑140 (MQ=255)                                                                                tCTGACCATCCGGTGGCGATGATGAAAGCCGTACAGGCGGTGGCTCACCATAACGAAACGGCTGATCGGGCATATGAACTCTATCTGAGGGCAAAACAGTAACTACCGATCTAAGGAAAAGAGTTATGCACATCACACTg                                     >  2:470085/1‑140 (MQ=255)                                                                                  tGACCATCCGGTGGCGATGATGAAAGCCGTACAGGCGGTGGCTCACCATAACGAAACGGCTGATCGGCCATATGAACTCTATCTGAGTGAAAAACAGTAACTGCGGATCTAAGGAGAAGAATTATGCACGTCACACTGGt                                   <  1:66165/140‑1 (MQ=255)                                                                                              ggCGATGATGAAAGCCGTACAGGCGGTGGCTCACCATAACGAAACGGCTGATCGGGCATATGAACTCTATCTGAGTGAAAAACAGTAACTGCGGATCTAAGGAGAAGAATTATGTACGTCACACTGGTTGAAATTAACGt                       >  2:915828/1‑140 (MQ=255)                                                                                                  atgatgAAAGCCGTACAGGCGGTGGCTCACCATAACGAAACGGCTGATCGGGCATATGAACTCTATCTGAGTGAAAAACAGTAACTGCGGATCTAAGGAGAAGAATTATGCACGTCACACTGGTTGAAATTAACGTTCAt                   <  1:514452/140‑1 (MQ=255)                                                                                                     atgaAAGCCGTACAGGCGGTGGCTCACCATAACGAAACGGCTGATCGGGCATATGAACTCTATCTGAGTGAAAAACAGTAACTGCGGATCTAAGGAGAAGAATTATGTACGTCACACTGGTTGAAATTAACGTTCATGaa                <  1:915828/140‑1 (MQ=255)                                                                                                         aaGCCGTACAGGCGGTGGCTCACCATAACGAAACGGCTGATCGGGCATATGAACTCTATCTGAGTGAAAAACAGTAACTGCGGATCTAAGGAGAAGAATTATGCACGTCACACTGGTTGAAATTAACGTTCATGAAGACa            <  2:24847/140‑1 (MQ=255)                                                                                                         aaGCCGTACAGGCGGAGGCTCAACACAACCAAACGGCCGATCGGGCATATGAACTCTATCTGAGTGAAAAACAGTAACTGCGGATCTAAGGAGAAGAATTATGCACGTCACACTGGTTGAAATTAACGTTCATGAAGACa            <  1:416227/140‑1 (MQ=255)                                                                                                               tACAGGCGGTGGCTCACCATAACGAAACGGCTGATCGGGCATATGAACTCTTTCTTAGTGAAAAACAGTAACTGCGGATCTAAGGAGAAGAATTATGCACGTCACACTGGGTGAAATTAACGTTCATGAAGACAAGGttt      >  2:749096/1‑139 (MQ=255)                                                                                                               tACAGGCGGTGGCTCACCATAACGAAACGGCTGATCGGGCATATGAACTCTATCTGAGTGAAAAACAGTAACTGCGGATCTAAGGAGAAGAATTATGCACGTCACACTGGTTGAAATTAACGTTCATGAAGACAAGGTTg      <  1:995330/140‑1 (MQ=255)                                                                                                                   ggCGGTGGCTCACCATAACGAAACGGCTGATCGGGCATATGAACTCTATCTGAGTGAAAAACAGTAACTGCGGATCTAAGGAGAAGAATTATGCACGTCACACTGGTTGAAATTAACGTTCATGAAGACAAGGTTGACGa  >  1:790472/1‑140 (MQ=255)                                                                                                                           |                                                                                                                                    GCGCGAGGCGCTGGAAATGTGCTGGCAGGCTATCGATCAGGGCGCTTCTGGTGTGGATATGGGGCGTAATATTTTCCAGTCTGACCATCCGGTGGCGATGATGAAAGCCGTACAGGCGGTGGTTCACCATAACGAAACGGCTGATCGGGCATATGAACTCTATCTGAGTGAAAAACAGTAACTGCGGATCTAAGGAGAAGAATTATGCACGTCACACTGGTTGAAATTAACGTTCATGAAGACAAGGTTGACGA  >  NC_000913/1606795‑1607048

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