Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,420,793	(G)6→7	coding (110/423 nt)	ydaT →	Rac prophage, uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,420,787	1	.	G	100.0%	151.8 / NA	43	Q35R (CAG→CGG)	ydaT	Rac prophage, uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (22/21);  total (22/21)

TCATACCTGAAACATCAAGAGGCAAATGATTCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCA‑GGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCAGCAATTCT  >  NC_000913/1420652‑1420916                                                                                                                                         |                                                                                                                                  tCATACCTGAAACATCAAGAGGCAAATGATTCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCA‑ggg                                                                                                                                >  1:495268/1‑139 (MQ=255)           aaCATCAAGAGGCAAATGATTCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGa                                                                                                                       >  2:581259/1‑139 (MQ=255)           aaCATCAAGAGGCAAATGATTCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCA‑gggggg                                                                                                                             <  2:791442/132‑1 (MQ=255)           aaCATCAAGAGGCAAATGATTCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCA‑gggggg                                                                                                                             >  1:791442/1‑132 (MQ=255)               tCAAGAGGCAAATGATTCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCA‑gggggg                                                                                                                             >  1:276292/1‑128 (MQ=255)               tCAAGAGGCAAATGATTCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCA‑gggggg                                                                                                                             <  2:276292/128‑1 (MQ=255)                    agGCAAATGATTCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCt                                                                                                              <  2:907158/139‑1 (MQ=255)                    agGCAAATGATTCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCt                                                                                                              <  1:224201/139‑1 (MQ=255)                       cAAATGATTCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGt                                                                                                           <  2:495268/139‑1 (MQ=255)                               tCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCa                                                                                                   <  1:106021/139‑1 (MQ=255)                               tCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCa                                                                                                   <  1:1068297/139‑1 (MQ=255)                                 aTGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAg                                                                                                 <  1:581259/139‑1 (MQ=255)                                      aaTCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTg                                                                                            <  2:134585/139‑1 (MQ=255)                                                 aGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATcagc                                                                                 >  2:49144/1‑139 (MQ=255)                                                    acaTCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGa                                                                              <  1:49144/139‑1 (MQ=255)                                                      aTCGAATCAGTGTTGTTTGCCCTAACAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAat                                                                            >  2:336743/1‑139 (MQ=255)                                                                tgttgtTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCAcc                                                                  >  2:153348/1‑139 (MQ=255)                                                                tgttgtTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCAcc                                                                  >  1:989757/1‑139 (MQ=255)                                                                  ttgttTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCgt                                                                >  1:788408/1‑139 (MQ=255)                                                                          ccTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGtt                                                                                                          >  1:705810/1‑89 (MQ=255)                                                                          ccTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGtt                                                                                                          <  2:705810/89‑1 (MQ=255)                                                                          ccTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGc                                                                                                    <  1:662026/95‑1 (MQ=255)                                                                          ccTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGc                                                                                                    >  2:662026/1‑95 (MQ=255)                                                                             agcagcCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAAt                                                                                     >  2:1000283/1‑107 (MQ=255)                                                                             agcagcCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAAt                                                                                     <  1:1000283/107‑1 (MQ=255)                                                                                agcCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGa                                                  <  2:788408/139‑1 (MQ=255)                                                                                                     ggTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATcag                                                                                  <  1:977683/86‑1 (MQ=255)                                                                                                     ggTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATcag                                                                                  >  2:977683/1‑86 (MQ=255)                                                                                                          ccAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATcag                                                                                  <  2:390436/81‑1 (MQ=255)                                                                                                          ccAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATcag                                                                                  >  1:390436/1‑81 (MQ=255)                                                                                                          ccAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTc                        >  2:940054/1‑139 (MQ=255)                                                                                                          ccAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTc                        >  1:363498/1‑139 (MQ=255)                                                                                                          ccAATGCAATTACTG‑AGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCAcc                                                                  <  1:692112/96‑1 (MQ=255)                                                                                                          ccAATGCAATTACTG‑AGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCAcc                                                                  >  2:692112/1‑96 (MQ=255)                                                                                                              tGCAATTACTGAAGTATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGaa                                                                             <  2:627541/82‑1 (MQ=255)                                                                                                              tGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGa                                                                                          >  2:320919/1‑69 (MQ=255)                                                                                                              tGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGa                                                                                          <  1:320919/69‑1 (MQ=255)                                                                                                              tGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATcagca                                                                                >  2:677915/1‑79 (MQ=255)                                                                                                              tGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATcagca                                                                                <  1:677915/79‑1 (MQ=255)                                                                                                              tGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGaa                                                                             >  1:627541/1‑82 (MQ=255)                                                                                                                     aCTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATc             >  1:948499/1‑139 (MQ=255)                                                                                                                        gaagaaTATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCa          >  2:70704/1‑139 (MQ=255)                                                                                                                        gaagaaTATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCa          <  2:948499/139‑1 (MQ=255)                                                                                                                           gaaTATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAg                                                                                                 >  2:177652/1‑49 (MQ=39)                                                                                                                           gaaTATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAg                                                                                                 <  1:177652/49‑1 (MQ=39)                                                                                                                                tCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCAGCAATTCt  >  1:814301/1‑139 (MQ=255)                                                                                                                                   gcgcCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCgt                                                                >  1:444577/1‑74 (MQ=255)                                                                                                                                   gcgcCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCgt                                                                <  2:444577/74‑1 (MQ=255)                                                                                                                                         |                                                                                                                                  TCATACCTGAAACATCAAGAGGCAAATGATTCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCA‑GGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCAGCAATTCT  >  NC_000913/1420652‑1420916

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