Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,351,696	T→C	K22R (AAG→AGG)	ycjD ←	DUF559 domain‑containing protein YcjD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,351,696	0	T	C	100.0%	110.2 / NA	33	K22R (AAG→AGG)	ycjD	DUF559 domain‑containing protein YcjD
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (14/19);  total (14/19)

ATGCTGCCCACCATCCAGCTCAACGACTACACGCGCCGAGCAGCAAGCAAAATCGAGAATGTAGCTCCCCACTGGATGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGC  >  NC_000913/1351563‑1351832                                                                                                                                      |                                                                                                                                         aTGCTGCCCACCATCCAGCTCAACGACTACACGCGCCGAGCAGCAAGCAAAATCGAGAATGTAGCTCCCCACTGGATGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCCTTCGtt                                                                                                                                    >  2:389822/1‑140 (MQ=255)    ctgcCCACCATCCAGCTCAACGACTACACGCGCCGAGCAGCAAGCAAAATCGAGAATGTAGCTCCCCACTGGATGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCCTTCGTTCCt                                                                                                                                 >  1:561600/1‑140 (MQ=255)       cccACCATCCAGCTCAACGACTACACGCGCCGAGCAGCAAGCAAAATCGAGAATGTAGCTCCCCACTGGATGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCCTTCGTTCCTGCa                                                                                                                              <  1:389822/140‑1 (MQ=255)                   cTCAACGACTACACGCGCCGAGCAGCAAGCAAAATCGAGAATGTAGCTCCCCACTGGATGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCCTTCGTTCCTGCAAAGTGAGAtt                                                                                                                    <  2:403514/138‑1 (MQ=255)                   cTCAACGACTACACGCGCCGAGCAGCAAGCAAAATCGAGAATGTAGCTCCCCACTGGATGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCCTTCGTTCCTGCAAAGTGAGAtt                                                                                                                    >  1:403514/1‑138 (MQ=255)                        cGACTACACGCGCCGAGCAGCAAGCAAAATCGAGAATGTAGCTCCCCACTGGATGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCCTTCGTTCCTGCAAAGTGAGAtt                                                                                                                    >  1:541179/1‑133 (MQ=255)                        cGACTACACGCGCCGAGCAGCAAGCAAAATCGAGAATGTAGCTCCCCACTGGATGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCCTTCGTTCCTGCAAAGTGAGAtt                                                                                                                    <  2:541179/133‑1 (MQ=255)                                  cgcCGAGCAGCAAGCAAAATCGAGAATGTAGCTCCCCACTGGATGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCCTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCAcgcg                                                                                                   <  1:478091/140‑1 (MQ=255)                                   gcCGAGCAGCAAGCAAAATCGAGAATGTAGCTCCCCACTGGATGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCCTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACgcgc                                                                                                  >  2:299971/1‑140 (MQ=255)                                    ccGAGCAGCAAGCAAAATCGAGAATGTAGCTCCCCACTGGATGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCCTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCa                                                                                                 <  1:299971/140‑1 (MQ=255)                                     cGAGCAGCAAGCAAAATCGAGAATGTAGCTCCCCACCGGATGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCCTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCAt                                                                                                >  1:60758/1‑140 (MQ=255)                                       agcagcAAGCAAAATCGAGAATGTAGCTCCCCACTGGATGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCCTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCAttt                                                                                              >  2:558580/1‑140 (MQ=255)                                         cagcaAGCAAAATCGAGAATGTAGCTCCCCACTGGATGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCCTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGa                                                                                            <  1:543021/140‑1 (MQ=255)                                            caagcaaAATCGAGAATGTAGCTCCCCACTGGATGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCCTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGAttt                                                                                         <  1:558580/140‑1 (MQ=255)                                                      cGAGAATGTAGCTCCCCACTGGATGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCCTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGAt                                                                                           >  1:140429/1‑128 (MQ=255)                                                      cGAGAATGTAGCTCCCCACTGGATGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCCTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGAt                                                                                           <  2:140429/128‑1 (MQ=255)                                                       gagaATGTAGCTCCCCACTGGATGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCCTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTtatcta                                                                              <  1:796189/140‑1 (MQ=255)                                                       gagaATGTAGCTCCCCACCGGATGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCCTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTtatcta                                                                              <  2:60758/140‑1 (MQ=255)                                                            tGTAGCTCCCCACTGGATGTTGACGGCAAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCCTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCAcc                                                                         <  2:407664/140‑1 (MQ=255)                                                               aGCTCCCCACTGGATGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCCTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAAttt                                                                                    >  2:570470/1‑126 (MQ=255)                                                               aGCTCCCCACTGGATGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCCTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAAttt                                                                                    <  1:570470/126‑1 (MQ=255)                                                                      cACTGGATGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCCTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAc                                                               <  2:675816/140‑1 (MQ=255)                                                                        cTGGATGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCCTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACaa                                                             >  2:588229/1‑140 (MQ=255)                                                                                tGACGGCGAAATTTTAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCCTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAAt                                                     >  2:55129/1‑140 (MQ=255)                                                                                 gACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCCTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATc                                                    >  2:345965/1‑140 (MQ=255)                                                                                   cGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCCTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCtt                                                  >  2:132719/1‑140 (MQ=255)                                                                                      cGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCCTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTg                                               <  2:181671/140‑1 (MQ=255)                                                                                      cGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCCTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTg                                               <  1:354490/140‑1 (MQ=255)                                                                                       gAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCCTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGt                                              <  2:561600/140‑1 (MQ=255)                                                                                                aaTCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCCTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAg                                     <  1:55129/140‑1 (MQ=255)                                                                                                     cTAAAACGTCGGCTGCGAAGATATCGCCAGAGCCTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGaa                                <  2:132029/140‑1 (MQ=255)                                                                                                              cGGCTGCGAAGATATCGCCAGAGCCTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGa                       >  1:614981/1‑140 (MQ=255)                                                                                                                                   agCCTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGc  <  1:284268/140‑1 (MQ=255)                                                                                                                                      |                                                                                                                                         ATGCTGCCCACCATCCAGCTCAACGACTACACGCGCCGAGCAGCAAGCAAAATCGAGAATGTAGCTCCCCACTGGATGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGC  >  NC_000913/1351563‑1351832

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