Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_017635	1,160,228	A→C	L182V (TTG→GTG)	yccM ←	4Fe‑4S binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_017635	1,160,228	0	A	C	100.0%	126.7 / NA	32	L182V (TTG→GTG)	yccM	4Fe‑4S binding protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (21/11);  total (21/11)

AAAGAAATCGAGCATTTTTACGTCCATCACCACGCTGTAGGGCGACAACATAAAATAGTGAATCGCCTGAGCGGGCATCAATAACGCGATA  >  NC_017635/1160181‑1160271                                                |                                            aaagaaaTCGAGCATTTTTACGTCCATCACCACGCTGTAGGGCGacacca                                           <  2:737345/50‑1 (MQ=255)  aagaaaTCGAGCATTTTTACGTCCATCACCACGCTGTAGGGCGACACCAt                                          <  2:29245/50‑1 (MQ=255)      aaTCGAGCATTTTTACGTCCATCACCACGCTGTAGGGCGACACCATaaaa                                      >  1:807049/1‑50 (MQ=255)      aaTCGAGCATTTTTACGTCCATCACCACGCTGTAGGGCGACACCATaaaa                                      >  1:325909/1‑50 (MQ=255)       aTCGAGCATTTTTACGTCCATCACCACGCTGTAGGGCGACACCATAAAAt                                     <  1:1288697/50‑1 (MQ=255)       aTCGAGCATTTTTACGTCCATCACCACGCTGTAGGGCGACACCATAAAAt                                     <  2:960954/50‑1 (MQ=255)       aTCGAGCATTTTTACGTCCATCACCACGCTGTAGGGCGACACCATAAAAt                                     <  2:616555/50‑1 (MQ=255)        tCGAGCATTTTTACGTCCATCACCACGCTGTAGGGCGACACCATAAAATa                                    <  1:1159833/50‑1 (MQ=255)         cGAGCATTTTTACGTCCATCACCACGCTGTAGGGCGACACCATAAAATAg                                   >  1:127094/1‑50 (MQ=255)         cGAGCATTTTTACGTCCATCACCACGCTGTAGGGCGACACCATAAAATAg                                   >  2:464825/1‑50 (MQ=255)          gAGCATTTTTACGTCCATCACCACGCTGTAGGGCGACACCATAAAATAGt                                  <  2:1752761/50‑1 (MQ=255)           aGCATTCTTACGTCCATCACCACGCTGTAGGGCGACACCATAAAATAGTg                                 <  1:80174/50‑1 (MQ=255)              aTTTTTACGTCCATCACCACGCTGTAGGGCGACACCATAAAATAGTGAAt                              >  2:414970/1‑50 (MQ=255)               tttttACGTCCATCACCACGCTGTAGGGCGACACCATAAAATAGTGAATc                             >  1:1986883/1‑50 (MQ=255)                ttttACGTCCATCACCACGCTGTAGGGCGACACCATAAAATAGTGAATCg                            >  2:845056/1‑50 (MQ=255)                 tttACGTCCATCACCACGCTGTAGGGCGACACCATAAAATAGTGAATCGc                           >  1:69008/1‑50 (MQ=255)                   tACGTCCATCACCACGCTGTAGGGCGACACCATAAAATAGTGAATCGCCt                         >  2:1389660/1‑50 (MQ=255)                     cGTCCATCACCACGCTGTAGGGCGACACCATAAAATAGTGAATCGCCt                         >  1:818131/1‑48 (MQ=255)                          aTCACCACGCTGTAGGGCGACACCATAAAATAGTGAATCGCCTGAGCggg                  <  1:998588/50‑1 (MQ=255)                           tCACCACGCTGTAGGGCGACACCATAAAATAGTGAATCGCCTGAGCggg                  >  2:1107845/1‑49 (MQ=255)                             accacGCTGTAGGGCGACACCATAAAATAGTGAATCGCCTGAGCGGGCAt               >  2:1590552/1‑50 (MQ=255)                               cacGCTGTAGGGCGACACCATAAAATAGTGAATCGCCTGAGCGGGCATc              >  2:244639/1‑49 (MQ=255)                                acGCTGTAGGGCGACACCATAAAATAGTGAATCGCCTGAGCGGGCATCaa            >  1:337397/1‑50 (MQ=255)                                acGCTGTAGGGCGACACCATAAAATAGTGAATCGCCTGAGCGGGCATCaa            >  1:2080481/1‑50 (MQ=255)                                 cGCTGTAGGGCGACACCATAAAATAGTGAATCGCCTGAGCGGGCATCAAt           <  2:2054719/50‑1 (MQ=255)                                   cTGTAGGGCGACACCATAAAATAGTGAATCGCCTGAGCGGGCATCAATaa         >  1:1719800/1‑50 (MQ=255)                                    tGTAGGGCGACACCATAAAATAGTGAATCGCCTGAGCGGGCATCAATAAc        >  1:781114/1‑50 (MQ=255)                                       aGGGCGACACCATAAAATAGTGAATCGCCTGAGCGGGCATCAATAAcgc      >  1:1838661/1‑49 (MQ=255)                                        gggCGACACCATAAAATAGTGAATCGCCTGAGCGGGCATCAATAACGCGa    <  1:453260/50‑1 (MQ=255)                                         ggCGACACCATAAAATAGTGAATCGCCTGAGCGGGCATCAATAACGCGa    >  1:207152/1‑49 (MQ=255)                                         ggCGACACCATAAAATAGTGAATCGCCTGAGCGGGCATCAATAACGCGa    >  1:1931328/1‑49 (MQ=255)                                          gCGACACCATAAAATAGTGAATCGCCTGAGCGGGCATCAATAACGCGata  >  2:1526971/1‑50 (MQ=255)                                          gCGACACCATAAAATAGTGAATCGCCTGAGCGGGCATCAATAACGCGata  <  2:408758/50‑1 (MQ=255)                                             acaccaTAAAATAGTGAATCGCCTGAGCGGGCATCAATAACGCGata  >  1:1669987/1‑47 (MQ=255)                                                |                                            AAAGAAATCGAGCATTTTTACGTCCATCACCACGCTGTAGGGCGACAACATAAAATAGTGAATCGCCTGAGCGGGCATCAATAACGCGATA  >  NC_017635/1160181‑1160271

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