Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,052,330	A→G	Y345H (TAC→CAC)	yccM ←	putative 4Fe‑4S membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,052,330	0	A	G	100.0%	59.4 / NA	18	Y345H (TAC→CAC)	yccM	putative 4Fe‑4S membrane protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (5/13);  total (5/13)

TGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAA  >  NC_000913/1052207‑1052441                                                                                                                            |                                                                                                                tGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTGCAATtcctc                                                                                                        <  1:385633/133‑1 (MQ=255)       ttATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTGCAATTCCTCCGGCAc                                                                                                  <  1:434872/133‑1 (MQ=255)        tATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTGCAATTCCTCCGGCAcc                                                                                                 <  1:296140/133‑1 (MQ=255)             tgTCAAAATCAAAATAATGTGATTTTGTGAACGTCACCCCGCGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTGCAATTCCTCCGGCACCGGGCt                                                                                            <  2:226629/133‑1 (MQ=255)              gTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTGCAATTCCTCCGGCACCGGGCtt                                                                                           <  1:309564/133‑1 (MQ=255)                caaaatcaaaatAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTGCAATTCCTCCGGCACCGGGCTTTg                                                                                         <  2:476517/133‑1 (MQ=255)                     tcaaaatAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTGCAATTCCTCCGGCACCGGGCTTTGCCATa                                                                                    <  2:6171/133‑1 (MQ=255)                              tgtgATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTGCAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCAt                                                                           >  2:353232/1‑133 (MQ=255)                                         gAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTGCAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTaa                                                                >  2:361224/1‑133 (MQ=255)                                                 ccccGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTGCAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCAc                                                        <  1:361224/133‑1 (MQ=255)                                                      tGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTGCAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCaa                                                   >  1:330163/1‑133 (MQ=255)                                                          aGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTGCAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGTATACATTGCGTAACCAATCACCGCAAACAt                                               >  1:440743/1‑133 (MQ=255)                                                           ggTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTGCAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATa                                              <  1:353232/133‑1 (MQ=255)                                                            gTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTGCAATTCCTCCGGCACCGGGCTTTGCCATACCCCCGCATACATTGCGTAACCAATCACCGCAAACATaa                                             <  2:338181/133‑1 (MQ=255)                                                            gTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATCGGCGCTTGTGGAATTAAGCGTCGGTGCAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATaa                                             <  2:351353/133‑1 (MQ=255)                                                                tGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTGCAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATccc                                         >  1:380384/1‑133 (MQ=255)                                                                         tgttgCTAGTTTAGTGACTAATCATTGGCGATTGTGGAATTAAGCGTCGGTGCAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAg                                <  1:242328/133‑1 (MQ=255)                                                                                                       cTTGTGGAATTAAGCGTCGGTGCAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGTATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCaaa  <  2:440743/133‑1 (MQ=255)                                                                                                                            |                                                                                                                TGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAA  >  NC_000913/1052207‑1052441

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: ‑