Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,052,605	T→G	D253A (GAC→GCC)	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,605	0	T	G	100.0%	175.6 / NA	46	D253A (GAC→GCC)	yccM	putative 4Fe‑4S membrane protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (28/18);  total (28/18)

TGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCACTGTTCGTACCTGAATTAATTTATCGACCGGGATCCGTGATGGGCAATTTTTTGCGCATTTGCCACAGTCGATACAACTTTCGGCATTGCGACGAATCTTAAACGGCGACAATAGCGAAACCACGCCCATCAGCGCGCCATATGGGCAAAGATAACGACACCAGGCATGGCGAATAAACAGGCTGGCAATCAGCAAAACGG  >  NC_000913/1052472‑1052737                                                                                                                                      |                                                                                                                                     tGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCACTGTTCGTACCTGAATTAATTTATCGACCGGGATCCGTGATGGGCAATTTTTTGCGCATTTGCCACAGGCGATAc                                                                                                                                >  1:536476/1‑140 (MQ=255)            aTGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCACTGTTCGTACCTGAATTAATTTATCGACCGGGATCCGTGATGGGCAATTTTTTGCGCATTTGCCACAGGCGATACAACTTTCGGCa                                                                                                                     <  2:650185/140‑1 (MQ=255)             tGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCACTGTTCGTACCTGAATTAATTTATCGACCGGGATCCGTGATGGGCAATTTTTTGCGCATTTGCCACAGGCGATACAACTTTCGGCAt                                                                                                                    >  2:11127/1‑140 (MQ=255)                tGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCACTGTTCGTACCTGAATTAATTTATCGACCGGGATCCGTGATGGGCAATTTTTTGCGCATTTGCCACAGGCGATACAACTTTCGGCATTGc                                                                                                                 >  1:870819/1‑140 (MQ=255)                              aTGACTCTACGCAAGTCATACAGCCGGTACATTCCACTGTTCGTACCTGAATTAATTTATCGACCGGGATCCGTGATGGGCAATTTTTTGCGCATTTGCCACAGGCGATACAACTTTCGGCATTGCGACGAATCTTAAAc                                                                                                   <  1:490415/140‑1 (MQ=255)                                           aGTCATACAGCCGGTACATTCCACTGTTCGTACCTGAATTAATTTATCGACCGGGATCCGTGATGGGCAATTTTTTGCGCATTTGCCACAGGCGATACAACTTTCGGCATTGCGACGAATCTTAAACGGCGACAATAGCg                                                                                      <  2:359305/140‑1 (MQ=255)                                            gTCATACAGCCGGTACATTCCACTGTTCGTACCTGAATTAATTTATCGACCGGGATCCGTGATGGGCAATTTTTTGCGCATTTGCCACAGGCGATACAACTTTCGGCATTGCGACGAATCTTAAACGGCGACAATAGCGa                                                                                     >  1:204583/1‑140 (MQ=255)                                            gTCATACAGCCGGTACATTCCACTGTTCGTACCTGAATTAATTTATCGACCGGGATCCGTGATGGGCAATTTTTTGCGCATTTGCCACAGGCGATACAACTTTCGGCATTGCGACGAATCTTAAACGGCGACAATAGCGa                                                                                     >  1:806430/1‑140 (MQ=255)                                            gTCATACAGCCGGTACATTCCACTGTTCGTACCTGAATTAATTTATCGACCGGGATCCGTGATGGGCAATTTTTTGCGCATTTGCCACAGGCGATACAACTTTCGGCATTGCGACGAATCTTAAACGGCGACAATAGCGa                                                                                     >  2:260316/1‑140 (MQ=255)                                              cATACAGCCGGTACATTCCACTGTTCGTACCTGAATTAATTTATCGACCGGGATCCGTGATGGGCAATTTTTTGCGCATTTGCCACAGGCGATACAACTTTCGGCATTGCGACGAATCTTAAACGGCGACAATAGCGaaa                                                                                   <  2:870819/140‑1 (MQ=255)                                                   aGCCGGTACATTCCACTGTTCGTACCTGAATTAATTTATCGACCGGGATCCGTGATGGGCAATTTTTTGCGCATTTGCCACAGGCGATACAACTTTCGGCATTGCGACGAATCTTAAACGGCGACAATAGCGAAACCACg                                                                              <  2:204583/140‑1 (MQ=255)                                                    gCCGGTACATTCCACTGTTCGTACCTGAATTAATTTATCGACCGGGATCCGTGATGGGCAATTTTTTGCGCATTTGCCACAGGCGATACAACTTTCGGCATTGCGACGAATCTTAAACGGCGACAATAGCGAAACCACGc                                                                             >  2:241704/1‑140 (MQ=255)                                                      cGGTACATTCCACTGTTCGTACCTGAATTAATTTATCGACCGGGATCCGTGATGGGCAATTTTTTGCGCATTTGCCACAGGCGATACAACTTTCGGCATTGCGACGAATCTTAAACGGCGACAATAGCGAAACCACGccc                                                                           >  2:126812/1‑140 (MQ=255)                                                      cGGTACATTCCACTGTTCGTACCTGAATTAATTTATCGACCGGGATCCGTGATGGGCAATTTTTTGCGCATTTGCCACAGGCGATACAACTTTCGGCATTGCGACGAATCTTAAACGGCGACAATAGCGAAACCACGccc                                                                           >  2:389829/1‑140 (MQ=255)                                                      cGGTACATTCCACTGTTCGTACCTGAATTAATTTATCGACCGGGATCCGTGATGGGCAATTTTTTGCGCATTTGCCACAGGCGATACAACTTTCGGCATTGCGACGAATCTTAAACGGCGACAATAGCGAAACCACGccc                                                                           >  1:535225/1‑140 (MQ=255)                                                      cGGTACATTCCACTGTTCGTACCTGAATTAATTTATCGACCGGGATCCGTGATGGGCAATTTTTTGCGCATTTGCCACAGGCGATACAACTTTCGGCATTGCGACGAATCTTAAACGGCGACAATAGCGAAACCACGccc                                                                           >  1:1006305/1‑140 (MQ=255)                                                           cATTCCACTGTTCGTACCTGAATTAATTTATCGACCGGGATCCGTGATGGGCAATTTTTTGCGCATTTGCCACAGGCGATACAACTTTCGGCATTGCGACGAATCTTAAACGGCGACAATAGCGAAACCACGCCCATCAg                                                                      <  1:357583/140‑1 (MQ=255)                                                             ttCCACTGTTCGTACCTGAATTAATTTATCGACCGGGATCCGTGATGGGCAATTTTTTGCGCATTTGCCACAGGCGATACAACTTTCGGCATTGCGACGAATCTTAAACGGCGACAATAGCGAAACCACGCCCATCAgcg                                                                    <  2:806430/140‑1 (MQ=255)                                                              tCCACTGTTCGTACCTGAATTAATTTATCGACCGGGATCCGTGATGGGCAATTTTTTGCGCATTTGCCACAGGCGATACAACTTTCGGCATTGCGACGAATCTTAAACGGCGACAATAGCGAAACCACGCCCATCAgcgc                                                                   <  2:419262/140‑1 (MQ=255)                                                              tCCACTGTTCGTACCTGAATTAATTTATCGACCGGGATCCGTGATGGGCAATTTTTTGCGCATTTGCCACAGGCGATACAACTTTCGGCATTGCGACGAATCTTAAACGGCGACAATAGCGAAACCACGCCCATCAgcgc                                                                   <  2:536476/140‑1 (MQ=255)                                                                cACTGTTCGTACCTGAATTAATTTATCGACCGGGATCCGTGATGGGCAATTTTTTGCGCATTTGCCACAGGCGATACAACTTTCGGCATTGCGACGAATCTTAAACGGCGACAATAGCGAAACCACGCCCATCAgcgcgc                                                                 >  2:93171/1‑140 (MQ=255)                                                                  cTGTTCGTACCTGAATTAATTTATCGACCGGGATCCGTGATGGGCAATTTTTTGCGCATTTGCCACAGGCGATACAACTTTCGGCATTGCGACGAATCTTAAACGGCGACAATAGCGAAACCACGCCCATCAGCGCGCCa                                                               <  1:11127/140‑1 (MQ=255)                                                                   tGTTCGTACCTGAATTAATTTATCGACCGGGATCCGTGATGGGCAATTTTTTGCGCATTTGCCACAGGCGATACAACTTTCGGCATTGCGACGAATCTTAAACGGCGACAATAGCGAAACCACGCCCATCAGCGCGCCat                                                              >  2:219146/1‑140 (MQ=255)                                                                   tGTTCGTACCTGAATTAATTTATCGACCGGGATCCGTGATGGGCAATTTTTTGCGCATTTGCCACAGGCGATACAACTTTCGGCATTGCGACGAATCTTAAACGGCGACAATAGCGAAACCACGCCCATCAGCGCGCCat                                                              >  2:118108/1‑140 (MQ=255)                                                                   tGTTCGTACCTGAATTAATTTATCGACCGGGATCCGTGATGGGCAATTTTTTGCGCATTTGCCACAGGCGATACAACTTTCGGCATTGCGACGAATCTTAAACGGCGACAATAGCGAAACCACGCCCATCAGCGCGCCat                                                              >  1:80903/1‑140 (MQ=255)                                                                   tGTTCGTACCTGAATTAATTTATCGACCGGGATCCGTGATGGGCAATTTTTTGCGCATTTGCCACAGGCGATACAACTTTCGGCATTGCGACGAATCTTAAACGGCGACAATAGCGAAACCACGCCCATCAGCGCGCCat                                                              >  1:724947/1‑140 (MQ=255)                                                                     ttCGTACCTGAATTAATTTATCGACCGGGATCCGTGATGGGCAATTTTTTGCGCATTTGCCACAGGCGATACAACTTTCGGCATTGCGACGAATCTTAAACGGCGACAATAGCGAAACCACGCCCATCAgcgcg                                                                  >  1:511018/1‑134 (MQ=255)                                                                     ttCGTACCTGAATTAATTTATCGACCGGGATCCGTGATGGGCAATTTTTTGCGCATTTGCCACAGGCGATACAACTTTCGGCATTGCGACGAATCTTAAACGGCGACAATAGCGAAACCACGCCCATCAgcgcg                                                                  <  2:511018/134‑1 (MQ=255)                                                                     ttCGTACCTGAATTAATTTATCGACCGGGATCCGTGATGGGCAATTTTTTGCGCATTTGCCACAGGCGATACAACTTTCGGCATTGCGACGAATCTTAAACGGCGACAATAGCGAAACCACGCCCATCAGCGCGCCatat                                                            >  1:629405/1‑140 (MQ=255)                                                                          aCCTGAATTAATTTATCGACCGGGATCCGTGATGGGCAATTTTTTGCGCATTTGCCACAGGCGATACAACTTTCGGCATTGCGACGAATCTTAAACGGCGACAATAGCGAAACCACGCCCATCa                                                                       <  1:961739/124‑1 (MQ=255)                                                                          aCCTGAATTAATTTATCGACCGGGATCCGTGATGGGCAATTTTTTGCGCATTTGCCACAGGCGATACAACTTTCGGCATTGCGACGAATCTTAAACGGCGACAATAGCGAAACCACGCCCATCa                                                                       >  2:961739/1‑124 (MQ=255)                                                                               aattaattTATCGACCGGGATCCGTGATGGGCAATTTTTTGCGCATTTGCCACAGGCGATACAACTTTCGGCATTGCGACGAATCTTAAACGGCGACAATAGCGAAACCACGCCCATCAGCGCGCCATATGGGCAAAGAt                                                  >  2:36465/1‑140 (MQ=255)                                                                                 ttaattTATCGACCGGGATCCGTGATGGGCAATTTTTTGCGCATTTGCCACAGGCGATACAACTTTCGGCATTGCGACGAATCTTAAACGGCGACAATAGCGAAACCACGCCCATCAGCGCGCCATATGGGCAAAGATaa                                                <  2:381360/140‑1 (MQ=255)                                                                                 ttaattTATCGACCGGGATCCGTGATGGGCAATTTTTTGCGCATTTGCCACAGGCGATACAACTTTCGGCATTGCGACGAATCTTAAACGGCGACAATAGCGAAACCACGCCCATCAGCGCGCCATATGGGCAAAGATaa                                                <  1:619090/140‑1 (MQ=255)                                                                                    attTATCGACCGGGATCCGTGATGGGCAATTTTTTGCGCATTTGCCACAGGCGATACAACTTTCGGCATTGCGACGAATCTTAAACGGCGACAATAGCGAAACCACGCCCATCAGCGCGCCATATGGGCAAAGATAACGa                                             <  1:693229/140‑1 (MQ=255)                                                                                              cGGGATCCGTGATGGGCAATTTTTTGCGCATTTGCCACAGGCGATACAACTTTCGGCATTGCGACGAATCTTAAACGGCGACAATAGCGAAACCACGCCCATCAGCGCGCCATATGGGCAAAGATAACGACACCAGGCAt                                   <  1:241704/140‑1 (MQ=255)                                                                                                       tGATGGGCAATTTTTTGCGCATTTGCCACAGGCGATACAACTTTCGGCATTGCGACGAATCTTAAACGGCGACAATAGCGAAACCACGCCCATCAGCGCGCCATATGGGCAAAGATAACGACACCAGGCATGGCGAATaa                          >  2:135733/1‑140 (MQ=255)                                                                                                         aTGGGCAATTTTTTGCGCATTTGCCACAGGCGATACAACTTTCGGCATTGCGACGAATCTTAAACGGCGACAATAGCGAAACCACGCCCATCAGCGCGCCATATGGGCAAAGATAACGACACCAGGCATGGCGAATAAAc                        >  1:611675/1‑140 (MQ=255)                                                                                                           gggCAATTTTTTGCGCATTTGCCACAGGCGATACAACTTTCGGCATTGCGACGAATCTTAAACGGCGACAATAGCGAAACCACGCCCATCAGCGCGCCATATGGGCAAAGATAACGACACCAGGCATGGCGAATAAACAg                      >  1:503334/1‑140 (MQ=255)                                                                                                            ggCAATTTTTTGCGCATTTGCCACAGGCGATACAACTTTCGGCATTGCGACGAATCTTAAACGGCGACAATAGCGAAACCACGCCCATCAGCGCGCCATATGGGCAAAGATAACGACACCAGGCATGGCGAATAAACAgg                     <  2:503334/140‑1 (MQ=255)                                                                                                            ggCAATTTTTTGCGCATTTGCCACAGGCGATACAACTTTCGGCATTGCGACGAATCTTAAACGGCGACAATAGCGAAACCACGCCCATCAGCGCGCCATATGGGCAAAGATAACGACACCAGGCATGGCGAATAAACAgg                     >  2:488491/1‑140 (MQ=255)                                                                                                            ggCAATTTTTTGCGCATTTGCCACAGGCGATACAACTTTCGGCATTGCGACGAATCTTAAACGGCGACAATAGCGAAACCACGCCCATCAGCGCGCCATATGGGCAAAGATAACGACACCAGGCATGGCGAATAAACAgg                     >  1:504066/1‑140 (MQ=255)                                                                                                             gCAATTTTTTGCGCATTTGCCACAGGCGATACAACTTTCGGCATTGCGACGAATCTTAAACGGCGACAATAGCGAAACCACGCCCATCAGCGCGCCATATGGGCAAAGATAACGACACCAGGCATGGCGAATAAACAGGc                    >  2:765533/1‑140 (MQ=255)                                                                                                                      tGCGCATTTGCCACAGGCGATACAACTTTCGGCATTGCGACGAATCTTAAACGGCGACAATAGCGAAACCACGCCCATCAGCGCGCCATATGGGCAAAGATAACGACACCAGGCATGGCGAATAAACAGGCTGGCAATCa           >  2:147160/1‑140 (MQ=255)                                                                                                                           aTTTGCCACAGGCGATACAACTTTCGGCATTGCGACGAATCTTAAACGGCGACAATAGCGAAACCACGCCCATCAGCGCGCCATATGGGCAAAGATAACGACACCAGGCATGGCGAATAAACAGGCTGGCAATCAGCaaa      >  1:254405/1‑140 (MQ=255)                                                                                                                               gCCACAGGCGATACAACTTTCGGCATTGCGACGAATCTTAAACGGCGACAATAGCGAAACCACGCCCATCAGCGCGCCATATGGGCAAAGATAACGACACCAGGCATGGCGAATAAACAGGCTGGCAATCAGCAAAACgg  <  1:93171/140‑1 (MQ=255)                                                                                                                                      |                                                                                                                                     TGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCACTGTTCGTACCTGAATTAATTTATCGACCGGGATCCGTGATGGGCAATTTTTTGCGCATTTGCCACAGTCGATACAACTTTCGGCATTGCGACGAATCTTAAACGGCGACAATAGCGAAACCACGCCCATCAGCGCGCCATATGGGCAAAGATAACGACACCAGGCATGGCGAATAAACAGGCTGGCAATCAGCAAAACGG  >  NC_000913/1052472‑1052737

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