Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,568,689	T→C	P45P (CCA→CCG)	yddW ←	putative lipoprotein YddW

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,568,689	0	T	C	100.0%	173.5 / NA	50	P45P (CCA→CCG)	yddW	putative lipoprotein YddW
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (23/27);  total (23/27)

TTGTTGTACACGGGCCCGGCTGGTGGGGTTACTAATGTTAACCGAGGAAACCGGTGGCCAGTCGAGCCGAGAAACCGTGGCCAGCCAGATGCCACGCATCGGTTGTGACGATTGTT‑‑GCGTCGTGGCTGGTGGCTTTGAACCTGCTGGTGGTGTCACCATGGACTCTGGAGGCGTGCTTTTACAACTACACAGTAAAAGTGCCAGTGCAACTAGTATCGCTGGTCTTCTAATCGTTAATTTCTTGTTTCGGGAGCAGATA  >  NC_000913/1568560‑1568818                                                                                                                                    |                                                                                                                                  ttgttgtACACGGGCCCGGCTGGTGGGGTTACTAATGTTAACCGAGGAAACCGGTGGCCAGTCGAGCCGAGAAACCGTGGCCAGCCAGATGGCACGCATCGGTTATCACGATTGTT‑‑GCGTCGTGGCTggcggc                                                                                                                                >  2:290146/1‑133 (MQ=255)           cGGGCCCGGCTGGTGGGGTTACTAATGTTAACCGAGGAAACCGGTGGCCAGTCGAGCCGAGAAACCGTGGCCAGCCAGATGCCACGCATCGGTTGTGACGATTGTT‑‑GCGTGGTGACTAGCGGCTTTGAACctg                                                                                                                      >  1:246135/1‑133 (MQ=255)           cGGGCCCGGCTGGTGGGGTTACTAATGTTAACCGAGGAAACCGGTGGCCAGTCGAGCCGAGAAACCGTGGCCAGCCAGATGCCACGCATCGGTTGTGACGATTGTT‑‑GCGTCGTGGCTGGCGGCTTTGAACctg                                                                                                                      <  1:447377/133‑1 (MQ=255)           cGGGCCCGGCTGGTGGGGTTACTAATGTTAACCGAGGAAACCGGTGGCCAGTCGAGCCGAGAAACCGTGGCCAGCCAGATGCCACGCATCGGTTGTGACGATAGTT‑‑GCGTTGTGACTCGCGGCTATGAACctg                                                                                                                      >  1:292203/1‑133 (MQ=255)            gggCCGGGCTGGTGGGGTTACTTATGTTAGCGGAGGAGACCGGTGGCGAGTCGAGCCGAGAAACCGTGGCCAGCCAGATGCCACGCATCGGTTGTGAGGATTGTT‑‑GCGTCGTGGCTGGCGGCTTTGAACctgc                                                                                                                     <  2:58434/133‑1 (MQ=255)            gggCCCGGCTGGTGGGGTTACTAATGTTAACCGAGGAAACCGGTGGCCAGTCGAGCCGAGAAACCGTGGCCAGCCAGATGCCACGCATCGGTTGTGACGATTGTT‑‑GCGTCGTGGCTGGCGGCTTTGAACctgc                                                                                                                     >  1:373107/1‑133 (MQ=255)              gCCCGGCTGGTGGGGTTACTAATGTTAACCGAGGAAACCGGTGGCCAGTCGAGCCGAGAAACCGTGGCCAGCCAGATGCCACGCATCGGTTGTGACGATTGTT‑‑GCGTCGTGGCTGGCGGCTTTGAACctgctg                                                                                                                   >  1:622623/1‑133 (MQ=255)              gCCCGGCTGGTGGGGTTACTAATGTTAACCGAGGAAACCGGTGGCCAGTCGAGCCGAGAAACCGTGGCCAGCCAGATGCCACGCATCGGTTGTGACGATTGTT‑‑GCGTCGTGGCTGGCGGCTTTGAACctgctg                                                                                                                   >  1:423422/1‑133 (MQ=255)              gCCCGGCTGGTGGGGTTACTAATGTTAACCGAGGAAACCGGTGGCCAGTCGAGCCGAGAAACCGTGGCCAGCCAGATGCCACGCATCGGTTGTGACGATTGTG‑‑GCGTCGTGGCTAGCGGCTTTGGACctgctg                                                                                                                   >  1:296682/1‑133 (MQ=255)                            gTTACTAATGTTAACCGAGGAAACCGGTGGCCAGTCGAGCCGAGAAACCGTGGCCAGCCAGATGCCACGCATCGGTTGTGACGATTGTT‑‑GCGTCGTGGCTGGCGGCTTTGAACCTGCTGGTGGTGTCACCATg                                                                                                     <  2:423422/133‑1 (MQ=255)                            gTTACTAATGTTAACCGAGGAAACCGGTGGCCAGTCGAGCCGAGAAACCGTGGCCAGCCAGATGCCACGCATCGGTTGTGACGATTGTT‑‑GCGTCGTGGCTGGCGGCTTTGAACCTGCTGGTGGTGTCACCATg                                                                                                     <  2:531559/133‑1 (MQ=255)                             ttACTAATGTTAACCGAGGAAACCGGTGGCCAGTCGAGCCGAGAAACCGTGGCCAGCCAGATGCCACGCATCGGTTGTGACGATTGTTGCGCGTCGTGGCTGGCGGCTTTGAACCTGCTGGTGGTGTCACCAt                                                                                                      >  2:142985/1‑133 (MQ=255)                              tACTAATGTTAACCGAGGAAACCGGTGGCCAGTCGAGCCGAGAAACCGTGGCCAGCCAGATGCCACGCATCGGTTGTGACGATTGTT‑‑GCGTCGTGGCTGGCGGCTTTGAACCTGATGGTGGTGTGACCATGGa                                                                                                   >  1:264292/1‑133 (MQ=255)                                     gTTAACCGAGGAAACCGGTGGCCAGTCGAGCCGAGAAACCGTGGCCAGCCAGATGCCACGCATCGGTTGTGACGATTGTT‑‑GCGTCGTGGCTGGCGGCTTTGAACCTGCTGGTGGTGTCACCATg                                                                                                     >  1:478107/1‑124 (MQ=255)                                     gTTAACCGAGGAAACCGGTGGCCAGTCGAGCCGAGAAACCGTGGCCAGCCAGATGCCACGCATCGGTTGTGACGATTGTT‑‑GCGTCGTGGCTGGCGGCTTTGAACCTGCTGGTGGTGTCACCATg                                                                                                     <  2:478107/124‑1 (MQ=255)                                       tAACCGAGGAAACCGGTGGCCAGTCGAGCCGAGAAACCGTGGCCAGCCAGATGCCACGCATCGGTTGTGACGATTGTT‑‑GCGTCGTGGCTGGCGGCTTTGAACCTGCCGGTGGTGTCACCATGGACTCTGGAgg                                                                                          >  1:311313/1‑133 (MQ=255)                                          ccGAGGAAACCGGTGGCCAGTCGAGCCGAGAAACCGTGGCCAGCCAGATGCCACGCATCGGTTGTGACGATTGTT‑‑GCGTCGTGGCTGGCGGCTTTGAACCTGCTGGTGGTGTCACCATGGACTCTGGAGGCGt                                                                                       <  1:404621/133‑1 (MQ=255)                                          ccGAGGAAACCGGTGGCCAGTCGAGCCGAGAAACCGTGGCCAGCCAGATGCCACGCATCGGTTGTGACGATTGTT‑‑GCGTCGTGGCTGGCGGCTTTGAACCTGCTGGTGGTGTCACCATGGACTCTGGAGGCGt                                                                                       >  2:367510/1‑133 (MQ=255)                                           cGAGGAAACCGGTGGCCAGTCGAGCCGAGAAACCGTGGCCAGCCAGATGCCACGCATCGGTTGTGACGATTGTT‑‑GCGTCGTGGCTGGCGGCTTTGAACCTGCTGGTGGTGTCACCATGGACTCTGGAGGCGTg                                                                                      >  2:411741/1‑133 (MQ=255)                                             aGGAAACCGGTGGCCAGTCGAGCCGAGAAACCGTGGCCAGCCAGATGCCACGCATCGGTTGTGACGATTGTT‑‑GCGTCGTGGCTGGCGGCTTTGAACCTGCTGGTGGTGTCACCATGGACTCTGGAGGCGt                                                                                       <  1:497031/130‑1 (MQ=255)                                             aGGAAACCGGTGGCCAGTCGAGCCGAGAAACCGTGGCCAGCCAGATGCCACGCATCGGTTGTGACGATTGTT‑‑GCGTCGTGGCTGGCGGCTTTGAACCTGCTGGTGGTGTCACCATGGACTCTGGAGGCGt                                                                                       >  2:497031/1‑130 (MQ=255)                                             aGGAAACCGGTGGCCAGTCGAGCCGAGAAACCGTGGCCAGCCAGATGCCACGCATCGGTTGTGACGATTGTT‑‑GCGTCGTGGCTGGCGGCTTTGAACCTGCTGGTGGTGTCACCATGGACTCTGGAGGCGTGCt                                                                                    >  1:79383/1‑133 (MQ=255)                                             aGGAAACCGGTGGCCAGTCGAGCCGAGAAACCGTGGCCAGCCAGATGCCACGCATCGGTTGTGACGATTGTA‑‑GCGTCGTGGCTGGCGGCTTTGAACCTGCTGGTGGTGTCACCATGGACTCTGGAGGCGTGCt                                                                                    >  1:35895/1‑133 (MQ=255)                                              ggAAACCGGTGGCCAGTCGAGCCGAGAAACCGTGGCCAGCCAGATGCCACGCATCGGTTGTGACGATTGTT‑‑GCGTCGTGGCTGGCGGCTTTGAACCTGCTGGTGGTGTCACCATGGACTCTGGAGGCGTGCtt                                                                                   >  1:84853/1‑133 (MQ=255)                                              ggAAACCGGTGGCCAGTCGAGCCGAGAAACCGTGGCCAGCCAGATGCCACGCATCGGTTGTGACGATTGTT‑‑GCGTCGTGGCTGGCGGCTTTGAACCTGCTGGTGGTGTCACCATCGACTCTGGACGCGTGCtt                                                                                   >  1:257778/1‑133 (MQ=255)                                               gAAACCGGTGGCCAGTCGAGCCGAGAAACCGTGGCCAGCCAGATGCCACGCATCGGTTGTGACGATTGTT‑‑GCGTCGTGGCTGGCGGCTTTGAACCTGCTGGTGGTGTCACCATGGACTCTGGAGGCGTGCttt                                                                                  <  1:364004/133‑1 (MQ=255)                                               gAAACCGGCGGCCAGACGTGCCAAGAAACCTTGGCCAGCCAGATGCCACGCATCGGTTGTGACGATTGTT‑‑GCGTCGTGGCTGGCGGCTTTGAACCTGCTGGTGGTGTCACCATGGACTCTGGAGGCGTGCttt                                                                                  <  1:290146/133‑1 (MQ=255)                                                    cGGTGGCCAGTCGAGCCGAGAAACCGTGGCCAGCCAGATGCCACGCATCGGTTGTGACGATTGTT‑‑GCGTCGTGGCTGGCGGCTTTGAACCTGCTGGTGGTGTCACCATGGACTCTGGAGGCGTGCTTTTACaa                                                                             <  1:411741/133‑1 (MQ=255)                                                    cGGTGGCCAGTCGAGCCGAGAAACCGTGGCCAGCCAGATGCCACGCATCGGTTGTGACGATTGTT‑‑GCGTCGTGGCTGGCGGCTTTGAACCTGCTGGTGGTGTCACCATGGACTCTGGAGGCGTGCTTTTACaa                                                                             <  1:566286/133‑1 (MQ=255)                                                     ggcggCCAGTCGAGCCGAGAAACCGTGGCCAGCCAGATGCCACGCATCGGTTGTGACGATTGTT‑‑GCGTCGTGGCTGGCGGCTTTGAACCTGCTGGTGGTGTCACCATGGACTCTGGAGGCGTGCTTTTACAAc                                                                            <  2:385995/130‑1 (MQ=255)                                                         gCCAGTCGAGCCGAGAAACCGTGGCCAGCCAGATGCCACGCATCGGTTGTGACGATTGTT‑‑GCGTCGTGGCTGGCGGCTTTGAACCTGCTGGTGGTGTCACCATGGACTCTGGAGGCGTGCTTTTACAACTaca                                                                        <  2:559113/133‑1 (MQ=255)                                                          ccAGTCGAGCCGAGAAACCGTGGCCAGCCAGATGCCACGCATCGGTTGTGACGATTGTT‑‑GCGTCGTGGCTGGCGGCTTTGAACCTGCTGGTGGTGTCACCATGGACTCTGGAGGCGTGCTTTTACAACTacac                                                                       <  2:414582/133‑1 (MQ=255)                                                            aGTCGAGCCGAGAAACCGTGGCCAGCCAGATGCCACGCATCGGTTGTGACGATTGTT‑‑GCGTCGTGGCTGGCGGCTTTGAACCTGCTGGTGGTGTCACCATGGACTCTGGAGGCGTGCTTTTACAACTACACAg                                                                     >  1:356641/1‑133 (MQ=255)                                                             gTCGAGCCGAGAAACCGTGGCCAGCCAGATGCCACGCATCGGTTGTGACGATTGTT‑‑GCGTCGTGGCTGGCGGCTTTGAACCTGCTGGTGGTGTCACCATGGACTCTGGAGGCGTGCTTTTACAACTACACAGt                                                                    >  2:575358/1‑133 (MQ=255)                                                              tCGAGCCGAGAAACCGTGGCCAGCCAGATGCCACGCATCGGTTGTGACGATTGTT‑‑GCGTCGTGGCTGGCGGCTTTGAACCTGCTGGTGGTGTCACCATGGACTCTGGAGGCGTGCTTTTACAAc                                                                            >  2:41874/1‑124 (MQ=255)                                                              tCGAGCCGAGAAACCGTGGCCAGCCAGATGCCACGCATCGGTTGTGACGATTGTT‑‑GCGTCGTGGCTGGCGGCTTTGAACCTGCTGGTGGTGTCACCATGGACTCTGGAGGCGTGCTTTTACAAc                                                                            <  1:41874/124‑1 (MQ=255)                                                                gAGCCGAGAAACCGTGGCCAGCCAGATGCCACGCATCGGTTGTGACGATTGTT‑‑GCGTCGTGGCTGGCGGCTTTGAACCTGCTGGTGGTGTCACCATGGACTCTGGAGGCGTGCTTTTACAACTACACAGTaaa                                                                 <  2:72682/133‑1 (MQ=255)                                                                 aGCCGAGAAACCGTGGCCAGCCAGATGCCACGCATCGGTTGTGACGATTGTT‑‑GCGTCGTGGCTGGCGGCTTTGAACCTGCTGGTGGTGTCACCATGGACTCTGGAGGCGTGCTTTTACAACTACACAGTaaaa                                                                >  2:32231/1‑133 (MQ=255)                                                                     gagaAACCGTTGCCATCCAGATGCCACGAATCTGTTGTGACGATTGTT‑‑GCGTCGTGGCTGGCGGCTTTGAACCTGCTGGTGGTGTCACCATGGACTCTGGAGGCGTGCTTTTACAACTACACAGTAAAAGTGc                                                            <  2:257778/133‑1 (MQ=255)                                                                      agaAACCGTGGCCAGCCAGATGCCACGCATCGGTTGTGACGATTGTTGCGCGTCGTGGCTGGCGGCTTTGAACCTGCTGGTGGTGTCACCATGGACTCTGGAGGCGTGCTTTTACAACTACACAGTAAAAGTg                                                             <  1:142985/133‑1 (MQ=255)                                                                         aaCCGTGGCCAGCCAGATGCCACGCATCGGTTGTGACGATTGTT‑‑GCGTCGTGGCTGGCGGCTTTGAACCTGCTGGTGGTGTCACCATGGACTCTGGAGGCGTGCTTTTACAACTGCACAGTAAAAGTGCCAGt                                                        <  2:79383/133‑1 (MQ=255)                                                                         aaCCGTGGCCAGCCAGATGCCACGCATCGGTTGTGACGATTGTT‑‑GCGTCGTGGCTGGCGGCTTTGAACCTGCTGGTGGTGTCACCATGGACTCTGGAGGCGTGCTTTTACAACTACACAGTAAAAGTGCCAGt                                                        <  1:575358/133‑1 (MQ=255)                                                                                       agATGCCACGCATCGGTTGTGACGTTTGAT‑‑GCATCGTGGCTGGCGGCTTTGAACCTGCTGGTGGTGTCACCATGGACTCTGGAGGCGTGCTTTTACAACTACACAGTAAAAGTGCCAGTGCAACTAGTATCGc                                          <  2:264292/133‑1 (MQ=255)                                                                                               cGCATCGGTTGTGACGATTGTT‑‑GCGTCGTGGCTGGCGGCTTTGAACCTGCTGGTGGTGTCACCATGGACTCTGGAGGCGTGCTTTTACAACAACACAGTAAAAGCGCCAGTGCAACTTGTATTGCTGGtcttc                                  >  2:272947/1‑133 (MQ=255)                                                                                                   tCGGTTGTGACGATTGTT‑‑GCGTCGTGGCTGGCGGCTTTGAACCTGCTGGTGGTGTCACCATGGACTCTGGAGGCGTGCTTTTACAACTACACAGTAAAAGTGCCAGTGCAACTAGTATCGCTGGTCTTCTAAt                              >  1:95636/1‑133 (MQ=255)                                                                                                                 tgtt‑‑gCGTCGTGGCTGGCGGCTTTGAACCTGCTGGTGGTGTCACCATGGACTCTGGAGGCGTGCTTTTACAACTACACAGTAAAAGTGCCAGTGCAACTAGTTTCGCTGGTCTTATAATCGTTACTTTCTTGt                >  1:305571/1‑133 (MQ=255)                                                                                                                        cgtcgtGGCTGGCGGTTTTGAACTTGCTCGTGGTGTCACCAAGGACTCTGGAGGCGTGCTTTTACAACTACACAGTAAAAGTGCCAGTGCAACTAGTATCGCTGGTCTTCTAATCGTTAATTTCTTGTTTCgg           <  2:292203/133‑1 (MQ=255)                                                                                                                         gtcgtGGCTGTCGGCGTTGAACCTGCTGCTGGTGTCACCATGGACTCTGGAGGCGTGCTTTTACAACTACACAGTAAAAGTGCCAGTGCAACTAGTATCGCTGGTCTTCTAATCGTTAATTTCTTGTTTCggg          <  2:246135/133‑1 (MQ=255)                                                                                                                         ggcgtGGCTCGCGGCTTTGAACCTGCCGTTGGTGTCACCATCGACTCTGGAGGCGTGCTTTTACAACTACACAGTAAAAGTGCCAGTGCAACTAGTATCGCTGGTCTTCTAATCGTTAATTTCTTGTTTCggg          <  2:294889/131‑1 (MQ=255)                                                                                                                           cgtGGCTGGCCGCTTCGAACCTGCGGGTAGTGCCACCATGGCCTCTGGAGGCGTGCTTTTACAACTACACAGTAAAAGTGCCAGTGCAACTAGTATCGCTGGTCTTCTAATCGTTAATTTCTTGTTTCGGGAg        <  2:164325/133‑1 (MQ=255)                                                                                                                            gtGGCTGGCTGCTTTTAACCTGCCAGTGGTGTAACCATGGAATCTGGAGGCGTGCTTTTACAACTACACAGTAAAAGTGCCAGTGCAACTAGTATCGCTGGTCTTCTAATCGTTAATTTCTTGTTTCGGGAGc       <  2:305571/133‑1 (MQ=255)                                                                                                                            gtGGCTGGCGGCTTTGAACCTGCTGGTGGTGTCACCATGGACTCTGGAGGCGTGCTTTTACAACTACACAGTAAAAGTGCCAGTGCAACTAGTATCGCTGGTCTTCTAATCGTTAATTTCTTGTTTCGGGAGc       <  1:396761/133‑1 (MQ=255)                                                                                                                                 tggcggCTTTGAACCTGCTGGTGGTGTCACCATGGACTCTGGAGGCGTGCTTTTACAACTACACAGTAAAAGTGCCAGTGCAACTAGTATCGCTGGTCTTCTAATCGTTAATTTCTTGTTTCGGGAGCAGata  <  2:35895/133‑1 (MQ=255)                                                                                                                                 tggcggCTTTGAACCTGCTGGTGGTGTCACCATGGACTCTGGAGGCGTGCTTTTACAACTACACAGTAAAAGTGCCAGTGCAACTAGTATCGCTGGTCTTCTAATCGTTAATTTCTTGTTTCGGGAGCAGata  <  1:367510/133‑1 (MQ=255)                                                                                                                                    |                                                                                                                                  TTGTTGTACACGGGCCCGGCTGGTGGGGTTACTAATGTTAACCGAGGAAACCGGTGGCCAGTCGAGCCGAGAAACCGTGGCCAGCCAGATGCCACGCATCGGTTGTGACGATTGTT‑‑GCGTCGTGGCTGGTGGCTTTGAACCTGCTGGTGGTGTCACCATGGACTCTGGAGGCGTGCTTTTACAACTACACAGTAAAAGTGCCAGTGCAACTAGTATCGCTGGTCTTCTAATCGTTAATTTCTTGTTTCGGGAGCAGATA  >  NC_000913/1568560‑1568818

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