Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,410,859	C→T	L103F (CTC→TTC)	yjfJ →	PspA family protein YjfJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,410,859	0	C	T	100.0%	93.9 / NA	29	L103F (CTC→TTC)	yjfJ	PspA family protein YjfJ
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (17/12);  total (17/12)

CCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAGCTGAAAGATTTACGTGAGCGCAAAGCCAGTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATCTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACG  >  NC_000913/4410627‑4411100                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                  ccAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAGCTGAAAGATTTACGTGAGCGCAAAGCCAGTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGc                                                                                                                                                                                                                                 <  1:133555/251‑1 (MQ=255) ccAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAGCTGAAAGATTTACGTGAGCGCAAAGCCAGTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGc                                                                                                                                                                                                                                 >  2:200506/1‑251 (MQ=255)   aGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAGCTGAAAGATTTACGTGAGCGCAAAGCCAGTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCa                                                                                                                                                                                                                                >  2:113709/1‑250 (MQ=255)   aGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCGGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAGCTGAAAGATTTACGTGAGCGCAAAGCCAGTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCaa                                                                                                                                                                                                                               <  2:152172/251‑1 (MQ=255)            gcATGCTGGAAGAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAGCTGAAAGATTTACGTGAGCGCAAAGCCAGTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCg                                                                                                                                                                                                                      >  1:200990/1‑251 (MQ=255)                  tGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATCAGCTGAAAGATTTACGTGAGCGCAAAGCCAGTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGCACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGCATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTg                                                                                                                                                                                                                >  2:202950/1‑251 (MQ=255)                   ggAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAGCTGAAAGATTTACGTGAGCGCAAAGCCAGTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTg                                                                                                                                                                                                                <  1:24399/250‑1 (MQ=255)                                       cTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAGCTGAAAGATTTACGTGAGCGCAAAGCCAGTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTg                                                                                                                                                                                           >  1:254911/1‑251 (MQ=255)                                                   tCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAGCTGAAAGATTTACGTGAGCGCAAAGCCAGTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGCTTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGACGTCTCCCGTGATGGCGTGGAAAAAGCAGtt                                                                                                                                                                               >  2:89293/1‑251 (MQ=255)                                                    cGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAGCTGAAAGATTTACGTGAGCGCAAAGCCAGTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGtt                                                                                                                                                                               >  2:19823/1‑250 (MQ=255)                                                     gATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAGCTGAAAGATTTACGTGAGCGCAAAGCCAGTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTTGAAAAAGCAGTTac                                                                                                                                                                             <  1:66999/251‑1 (MQ=255)                                                            ccGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAGCTGAAAGATTTACGTGAGCGCAAAGCCAGTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTacagcgac                                                                                                                                                                       >  2:158061/1‑250 (MQ=255)                                                                                         cgGGTGAAATTAAGTCACGATAAGCTGAAAGCTTTACGTGAGCGCAAAGCAAGTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGagcagc                                                                                                                                          >  1:74200/1‑250 (MQ=255)                                                                                                       tCACGATAAGCTGAAAGATTTACGTGAGCGCAAAGCCAGTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGtt                                                                                                                            >  2:140820/1‑250 (MQ=255)                                                                                                                                     cAAATCCAGTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAAAGGTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAAGCGCGCGCCTTGAGAATTTCATTACCGATGAAGAGCAAGTGGTGTCGAATCTGGAAGTCTCCCGTGAGGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGGGCacaac                                                                                             <  2:237054/251‑1 (MQ=255)                                                                                                                                          ccAGTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCt                                                                                        >  2:133277/1‑251 (MQ=255)                                                                                                                                                    aGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGttt                                                                                                                                                 <  1:87335/184‑1 (MQ=255)                                                                                                                                                    aGCCCGAGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGCAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGACGTCTCCCGTGATGGCGTGGAACAAGCAGTTACAGAGACAGCGCAGCGTATCGCTCAGttt                                                                                                                                                 >  2:87335/1‑184 (MQ=255)                                                                                                                                                                     cGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCt                                                             <  1:198723/251‑1 (MQ=255)                                                                                                                                                                      gTTAAGCAAGACCGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCtt                                                            >  2:26773/1‑251 (MQ=255)                                                                                                                                                                             aaGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGt                                                     >  1:126727/1‑251 (MQ=255)                                                                                                                                                                                    ttAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGtt                                                    >  2:10349/1‑245 (MQ=255)                                                                                                                                                                                    ttAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGtt                                                    <  1:10349/245‑1 (MQ=255)                                                                                                                                                                                       aTCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCg                                           <  1:19823/251‑1 (MQ=255)                                                                                                                                                                                         ccGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGc                                         <  1:113709/251‑1 (MQ=255)                                                                                                                                                                                                  aTTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCtt                                <  1:140820/251‑1 (MQ=255)                                                                                                                                                                                                                              gcCTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGa     >  1:26448/1‑250 (MQ=255)                                                                                                                                                                                                                               ccTTGAGAATTTCATTACCGCTGAAGAGCAAGTGTTGTTGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTTCCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAAc   <  2:51839/251‑1 (MQ=255)                                                                                                                                                                                                                                cTTGAGAATTTCATTACAGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCAAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCACACGCTACAAGCCGAACg  >  1:554/1‑251 (MQ=255)                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                  CCAGGGCGTGCGCATGCTGGAACAGCATATTCGTGACGCTAAAGCTGAACTCGATAAAGCCGGAAAATCTCGCGTTGATCTGCTGGCGCGGGTGAAATTAAGTCACGATAAGCTGAAAGATTTACGTGAGCGCAAAGCCAGTCTGGAAGCCCGTGCGCTGGAAGCGTTAAGCAAGAACGTTAATCCGTCGTTGATTAACGAAGTTGCTGAAGAAATCGCGCGCCTTGAGAATCTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTAAAACGCCTGCAAACGCGTCAGGCCGAACG  >  NC_000913/4410627‑4411100

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