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%	153.3 / NA	42	Q35R (CAG→CGG)	ydaT	Rac prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (22/20);  total (22/20)

TCCTCATTTGATCATACCTGAAACATCAAGAGGCAAATGATTCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCA‑‑GGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCAGCAATTCTGGCAATCCTTCCGCGC  >  NC_000913/1420641‑1420932                                                                                                                                                    |                                                                                                                                                   tCCTCATTTGATCATACCTGAAACATCAAGAGGCAAATGATTCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCA‑‑gggg                                                                                                                                               >  1:390928/1‑151 (MQ=255)     cATTTGATCATACCTGAAACATCAAGAGGCAAATGATTCATGAAAATCAAGCATGAGCGCATCGAATCATTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCA‑‑gggggggc                                                                                                                                           <  2:42967/151‑3 (MQ=255)     cATTTGATCATACCTGAAACATCAAGAGGCAAATGATTCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTTCGCCA‑‑gggggggc                                                                                                                                           <  1:604963/151‑3 (MQ=255)     cATTTGATCATACCTGAAACATCAAGAGGCAAATGATTCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCA‑‑gggggg                                                                                                                                             <  1:396718/149‑1 (MQ=255)       tttGATCATACCTGAAACATCAAGAGGCAAATGATTCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCA‑‑gggggggcgg                                                                                                                                         >  1:502717/1‑147 (MQ=255)              aTACCTGAAACATCAAGAGGC‑AATGATTCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAG‑GGGGGGCGGCGAATTGc                                                                                                                                  <  2:11261/151‑1 (MQ=255)                aCCTGAAACATCAAGAGGCAAATGATTCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAG‑GGGGGGCGGCGAATTGcc                                                                                                                                 >  2:215569/1‑151 (MQ=255)                       aCATCAAGAGGCAAATGATTCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAG‑GGGGGGCGGCGAATTGCCCCTGGtt                                                                                                                          >  2:484206/1‑151 (MQ=255)                       aCATCAAGAGGCAAATGATTCATGAAAATCAAGCATGAGCACATCGAATCAGCGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAG‑GGGGGGCGGCGAATTGCCCCTGGtt                                                                                                                          <  1:523102/151‑1 (MQ=255)                         aTCAAGAGGCAAATGATTCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAG‑GGGGGGCGGCGAATTGCCCCTGGTTcc                                                                                                                        >  1:191079/1‑151 (MQ=255)                         aTCAAGAGGCAAATGATTCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAG‑GGGGGGCGGCGAATTGCCCCTGGTTcc                                                                                                                        >  1:674309/1‑151 (MQ=255)                          tCAAGAGGCAAATGATTCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAG‑GGGGGGCGGCGAATTGCCCCTGGTTcc                                                                                                                        >  1:499199/1‑150 (MQ=255)                                   aaaTGATTCATGAAAATCAAGCATGAGCACATCGAATCAGGGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAG‑GGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGAc                                                                                                              >  1:157701/1‑151 (MQ=255)                                     aTGATTCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAG‑GGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTg                                                                                                            >  2:633098/1‑151 (MQ=255)                                          tCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAG‑GGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAAc                                                                                                        >  2:577007/1‑150 (MQ=255)                                           cATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAG‑GGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACa                                                                                                       >  1:122137/1‑150 (MQ=255)                                              gAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAG‑GGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATc                                                                                                    >  1:166633/1‑150 (MQ=255)                                               aaaaTCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAG‑GGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATcag                                                                                                  <  1:638429/151‑1 (MQ=255)                                                aaaTCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAG‑GGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATcagc                                                                                                 >  1:364547/1‑151 (MQ=255)                                                  aTCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAG‑GGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATcagcag                                                                                               <  1:504159/151‑1 (MQ=255)                                                    cAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAG‑GGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGa                                                                                              <  1:62515/150‑1 (MQ=255)                                                    cAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAG‑GGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGa                                                                                              >  2:339020/1‑150 (MQ=255)                                                        cATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAG‑GGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAatat                                                                                          <  2:45970/150‑1 (MQ=255)                                                               acaTCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAG‑GGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCAcc                                                                                  >  1:549416/1‑151 (MQ=255)                                                                  tCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAG‑GGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCgtt                                                                               >  2:63428/1‑151 (MQ=255)                                                                    gAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAG‑GGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCgttgg                                                                             >  1:169951/1‑151 (MQ=255)                                                                       tCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAG‑GGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCgttggttg                                                                          >  2:66683/1‑151 (MQ=255)                                                                        cAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAG‑GGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGa                                                                         <  1:498964/151‑1 (MQ=255)                                                                         aGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAG‑GGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGa                                                                         <  1:562828/150‑1 (MQ=255)                                                                              tgttTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGt                                                                    >  1:433449/1‑151 (MQ=255)                                                                               gttTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAG‑GGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTg                                                                   >  2:216861/1‑150 (MQ=255)                                                                                     ccTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAG‑GGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGaa                                                            <  2:446978/151‑1 (MQ=255)                                                                                        cgcagcCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAG‑GGGGGGGGGCGAATTGCCCCGGGTTCCAGGCAAGGACGGGAACAATCAGCAGAATATCTACCACCGTTGGTGGAAGGGGGAAACGaaaa                                                          >  1:102502/2‑150 (MQ=255)                                                                                         gcagcCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAG‑GGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACg                                                        <  2:674288/151‑1 (MQ=255)                                                                                          cagcCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAG‑GGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAAc                                                         <  2:611305/149‑1 (MQ=255)                                                                                                     gggCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAG‑GGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGaaaaa                                            <  1:213426/151‑1 (MQ=255)                                                                                                              tGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAG‑GGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGa                                     >  2:164432/1‑149 (MQ=255)                                                                                                                      cAATGCAATTACTGAAGAATATCTGCGCCAG‑GGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATccc                           <  1:680490/151‑1 (MQ=255)                                                                                                                       aaTGCAATTACTGAAGAATATCTGCGCCAG‑GGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCa                          <  2:673544/151‑1 (MQ=255)                                                                                                                        aTGCAATTACTGAAGAATATCTGCGCCAG‑GGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCAg                         <  2:126915/151‑1 (MQ=255)                                                                                                                               tACTGAAGAATATCTGCGCCAG‑GGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCAGCAATTCt                  >  1:307280/1‑151 (MQ=255)                                                                                                                                   gaagaaTATCTGCGCCAG‑GGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCAGCAATTCTGGc               <  1:417287/150‑1 (MQ=255)                                                                                                                                        atatCTGCGCCAG‑GGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCAGCAATTCTGGCAATCCt         <  2:533891/151‑1 (MQ=255)                                                                                                                                        atatCTGCGCCAG‑GGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCAGCAATTCTGGCAATCCt         >  1:561988/1‑151 (MQ=255)                                                                                                                                        atatCTGCGCCAG‑GGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCAGCAATTCTGGCAATCCt         <  1:633117/151‑1 (MQ=255)                                                                                                                                              gcgcCAG‑GGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCAGCAATTCTGGCAATCCTTc       <  2:593315/147‑1 (MQ=255)                                                                                                                                               cgcCAG‑GGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCAGCAATTCTGGCAATCCTTCcgcgc  <  1:567743/151‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                   TCCTCATTTGATCATACCTGAAACATCAAGAGGCAAATGATTCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCA‑‑GGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCAGCAATTCTGGCAATCCTTCCGCGC  >  NC_000913/1420641‑1420932

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