Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,134,823	G→A	*231* (TAG→TAA)	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,823	0	G	A	100.0%	102.5 / NA	32	*231* (TAG→TAA)	yghT	putative ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (13/19);  total (13/19)

TTGCCATAACGCCACAGTTAACCTTTGGTGGTGCACAACTGGTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGTGGTTATCACTAGTCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGGTAAAGTAAATACCCACGCCATCA  >  NC_000913/3134690‑3134951                                                                                                                                      |                                                                                                                                 ttGCCATAACGCCACAGTTAACCTTTGGTGGTGCACAACTGGTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCaaaa                                                                                                                            <  2:20540/139‑1 (MQ=255) ttGCCATAACGCCACAGTTAACCTTTGGTGGTGCACAACTGGTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCaaaa                                                                                                                            <  1:558021/139‑1 (MQ=255)  tGCCATAACGCCACAGTTAACCTTTGGTGGTGCACAACTGGTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCaaaaa                                                                                                                           >  2:40458/1‑139 (MQ=255)       tAACGCCACAGTTAACCTTTGGTGGTGCACAACTGGTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGaa                                                                                                                      >  1:451904/1‑139 (MQ=255)              acaGTTAACCTTTGGTGGTGCACAACTGGTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCg                                                                                                               <  2:417006/139‑1 (MQ=255)              acaGTTAACCTTTGGTGGTGCACAACTGGTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCg                                                                                                               <  2:363120/139‑1 (MQ=255)                    aaCCTTTGGTGGTGCACAACTGGTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCGATCGcc                                                                                                         >  2:187250/1‑139 (MQ=255)                     aCCTTTGGTGGTGCACAACTGGTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCa                                                                                                        <  2:334424/139‑1 (MQ=255)                          tggtggTGCACAACTGGTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGAcc                                                                                                   >  1:561894/1‑139 (MQ=255)                                       cTGGTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAgg                                                                                      <  2:451904/139‑1 (MQ=255)                                          gTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAg                                                                                   >  2:484696/1‑139 (MQ=255)                                                       ccAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTg                                                                      <  2:741472/139‑1 (MQ=255)                                                            cAGCCACTGTAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATg                                                                 <  1:660398/139‑1 (MQ=255)                                                            cAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATg                                                                 <  2:461377/139‑1 (MQ=255)                                                            cAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATg                                                                 <  1:187250/139‑1 (MQ=255)                                                               ccACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGAc                                                                                                    >  1:242978/1‑101 (MQ=255)                                                               ccACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGAc                                                                                                    <  2:242978/101‑1 (MQ=255)                                                                      aaCAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTa                                                       <  1:703186/139‑1 (MQ=255)                                                                             gttggttgATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCcca                                                >  2:214563/1‑139 (MQ=255)                                                                                ggttgATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATgg                                                                                                                        <  2:319306/64‑1 (MQ=255)                                                                                ggttgATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATgg                                                                                                                        >  1:319306/1‑64 (MQ=255)                                                                                     aTGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCagaaag                                        >  2:246022/1‑139 (MQ=255)                                                                                     aTGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCagaaag                                        >  2:782730/1‑139 (MQ=255)                                                                                     aTGCAGAAAAAGCGATTACGCCTTTTATGTCCGCG‑GGGGTTATCACTAATCAAAAATGGAAATGATCGATCGCCAGGACCGGGCATTTTCAGGAAGCTTAAATCAACTGCAATGCTATCCACTACAGCACCCCagaaag                                        >  1:36818/1‑139 (MQ=255)                                                                                      tGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCagaaaga                                       >  1:487133/1‑139 (MQ=255)                                                                                               aGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTgcc                              <  2:487133/139‑1 (MQ=255)                                                                                                 cGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTgccgc                            <  1:782730/139‑1 (MQ=255)                                                                                                     tACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGgt                        <  2:196232/139‑1 (MQ=255)                                                                                                      aCGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGgta                       <  1:484696/139‑1 (MQ=255)                                                                                                                   ccGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATc                                                            <  1:529464/89‑1 (MQ=255)                                                                                                                   ccGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATc                                                            >  2:529464/1‑89 (MQ=255)                                                                                                                            gTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGGTAAAGTAAATACCCACGCCATCa  <  1:246022/139‑1 (MQ=255)                                                                                                                                      |                                                                                                                                 TTGCCATAACGCCACAGTTAACCTTTGGTGGTGCACAACTGGTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGTGGTTATCACTAGTCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGGTAAAGTAAATACCCACGCCATCA  >  NC_000913/3134690‑3134951

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