Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,134,811	Δ1 bp	coding (681/693 nt)	yghT →	putative ATP‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,134,811	0	T	.	100.0%	110.1 / NA	27	coding (681/693 nt)	yghT	putative ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (13/14);  total (13/14)

TTGCCATAACGCCACAGTTAACCTTTGGTGGTGCACAACTGGTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGTGGTTATCACTAGTCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGGTAAAGTAAATACCCA  >  NC_000913/3134690‑3134943                                                                                                                          |                                                                                                                                     ttGCCATAACGCCACAGTTAACCTTTGGTGGTGCACAACTGGTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCaaaa                                                                                                                    >  2:743012/1‑139 (MQ=255)      aTAACGCCACAGTTAACCTTTGGTGGTGCACAACTGGTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGa                                                                                                               >  1:310532/1‑139 (MQ=255)              acaGTTAACCTTTGGTGGTGCACAACTGGTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCg                                                                                                       <  2:764272/139‑1 (MQ=255)              acaGTTAACCTTTGGTGGTGCACAACTGGTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCg                                                                                                       <  1:743012/139‑1 (MQ=255)                    aaCCTTTGGTGGTGCACAACTGGTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCGATCGcc                                                                                                 >  1:649490/1‑139 (MQ=255)                          tggtggTGCACAACTGGTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCg                                                                                                       <  2:47482/127‑1 (MQ=255)                          tggtggTGCACAACTGGTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCg                                                                                                       >  1:47482/1‑127 (MQ=255)                            gtggtgCACAACTGGTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCgg                                                                                         >  1:466077/1‑139 (MQ=255)                                tgCACAACTGGTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCAt                                                                                     >  1:220558/1‑139 (MQ=255)                                       cTGGTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAgg                                                                              >  2:433975/1‑139 (MQ=255)                                       cTGGTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATAGAAATGCCCGATCGCCAGGACCGGGCATTTTCAgg                                                                              >  1:464903/1‑139 (MQ=255)                                                    aTGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCa                                                                                                                              >  1:283790/1‑78 (MQ=255)                                                    aTGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCa                                                                                                                              <  2:283790/78‑1 (MQ=255)                                                    aTGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCaa                                                                 >  1:28988/1‑139 (MQ=255)                                                      gCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACt                                                               <  1:647188/139‑1 (MQ=255)                                                       ccAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTg                                                              <  2:601789/139‑1 (MQ=255)                                                       ccAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTg                                                              <  1:271419/139‑1 (MQ=255)                                                            cAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATg                                                         <  2:10637/139‑1 (MQ=255)                                                            cAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATg                                                         <  2:28988/139‑1 (MQ=255)                                                            cAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATg                                                         <  2:649490/139‑1 (MQ=255)                                                                             gttggttgATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCcca                                        >  2:665727/1‑139 (MQ=255)                                                                               tggttgATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCccacc                                      <  2:770476/139‑1 (MQ=255)                                                                                     aTGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGa                                                                             >  2:592277/1‑94 (MQ=255)                                                                                     aTGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGa                                                                             <  1:592277/94‑1 (MQ=255)                                                                                     aTGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCagaaag                                >  1:559513/1‑139 (MQ=255)                                                                                                      aCGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGgta               <  1:665727/139‑1 (MQ=255)                                                                                                                   ccGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGGTAAAGTAAATACCCa  <  1:732235/139‑1 (MQ=255)                                                                                                                          |                                                                                                                                     TTGCCATAACGCCACAGTTAACCTTTGGTGGTGCACAACTGGTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGTGGTTATCACTAGTCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGGTAAAGTAAATACCCA  >  NC_000913/3134690‑3134943

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