Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	173,839	A→G	V348V (GTT→GTC)	hemL ←	glutamate‑1‑semialdehyde aminotransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	173,839	0	A	G	100.0%	124.6 / NA	36	V348V (GTT→GTC)	hemL	glutamate‑1‑semialdehyde aminotransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (22/14);  total (22/14)

TCGTCCAGCATCATATGGAAGAAACGCTTAAAGCGTTCCACGTCACAGGCCATCACATCCTGATAGCACGTCACGGACTCGGCGTCGGTAAAGAAAATACCGAACATGCCGCCAACGTGGTTAACGACCAGCGGAATTCCGGCTTCTTCTGCCGCTTCCAGCAGACCTTCTGCCAGACGTGTTGTCAGCTCATCCAGCGTTTCGTGAACGCCCGGCTGCGCGACTTCATTCAGACAGGCGAAACCCGC  >  NC_000913/173717‑173964                                                                                                                           |                                                                                                                              tCGTCCAGCATCATATGGAAGAAACGCTTAAAGCGTTCCACGTCACAGGCCATCACATCCTGATAGCACGTCACGGACTCGGCGTCGGTAAAGAAAATACCGAACATGCCGCCAACGTGGTTGACGACCAGCGGAATTcc                                                                                                              >  1:695174/1‑140 (MQ=255) tCGTCCAGCATCATATGGAAGAAACGCTTAAAGCGTTCCACGTCACAGGCCATCACATCCTGATAGCACGTCACGGACTCGGCGTCGGTAAAGAAAATACCGAACATGCCGCCAACGTGGTTGACGACCAGCGGAATTcc                                                                                                              >  2:457619/1‑140 (MQ=255)                  gaagaaACGCTTAAAGCGTTCCACGTCACAGGCCATCACATCCTGATAGCACGTCACGGACTCGGCGTCGGTAAAGAAAATACCGAACATGCCGCCAACGTGGTTGACGACCAGCGGAATTCCGGCTTCTTCTGCCGCtt                                                                                             <  1:470649/140‑1 (MQ=255)                  gaagaaACGCTTAAAGCGTTCCACGTCACAGGCCATCACATCCTGATAGCACGTCACGGACTCGGCGTCGGTAAAGAAAATACCGAACATGCCGCCAACGTGGTTGACGACCAGCGGAATTCCGGCTTCTTCTGCCGCtt                                                                                             <  2:512446/140‑1 (MQ=255)                   aagaaACGCTTAAAGCGTTCCACGTCACAGGCCATCACATCCTGATAGCACGTCACGGACTCGGCGTCGGTAAAGAAAATACCGAACATGCCGCCAACGTGGTTGACGACCAGCGGAATTCCGGCTTCTTCTGCCGCTTc                                                                                            >  1:282442/1‑140 (MQ=255)                            ttAAAGCGTTCCACGTCACAGGCCATCACATCCTGATAGCACGTCACGGACTCGGCGTCGGTAAAGAAAATACCGAACATGCCGCCAACGTGGTTGACGACCAGCGGAATTCCGGCTTCTTCTGCCGCTTCCAGCAGAcc                                                                                   <  1:640793/140‑1 (MQ=255)                                      ccACGTCACAGGCCATCACATCCTGATAGCACGTCACGGACTCGGCGTCGGTAAAGAAAATACCGAACATGCCGCCAACGTGGTTGACGACCAGCGGAATTCCGGCTTCTTCTGCCGCTTCCAGCAGACCTTCTGCCAGa                                                                         >  2:397117/1‑140 (MQ=255)                                            cacaGGCCATCACATCCTGATAGCACGTCACGGACTCGGCGTCGGTAAAGAAAATACCGAACATGCCGCCAACGTGGTTGACGACCAGCGGAATTCCGGCTTCTTCTGCCGCTTCCAGCAGACCTTCTGCCAGACGtgtt                                                                   >  1:655188/1‑140 (MQ=255)                                             acaGGCCATCACATCCTGATAGCACGTCACGGACTCGGCGTCGGTAAAGAAAATACCGAACATGCCGCCAACGTGGTTGACGACCAGCGGAATTCCGGCTTCTTCTGCCGCTTCCAGCAGACCTTCTGCCAGACGtgttg                                                                  <  1:323610/140‑1 (MQ=255)                                                      cacaTCCTGATAGCACGTCACGGACTCGGCGTCGGTAAAGAAAATACCGAACATGCCGCCAACGTGGTTGACGACCAGCGGAATTCCGGCTTCTTCTGCCGCTTCCAGCAGACCTTCTGCCAGACGTGTTGTCAGCTCAt                                                         >  1:313284/1‑140 (MQ=255)                                                      cacaTCCTGATAGCACGTCACGGACTCGGCGTCGGTAAAGAAAATACCGAACATGCCGCCAACGTGGTTGACGACCAGCGGAATTCCGGCTTCTTCTGCCGCTTCCAGCAGACCTTCTGCCAGACGTGTTGTCAGCTCAt                                                         >  2:369553/1‑140 (MQ=255)                                                               aTAGCACGTCACGGACTCGGCGTCGGTAAAGAAAATACCGAACATGCCGCCAACGTGGTTGACGACCAGCGGAATTCCGGCTTCTTCTGCCGCTTCCAGCAGACCTTCTGCCAGACGTGTTGTCAGCTCATCCAGCGttt                                                >  2:557515/1‑140 (MQ=255)                                                               aTAGCACGTCACGGACTCGGCGTCGGTAAAGAAAATACCGAACATGCCGCCAACGTGGTTGACGACCAGCGGAATTCCGGCTTCTTCTGCCGCTTCCAGCAGACCTTCTGCCAGACGTGTTGTCAGCTCATCCAGCGttt                                                >  2:398704/1‑140 (MQ=255)                                                                  gCACGTCACGGACTCGGCGTCGGTAAAGAAAATACCGAACATGCCGCCAACGTGGTTGACGACCAGCGGAATTCCGGCTTCTTCTGCCGCTTCCAGCAGACCTTCTGCCAGACGTGTTGTCAGCTCATCCAGCGTTTCGt                                             <  2:282442/140‑1 (MQ=255)                                                                   cACGTCACGGACTCGGCGTCGGTAAAGAAAATACCGAACATGCCGCCAACGTGGTTGACGACCAGCGGAATTCCGGCTTCTTCTGCCGCTTCCAGCAGACCTTCTGCCAGACGTGTTGTCAGCTCATCCAGCGtt                                                 <  2:682952/135‑1 (MQ=255)                                                                   cACGTCACGGACTCGGCGTCGGTAAAGAAAATACCGAACATGCCGCCAACGTGGTTGACGACCAGCGGAATTCCGGCTTCTTCTGCCGCTTCCAGCAGACCTTCTGCCAGACGTGTTGTCAGCTCATCCAGCGtt                                                 >  1:682952/1‑135 (MQ=255)                                                                   cACGTCACGGACTCGGCGTCGGTAAAGAAAATACCGAACATGCCGCCAACGTGGTTGACGACCAGCGGAATTCCGGCTTCTTCTGCCGCTTCCAGCAGACCTTCTGCCAGACGTGTTGTCAGCTCATCCAGCGTTTCGTg                                            >  1:591061/1‑140 (MQ=255)                                                                      gTCACGGACTCGGCGTCGGTAAAGAAAATACCGAACATGCCGCCAACGTGGTTGACGACCAGCGGAATTCCGGCTTCTTCTGCCGCTTCCAGCAGACCTTCTGCCAGACGTGTTGTCAGCTCATCCAGCGTTTCGTGAAc                                         >  2:270170/1‑140 (MQ=255)                                                                       tCACGGACTCGGCGTCGGTAAAGAAAATACCGAACATGCCGCCAACGTGGTTGACGACCAGCGGAATTCCGGCTTCTTCTGCCGCTTCCAGCAGACCTTCTGCCAGACGTGTTGTCAGCTCATCCAGCGTTTCGt                                             <  2:96256/135‑1 (MQ=255)                                                                       tCACGGACTCGGCGTCGGTAAAGAAAATACCGAACATGCCGCCAACGTGGTTGACGACCAGCGGAATTCCGGCTTCTTCTGCCGCTTCCAGCAGACCTTCTGCCAGACGTGTTGTCAGCTCATCCAGCGTTTCGt                                             >  1:96256/1‑135 (MQ=255)                                                                         aCGGACTCGGCGTCGGTAAAGAAAATACCGAACATGCCGCCAACGTGGTTGACGACCAGCGGAATTCCGGCTTCTTCTGCCGCTTCCAGCAGACCTTCTGCCAGACGTGTTGTCAGCTCATCCAGCGTTTCGTGAACGcc                                      >  1:464146/1‑140 (MQ=255)                                                                            gACTCGGCGTCGGTAAAGAAAATACCGAACATGCCGCCAACGTGGTTGACGACCAGCGGAATTCCGGCTTCTTCTGCCGCTTCCAGCAGACCTTCTGCCAGACGTGTTGTCAGCTCATCCAGCGTTTCGTGAACGCCCgg                                   <  2:655188/140‑1 (MQ=255)                                                                                cGGCGTCGGTAAAGAAAATACCGAACATGCCGCCAACGTGGTTGACGACCAGCGGAATTCCGGCTTCTTCTGCCGCTTCCAGCAGACCTTCTGCCAGACGTGTTGTCAGCTCATCCAGCGTTTCGTGAACGCCCGGCTgc                               >  2:406871/1‑140 (MQ=255)                                                                                 ggCGTCGGTAAAGAAAATACCGAACATGCCGCCAACGTGGTTGACGACCAGCGGAATTCCGGCTTCTTCTGCCGCTTCCAGCAGACCTTCTGCCAGACGTGTTGTCAGCTCATCCAGCGTTTCGTGAACGCCCGGCTgcg                              >  2:343127/1‑140 (MQ=255)                                                                                  gCGTCGGTAAAGAAAATACCGAACATGCCGCCAACGTGGTTGACGACCAGCGGAATTCCGGCTTCTTCTGCCGCTTCCAGCAGACCTTCTGCCAGACGTGTTGTCAGCTCATCCAGCGTTTCGTGa                                           <  1:171831/126‑1 (MQ=255)                                                                                  gCGTCGGTAAAGAAAATACCGAACATGCCGCCAACGTGGTTGACGACCAGCGGAATTCCGGCTTCTTCTGCCGCTTCCAGCAGACCTTCTGCCAGACGTGTTGTCAGCTCATCCAGCGTTTCGTGa                                           >  2:171831/1‑126 (MQ=255)                                                                                     tCGGTAAAGAAAATACCGAACATGCCGCCAACGTGGTTGACGACCAGCGGAATTCCGGCTTCTTCTGCCGCTTCCAGCAGACCTTCTGCCAGACGTGTTGTCAGCTCATCCAGCGTTTCGTGAACGCCCGGCTGCGCGAc                          <  2:221205/140‑1 (MQ=255)                                                                                      cGGTAAAGAAAATACCGAACATGCCGCCAACGTGGTTGACGACCAGCGGAATTCCGGCTTCTTCTGCCGCTTCCAGCAGACCTTCTGCCAGACGTGTTGTCAGCTCATCCAGCGTTTCGTGAACGCCCGGCTGCGCGACt                         <  1:343127/140‑1 (MQ=255)                                                                                       ggTAAAGAAAATACCGAACATGCCGCCAACGTGGTTGACGACCAGCGGAATTCCGGCTTCTTCTGCCGCTTCCAGCAGACCTTCTGCCAGACGTGTTGTCAGCTCATCCAGCGTTTCGTGAACGCCCGGCTGCGCGACtt                        >  1:395938/1‑140 (MQ=255)                                                                                       ggTAAAGAAAATACCGAACATGCCGCCAACGTGGTTGACGACCAGCGGAATT‑CGGCTTCTTCTGCCGCTTCCAGCAGACCTTCTGCCAGACGTGTTGTCAGCTCATCCAGCGTTTCGTGAACGCCCGGCTGCGCGAc                          <  2:32771/137‑1 (MQ=255)                                                                                       ggTAAAGAAAATACCGAACATGCCGCCAACGTGGTTGACGACCAGCGGAATT‑CGGCTTCTTCTGCCGCTTCCAGCAGACCTTCTGCCAGACGTGTTGTCAGCTCATCCAGCGTTTCGTGAACGCCCGGCTGCGCGAc                          >  1:32771/1‑137 (MQ=255)                                                                                               aaaTACCGAACATGCCGCCAACGTGGTTGACGACCAGCGGAATTCCGGCTTCTTCTGCCGCTTCCAGCAGACCTTCTGCCAGACGTGTTGTCAGCTCATCCAGCGTTTCGTGAACGCCCGGCTGCGCGACTTCATTCAGa                <  1:406871/140‑1 (MQ=255)                                                                                                   aCCGAACATGCCGCCAACGTGGTTGACGACCAGCGGAATTCCGGCTTCTTCTGCCGCTTCCAGCAGACCTTCTGCCAGACGTGTTGTCAGCTCATCCAGCGTTTCGTGAACGCCCGGCTGCGCGACTTCATTCAGACAgg            >  2:662214/1‑140 (MQ=255)                                                                                                        aCATGCCGCCAACGTGGTTGACGACCAGCGGAATTCCGGCTTCTTCTGCCGCTTCCAGCAGACCTTCTGCCAGACGTGTTGTCAGCTCATCCAGCGTTTCGTGAACGCCCGGCTGCGCGACTTCATTCAGACAGGCGa         <  1:695115/138‑1 (MQ=255)                                                                                                        aCATGCCGCCAACGTGGTTGACGACCAGCGGAATTCCGGCTTCTTCTGCCGCTTCCAGCAGACCTTCTGCCAGACGTGTTGTCAGCTCATCCAGCGTTTCGTGAACGCCCGGCTGCGCGACTTCATTCAGACAGGCGa         >  2:695115/1‑138 (MQ=255)                                                                                                          aTGCCGCCAACGTGGTTGACGACCAGCGGAATTCCGGCTTCTTCTGCCGCTTCCAGCAGACCTTCTGCCAGACGTGTTGTCAGCTCATCCAGCGTTTCGTGAACGCCCGGCTGCGCGACTTCATTCAGACAGGCGAAAcc     <  1:662214/140‑1 (MQ=255)                                                                                                            gccgccAACGTGG‑TGACGACCAGCGGAATTCCGGCTTCTTCTGCCGCTTCCAGCAGACCTTCTGCCAGACGTGTTGTCAGCTCATCCAGCGTTTCGTGAACGCCCGGCTGCGCGACTTCATTCAGACAGGCGAAACCCGc  >  1:489694/1‑140 (MQ=255)                                                                                                                           |                                                                                                                              TCGTCCAGCATCATATGGAAGAAACGCTTAAAGCGTTCCACGTCACAGGCCATCACATCCTGATAGCACGTCACGGACTCGGCGTCGGTAAAGAAAATACCGAACATGCCGCCAACGTGGTTAACGACCAGCGGAATTCCGGCTTCTTCTGCCGCTTCCAGCAGACCTTCTGCCAGACGTGTTGTCAGCTCATCCAGCGTTTCGTGAACGCCCGGCTGCGCGACTTCATTCAGACAGGCGAAACCCGC  >  NC_000913/173717‑173964

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