Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,624,281	A→G	L170P (CTG→CCG)	ytjB ←	protein Smp

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,624,281	0	A	G	100.0%	127.2 / NA	39	L170P (CTG→CCG)	ytjB	protein Smp
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (20/19);  total (20/19)

TCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTGCGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAGAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTT  >  NC_000913/4624145‑4624417                                                                                                                                         |                                                                                                                                         tCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTGCGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAGagcggc                                                                                                                                        >  1:319551/1‑139 (MQ=255)    ctctTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTGCGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAGAGCGGCAACa                                                                                                                                    >  1:428932/1‑140 (MQ=255)     tctTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTGCGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAGAGCGGCAACAt                                                                                                                                   <  2:744/140‑1 (MQ=255)         ttttCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTGCGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAGAGCGGCAACATCAgg                                                                                                                               >  1:228870/1‑140 (MQ=255)               cGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTGCGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAGAGCGGCAACATCAGGCGTaaa                                                                                                                         >  2:517646/1‑140 (MQ=255)                 cTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTGCGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAGAGCGGCAACATCAGGCGTAAAAt                                                                                                                       >  1:164713/1‑140 (MQ=255)                      ttcctcttccGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTGCGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAGAGCGGCAACATCAGGCGTAAAATGTTAg                                                                                                                  <  1:442858/140‑1 (MQ=255)                        cctcttccGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTGCGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAGAGCGGCAACATCAGGCGTAAAATGTTAgtg                                                                                                                <  2:164713/140‑1 (MQ=255)                            ttccGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGGGCGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAGAGCGGCAACATCAGGCGTAAAATGTTAGTGgtgt                                                                                                            >  1:546137/1‑140 (MQ=255)                              ccGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTGCGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAGAGCGGCAACATCAGGCGTAAAATGTTAGTGGTGTTa                                                                                                          >  2:565264/1‑140 (MQ=255)                              ccGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTGCGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACACCGATTGCCAGAGAGAGCGGCAATATCAGGCGTAAAATGTTAGTGGTGTTa                                                                                                          >  1:391081/1‑140 (MQ=255)                                ggCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGGGCGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAGAGCGGCAACATCAGGCGTAAAATGTTAGTGGTGTTATc                                                                                                        >  1:807741/1‑140 (MQ=255)                                               gCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTGCGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAGAGCGGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTc                                                                                         <  2:282189/140‑1 (MQ=255)                                               gCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTGCGTTTACCCTGTAGCAGCGTGCGGGGCAGCACTACACCGATTTCCAGTGAGAGCGGCAACCTCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTc                                                                                         >  1:564429/1‑140 (MQ=255)                                                cGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTGCGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAGAGCGGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCg                                                                                        <  1:565264/140‑1 (MQ=255)                                                     aaCAGGAAGGGCGATTGCTGCCAGCGGGTGCGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAGAGCGGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGc                                                                                   <  2:428932/140‑1 (MQ=255)                                                        aGGAAGGGCGATTGCTGCCAGCGGGTGCGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAGAGCGGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAg                                                                                <  1:529389/140‑1 (MQ=255)                                                             gggCGATTGCTGCCAGCGGGTGCGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAGAGCGGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTAt                                                                           >  1:709095/1‑140 (MQ=255)                                                                cGATTGCTGCCAGCGGGTGCGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAGAGCGGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATggg                                                                        <  2:228870/140‑1 (MQ=255)                                                                    tgctgcCAGCGGGTGCGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAGAGCGGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGc                                                                                   >  1:247334/1‑125 (MQ=255)                                                                    tgctgcCAGCGGGTGCGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAGAGCGGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGc                                                                                   <  2:247334/125‑1 (MQ=255)                                                                             cGGGTGCGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAGAGCGGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAg                                                           >  1:644451/1‑140 (MQ=255)                                                                                           cTGTAGCAGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAGAGCGGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGa                                             <  2:546137/140‑1 (MQ=255)                                                                                              tAGCAGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAGAGCGGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCg                                          >  2:401643/1‑140 (MQ=255)                                                                                                gcagcGTGCGGGTCAGCACTACACCGATTGCCAGAGAGAGCGGCAATATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGt                                        <  2:391081/140‑1 (MQ=255)                                                                                                        cGGGTCAGCACTACACCGATTGCCAGTGAGAGCGGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGttttt                                >  1:598436/1‑140 (MQ=255)                                                                                                          ggTCAGCACTACACCGATTGCCAGTGAGAGCGGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTAc                              <  1:698875/140‑1 (MQ=255)                                                                                                               gCACTACACCGATTGCCAGTGAGAGCGGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCGCa                         <  2:807741/140‑1 (MQ=255)                                                                                                                  cTACACCGATTGCCAGTGAGAGCGGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCGCAAtt                      <  2:644451/140‑1 (MQ=255)                                                                                                                     cacCGATTGCCAGTGAGAGCGGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGc                   >  1:552245/1‑140 (MQ=255)                                                                                                                     cacCGATTGCCAGTGAGAGCGGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGc                   >  1:318083/1‑140 (MQ=255)                                                                                                                          aTTGCCAGTGAGAGCGGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGCTCGAc              >  1:782547/1‑140 (MQ=255)                                                                                                                          aTTGCCAGTGAGAGCGGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGCTCGAc              >  2:564406/1‑140 (MQ=255)                                                                                                                               cAGTGAGAGCGGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATct         >  2:623552/1‑140 (MQ=255)                                                                                                                                aGTGAGAGCGGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATctg        >  2:20916/1‑140 (MQ=255)                                                                                                                                aGTGAGAGCGGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATctg        <  2:645376/140‑1 (MQ=255)                                                                                                                                aGTGAGAGCGGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATctg        <  2:709095/140‑1 (MQ=255)                                                                                                                                aGTGAGAGCGGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATctg        <  1:20916/140‑1 (MQ=255)                                                                                                                                aGTGAGAGCGGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATAGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATctg        <  1:517646/140‑1 (MQ=255)                                                                                                                                   gagaGCGGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATctgctg     <  1:401643/140‑1 (MQ=255)                                                                                                                                      agcggcAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATCTGCTGGtt  >  2:772466/1‑140 (MQ=255)                                                                                                                                         |                                                                                                                                         TCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTGCGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAGAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTT  >  NC_000913/4624145‑4624417

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