Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,949,786	C→A	A62S (GCG→TCG)	yifB ←	putative magnesium chelatase YifB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,949,786	0	C	A	100.0%	175.4 / NA	51	A62S (GCG→TCG)	yifB	putative magnesium chelatase YifB
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (23/28);  total (23/28)

ACTAATTCATATTCATCTAACTTATTGGCTGTAAGCTGTTCTGAGGCCGCCAGCAACGCAATGGCGATAGGTAAATCATATCGTCCCCCTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGATTTTTTTCGCCGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGCTTCTTTTACCGTTGTTTCTGGTAAGCCCACCATCGTTAAGCCGGGTAGACCTTTACTGATATGTACCTCAACAGTGATCGGGG  >  NC_000913/3949650‑3949920                                                                                                                                         |                                                                                                                                       aCTAATTCATATTCATCTAACTTATTGGCTGTAAGCTGTTCTGAGGCCGCCAGCAACGCAATGGCGATAGGTAAATCATATCGTCCCCCTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGATTTTTTTCGACgg                                                                                                                                     <  2:652082/140‑1 (MQ=255) aCTAATTCATATTCATCTAACTTATTGGCTGTAAGCTGTTCTGAGGCCGCCAGCAACGCAATGGCGATAGGTAAATCATATCGTCCCCCTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGATTTTTTTCGACgg                                                                                                                                     <  2:49851/140‑1 (MQ=255) aCTAATTCATATTCATCTAACTTATTGGCTGTAAGCTGTTCTGAGGCCGCCAGCAACGCAATGGCGATAGGTAAATCATATCGTCCCCCTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGATTTTTTTCGACgg                                                                                                                                     <  2:703898/140‑1 (MQ=255)  cTAATTCATATTCATCTAACTTATTGGCTGTAAGCTGTTCTGAGGCCGCCAGCAACGCAATGGCGATAGGTAAATCATATCGTCCCCCTTATTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGATTTTTTTCGACGGa                                                                                                                                    <  1:553508/140‑1 (MQ=255)    aaTTCATATTCATCTAACTTATTGGCTGTAAGCTGTTCTGAGGCCGCCAGCAACGCAATGGCGATAGGTAAATCATATCGTCCCCCTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGATTTTTTTCGACGGata                                                                                                                                  <  2:38414/140‑1 (MQ=255)              catCTAACTTATTGGCTGTAAGCTGTTCTGAGGCCGCCAGCAACGCAATGGCGATAGGTAAATCATATCGTCCCCCTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGATTTTTTTCGACGGATATTCATATCCg                                                                                                                        <  1:212430/140‑1 (MQ=255)              catCTAACTTATTGGCTGTAAGCTGTTCTGAGGCCGCCAGCAACGCAATGGCGATAGGTAAATCATATCGTCCCCCTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGATTTTTTTCGACGGATATTCATATCCg                                                                                                                        <  2:462808/140‑1 (MQ=255)                tCTAACTTATTGGCTGTAAGCTGTTCTGAGGCCGCCAGCAACGCAATGGCGATAGGTAAATCATATCGTCCCCCTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGATTTTTTTCGACGGATATTCATATCCGCt                                                                                                                      >  1:215103/1‑140 (MQ=255)                         ttGGCTGTAAGCTGTTCTGAGGCCGCCAGCAACGCAATGGCGATAGGTAAATCATATCGTCCCCCTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGATTTTTTTCGACGGATATTCATATCCg                                                                                                                        >  1:80342/1‑129 (MQ=255)                         ttGGCTGTAAGCTGTTCTGAGGCCGCCAGCAACGCAATGGCGATAGGTAAATCATATCGTCCCCCTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGATTTTTTTCGACGGATATTCATATCCg                                                                                                                        <  2:80342/129‑1 (MQ=255)                         ttGGCTGTAAGCTGTTCTGAGGCCGCCAGCAACGCAATGGCGATAGGTAAATCATATCGTCCCCCTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGATTTTTTTCGACGGATATTCATATCCGCTATTGATAAt                                                                                                             >  1:848409/1‑140 (MQ=255)                           ggCTGTAAGCTGTTCTGAGGCCGCCAGCAACGCAATGGCGATAGGTAAATCATATCGTCCCCCTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGATTTTTTTCGACGGATATTCATATCCGCTATTGATAATgg                                                                                                           >  2:451886/1‑140 (MQ=255)                               gTAAGCTGTTCTGAGGCCGCCAGCAACGCAATGGCGATAGGTAAATCATATCGTCCCCCTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGATTTTTTTCGACGGTTATTCATATCCGCTATTGATAATGGCGCt                                                                                                       <  2:494359/140‑1 (MQ=255)                                       ttCTGAGGCCGCCAGCAACGCAATGGCGATAGGTAAATCATATCGTCCCCCTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGATTTTTTTCGACGGATATTCATATCCGCTATTGATAATGGCGCTGCGCAcgc                                                                                               >  1:552155/1‑140 (MQ=255)                                                 gccAGCAACGCAATGGCGATAGGTAAATCATATCGTCCCCCTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGATTTTTTTCGACGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGc                                                                                     >  1:319809/1‑140 (MQ=255)                                                 gccAGCAACGCAATGGCGATAGGTAAATCATATCGTCCCCCTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGATTTTTTTCGACGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGc                                                                                     <  2:319809/140‑1 (MQ=255)                                                 gccAGCAACGCAATGGCGATAGGTAAATCATATCGTCCCCCTTCTTTTGGCAGA‑CAGCTGGAGCCAGGTTGATGGTGATTTTTTTCGACGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGct                                                                                    >  2:147432/1‑140 (MQ=255)                                                   cAGCAACGCAATGGCGATAGGTAAATCATATCGTCCCCCTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGATTTTTTTCGACGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGctt                                                                                   >  2:708966/1‑140 (MQ=255)                                                       aaCGCAATGGCGATAGGTAAATCATATCGTCCCCCTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGATTTTTTTCGACGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGCTTCttt                                                                               <  1:451886/140‑1 (MQ=255)                                                           cAATGGCGATAGGTAAATCATATCGTCCCCCTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGATTTTTTTCGACGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGCTTCTTTTAcc                                                                           <  1:210717/140‑1 (MQ=255)                                                            aaTGGCGATAGGTAAATCATATCGTCCCCCTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGATTTTTTTCGACGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGCTTCTTTTACCg                                                                          <  2:357293/140‑1 (MQ=255)                                                                                atCGTCCCCCTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGATTTTTTTCGACGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGCTTCTTTTACCGTTGTTTCTGGTAAGCccacc                                                      >  1:22763/1‑140 (MQ=255)                                                                                atCGTCCCCCTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGATTTTTTTCGACGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGCTTCTTTTACCGTTGTTTCTGGTAAGCccacc                                                      <  2:188642/140‑1 (MQ=255)                                                                                  cGTCCCCCTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGATTTTTTTCGACGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGCTTCTTTTACCGTTGTTTCTGGTAAGCCCACCAt                                                    >  1:339438/1‑140 (MQ=255)                                                                                    tCCCCCTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGATTTTTTTCGACGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGCTTCTTTTACCGTTGTTTCTGGTAAGcc                                                          >  2:211316/1‑132 (MQ=255)                                                                                    tCCCCCTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGATTTTTTTCGACGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGCTTCTTTTACCGTTGTTTCTGGTAAGcc                                                          <  1:211316/132‑1 (MQ=255)                                                                                      ccccTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGATTTTTTTCGACGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGCTTCTTTTACCGTTGTTTCTGGTAAGCCCACCATctgt                                                >  2:841772/1‑137 (MQ=255)                                                                                         cttcttTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGATTTTTTTCGACGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGCTTCTTTTACCGTTGTTTCTGGTAAGCCCACCATCGTTAAg                                             <  2:660931/140‑1 (MQ=255)                                                                                              tttGGCAGATCAGCTGGAGCCAGGTTGATGGTGATTTTTTTCGACGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGCTTCTTTTACCGTTGTTTCTGGTAAGCCCACCATCGTTAAGCCggg                                        >  2:249482/1‑140 (MQ=255)                                                                                               ttGGCAGATCAGCTGGAGCCAGGTTGATGGTGATTTTTTTCGACGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGCTTCTTTTACCGTTGTTTCTGGTAAGCCCACCATCGTTAAGCCGGGt                                       >  1:695567/1‑140 (MQ=255)                                                                                                tGGCAGATCAGCTGGAGCCAGGTTGATGGTGATTTTTTTCGACGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGCTTCTTTTACCGTTGTTTCTGGTAAGCCCACCATCGTTAAGCCGGGTa                                      <  1:249482/140‑1 (MQ=255)                                                                                                   cAGATCAGCTGGAGCCAGGTTGATGGTGATTTTTTTCGACGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGCTTCTTTTACCGTTGTTTCTGGTAAGCCCACCATCGTTAAGCCGGGTAGAc                                   >  2:155057/1‑140 (MQ=255)                                                                                                   cAGATCAGCTGGAGCCAGGTTGATGGTGATTTTTTTCGACGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGCTTCTTTTACCGTTGTTTCTGGTAAGCCCACCATCGTTAAGCCGGGTAGAc                                   <  2:215103/140‑1 (MQ=255)                                                                                                   cAGATCAGCTGGAGCCAGGTTGATGGTGATTTTTTTCGACGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGCTTCTTTTACCGTTGTTTCTGGTAAGCCCACCATCGTTAAGCCGGGTAGAc                                   <  2:695567/140‑1 (MQ=255)                                                                                                   cAGATCAGCTGGAGCCAGGTTGATGGTGATTTTTTTCGACGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGCTTCTTTTACCGTTGTTTCTGGTAAGCCCACCATCGTTAAGCCGGGTAGAc                                   >  1:275929/1‑140 (MQ=255)                                                                                                     gATCAGCTGGAGCCAGGTTGATGGTGATTTTTTTCGACGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGCTTCTTTTACCGTTGTTTCTGGTAAGCCCACCATCGTTAAGCCGGGTAGACCt                                 <  1:792750/140‑1 (MQ=255)                                                                                                     gATCAGCTGGAGCCAGGTTGATGGTGATTTTTTTCGACGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGCTTCTTTTACCGTTGTTTCTGGTAAGCCCACCATCGTTAAGCCGGGTAGACCt                                 <  1:445083/140‑1 (MQ=255)                                                                                                        cAGCTGGAGCCAGGTTGATGGTGATTTTTTTCGACGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGCTTCTTTTACCGTTGTTTCTGGTAAGCCCACCATCGTTAAGCCGGGTAGACCTTTa                              <  1:495360/140‑1 (MQ=255)                                                                                                         aGCTGGAGCCAGGTTGATGGTGATTTTTTTCGACGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGCTTCTTTTACCGTTGTTTCTGGTAAGCCCACCATCGTTAAGCCGGGTAGACCTTTAc                             <  1:421574/140‑1 (MQ=255)                                                                                                           cTGGAGCCAGGTTGATGGTGATTTTTTTCGACGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGCTTCTTTTACCGTTGTTTCTGGTAAGCCCACCATCGTTAAGCCGGGTAGACCTTTACTg                           >  2:550811/1‑140 (MQ=255)                                                                                                              gAGCCAGGTTGATGGTGATTTTTTTCGACGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGCTTCTTTTACCGTTGTTTCTGGTAAGCCCACCATCGTTAAGCCGGGTAGACCTTTACTGata                        <  1:155057/140‑1 (MQ=255)                                                                                                                gCCAGGTTGATGGTGATTTTTTTCGACGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGCTTCTTTTACCGTTGTTTCTGGTAAGCCCACCATCGTTAAGCCGGGTAGAc                                   <  2:527407/127‑1 (MQ=255)                                                                                                                gCCAGGTTGATGGTGATTTTTTTCGACGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGCTTCTTTTACCGTTGTTTCTGGTAAGCCCACCATCGTTAAGCCGGGTAGAc                                   >  1:527407/1‑127 (MQ=255)                                                                                                                gCCAGGTTGATGGTGATTTTTTTCGACGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGCTTCTTTTACCGTTGTTTCTGGTAAGCCCACCATCGTTAAGCCGGGTAGACCTTTACTGATATg                      <  2:848409/140‑1 (MQ=255)                                                                                                                     gTTGATGGTGATTTTTTTCGACGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGCTTCTTTTACCGTTGTTTCTGGTAAGCCCACCATCGTTAAGCCGGGTAGACCTTTACTGATATGTACCt                 >  1:629271/1‑140 (MQ=255)                                                                                                                     gTTGATGGTGATTTTTTTCGACGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGCTTCTTTTACCGTTGTTTCTGGTAAGCCCACCATCGTTAAGCCGGGTAGACCTTTACTGATATGTACCt                 >  1:616413/1‑140 (MQ=255)                                                                                                                     gTTGATGGTGATTTTTTTCGACGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGCTTCTTTTACCGTTGTTTCTGGTAAGCCCACCATCGTTAAGCCGGGTAGACCTTTACTGATATGTACCt                 >  2:87369/1‑140 (MQ=255)                                                                                                                      ttGATGGTGATTTTTTTCGACGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGCTTCTTTTACCGTTGTTTCTGGTAAGCCCACCATCGTTAAGCCGGGTAGACCTTTACTGATATGTACCTc                <  2:629271/140‑1 (MQ=255)                                                                                                                               aTTTTTTTCGACGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGCTTCTTTTACCGTTGTTTCTGGTAAGCCCACCATCGTTAAGCCGGGTAGACCTTTACTGATATGTACCTCAACAGTGAt       >  1:129187/1‑140 (MQ=255)                                                                                                                                   ttttCGACGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGCTTCTTTTACCGTTGTTTCTGGTAAGCCCACCATCGTTAAGCCGGGTAGACCTTTACTGATATGTACCTCAACAGTGATCggg   <  2:714753/140‑1 (MQ=255)                                                                                                                                    tttCGACGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGCTTCTTTTACCGTTGTTTCTGGTAAGCCCACCATCGTTAAGCCGGGTAGACCTTTACTGATATGTACCTCAACAGTGATCgggg  >  1:565727/1‑140 (MQ=255)                                                                                                                                         |                                                                                                                                       ACTAATTCATATTCATCTAACTTATTGGCTGTAAGCTGTTCTGAGGCCGCCAGCAACGCAATGGCGATAGGTAAATCATATCGTCCCCCTTCTTTTGGCAGATCAGCTGGAGCCAGGTTGATGGTGATTTTTTTCGCCGGATATTCATATCCGCTATTGATAATGGCGCTGCGCACGCGATCGCGAGCTTCTTTTACCGTTGTTTCTGGTAAGCCCACCATCGTTAAGCCGGGTAGACCTTTACTGATATGTACCTCAACAGTGATCGGGG  >  NC_000913/3949650‑3949920

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