Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,497,153	G→A	D381D (GAC→GAT)	alaC ←	glutamate‑‑pyruvate aminotransferase AlaC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,497,153	0	G	A	100.0%	88.2 / NA	28	D381D (GAC→GAT)	alaC	glutamate‑‑pyruvate aminotransferase AlaC
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (12/16);  total (12/16)

ATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACGAATACGGTCGCGGTTTTCAATCAGTGCAAAGCGAACATGCGTATCACCATAGTCGCCAAAGCCAATCCCCGGCGAGACACAAACCTTCGCTTCGTTAAGGAGCTTCTTGGCAAATTCCAGCGATCCCATGGCCGCATAT  >  NC_000913/2497023‑2497285                                                                                                                                   |                                                                                                                                     aTCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACGAATACGATCGCGGgtt                                                                                                                             <  2:142799/140‑1 (MQ=255) aTCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGGAGTAAACCGGTGGGCCGGAACATCGCTTTAATCCCCCTGATGGCCTGACGAATACGATCGCGGttt                                                                                                                             >  2:712077/1‑140 (MQ=255)      gCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACGAATACGATCGCGGTTTTCAAt                                                                                                                        >  1:524654/1‑140 (MQ=255)       cTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACGAATACGATCGCGGTTTTCAATc                                                                                                                       >  1:371751/1‑140 (MQ=255)          cTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACGAATACGATCGCGGTTTTCAATCAGt                                                                                                                    >  1:716140/1‑140 (MQ=255)             tttGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACGAATACGATCGCGGTTTTCAATCAGTGCa                                                                                                                 >  1:37899/1‑140 (MQ=255)                       tGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGCTGGCCTGAAGAAGACGTTCGCGGTTATCAAGCAGTGCAAAGCGAACAt                                                                                                       >  2:214692/1‑140 (MQ=255)                          cgacgGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACGAATACGATCGCGGTTTTCAATCAGTGCAAAGCGAACATGCg                                                                                                    <  2:638017/140‑1 (MQ=255)                          cgacgGTGTTTATTCCGCGTTTCCGCGAACATGTTTGTTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACGAATACGATCGCGGTTTTCAATCAGTGCAAAGCGAACATGCg                                                                                                    <  1:214692/140‑1 (MQ=255)                                    tATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACGAATACGATCGCGGTTTTCAATCAGTGCAAAGCGAACATGCGTATCACCATa                                                                                          <  1:633371/140‑1 (MQ=255)                                     aTTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACGAATACGATCGCGGTTTTCAATCAGTGCAAAGCGAACATGCGTATCACCATAg                                                                                         <  2:716140/140‑1 (MQ=255)                                            gTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACGAATACGATCGCGGTTTTCAATCAGTGCAAAGCGAACATGCGTATCACCATAGTCGCCaa                                                                                  >  1:167626/1‑140 (MQ=255)                                                           tttGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACGAATACGATCGCGGTTTTCAATCAGTGCAAAGCGAACATGCGTATCACCATAGTCGCCAAAGCCAATCCCCGGCg                                                                   <  2:400994/140‑1 (MQ=255)                                                                    ggCGGGGAGTAAACCGTCGGCACGGGACATCGCTTTAATCCCTCTGATGGCCTGACGAATACGATCGCGGTTTTCAATCAGTGCAAAGCGAACATGCGTATCACCATAGTCGCCAAAGCCAATCCCCGGCGAGACACaaa                                                          >  2:652259/1‑140 (MQ=255)                                                                          tagtaAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACGAATGCGATCGCCGTTTTCAAACAGGGCAAAGCGAACATGCGTATCACCATAGTCGTCAAAGCAAATCCCCGGCGAGACACAAACcttcg                                                    <  2:601821/140‑1 (MQ=255)                                                                          tagtaAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACGAATACGATCGCGGTTTTCAATCAGTGCAAAGCGAACATGCGTATCACCATAGTCGCCAAAGCCAATCCCCGGCGAGACACAAACcttcg                                                    <  2:167145/140‑1 (MQ=255)                                                                             taAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACGAATACGATCGCGGTTTTCAATCAGTGCAAAGCGAACATGCGTATCACCATAGTCGCCAAAGCCAATCCCCGGCGAGACACAAACcttcgctt                                                 <  2:540377/140‑1 (MQ=255)                                                                             taAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACGAATACGATCGCGGTTTTCAATCAGTGCAAAGCGAACATGCGTATCACCATAGTCGCCAAAGCCAATCCCCGGCGAGACACAAACcttcgctt                                                 <  2:167626/140‑1 (MQ=255)                                                                                   gTCGGCACGGAACATTGCTTTACTCCCTCAGATGGCCTGACGAATACGATCGCGGTTTTCAATCAGTGCAAAGCGAACATGCGTATCACCATAGTCGCCAAAGCCAATCCCCGGCGAGACACAAACCTTCGCTTCGTTaa                                           <  1:356654/140‑1 (MQ=255)                                                                                     cGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACGAATACGATCGCGGTTTTCAATCAGTGCAAAGCGAACATGCGTATCACCATAGTCGCCAAATCCAATCCCCGGCGAGAGACAATGCTTCGCTCCGTTAAgg                                         >  2:356654/1‑140 (MQ=255)                                                                                     cGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACGAATACGATCGCGGTTTTCAATCAGTGCAAAGCGAACATGCGTATCACCATAGTCGCCAAAGCCAATCCCCGGCGAGACACAAACCTTCGCTTCGTTAAgg                                         >  1:447843/1‑140 (MQ=255)                                                                                                 tCGCTTTAATCCCTCTGATGGCCTGACGAATACGATCGCGGTTTTCAATCAGTGCAAAGCGAACATGCGTATCACCATAGTCGCCAAAGCCAATCCCCGGCGAGACACAAACCTTCGCTTCGTTAAGGAGCTTCTTGGCa                             <  2:386188/140‑1 (MQ=255)                                                                                                       tAATCCCTCTGATGGCCTGACGAATACGATCGCGGTTTTCAATCAGTGCAAAGCGAACATGCGTATCACCATAGTCGCCAAAGCCAATCCCCGGCGAGACACAAACCTTCGCTTCGTTAAGGAGCTTCTTGGCAAATTcc                       >  1:594076/1‑140 (MQ=255)                                                                                                               ctGATGGCCTGACGAATACGATCGCGGTTTTCAATCAGTGCAAAGCGAACATGCGTATCACCATAGTCGCCAAAGCCAATCCCCGGCGAGACACAAACCTTCGCTTCGTTAAGGAGCTTCTTGGCAAAt                          <  1:514400/129‑1 (MQ=255)                                                                                                               ctGATGGCCTGACGAATACGATCGCGGTTTTCAATCAGTGCAAAGCGAACATGCGTATCACCATAGTCGCCAAAGCCAATCCCCGGCGAGACACAAACCTTCGCTTCGTTAAGGAGCTTCTTGGCAAAt                          >  2:514400/1‑129 (MQ=255)                                                                                                               ctGATGGCCTGACGAATACGATCGCGGTTTTCAATCAGTGCAAAGCGAACATGCGTATCACCATAGTCGCCAAAGCCAATCCCCGGCGAGACACAAACCTTCACTTCGTTAAGGAGCTTCTTGGGAAATTCCAGCGATcc               >  1:355411/1‑140 (MQ=255)                                                                                                                 gATGGCCTGACGAATACGATCGCGGTTTTCAATCAGTGCAAAGCGAACATGCGTATCACCATAGTCGCCAAAGCCAATCCCCGGCGAGACACAAACCTTCGCTTCGTTAAGGAGCTTCTTGGCAAATTCCAGCGATCCCa             <  1:577203/140‑1 (MQ=255)                                                                                                                    ggCCTGACGAATACGATCGCGGTTTTCAATCAGTGCAAAGCGAACATGCGTATCACCATAGTCGCCAAAGCCAATCCCCGGCGAGACACAAACCTTCGCTTCGTTAAGGAGCTTCTTGGCAAATTCCAGCGATCCCATgg          <  1:474197/140‑1 (MQ=255)                                                                                                                     gCCTGACGAATACGATCGCGGTTTTCAATCAGTGCAAAGCGAACATGCGTATCACCATAGTCGCCAAAGCCAATCCCCGGCGAGACACAAACCTTCGCTTCGTTAAGGAGCTTCTTGGCAAATTCCAGCGATCCCATGGc         <  1:697600/140‑1 (MQ=255)                                                                                                                       cTGACGAATACGATCGCGGTTTTCAATCAGTGCAAAGCGAACATGCGTATCACCATAGTCGCCAAAGCCAATCCCCGGCGAGACACAAACCTTCGCTTCGTTAAGGAGCTTCTTGGCAAATTCCAGCGATCCCATGGCCg       >  1:661719/1‑140 (MQ=255)                                                                                                                           cGAATACGATCGCGGTTTTCAATCAGTGCAAAGCGAACATGCGTATCACCATAGTCGCCAAAGCCAATCCCCGGCGAGACACAAACCTTCGCTTCGTTAAGGAGCTTCTTGGCAAATTCCAGCGATCCCATGGCCGCata   <  2:524654/140‑1 (MQ=255)                                                                                                                            gAATACGAACGCGGTTTTCAATCAGTGCAAAGCGAACATGCGTATCACCATAGTCGCCAAAGCCAATCCCCGGCGAGACACAAACCTTCGCTTCGTTAAGGAGCTTCTTGGCAAATTCCAGCGATCCCATGGCCGCatat  >  2:177891/1‑140 (MQ=255)                                                                                                                                   |                                                                                                                                     ATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACGAATACGGTCGCGGTTTTCAATCAGTGCAAAGCGAACATGCGTATCACCATAGTCGCCAAAGCCAATCCCCGGCGAGACACAAACCTTCGCTTCGTTAAGGAGCTTCTTGGCAAATTCCAGCGATCCCATGGCCGCATAT  >  NC_000913/2497023‑2497285

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: ‑