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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,420,788	0	G	.	100.0%	188.8 / NA	44	coding (105/423 nt)	ydaT	Rac prophage; protein YdaT
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (23/21);  total (23/21)

TCATACCTGAAACATCAAGAGGCAAATGATTCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCAGCAATTCTGGC  >  NC_000913/1420652‑1420919                                                                                                                                         |                                                                                                                                    tCATACCTGAAACATCAAGAGGCAAATGATTCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAgtgg                                                                                                                                  <  2:995182/140‑4 (MQ=255)    tACCTGAAACATCAAGAGGCAAATGATTCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAgggggcg                                                                                                                               >  2:50577/1‑138 (MQ=255)          aaaCATCAAGAGGCAAATGATTCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCA‑GGGGGCGGCGAAt                                                                                                                        >  2:164417/1‑140 (MQ=255)          aaaCATCAAGAGGCAAATGATTCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCA‑GGGGGCGGCGAAt                                                                                                                        >  1:910874/1‑140 (MQ=255)               tCAAGAGGCAAATGATTCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCA‑GGGGGCGGCGAATTGccc                                                                                                                   >  2:297247/1‑140 (MQ=255)                cAAGAGGCAAATGATTCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCA‑GGGGGCGGCGAATTGcccc                                                                                                                  <  1:164417/140‑1 (MQ=255)                       cAAATGATTCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCA‑GGGGGCGGCGAATTGCCCCTGGTTcc                                                                                                           >  1:627573/1‑140 (MQ=255)                            gATTCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCA‑GGGGGCGGCGAATTGCCCCTGGTTCCAGGCa                                                                                                      <  2:910874/140‑1 (MQ=255)                                        tCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCA‑GGGGGCGGCGAATTGCCCCTGGTTCCAGGCa                                                                                                      >  2:54888/1‑128 (MQ=255)                                        tCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCA‑GGGGGCGGCGAATTGCCCCTGGTTCCAGGCa                                                                                                      <  1:54888/128‑1 (MQ=255)                                           aGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCA‑GGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATc                                                                                       <  1:21581/140‑1 (MQ=255)                                             cATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCA‑GGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATcag                                                                                     <  1:297247/140‑1 (MQ=255)                                               tGAGCACATCGAATCAGTGTTGTTTGCGCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCA‑GGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATcagca                                                                                   <  2:746103/140‑1 (MQ=255)                                                 aGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCA‑GGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGa                                                                                 <  1:656140/140‑1 (MQ=255)                                                   cacaTCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCA‑GGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAat                                                                               <  1:814121/140‑1 (MQ=255)                                                    acaTCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCA‑GGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAata                                                                              <  1:865537/140‑1 (MQ=255)                                                      aTCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCA‑GGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATcagcag                                                                                  <  2:385610/134‑1 (MQ=255)                                                      aTCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCA‑GGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATcagcag                                                                                  >  1:385610/1‑134 (MQ=255)                                                      aTCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCA‑GGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAatatt                                                                            >  1:239461/1‑139 (MQ=255)                                                      aTCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCA‑GGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAatatt                                                                            >  1:222185/1‑139 (MQ=255)                                                           aTCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCA‑GGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCa                                                                       >  1:666477/1‑140 (MQ=255)                                                           aTCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCA‑GGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCa                                                                       >  2:989233/1‑140 (MQ=255)                                                            tCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCA‑GGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCAc                                                                      <  2:842359/140‑1 (MQ=255)                                                                tgttgtTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCA‑GGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATcagca                                                                                   <  2:899821/123‑1 (MQ=255)                                                                tgttgtTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCA‑GGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATcagca                                                                                   >  1:899821/1‑123 (MQ=255)                                                                tgttgtTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCA‑GGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCg                                                                    >  1:418644/1‑138 (MQ=255)                                                                tgttgtTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCA‑GGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCg                                                                    <  2:418644/138‑1 (MQ=255)                                                                    gttTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCA‑GGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCgt                                                                   <  1:482748/135‑1 (MQ=255)                                                                    gttTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCA‑GGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCgt                                                                   >  2:482748/1‑135 (MQ=255)                                                                     tttGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCA‑GGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCgttggttg                                                             >  1:752404/1‑140 (MQ=255)                                                                       tGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCA‑GGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGaa                                                           <  1:989233/140‑1 (MQ=255)                                                                          ccTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCA‑GGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAgg                                                        >  1:570082/1‑140 (MQ=255)                                                                          ccTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCA‑GGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAgg                                                        >  1:860670/1‑140 (MQ=255)                                                                            tAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCA‑GGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTg                                                      >  1:537309/1‑140 (MQ=255)                                                                             agcagcCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCA‑GGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGa                                                     >  2:1045877/1‑140 (MQ=255)                                                                                agcCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCA‑GGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATAT‑TATCACCGTTGGTTGAAAGGTGAAACg                                                 <  2:239461/140‑1 (MQ=255)                                                                                agcCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCA‑GGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATAT‑TATCACCGTTGGTTGAAAGGTGAAACg                                                 <  2:222185/140‑1 (MQ=255)                                                                                   cGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCA‑GGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGaa                                               >  2:331674/1‑140 (MQ=255)                                                                                               ggcaTGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCA‑GGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAgaga                                   <  2:860670/140‑1 (MQ=255)                                                                                                gcaTGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCA‑GGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGaa                                  >  1:23077/1‑140 (MQ=255)                                                                                                              tGCAATTACTGAAGAATATCTGCGCCA‑GGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCt                    >  1:604673/1‑140 (MQ=255)                                                                                                                 aaTTACTGAAGAATATCTGCGCCA‑GGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGAt                 <  2:383726/140‑1 (MQ=255)                                                                                                                    tACTGAAGAATATCTGCGCCA‑GGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATccc              >  1:807883/1‑140 (MQ=255)                                                                                                                          agaaTATCTGCGCCA‑GGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCAGCAAt        <  1:904894/140‑1 (MQ=255)                                                                                                                           gaaTATCTGCGCCA‑GGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCAGCAAtt       <  2:511777/140‑1 (MQ=255)                                                                                                                                tCTGCGCCA‑GGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCAGCAATTCTGGc  >  1:433473/1‑140 (MQ=255)                                                                                                                                         |                                                                                                                                    TCATACCTGAAACATCAAGAGGCAAATGATTCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCAGCAATTCTGGC  >  NC_000913/1420652‑1420919

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